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 261 for trunk/NEMO/TOP_SRC – NEMO

Changeset 261 for trunk/NEMO/TOP_SRC


Ignore:
Timestamp:
2005-09-08T15:18:32+02:00 (19 years ago)
Author:
opalod
Message:

nemo_v1_update_05:RB+OA:Update and rewritting of (part of) the TOP component

Location:
trunk/NEMO/TOP_SRC/TRP
Files:
27 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/TOP_SRC/TRP/trcadv_cen2.F90

    r247 r261  
    2424#  include "passivetrc_substitute.h90" 
    2525   !!---------------------------------------------------------------------- 
    26    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    27    !! $Header$  
    28    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     26   !!   OPA 9.0 , LODYC-IPSL (2003) 
    2927   !!---------------------------------------------------------------------- 
    3028 
    … …  
    257255         !                                                ! =============== 
    258256 
    259          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    260             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     257         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     258            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    261259            WRITE(numout,*) ' trc/had - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)  
    262260            tra_ctl(jn) = ztra  
    … …  
    320318         END DO 
    321319 
    322          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    323             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     320         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     321            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    324322            WRITE(numout,*) ' trc/zad - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' centered2'  
    325323            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcadv_muscl.F90

    r247 r261  
    1414   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    1515   USE trcbbl          ! advective passive tracers in the BBL 
     16   USE lib_mpp 
    1617 
    1718   IMPLICIT NONE 
    … …  
    2122   PUBLIC trc_adv_muscl  ! routine called by trcstp.F90 
    2223 
    23    !! * Module variable 
    24    REAL(wp), DIMENSION(jpk) ::   & 
    25       rdttrc                     ! vertical profile of tracer time-step 
    26  
    2724   !! * Substitutions 
    2825#  include "passivetrc_substitute.h90" 
    2926   !!---------------------------------------------------------------------- 
    30    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    31    !! $Header$  
    32    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     27   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3328   !!---------------------------------------------------------------------- 
    3429 
    … …  
    7772      REAL(wp) ::   zu, zv, zw, zeu, zev, zew, zbtr, ztra 
    7873      REAL(wp) ::   z0u, z0v, z0w 
    79       REAL(wp) ::   zzt1, zzt2, zalpha 
     74      REAL(wp) ::   zzt1, zzt2, zalpha, z2dtt 
    8075#if defined key_trc_diatrd 
    8176      REAL(wp) ::   ztai, ztaj 
    … …  
    9085         WRITE(numout,*) 'trc_adv : MUSCL advection scheme' 
    9186         WRITE(numout,*) '~~~~~~~' 
    92          rdttrc(:) =  rdttra(:) * FLOAT(ndttrc) 
    9387      ENDIF 
    9488 
    … …  
    145139         ztp2(:,:,jpk) = 0.e0 
    146140 
    147  
    148141         ! Slopes limitation 
    149142         DO jk = 1, jpkm1 
    150143            DO jj = 2, jpj 
    151                DO ji = fs_2, fs_jpim1   ! vector opt. 
     144               DO ji = fs_2, jpi   ! vector opt. 
    152145                  ztp1(ji,jj,jk) = SIGN( 1., ztp1(ji,jj,jk) )   & 
    153146                     &           * MIN(    ABS( ztp1(ji  ,jj,jk) ),   & 
    … …  
    163156            END DO 
    164157         END DO 
    165  
    166158 
    167159         ! Advection terms 
    … …  
    179171#endif 
    180172                  ! MUSCL fluxes 
     173                  z2dtt = rdttra(jk) * FLOAT(ndttrc) 
    181174                  z0u = SIGN( 0.5, zun(ji,jj,jk) )             
    182175                  zalpha = 0.5 - z0u 
    183                   zu  = z0u - 0.5 * zun(ji,jj,jk) * rdttrc(jk) / e1u(ji,jj) 
     176                  zu  = z0u - 0.5 * zun(ji,jj,jk) * z2dtt / e1u(ji,jj) 
    184177                  zzt1 = trb(ji+1,jj,jk,jn) + zu*ztp1(ji+1,jj,jk) 
    185178                  zzt2 = trb(ji  ,jj,jk,jn) + zu*ztp1(ji  ,jj,jk) 
    186179                  zt1(ji,jj,jk) = zeu * ( zalpha * zzt1 + (1.-zalpha) * zzt2 ) 
    187  
    188180                  z0v = SIGN( 0.5, zvn(ji,jj,jk) )             
    189181                  zalpha = 0.5 - z0v 
    190                   zv  = z0v - 0.5 * zvn(ji,jj,jk) * rdttrc(jk) / e2v(ji,jj) 
     182                  zv  = z0v - 0.5 * zvn(ji,jj,jk) * z2dtt / e2v(ji,jj) 
    191183                  zzt1 = trb(ji,jj+1,jk,jn) + zv*ztp2(ji,jj+1,jk) 
    192184                  zzt2 = trb(ji,jj  ,jk,jn) + zv*ztp2(ji,jj  ,jk) 
    … …  
    238230         END DO 
    239231 
    240          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    241             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     232         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     233            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    242234            WRITE(numout,*) ' trc/had  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl'  
    243235            tra_ctl(jn) = ztra  
    … …  
    284276            DO jj = 2, jpjm1       
    285277               DO ji = fs_2, fs_jpim1   ! vector opt. 
     278                  z2dtt = rdttra(jk) * FLOAT(ndttrc) 
    286279                  zew = zwn(ji,jj,jk+1) 
    287280                  z0w = SIGN( 0.5, zwn(ji,jj,jk+1) ) 
    288281                  zalpha = 0.5 + z0w 
    289                   zw  = z0w - 0.5 * zwn(ji,jj,jk+1)*rdttrc(jk) / fse3w(ji,jj,jk+1) 
     282                  zw  = z0w - 0.5 * zwn(ji,jj,jk+1)*z2dtt / fse3w(ji,jj,jk+1) 
    290283                  zzt1 = trb(ji,jj,jk+1,jn) + zw*ztp1(ji,jj,jk+1) 
    291284                  zzt2 = trb(ji,jj,jk  ,jn) + zw*ztp1(ji,jj,jk  ) 
    … …  
    322315         END DO 
    323316 
    324          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    325             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     317         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     318            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    326319            WRITE(numout,*) ' trc/zad  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl'  
    327320            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcadv_muscl2.F90

    r247 r261  
    2828#  include "passivetrc_substitute.h90" 
    2929   !!---------------------------------------------------------------------- 
    30    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    31    !! $Header$  
    32    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     30   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3331   !!---------------------------------------------------------------------- 
    3432 
    … …  
    280278         END DO 
    281279 
    282          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    283             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     280         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     281            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    284282            WRITE(numout,*) ' trc/had  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl'  
    285283            tra_ctl(jn) = ztra  
    … …  
    379377         END DO 
    380378 
    381          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    382             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     379         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     380            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    383381            WRITE(numout,*) ' trc/zad  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl2'  
    384382            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcadv_smolar.F90

    r247 r261  
    2929#  include "passivetrc_substitute.h90" 
    3030   !!---------------------------------------------------------------------- 
    31    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    32    !! $Header$  
    33    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     31   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3432   !!---------------------------------------------------------------------- 
    3533CONTAINS 
    … …  
    501499#endif 
    502500 
    503         IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    504            ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     501        IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     502           ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    505503           WRITE(numout,*) ' trc/zad  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' smolar'  
    506504           tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcadv_tvd.F90

    r247 r261  
    2626#  include "passivetrc_substitute.h90" 
    2727   !!---------------------------------------------------------------------- 
    28    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    29    !! $Header$  
    30    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     28   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3129   !!---------------------------------------------------------------------- 
    3230 
    … …  
    234232         END DO 
    235233 
    236          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    237             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     234         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     235            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    238236            WRITE(numout,*) ' trc/zad  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' tvd'  
    239237            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcbbc.F90

    r247 r261  
    3939#  include "passivetrc_substitute.h90" 
    4040   !!---------------------------------------------------------------------- 
    41    !!  TOP 1.0 , LOCEAN-IPSL (2005)  
    42    !! $Header$  
    43    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     41   !!  OPA 9.0 , LODYC-IPSL (2003) 
    4442   !!---------------------------------------------------------------------- 
    4543 
    … …  
    106104            END DO 
    107105#endif 
    108             IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    109                ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     106            IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     107               ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    110108               WRITE(numout,*) ' trc/bbc  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)  
    111109               tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcbbl.F90

    r247 r261  
    4848#  include "passivetrc_substitute.h90" 
    4949   !!---------------------------------------------------------------------- 
    50    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    51    !! $Header$  
    52    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     50   !!   OPA 9.0 , LODYC-IPSL (2003) 
    5351   !!---------------------------------------------------------------------- 
    5452 
    … …  
    341339#  endif 
    342340         END DO 
    343          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    344             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     341         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     342            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    345343            WRITE(numout,*) ' trc/bbl  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)  
    346344            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcbbl_adv.h90

    r247 r261  
    44 
    55   !!---------------------------------------------------------------------- 
    6    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    7    !! $Header$  
    8    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     6   !!   OPA 9.0 , LODYC-IPSL  (2003) 
    97   !!---------------------------------------------------------------------- 
    108 
    … …  
    384382#endif 
    385383          END DO 
    386           IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    387              ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     384          IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     385             ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    388386             WRITE(numout,*) ' trc/bbl  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)  
    389387             tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcdmp.F90

    r247 r261  
    1818   USE trc             ! ocean passive tracers variables 
    1919   USE trctrp_lec      ! passive tracers transport 
     20   USE trcdta 
    2021 
    2122   IMPLICIT NONE 
    … …  
    3435#  include "passivetrc_substitute.h90" 
    3536   !!---------------------------------------------------------------------- 
    36    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    37    !! $Header$  
    38    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     37   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3938   !!---------------------------------------------------------------------- 
    4039 
    … …  
    8685 
    8786      ! Initialize the input fields for newtonian damping 
    88       CALL trcdta( kt ) 
     87      CALL dta_trc( kt ) 
    8988 
    9089      DO jn = 1, jptra 
    … …  
    156155         ENDIF 
    157156 
    158          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    159             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     157         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     158            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    160159            WRITE(numout,*) ' trc/dmp  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)  
    161160            tra_ctl(jn) = ztra  
    … …  
    247246      !!---------------------------------------------------------------------- 
    248247      !! * Local declarations 
    249       INTEGER ::   ji, jj, jk, jn      ! dummy loop indices 
     248      INTEGER ::   ji, jj, jk, jn       ! dummy loop indices 
    250249      REAL(wp) ::   & 
    251250         zlat, zlat0, zlat1, zlat2     ! temporary scalar 
    … …  
    266265      ! damping along the forced closed boundary over 6 grid-points 
    267266      DO jn = 1, 6 
    268          IF( lzoom_w )   resto( mi0(jn+jpizoom):mi1(jn+jpizoom), : , : ) = zfact(jn) ! west  closed 
    269          IF( lzoom_s )   resto( : , mj0(jn+jpjzoom):mj1(jn+jpjzoom), : ) = zfact(jn) ! south closed  
    270          IF( lzoom_e )   resto( mi0(jpiglo+jpizoom-1-jn):mi1(jpiglo+jpizoom-1-jn) , : , : ) & 
     267         IF( lzoom_w )   restotr( mi0(jn+jpizoom):mi1(jn+jpizoom), : , : , : ) = zfact(jn) ! west  closed 
     268         IF( lzoom_s )   restotr( : , mj0(jn+jpjzoom):mj1(jn+jpjzoom), : , : ) = zfact(jn) ! south closed  
     269         IF( lzoom_e )   restotr( mi0(jpiglo+jpizoom-1-jn):mi1(jpiglo+jpizoom-1-jn) , : , : , : ) & 
    271270                       &              = zfact(jn)                                 ! east  closed  
    272          IF( lzoom_n )   resto( : , mj0(jpjglo+jpjzoom-1-jn):mj1(jpjglo+jpjzoom-1-jn) , : ) & 
     271         IF( lzoom_n )   restotr( : , mj0(jpjglo+jpjzoom-1-jn):mj1(jpjglo+jpjzoom-1-jn) , : , : ) & 
    273272                       &              = zfact(jn)                                 ! north closed 
    274273      END DO 

  • trunk/NEMO/TOP_SRC/TRP/trcldf_bilap.F90

    r247 r261  
    2323#  include "passivetrc_substitute.h90" 
    2424   !!---------------------------------------------------------------------- 
    25    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    26    !! $Header$  
    27    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     25   !!   OPA 9.0 , LODYC-IPSL (2003) 
    2826   !!---------------------------------------------------------------------- 
    2927 
    … …  
    205203         CALL lbc_lnk( trtrd(1,1,1,jn,5), 'T', 1. )   
    206204#endif 
    207          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    208             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     205         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     206            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    209207            WRITE(numout,*) ' trc/ldf  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    210208            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcldf_bilapg.F90

    r247 r261  
    2626#  include "passivetrc_substitute.h90" 
    2727   !!---------------------------------------------------------------------- 
    28    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    29    !! $Header$  
    30    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     28   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3129   !!---------------------------------------------------------------------- 
    3230    
    … …  
    118116         !                                                ! =============== 
    119117 
    120          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    121             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     118         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     119            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    122120            WRITE(numout,*) ' trc/ldf  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    123121            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcldf_iso.F90

    r247 r261  
    2525#  include "passivetrc_substitute.h90" 
    2626   !!---------------------------------------------------------------------- 
    27    !!---------------------------------------------------------------------- 
    28    !!  TOP 1.0 , LOCEAN-IPSL (2005) 
    29    !! $Header$ 
    30    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    31    !!---------------------------------------------------------------------- 
     27 
    3228CONTAINS 
    3329 
    … …  
    203199         END DO                                        !   End of slab   
    204200         !                                             ! =============== 
    205          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    206             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     201         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     202            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    207203            WRITE(numout,*) ' trc/ldf  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    208204            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90

    r247 r261  
    2424#  include "passivetrc_substitute.h90" 
    2525   !!---------------------------------------------------------------------- 
    26    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    27    !! $Header$  
    28    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     26   !!   OPA 9.0 , LODYC-IPSL (2003) 
    2927   !!---------------------------------------------------------------------- 
    3028 
    … …  
    240238         END DO                                        !   End of slab   
    241239         !                                             ! =============== 
    242          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    243             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     240         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     241            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    244242            WRITE(numout,*) ' trc/ldf  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    245243            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcldf_lap.F90

    r247 r261  
    2222#  include "passivetrc_substitute.h90" 
    2323   !!---------------------------------------------------------------------- 
    24    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    25    !! $Header$  
    26    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     24   !!   OPA 9.0 , LODYC-IPSL (2003) 
    2725   !!---------------------------------------------------------------------- 
    2826    
    … …  
    139137         !                                                  ! ============= 
    140138 
    141          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    142             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     139         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     140            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    143141            WRITE(numout,*) ' trc/ldf  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    144142            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcnxt.F90

    r247 r261  
    2020   PUBLIC trc_nxt          ! routine called by step.F90 
    2121   !!---------------------------------------------------------------------- 
    22    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    23    !! $Header$  
    24    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     22   !!   OPA 9.0 , LODYC-IPSL   (2003) 
    2523   !!---------------------------------------------------------------------- 
    2624 
    … …  
    6967      !!---------------------------------------------------------------------- 
    7068 
    71       IF( kt == nittrc000 ) THEN 
     69      IF( kt == nittrc000 .AND. lwp ) THEN 
    7270         WRITE(numout,*) 
    7371         WRITE(numout,*) 'trc_nxt : time stepping on passive tracers' 
    … …  
    139137         END DO                                           !   End of slab 
    140138         !                                                ! =============== 
    141          IF(l_ctl) THEN         ! print mean field (used for debugging) 
    142             ztra = SUM( trn(2:nictl,2:njctl,1:jpkm1,jn)*tmask(2:nictl,2:njctl,1:jpkm1) )  
     139         IF(ln_ctl) THEN         ! print mean field (used for debugging) 
     140            ztra = SUM( trn(2:nictle,2:njctle,1:jpkm1,jn)*tmask(2:nictle,2:njctle,1:jpkm1) )  
    143141            WRITE(numout,*) ' trc/nxt  - ',ctrcnm(jn),' : ', ztra 
    144142         ENDIF 

  • trunk/NEMO/TOP_SRC/TRP/trcrad.F90

    r247 r261  
    1111   USE oce_trc             ! ocean dynamics and tracers variables 
    1212   USE trc                 ! ocean passive tracers variables 
     13   USE lib_mpp 
    1314 
    1415   IMPLICIT NONE 
    … …  
    2021#  include "passivetrc_substitute.h90" 
    2122   !!---------------------------------------------------------------------- 
    22    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    23    !! $Header$  
    24    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     23   !!   OPA 9.0 , LODYC-IPSL   (2003) 
    2524   !!---------------------------------------------------------------------- 
    2625CONTAINS 
    … …  
    8483            END DO 
    8584         END DO 
     85 
     86         IF( lk_mpp ) THEN 
     87           CALL mpp_sum( trcorb )   ! sum over the global domain 
     88           CALL mpp_sum( trcorn )   ! sum over the global domain 
     89           CALL mpp_sum( trmasb )   ! sum over the global domain 
     90           CALL mpp_sum( trmasn )   ! sum over the global domain 
     91         ENDIF 
     92 
    8693         DO jk = 1, jpkm1 
    8794            DO jj = 1, jpj 
    … …  
    135142         END DO 
    136143 
     144         IF( lk_mpp ) THEN 
     145           CALL mpp_sum( trcorb )   ! sum over the global domain 
     146           CALL mpp_sum( trcorn )   ! sum over the global domain 
     147           CALL mpp_sum( trmasb )   ! sum over the global domain 
     148           CALL mpp_sum( trmasn )   ! sum over the global domain 
     149         ENDIF 
     150 
    137151         IF( trcorb /= 0) THEN 
    138152            DO jk = 1, jpkm1 
    … …  
    161175#endif 
    162176      DO jn = 1, jptra 
    163          IF(l_ctl) THEN         ! print mean field (used for debugging) 
    164             ztra = SUM( trn(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) )  
     177         IF(ln_ctl) THEN         ! print mean field (used for debugging) 
     178            ztra = SUM( trn(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )  
    165179            WRITE(numout,*) ' trc/rad  - ',ctrcnm(jn),' : ', ztra 
    166180         ENDIF 

  • trunk/NEMO/TOP_SRC/TRP/trcsbc.F90

    r247 r261  
    2121#  include "passivetrc_substitute.h90" 
    2222   !!---------------------------------------------------------------------- 
    23    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    24    !! $Header$  
    25    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     23   !!   OPA 9.0 , LODYC-IPSL (2003) 
    2624   !!---------------------------------------------------------------------- 
    2725 
    … …  
    8987         END DO 
    9088          
    91          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    92             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     89         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     90            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    9391            WRITE(numout,*) ' trc/sbc  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    9492            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trcstp.F90

    r247 r261  
    1212   USE trc              ! ocean passive tracers variables  
    1313   USE trctrp           ! passive tracers transport 
     14   USE trcsms           ! passive tracers sources and sinks 
    1415 
    1516   IMPLICIT NONE 
    … …  
    1920   PUBLIC trc_stp           ! called by step 
    2021   !!---------------------------------------------------------------------- 
    21    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    22    !! $Header$  
    23    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     22   !!   OPA 9.0 , LODYC-IPSL (2003) 
    2423   !!---------------------------------------------------------------------- 
    2524 
    … …  
    4140      !! * Arguments 
    4241      INTEGER, INTENT( in ) ::  kt  ! ocean time-step index 
     42      INTEGER :: jn 
     43      REAL(wp) :: ztra 
    4344 
    4445 
    … …  
    4950 
    5051 
    51       CALL trcsms( kt ) 
     52      CALL trc_sms( kt ) 
     53 
     54       DO jn = 1, jptra 
     55         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     56            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
     57            WRITE(numout,*) ' trc/sms  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
     58            tra_ctl(jn) = ztra 
     59         ENDIF 
     60      END DO 
     61 
    5262  
    5363 
    … …  
    5565 
    5666      CALL trc_trp( kt ) 
     67 
     68       DO jn = 1, jptra 
     69         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     70            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
     71            WRITE(numout,*) ' trc/trp  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
     72            tra_ctl(jn) = ztra 
     73         ENDIF 
     74      END DO 
     75 
    5776 
    5877 

  • trunk/NEMO/TOP_SRC/TRP/trctrp.F90

    r247 r261  
    2424  
    2525   USE trcnxt          ! time-stepping                       (trc_nxt routine) 
    26    USE trcrad          ! positivity                          (trc_rad routine) 
     26!   USE trcrad          ! positivity                          (trc_rad routine) 
    2727 
    2828   USE trcadv_cen2     ! 2nd order centered advection   (trc_adv_cen2 routine) 
    … …  
    5050#  include "domzgr_substitute.h90" 
    5151   !!---------------------------------------------------------------------- 
    52    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    53    !! $Header$  
    54    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     52   !!   OPA 9.0 , LODYC-IPSL (2003) 
    5553   !!---------------------------------------------------------------------- 
    5654 
    … …  
    9391 
    9492      IF( lk_trcdmp        )   CALL trc_dmp( kt )            ! internal damping trends 
    95   
     93 
    9694      !                                                      ! horizontal & vertical advection 
    9795      IF( ln_trcadv_cen2   )   CALL trc_adv_cen2  ( kt )             ! 2nd order centered scheme 
    … …  
    10098      IF( ln_trcadv_tvd    )   CALL trc_adv_tvd   ( kt )             ! TVD scheme 
    10199      IF( ln_trcadv_smolar )   CALL trc_adv_smolar( kt )             ! SMOLARKIEWICZ scheme 
     100 
    102101  
    103102      IF( n_cla == 1   )   THEN  
    … …  
    124123                               CALL trc_nxt( kt )            ! tracer fields at next time step 
    125124  
    126                                CALL trc_rad( kt )            ! Correct artificial negative concentrations  
     125       !!                       CALL trc_rad( kt )            ! Correct artificial negative concentrations  
    127126      !                                                      ! for isopycnal scheme 
    128127 
    129       IF( lk_zps           )   CALL zps_hde_trc( kt, trb, gtru, gtrv )  ! Partial steps: now horizontal gradient 
     128      ! IF( lk_zps           )   CALL zps_hde_trc( kt, trb, gtru, gtrv )  ! Partial steps: now horizontal gradient 
    130129      !                                                                 ! of passive tracers at the bottom ocean level 
    131130 

  • trunk/NEMO/TOP_SRC/TRP/trctrp_ctl.F90

    r247 r261  
    2727 
    2828   !!---------------------------------------------------------------------- 
    29    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    30    !! $Header$  
    31    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     29   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3230   !!---------------------------------------------------------------------- 
    3331 

  • trunk/NEMO/TOP_SRC/TRP/trctrp_lec.F90

    r247 r261  
    7070#endif 
    7171   !!---------------------------------------------------------------------- 
    72    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    73    !! $Header$  
    74    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     72   !!   OPA 9.0 , LODYC-IPSL (2003) 
    7573   !!---------------------------------------------------------------------- 
    7674 

  • trunk/NEMO/TOP_SRC/TRP/trczdf_exp.F90

    r247 r261  
    2828#  include "passivetrc_substitute.h90" 
    2929   !!---------------------------------------------------------------------- 
    30    !!   TOP 1.0 , LOCEAN-IPSL (2005) 
    31    !! $Header$ 
    32    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
     30   !!  OPA 9.0, LODYC-IPSL (2003) 
    3331   !!---------------------------------------------------------------------- 
    3432 

  • trunk/NEMO/TOP_SRC/TRP/trczdf_imp.F90

    r247 r261  
    2727#  include "passivetrc_substitute.h90" 
    2828   !!---------------------------------------------------------------------- 
    29    !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    30    !! $Header$  
    31    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     29   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3230   !!---------------------------------------------------------------------- 
    3331 
    … …  
    183181         END DO 
    184182 
    185          IF(l_ctl) THEN  
     183         IF(ln_ctl) THEN  
    186184            DO jk = 1, jpkm1 
    187185               DO jj = 2, jpjm1 
    … …  
    193191         ENDIF 
    194192       
    195          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    196             ztra = SUM( ztrd(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     193         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     194            ztra = SUM( ztrd(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    197195            WRITE(numout,*) ' trc/zdf  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    198196            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trczdf_imp.F90_vopt

    r247 r261  
    2727#  include "passivetrc_substitute.h90" 
    2828   !!---------------------------------------------------------------------- 
    29    !!  TOP 1.0 , LOCEAN-IPSL (2005) 
    30    !! $Header$ 
    31    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
     29   !!   OPA 9.0 , LODYC-IPSL (2003) 
    3230   !!---------------------------------------------------------------------- 
    3331 
    … …  
    183181         END DO 
    184182 
    185          IF( l_ctl .AND. lwp ) THEN  
     183         IF( ln_ctl .AND. lwp ) THEN  
    186184            DO jk = 1, jpkm1 
    187185               DO jj = 2, jpjm1 
    … …  
    193191         ENDIF 
    194192       
    195          IF( l_ctl .AND. lwp ) THEN         ! print mean trends (used for debugging) 
     193         IF( ln_ctl .AND. lwp ) THEN         ! print mean trends (used for debugging) 
    196194            ztra = SUM( ztrd(2:jpim1,2:jpjm1,1:jpkm1,jn) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    197195            WRITE(numout,*) ' trc/zdf  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 

  • trunk/NEMO/TOP_SRC/TRP/trczdf_iso.F90

    r247 r261  
    3131#  include "passivetrc_substitute.h90" 
    3232   !!---------------------------------------------------------------------- 
    33    !!---------------------------------------------------------------------- 
    34    !!  TOP 1.0 , LOCEAN-IPSL (2005) 
    35    !! $Header$ 
    36    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    37    !!---------------------------------------------------------------------- 
     33 
    3834CONTAINS 
    3935 
    … …  
    114110      !!   9.0  !  04-03  (C. Ethe )  adapted for passive tracers 
    115111      !!--------------------------------------------------------------------- 
     112      !! * Modules used 
     113      USE oce_trc               ,   & 
     114         zavs => va 
     115 
    116116      !! * Arguments 
    117117      INTEGER, INTENT( in ) ::   kt           ! ocean time-step index 
    … …  
    134134         zwx, zwy, zwz, zwt        ! ??? 
    135135      REAL(wp), DIMENSION(jpi,jpk) ::   & 
    136          ztfw, zdit, zdjt, zdj1t,   & 
    137          zavs 
     136         ztfw, zdit, zdjt, zdj1t 
    138137#if defined key_trcldf_eiv   ||   defined key_esopa 
    139138      REAL(wp), DIMENSION(jpi,jpk) ::   & 
    … …  
    145144         zuwki, zvwki 
    146145#endif 
     146      !!--------------------------------------------------------------------- 
     147      !!  OPA 8.5, LODYC-IPSL (2002) 
     148      !!--------------------------------------------------------------------- 
    147149 
    148150      IF( kt == nittrc000 ) THEN 
    … …  
    257259 
    258260                  ! save avs in zavs to recover avs in output files 
    259                   zavs(ji,jk) = fstravs(ji,jj,jk) 
     261                  zavs(ji,jj,jk) = fstravs(ji,jj,jk) 
    260262                  ! add isopycnal vertical coeff. to avs 
    261263                  fstravs(ji,jj,jk) = fstravs(ji,jj,jk) + zavi 
    … …  
    343345         !                                                ! =============== 
    344346 
    345          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    346             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     347         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     348            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    347349            WRITE(numout,*) ' trc/zdf 1  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    348350            tra_ctl(jn) = ztra  
    … …  
    425427            DO jk = 2, jpkm1 
    426428               DO ji = 2, jpim1 
    427                   fstravs(ji,jj,jk) = zavs(ji,jk) 
     429                  fstravs(ji,jj,jk) = zavs(ji,jj,jk) 
    428430               END DO 
    429431            END DO 
    … …  
    434436 
    435437 
    436          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    437             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     438         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     439            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    438440            WRITE(numout,*) ' trc/zdf 2  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    439441            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/trczdf_iso_vopt.F90

    r247 r261  
    3434#  include "passivetrc_substitute.h90" 
    3535   !!---------------------------------------------------------------------- 
    36    !!---------------------------------------------------------------------- 
    37    !!  TOP 1.0 , LOCEAN-IPSL (2005) 
    38    !! $Header$ 
    39    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    40    !!---------------------------------------------------------------------- 
     36 
    4137CONTAINS 
    4238    
    … …  
    5551      !! * Arguments 
    5652      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
     53      !!--------------------------------------------------------------------- 
     54      !!  OPA 8.5, LODYC-IPSL (2002) 
     55      !!--------------------------------------------------------------------- 
    5756 
    5857      IF( kt == nittrc000 ) THEN 
    … …  
    151150         zwi, zwt, zavsi                     ! temporary workspace arrays 
    152151      REAL(wp) ::    ztra              !temporary scalars 
     152      !!--------------------------------------------------------------------- 
     153      !!  OPA 8.5, LODYC-IPSL (2002) 
     154      !!--------------------------------------------------------------------- 
    153155 
    154156 
    … …  
    312314#  endif 
    313315#endif 
    314          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    315             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     316         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     317            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    316318            WRITE(numout,*) ' trc/zdf 1 - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    317319            tra_ctl(jn) = ztra  
    … …  
    400402         zwz, zwt, ztfw             ! temporary workspace arrays 
    401403      REAL(wp) ::    ztra              !temporary scalars 
    402  
     404      !!--------------------------------------------------------------------- 
     405      !!  OPA 8.5, LODYC-IPSL (2002) 
     406      !!--------------------------------------------------------------------- 
    403407 
    404408      DO jn = 1, jptra 
    … …  
    544548         END DO 
    545549 
    546          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    547             ztra = SUM( tra(2:nictl,2:njctl,1:jpkm1,jn) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     550         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     551            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 
    548552            WRITE(numout,*) ' trc/zdf 2 - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 
    549553            tra_ctl(jn) = ztra  

  • trunk/NEMO/TOP_SRC/TRP/zpshde_trc.F90

    r247 r261  
    2929#  include "passivetrc_substitute.h90" 
    3030   !!---------------------------------------------------------------------- 
    31    !!---------------------------------------------------------------------- 
    32    !!  TOP 1.0 , LOCEAN-IPSL (2005) 
    33    !! $Header$ 
    34    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    35    !!---------------------------------------------------------------------- 
     31 
    3632CONTAINS 
    3733 
    … …  
    9995         zmaxu1, zmaxu2,         &  !    "         " 
    10096         zmaxv1, zmaxv2             !    "         " 
     97      !!---------------------------------------------------------------------- 
     98      !!  OPA 8.5, LODYC-IPSL (2002) 
     99      !!---------------------------------------------------------------------- 
    101100 
    102101      ! Initialization (first time-step only): compute mbatu and mbatv 
    … …  
    187186CONTAINS 
    188187   SUBROUTINE zps_hde_trc ( kt, ptra, pgtru, pgtrv ) ! Empty routine 
    189       REAL(wp) :: kt,ptra,pgtru, pgtrv 
     188      INTEGER, INTENT( in) :: kt 
     189      REAL(wp), DIMENSION(:,:,:,:) :: ptra 
     190      REAL(wp) :: pgtru, pgtrv 
    190191      WRITE(*,*) 'zps_hde_trc: You should not have seen this print! error?',   & 
    191192         kt, ptra, pgtru, pgtrv 
Note: See TracChangeset for help on using the changeset viewer.