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 6225 for branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z – NEMO

Ignore:
Timestamp:
2016-01-08T10:35:19+01:00 (8 years ago)
Author:
jamesharle
Message:

Update MPP_BDY_UPDATE branch to be consistent with head of trunk

Location:
branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90

    • Property svn:keywords set to Id
    r4624 r6225  
    2121   USE lbclnk          !  
    2222   USE prtctl_trc      ! Print control for debbuging 
    23    USE trdmod_oce 
    24    USE trdmod_trc 
     23   USE trd_oce 
     24   USE trdtrc 
    2525   USE iom 
    2626    
    … …  
    5959   REAL(wp) ::   fdbod      ! zooplankton mortality fraction that goes to detritus 
    6060 
    61    !!* Substitution 
    62 #  include "top_substitute.h90" 
     61   !! * Substitutions 
     62#  include "vectopt_loop_substitute.h90" 
    6363   !!---------------------------------------------------------------------- 
    6464   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    65    !! $Id: p2zbio.F90 3294 2012-01-28 16:44:18Z rblod $  
     65   !! $Id$  
    6666   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    6767   !!---------------------------------------------------------------------- 
    68  
    6968CONTAINS 
    7069 
    … …  
    110109      IF( nn_timing == 1 )  CALL timing_start('p2z_bio') 
    111110      ! 
    112       IF( ln_diatrc ) THEN 
     111      IF( ln_diatrc .OR. lk_iomput ) THEN 
    113112         CALL wrk_alloc( jpi, jpj,     17, zw2d ) 
    114113         CALL wrk_alloc( jpi, jpj, jpk, 3, zw3d ) 
    … …  
    122121 
    123122      xksi(:,:) = 0.e0        ! zooplakton closure ( fbod) 
    124       IF( ln_diatrc ) THEN 
     123      IF( ln_diatrc .OR. lk_iomput ) THEN 
    125124         zw2d  (:,:,:) = 0.e0 
    126125         zw3d(:,:,:,:) = 0.e0 
    … …  
    186185               !    closure : flux grazing is redistributed below level jpkbio 
    187186               zzoobod = tmminz * zzoo * zzoo 
    188                xksi(ji,jj) = xksi(ji,jj) + (1-fdbod) * zzoobod * fse3t(ji,jj,jk) 
     187               xksi(ji,jj) = xksi(ji,jj) + (1-fdbod) * zzoobod * e3t_n(ji,jj,jk) 
    189188               zboddet = fdbod * zzoobod 
    190189 
    … …  
    239238                  !  trend number 17 in p2zexp 
    240239                ENDIF 
    241                 IF( ln_diatrc ) THEN 
     240                IF( ln_diatrc .OR. lk_iomput ) THEN 
    242241                  ! convert fluxes in per day 
    243                   ze3t = fse3t(ji,jj,jk) * 86400. 
     242                  ze3t = e3t_n(ji,jj,jk) * 86400._wp 
    244243                  zw2d(ji,jj,1)  = zw2d(ji,jj,1)  + zno3phy * ze3t 
    245244                  zw2d(ji,jj,2)  = zw2d(ji,jj,2)  + znh4phy * ze3t 
    … …  
    360359                  !  trend number 17 in p2zexp  
    361360                ENDIF 
    362                 IF( ln_diatrc ) THEN 
     361                IF( ln_diatrc .OR. lk_iomput ) THEN 
    363362                  ! convert fluxes in per day 
    364                   ze3t = fse3t(ji,jj,jk) * 86400. 
     363                  ze3t = e3t_n(ji,jj,jk) * 86400._wp 
    365364                  zw2d(ji,jj,1)  = zw2d(ji,jj,1)  + zno3phy * ze3t 
    366365                  zw2d(ji,jj,2)  = zw2d(ji,jj,2)  + znh4phy * ze3t 
    … …  
    381380                  zw2d(ji,jj,17) = zw2d(ji,jj,17) + zdetdom * ze3t 
    382381                  !    
    383                   zw3d(ji,jj,jk,1) = zno3phy * 86400 
    384                   zw3d(ji,jj,jk,2) = znh4phy * 86400 
    385                   zw3d(ji,jj,jk,3) = znh4no3 * 86400 
     382                  zw3d(ji,jj,jk,1) = zno3phy * 86400._wp 
     383                  zw3d(ji,jj,jk,2) = znh4phy * 86400._wp 
     384                  zw3d(ji,jj,jk,3) = znh4no3 * 86400._wp 
    386385                   ! 
    387386                ENDIF 
    … …  
    390389      END DO 
    391390 
    392       IF( ln_diatrc ) THEN 
    393          ! 
     391      IF( ln_diatrc .OR. lk_iomput ) THEN 
    394392         DO jl = 1, 17  
    395393            CALL lbc_lnk( zw2d(:,:,jl),'T', 1. ) 
    … …  
    398396            CALL lbc_lnk( zw3d(:,:,:,jl),'T', 1. ) 
    399397         END DO 
    400          IF( lk_iomput ) THEN 
     398      ENDIF 
     399      IF( lk_iomput ) THEN 
    401400            ! Save diagnostics 
    402             CALL iom_put( "TNO3PHY", zw2d(:,:,1) ) 
    403             CALL iom_put( "TNH4PHY", zw2d(:,:,2) ) 
    404             CALL iom_put( "TPHYDOM", zw2d(:,:,3) ) 
    405             CALL iom_put( "TPHYNH4", zw2d(:,:,4) ) 
    406             CALL iom_put( "TPHYZOO", zw2d(:,:,5) ) 
    407             CALL iom_put( "TPHYDET", zw2d(:,:,6) ) 
    408             CALL iom_put( "TDETZOO", zw2d(:,:,7) ) 
    409             CALL iom_put( "TZOODET", zw2d(:,:,8) ) 
    410             CALL iom_put( "TZOOBOD", zw2d(:,:,9) ) 
    411             CALL iom_put( "TZOONH4", zw2d(:,:,10) ) 
    412             CALL iom_put( "TZOODOM", zw2d(:,:,11) ) 
    413             CALL iom_put( "TNH4NO3", zw2d(:,:,12) ) 
    414             CALL iom_put( "TDOMNH4", zw2d(:,:,13) ) 
    415             CALL iom_put( "TDETNH4", zw2d(:,:,14) ) 
    416             CALL iom_put( "TPHYTOT", zw2d(:,:,15) ) 
    417             CALL iom_put( "TZOOTOT", zw2d(:,:,16) ) 
    418             !  
    419             CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) ) 
    420             CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) ) 
    421             CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) ) 
    422             ! 
    423          ELSE 
     401        CALL iom_put( "TNO3PHY", zw2d(:,:,1) ) 
     402        CALL iom_put( "TNH4PHY", zw2d(:,:,2) ) 
     403        CALL iom_put( "TPHYDOM", zw2d(:,:,3) ) 
     404        CALL iom_put( "TPHYNH4", zw2d(:,:,4) ) 
     405        CALL iom_put( "TPHYZOO", zw2d(:,:,5) ) 
     406        CALL iom_put( "TPHYDET", zw2d(:,:,6) ) 
     407        CALL iom_put( "TDETZOO", zw2d(:,:,7) ) 
     408        CALL iom_put( "TZOODET", zw2d(:,:,8) ) 
     409        CALL iom_put( "TZOOBOD", zw2d(:,:,9) ) 
     410        CALL iom_put( "TZOONH4", zw2d(:,:,10) ) 
     411        CALL iom_put( "TZOODOM", zw2d(:,:,11) ) 
     412        CALL iom_put( "TNH4NO3", zw2d(:,:,12) ) 
     413        CALL iom_put( "TDOMNH4", zw2d(:,:,13) ) 
     414        CALL iom_put( "TDETNH4", zw2d(:,:,14) ) 
     415        CALL iom_put( "TPHYTOT", zw2d(:,:,15) ) 
     416        CALL iom_put( "TZOOTOT", zw2d(:,:,16) ) 
     417         !  
     418        CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) ) 
     419        CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) ) 
     420        CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) ) 
     421         ! 
     422       ELSE 
     423          IF( ln_diatrc ) THEN 
    424424            ! 
    425425            trc2d(:,:,jp_pcs0_2d    ) = zw2d(:,:,1)  
    … …  
    457457      IF( l_trdtrc ) THEN 
    458458         DO jl = jp_pcs0_trd, jp_pcs1_trd 
    459             CALL trd_mod_trc( trbio(:,:,:,jl), jl, kt )   ! handle the trend 
     459            CALL trd_trc( trbio(:,:,:,jl), jl, kt )   ! handle the trend 
    460460         END DO 
    461461      ENDIF 
    … …  
    467467      ENDIF 
    468468      ! 
    469       IF( ln_diatrc ) THEN 
     469      IF( ln_diatrc .OR. lk_iomput ) THEN 
    470470         CALL wrk_dealloc( jpi, jpj,     17, zw2d ) 
    471471         CALL wrk_dealloc( jpi, jpj, jpk, 3, zw3d ) 
    … …  
    598598 
    599599   !!====================================================================== 
    600 END MODULE  p2zbio 
     600END MODULE p2zbio 

  • branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zexp.F90

    • Property svn:keywords set to Id
    r3446 r6225  
    2222   USE lbclnk 
    2323   USE prtctl_trc      ! Print control for debbuging 
    24    USE trdmod_oce 
    25    USE trdmod_trc 
     24   USE trd_oce 
     25   USE trdtrc 
    2626   USE iom 
    2727 
    … …  
    4141   REAL(wp)                                ::   areacot   !: surface coastal area 
    4242 
    43    !!* Substitution 
    44 #  include "top_substitute.h90" 
     43   !! * Substitutions 
     44#  include "vectopt_loop_substitute.h90" 
    4545   !!---------------------------------------------------------------------- 
    4646   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    47    !! $Id: trcexp.F90 3294 2012-01-28 16:44:18Z rblod $  
     47   !! $Id$  
    4848   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4949   !!---------------------------------------------------------------------- 
    50  
    5150CONTAINS 
    5251 
    … …  
    9493         DO jj = 2, jpjm1 
    9594            DO ji = fs_2, fs_jpim1 
    96                ze3t = 1. / fse3t(ji,jj,jk) 
     95               ze3t = 1. / e3t_n(ji,jj,jk) 
    9796               tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + ze3t * dmin3(ji,jj,jk) * xksi(ji,jj) 
    9897            END DO 
    … …  
    109108         DO ji = fs_2, fs_jpim1 
    110109            ikt = mbkt(ji,jj)  
    111             tra(ji,jj,ikt,jpno3) = tra(ji,jj,ikt,jpno3) + sedlam * sedpocn(ji,jj) / fse3t(ji,jj,ikt)  
     110            tra(ji,jj,ikt,jpno3) = tra(ji,jj,ikt,jpno3) + sedlam * sedpocn(ji,jj) / e3t_n(ji,jj,ikt)  
    112111            ! Deposition of organic matter in the sediment 
    113112            zwork = vsed * trn(ji,jj,ikt,jpdet) 
    … …  
    120119      DO jj = 2, jpjm1 
    121120         DO ji = fs_2, fs_jpim1 
    122             tra(ji,jj,1,jpno3) = tra(ji,jj,1,jpno3) + zgeolpoc * cmask(ji,jj) / areacot / fse3t(ji,jj,1) 
     121            tra(ji,jj,1,jpno3) = tra(ji,jj,1,jpno3) + zgeolpoc * cmask(ji,jj) / areacot / e3t_n(ji,jj,1) 
    123122         END DO 
    124123      END DO 
    … …  
    127126  
    128127      ! Oa & Ek: diagnostics depending on jpdia2d !          left as example 
    129       IF( ln_diatrc ) THEN 
    130          IF( lk_iomput ) THEN   ;   CALL iom_put( "SEDPOC" , sedpocn ) 
    131          ELSE                   ;   trc2d(:,:,jp_pcs0_2d + 18) = sedpocn(:,:) 
    132          ENDIF 
     128      IF( lk_iomput ) THEN   
     129         CALL iom_put( "SEDPOC" , sedpocn ) 
     130      ELSE 
     131         IF( ln_diatrc )           trc2d(:,:,jp_pcs0_2d + 18) = sedpocn(:,:) 
    133132      ENDIF 
    134133 
    … …  
    164163         ztrbio(:,:,:) = tra(:,:,:,jpno3) - ztrbio(:,:,:) 
    165164         jl = jp_pcs0_trd + 16 
    166          CALL trd_mod_trc( ztrbio, jl, kt )   ! handle the trend 
     165         CALL trd_trc( ztrbio, jl, kt )   ! handle the trend 
    167166         CALL wrk_dealloc( jpi, jpj, jpk, ztrbio )   ! temporary save of trends 
    168167      ENDIF 
    … …  
    211210         DO jj = 1, jpj 
    212211            DO ji = 1, jpi 
    213                zfluo = ( fsdepw(ji,jj,jk  ) / fsdepw(ji,jj,jpkb) )**xhr 
    214                zfluu = ( fsdepw(ji,jj,jk+1) / fsdepw(ji,jj,jpkb) )**xhr 
     212               zfluo = ( gdepw_n(ji,jj,jk  ) / gdepw_n(ji,jj,jpkb) )**xhr 
     213               zfluu = ( gdepw_n(ji,jj,jk+1) / gdepw_n(ji,jj,jpkb) )**xhr 
    215214               IF( zfluo.GT.1. )   zfluo = 1._wp 
    216215               zdm0(ji,jj,jk) = zfluo - zfluu 

  • branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zopt.F90

    • Property svn:keywords set to Id
    r4624 r6225  
    4040   REAL(wp), PUBLIC ::  reddom    ! redfield ratio (C:N) for DOM 
    4141 
    42    !!* Substitution 
    43 #  include "top_substitute.h90" 
    4442   !!---------------------------------------------------------------------- 
    4543   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    46    !! $Id: trcopt.F90 3294 2012-01-28 16:44:18Z rblod $  
     44   !! $Id$  
    4745   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4846   !!---------------------------------------------------------------------- 
    49  
    5047CONTAINS 
    5148 
    … …  
    8986 
    9087      !                                          ! surface irradiance 
    91       zpar0m (:,:)   = qsr   (:,:) * 0.43        ! ------------------ 
     88      !                                          ! ------------------ 
     89      IF( ln_dm2dc ) THEN   ;   zpar0m(:,:) = qsr_mean(:,:) * 0.43 
     90      ELSE                  ;   zpar0m(:,:) = qsr     (:,:) * 0.43 
     91      ENDIF 
    9292      zpar100(:,:)   = zpar0m(:,:) * 0.01 
    9393      zparr  (:,:,1) = zpar0m(:,:) * 0.5 
    … …  
    102102               zkr  = xkr0 + xkrp * EXP( xlr * zpig ) 
    103103               zkg  = xkg0 + xkgp * EXP( xlg * zpig ) 
    104                zparr(ji,jj,jk) = zparr(ji,jj,jk-1) * EXP( -zkr * fse3t(ji,jj,jk-1) ) 
    105                zparg(ji,jj,jk) = zparg(ji,jj,jk-1) * EXP( -zkg * fse3t(ji,jj,jk-1) ) 
     104               zparr(ji,jj,jk) = zparr(ji,jj,jk-1) * EXP( -zkr * e3t_n(ji,jj,jk-1) ) 
     105               zparg(ji,jj,jk) = zparg(ji,jj,jk-1) * EXP( -zkg * e3t_n(ji,jj,jk-1) ) 
    106106            END DO 
    107107        END DO 
    … …  
    113113               zkr  = xkr0 + xkrp * EXP( xlr * zpig ) 
    114114               zkg  = xkg0 + xkgp * EXP( xlg * zpig ) 
    115                zparr(ji,jj,jk) = zparr(ji,jj,jk) / ( zkr * fse3t(ji,jj,jk) ) * ( 1 - EXP( -zkr * fse3t(ji,jj,jk) ) ) 
    116                zparg(ji,jj,jk) = zparg(ji,jj,jk) / ( zkg * fse3t(ji,jj,jk) ) * ( 1 - EXP( -zkg * fse3t(ji,jj,jk) ) ) 
     115               zparr(ji,jj,jk) = zparr(ji,jj,jk) / ( zkr * e3t_n(ji,jj,jk) ) * ( 1 - EXP( -zkr * e3t_n(ji,jj,jk) ) ) 
     116               zparg(ji,jj,jk) = zparg(ji,jj,jk) / ( zkg * e3t_n(ji,jj,jk) ) * ( 1 - EXP( -zkg * e3t_n(ji,jj,jk) ) ) 
    117117               etot (ji,jj,jk) = MAX( zparr(ji,jj,jk) + zparg(ji,jj,jk), 1.e-15 ) 
    118118            END DO 
    … …  
    128128              IF( etot(ji,jj,jk) >= zpar100(ji,jj) )   neln(ji,jj) = jk + 1  
    129129              !                                       ! nb. this is to ensure compatibility with 
    130               !                                       ! nmld_trc definition in trd_mld_trc_zint 
     130              !                                       ! nmld_trc definition in trd_mxl_trc_zint 
    131131           END DO 
    132132         END DO 
    … …  
    135135      DO jj = 1, jpj 
    136136         DO ji = 1, jpi 
    137             heup(ji,jj) = fsdepw(ji,jj,neln(ji,jj)) 
     137            heup(ji,jj) = gdepw_n(ji,jj,neln(ji,jj)) 
    138138         END DO 
    139139      END DO  

  • branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90

    • Property svn:keywords set to Id
    r4624 r6225  
    1818   USE sms_pisces 
    1919   USE lbclnk 
    20    USE trdmod_oce 
    21    USE trdmod_trc 
     20   USE trd_oce 
     21   USE trdtrc 
    2222   USE iom 
    2323   USE prtctl_trc      ! Print control for debbuging 
    … …  
    3434   REAL(wp), PUBLIC ::   xhr         ! coeff for martin''s remineralisation profile 
    3535 
    36    !!* Substitution 
    37 #  include "top_substitute.h90" 
    3836   !!---------------------------------------------------------------------- 
    3937   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    40    !! $Id: p2z_sed.F90 3294 2012-01-28 16:44:18Z rblod $  
     38   !! $Id$  
    4139   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4240   !!---------------------------------------------------------------------- 
    … …  
    102100         DO jj = 1, jpj 
    103101            DO ji = 1, jpi 
    104                ztra(ji,jj,jk)  = - ( zwork(ji,jj,jk) - zwork(ji,jj,jk+1) ) / fse3t(ji,jj,jk) 
     102               ztra(ji,jj,jk)  = - ( zwork(ji,jj,jk) - zwork(ji,jj,jk+1) ) / e3t_n(ji,jj,jk) 
    105103               tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + ztra(ji,jj,jk)  
    106104            END DO 
    … …  
    108106      END DO 
    109107 
    110       IF( ln_diatrc ) THEN  
    111          CALL wrk_alloc( jpi, jpj, zw2d ) 
    112          zw2d(:,:) =  ztra(:,:,1) * fse3t(:,:,1) * 86400. 
    113          DO jk = 2, jpkm1 
    114             zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400. 
    115          END DO 
    116          IF( lk_iomput )  THEN 
    117            CALL iom_put( "TDETSED", zw2d ) 
    118          ELSE 
    119            trc2d(:,:,jp_pcs0_2d + 7) = zw2d(:,:) 
     108      IF( lk_iomput )  THEN 
     109         IF( iom_use( "TDETSED" ) ) THEN 
     110            CALL wrk_alloc( jpi, jpj, zw2d ) 
     111            zw2d(:,:) =  ztra(:,:,1) * e3t_n(:,:,1) * 86400._wp 
     112            DO jk = 2, jpkm1 
     113               zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * e3t_n(:,:,jk) * 86400._wp 
     114            END DO 
     115            CALL iom_put( "TDETSED", zw2d ) 
     116            CALL wrk_dealloc( jpi, jpj, zw2d ) 
    120117         ENDIF 
    121          CALL wrk_dealloc( jpi, jpj, zw2d ) 
     118      ELSE 
     119         IF( ln_diatrc ) THEN  
     120            CALL wrk_alloc( jpi, jpj, zw2d ) 
     121            zw2d(:,:) =  ztra(:,:,1) * e3t_n(:,:,1) * 86400._wp 
     122            DO jk = 2, jpkm1 
     123               zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * e3t_n(:,:,jk) * 86400._wp 
     124            END DO 
     125            trc2d(:,:,jp_pcs0_2d + 7) = zw2d(:,:) 
     126            CALL wrk_dealloc( jpi, jpj, zw2d ) 
     127         ENDIF 
    122128      ENDIF 
    123129      ! 
    … …  
    128134         ztrbio(:,:,:) = tra(:,:,:,jpdet) - ztrbio(:,:,:) 
    129135         jl = jp_pcs0_trd + 7 
    130          CALL trd_mod_trc( ztrbio, jl, kt )   ! handle the trend 
     136         CALL trd_trc( ztrbio, jl, kt )   ! handle the trend 
    131137         CALL wrk_dealloc( jpi, jpj, jpk, ztrbio ) 
    132138      ENDIF 

  • branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsms.F90

    • Property svn:keywords set to Id
    r4624 r6225  
    2020   USE p2zsed 
    2121   USE p2zexp 
    22    USE trdmod_oce 
    23    USE trdmod_trc_oce 
    24    USE trdmod_trc 
    25    USE trdmld_trc 
     22   USE trd_oce 
     23   USE trdtrc_oce 
     24   USE trdtrc 
     25   USE trdmxl_trc 
    2626 
    2727   IMPLICIT NONE 
    … …  
    3232   !!---------------------------------------------------------------------- 
    3333   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    34    !! $Id: p2zsms.F90 3294 2012-01-28 16:44:18Z rblod $  
     34   !! $Id$  
    3535   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
    … …  
    6161      IF( l_trdtrc ) THEN 
    6262         DO jn = jp_pcs0, jp_pcs1 
    63            CALL trd_mod_trc( tra(:,:,:,jn), jn, jptra_trd_sms, kt )   ! save trends 
     63           CALL trd_trc( tra(:,:,:,jn), jn, jptra_sms, kt )   ! save trends 
    6464         END DO 
    6565      END IF 
    6666 
    67       IF( lk_trdmld_trc )  CALL trd_mld_bio( kt )   ! trends: Mixed-layer 
     67      IF( lk_trdmxl_trc )  CALL trd_mxl_bio( kt )   ! trends: Mixed-layer 
    6868      ! 
    6969      IF ( lwm .AND. kt == nittrc000 ) CALL FLUSH    ( numonp )     ! flush output namelist PISCES 
    … …  
    8484 
    8585   !!====================================================================== 
    86 END MODULE  p2zsms 
     86END MODULE p2zsms 
Note: See TracChangeset for help on using the changeset viewer.