New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 6639 for branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC – NEMO

Ignore:
Timestamp:
2016-05-27T14:58:40+02:00 (8 years ago)
Author:
frrh
Message:

Merge NERC/dev_r5518_NOC_MEDUSA_Stable rev 5736 to 6509 inclusive

Location:
branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC/par_cfc.F90

    r6636 r6639  
    1515   USE par_pisces , ONLY : jp_pisces_trd   !: number of biological diag in PISCES 
    1616 
     17   USE par_medusa , ONLY : jp_medusa       !: number of tracers in MEDUSA 
     18   USE par_medusa , ONLY : jp_medusa_2d    !: number of 2D diag in MEDUSA 
     19   USE par_medusa , ONLY : jp_medusa_3d    !: number of 3D diag in MEDUSA 
     20   USE par_medusa , ONLY : jp_medusa_trd   !: number of biological diag in MEDUSA 
     21 
     22   USE par_idtra  , ONLY : jp_idtra        !: number of tracers in MEDUSA 
     23   USE par_idtra  , ONLY : jp_idtra_2d     !: number of tracers in MEDUSA 
     24   USE par_idtra  , ONLY : jp_idtra_3d     !: number of tracers in MEDUSA 
     25   USE par_idtra  , ONLY : jp_idtra_trd    !: number of tracers in MEDUSA 
     26 
    1727   IMPLICIT NONE 
    1828 
    19    INTEGER, PARAMETER ::   jp_lc      =  jp_pisces     !: cumulative number of passive tracers 
    20    INTEGER, PARAMETER ::   jp_lc_2d   =  jp_pisces_2d  !: 
    21    INTEGER, PARAMETER ::   jp_lc_3d   =  jp_pisces_3d  !: 
    22    INTEGER, PARAMETER ::   jp_lc_trd  =  jp_pisces_trd !: 
     29   INTEGER, PARAMETER ::   jp_lc      =  jp_pisces     + jp_medusa     + & 
     30                      jp_idtra     !: cumulative number of passive tracers 
     31   INTEGER, PARAMETER ::   jp_lc_2d   =  jp_pisces_2d  + jp_medusa_2d  + & 
     32                      jp_idtra_2d !: 
     33   INTEGER, PARAMETER ::   jp_lc_3d   =  jp_pisces_3d  + jp_medusa_3d  + & 
     34                      jp_idtra_3d !: 
     35   INTEGER, PARAMETER ::   jp_lc_trd  =  jp_pisces_trd + jp_medusa_trd + & 
     36                      jp_idtra_trd !: 
    2337    
    2438#if defined key_cfc 
  • branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

    r6636 r6639  
    1515   !!   cfc_init     :  sets constants for CFC surface forcing computation 
    1616   !!---------------------------------------------------------------------- 
     17   USE dom_oce       ! ocean space and time domain 
    1718   USE oce_trc       ! Ocean variables 
    1819   USE par_trc       ! TOP parameters 
     
    176177         !                                                  !----------------! 
    177178      END DO                                                !  end CFC loop  ! 
    178       ! 
    179       IF( lrst_trc ) THEN 
    180          IF(lwp) WRITE(numout,*) 
    181          IF(lwp) WRITE(numout,*) 'trc_sms_cfc : cumulated input function fields written in ocean restart file ',   & 
    182             &                    'at it= ', kt,' date= ', ndastp 
    183          IF(lwp) WRITE(numout,*) '~~~~' 
    184          DO jn = jp_cfc0, jp_cfc1 
    185             CALL iom_rstput( kt, nitrst, numrtw, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jn) ) 
    186          END DO 
    187       ENDIF                                             
     179         ! 
     180      IF( kt == nittrc000 ) THEN 
     181         DO jl = 1, jp_cfc    
     182             WRITE(NUMOUT,*) ' ' 
     183             WRITE(NUMOUT,*) 'CFC interpolation verification '  !! Jpalm   
     184             WRITE(NUMOUT,*) '################################## ' 
     185             WRITE(NUMOUT,*) ' ' 
     186               if (jl.EQ.1) then 
     187                   WRITE(NUMOUT,*) 'Traceur = CFC11: ' 
     188               elseif (jl.EQ.2) then 
     189                   WRITE(NUMOUT,*) 'Traceur = CFC12: ' 
     190               endif 
     191             WRITE(NUMOUT,*) 'nyear    = ', nyear 
     192             WRITE(NUMOUT,*) 'nmonth   = ', nmonth 
     193             WRITE(NUMOUT,*) 'iyear_beg= ', iyear_beg 
     194             WRITE(NUMOUT,*) 'iyear_end= ', iyear_end 
     195             WRITE(NUMOUT,*) 'p_cfc(iyear_beg)= ',p_cfc(iyear_beg, 1, jl) 
     196             WRITE(NUMOUT,*) 'p_cfc(iyear_end)= ',p_cfc(iyear_end, 1, jl) 
     197             WRITE(NUMOUT,*) 'Im1= ',im1 
     198             WRITE(NUMOUT,*) 'Im2= ',im2 
     199             WRITE(NUMOUT,*) 'zpp_cfc = ',zpp_cfc 
     200             WRITE(NUMOUT,*) ' ' 
     201         END DO   
     202# if defined key_debug_medusa 
     203         CALL flush(numout) 
     204# endif 
     205      ENDIF 
     206        ! 
     207      !IF( lrst_trc ) THEN 
     208      !   IF(lwp) WRITE(numout,*) 
     209      !   IF(lwp) WRITE(numout,*) 'trc_sms_cfc : cumulated input function fields written in ocean restart file ',   & 
     210      !      &                    'at it= ', kt,' date= ', ndastp 
     211      !   IF(lwp) WRITE(numout,*) '~~~~' 
     212      !   DO jn = jp_cfc0, jp_cfc1 
     213      !      CALL iom_rstput( kt, nitrst, numrtw, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jn) ) 
     214      !   END DO 
     215      !ENDIF                                             
    188216      ! 
    189217      IF( lk_iomput ) THEN 
     
    203231      END IF 
    204232      ! 
     233# if defined key_debug_medusa 
     234      IF(lwp) WRITE(numout,*) '   CFC - Check: nn_timing = ', nn_timing 
     235      CALL flush(numout) 
     236# endif 
    205237      IF( nn_timing == 1 )  CALL timing_stop('trc_sms_cfc') 
    206238      ! 
     
    254286      sca(4,2) =  -0.067430 
    255287 
    256       IF( ln_rsttr ) THEN 
    257          IF(lwp) WRITE(numout,*) 
    258          IF(lwp) WRITE(numout,*) ' Read specific variables from CFC model ' 
    259          IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 
    260          ! 
    261          DO jn = jp_cfc0, jp_cfc1 
    262             CALL iom_get( numrtr, jpdom_autoglo, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jn) )  
    263          END DO 
    264       ENDIF 
     288      !IF( ln_rsttr ) THEN 
     289      !   IF(lwp) WRITE(numout,*) 
     290      !   IF(lwp) WRITE(numout,*) ' Read specific variables from CFC model ' 
     291      !   IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 
     292      !   ! 
     293      !   DO jn = jp_cfc0, jp_cfc1 
     294      !      CALL iom_get( numrtr, jpdom_autoglo, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jn) )  
     295      !   END DO 
     296      !ENDIF 
    265297      IF(lwp) WRITE(numout,*) 
    266298      ! 
Note: See TracChangeset for help on using the changeset viewer.