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

Changeset 132 for trunk/NEMO


Ignore:
Timestamp:
2004-07-07T14:51:53+02:00 (20 years ago)
Author:
opalod
Message:

CT : UPDATE082 : Finalization of the poleward transport diagnostics

Location:
trunk/NEMO/OPA_SRC
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/TRA/traadv_cen2.F90

    r106 r132  
    2020   USE lbclnk          ! ocean lateral boundary condition (or mpp link) 
    2121   USE in_out_manager  ! I/O manager 
    22    USE ptr             ! poleward transport diagnostics 
     22   USE diaptr          ! poleward transport diagnostics 
    2323 
    2424   IMPLICIT NONE 
    … …  
    307307      ENDIF 
    308308 
    309 #if defined key_diaptr 
    310309      ! "zonal" mean advective heat and salt transport  
    311       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     310      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    312311# if defined key_s_coord || defined key_partial_steps 
    313          pht_adv(:) = prt_vj( zwy(:,:,:) ) 
    314          pst_adv(:) = prt_vj( zwz(:,:,:) ) 
     312         pht_adv(:) = ptr_vj( zwy(:,:,:) ) 
     313         pst_adv(:) = ptr_vj( zwz(:,:,:) ) 
    315314# else 
    316315         DO jk = 1, jpkm1 
    … …  
    322321            END DO 
    323322         END DO 
    324          pht_adv(:) = prt_vj( zwy(:,:,:) ) 
    325          pst_adv(:) = prt_vj( zwz(:,:,:) ) 
     323         pht_adv(:) = ptr_vj( zwy(:,:,:) ) 
     324         pst_adv(:) = ptr_vj( zwz(:,:,:) ) 
    326325# endif 
    327326      ENDIF 
    328 #endif 
    329327 
    330328 

  • trunk/NEMO/OPA_SRC/TRA/traadv_cen2_atsk.h90

    r106 r132  
    256256      ENDIF 
    257257 
    258 #if defined key_diaptr 
    259258      ! "zonal" mean advective heat and salt transport 
    260       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     259      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    261260# if defined key_s_coord || defined key_partial_steps 
    262          pht_adv(:) = prt_vj( zwy(:,:,:) ) 
    263          pst_adv(:) = prt_vj( zwz(:,:,:) ) 
     261         pht_adv(:) = ptr_vj( zwy(:,:,:) ) 
     262         pst_adv(:) = ptr_vj( zwz(:,:,:) ) 
    264263# else 
    265264         DO jk = 1, jpkm1 
    … …  
    271270            END DO 
    272271         END DO 
    273          pht_adv(:) = prt_vj( zwy(:,:,:) ) 
    274          pst_adv(:) = prt_vj( zwz(:,:,:) ) 
     272         pht_adv(:) = ptr_vj( zwy(:,:,:) ) 
     273         pst_adv(:) = ptr_vj( zwz(:,:,:) ) 
    275274# endif 
    276275      ENDIF 
    277 #endif 
    278276 
    279277      ! II. Vertical advection 

  • trunk/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r106 r132  
    1919   USE lib_mpp 
    2020   USE lbclnk          ! ocean lateral boundary condition (or mpp link)  
    21    USE ptr             ! poleward transport diagnostics 
     21   USE diaptr          ! poleward transport diagnostics 
    2222 
    2323   IMPLICIT NONE 
    … …  
    274274      ENDIF 
    275275 
    276 #if defined key_diaptr  
    277276      ! "zonal" mean advective heat and salt transport  
    278       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     277      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    279278# if defined key_s_coord || defined key_partial_steps 
    280          pht_adv(:) = prt_vj( zt2(:,:,:) ) 
    281          pst_adv(:) = prt_vj( zs2(:,:,:) ) 
     279         pht_adv(:) = ptr_vj( zt2(:,:,:) ) 
     280         pst_adv(:) = ptr_vj( zs2(:,:,:) ) 
    282281# else 
    283282         DO jk = 1, jpkm1 
    … …  
    289288            END DO 
    290289         END DO 
    291          pht_adv(:) = prt_vj( zt2(:,:,:) ) 
    292          pst_adv(:) = prt_vj( zs2(:,:,:) ) 
     290         pht_adv(:) = ptr_vj( zt2(:,:,:) ) 
     291         pst_adv(:) = ptr_vj( zs2(:,:,:) ) 
    293292# endif 
    294293      ENDIF 
    295 #endif 
    296  
    297294 
    298295      ! II. Vertical advective fluxes 

  • trunk/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r106 r132  
    1919   USE lib_mpp 
    2020   USE lbclnk          ! ocean lateral boundary condition (or mpp link)  
    21    USE ptr             ! poleward transport diagnostics 
     21   USE diaptr          ! poleward transport diagnostics 
    2222 
    2323 
    … …  
    335335      ENDIF 
    336336 
    337 #if defined key_diaptr 
    338337      ! "zonal" mean advective heat and salt transport 
    339       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     338      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    340339# if defined key_s_coord || defined key_partial_steps 
    341          pht_adv(:) = prt_vj( zt2(:,:,:) ) 
    342          pst_adv(:) = prt_vj( zs2(:,:,:) ) 
     340         pht_adv(:) = ptr_vj( zt2(:,:,:) ) 
     341         pst_adv(:) = ptr_vj( zs2(:,:,:) ) 
    343342# else 
    344343         DO jk = 1, jpkm1 
    … …  
    350349            END DO 
    351350         END DO 
    352          pht_adv(:) = prt_vj( zt2(:,:,:) ) 
    353          pst_adv(:) = prt_vj( zs2(:,:,:) ) 
     351         pht_adv(:) = ptr_vj( zt2(:,:,:) ) 
     352         pst_adv(:) = ptr_vj( zs2(:,:,:) ) 
    354353# endif 
    355354      ENDIF 
    356 #endif 
    357355 
    358356      ! II. Vertical advective fluxes 

  • trunk/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r106 r132  
    2121   USE lib_mpp 
    2222   USE lbclnk          ! ocean lateral boundary condition (or mpp link)  
    23    USE ptr             ! poleward transport diagnostics 
     23   USE diaptr          ! poleward transport diagnostics 
    2424 
    2525 
    … …  
    267267      ENDIF 
    268268 
    269 #if defined key_diaptr 
    270       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
    271          ! "zonal" mean advective heat and salt transport 
    272          pht_adv(:) = prt_vj( ztv(:,:,:) ) 
    273          pst_adv(:) = prt_vj( zsv(:,:,:) ) 
     269      ! "zonal" mean advective heat and salt transport 
     270      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
     271         pht_adv(:) = ptr_vj( ztv(:,:,:) ) 
     272         pst_adv(:) = ptr_vj( zsv(:,:,:) ) 
    274273      ENDIF 
    275 #endif 
    276274 
    277275   END SUBROUTINE tra_adv_tvd 

  • trunk/NEMO/OPA_SRC/TRA/traldf_bilap.F90

    r74 r132  
    1717   USE ldfslp          ! iso-neutral slopes  
    1818   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    19    USE ptr             ! poleward transport diagnostics 
     19   USE diaptr          ! poleward transport diagnostics 
    2020 
    2121   IMPLICIT NONE 
    … …  
    8484      !!   8.5  !  02-08  (G. Madec)  F90: Free form and module 
    8585      !!---------------------------------------------------------------------- 
     86      !! * Modules used 
     87      USE oce           , ztu => ua,  &  ! use ua as workspace 
     88         &                ztv => va      ! use va as workspace 
     89 
    8690      !! * Arguments 
    8791      INTEGER, INTENT( in ) ::   kt       ! ocean time-step index 
    … …  
    9599      REAL(wp), DIMENSION(jpi,jpj) ::   &  
    96100         zeeu, zeev, zbtr,             &  ! workspace 
    97          zlt, ztu, ztv,                & 
    98          zls, zsu, zsv 
     101         zlt, zls 
     102      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   &  
     103         zsu, zsv 
    99104      !!---------------------------------------------------------------------- 
    100105 
    … …  
    134139         DO jj = 1, jpjm1 
    135140            DO ji = 1, fs_jpim1   ! vector opt. 
    136                ztu(ji,jj) = zeeu(ji,jj) * ( tb(ji+1,jj  ,jk) - tb(ji,jj,jk) ) 
    137                zsu(ji,jj) = zeeu(ji,jj) * ( sb(ji+1,jj  ,jk) - sb(ji,jj,jk) ) 
    138                ztv(ji,jj) = zeev(ji,jj) * ( tb(ji  ,jj+1,jk) - tb(ji,jj,jk) ) 
    139                zsv(ji,jj) = zeev(ji,jj) * ( sb(ji  ,jj+1,jk) - sb(ji,jj,jk) ) 
     141               ztu(ji,jj,jk) = zeeu(ji,jj) * ( tb(ji+1,jj  ,jk) - tb(ji,jj,jk) ) 
     142               zsu(ji,jj,jk) = zeeu(ji,jj) * ( sb(ji+1,jj  ,jk) - sb(ji,jj,jk) ) 
     143               ztv(ji,jj,jk) = zeev(ji,jj) * ( tb(ji  ,jj+1,jk) - tb(ji,jj,jk) ) 
     144               zsv(ji,jj,jk) = zeev(ji,jj) * ( sb(ji  ,jj+1,jk) - sb(ji,jj,jk) ) 
    140145            END DO 
    141146         END DO 
    … …  
    147152               ikv = MIN ( mbathy(ji,jj), mbathy(ji  ,jj+1) ) - 1 
    148153               IF( iku == jk ) THEN 
    149                   ztu(ji,jj) = zeeu(ji,jj) * gtu(ji,jj) 
    150                   zsu(ji,jj) = zeeu(ji,jj) * gsu(ji,jj) 
     154                  ztu(ji,jj,jk) = zeeu(ji,jj) * gtu(ji,jj) 
     155                  zsu(ji,jj,jk) = zeeu(ji,jj) * gsu(ji,jj) 
    151156               ENDIF 
    152157               IF( ikv == jk ) THEN 
    153                   ztv(ji,jj) = zeev(ji,jj) * gtv(ji,jj) 
    154                   zsv(ji,jj) = zeev(ji,jj) * gsv(ji,jj) 
     158                  ztv(ji,jj,jk) = zeev(ji,jj) * gtv(ji,jj) 
     159                  zsv(ji,jj,jk) = zeev(ji,jj) * gsv(ji,jj) 
    155160               ENDIF 
    156161            END DO 
    … …  
    161166         DO jj = 2, jpjm1 
    162167            DO ji = fs_2, fs_jpim1   ! vector opt. 
    163                zlt(ji,jj) = zbtr(ji,jj) * (  ztu(ji,jj) - ztu(ji-1,jj) + ztv(ji,jj) - ztv(ji,jj-1)  ) 
    164                zls(ji,jj) = zbtr(ji,jj) * (  zsu(ji,jj) - zsu(ji-1,jj) + zsv(ji,jj) - zsv(ji,jj-1)  ) 
     168               zlt(ji,jj) = zbtr(ji,jj) * (  ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk)  ) 
     169               zls(ji,jj) = zbtr(ji,jj) * (  zsu(ji,jj,jk) - zsu(ji-1,jj,jk) + zsv(ji,jj,jk) - zsv(ji,jj-1,jk)  ) 
    165170            END DO 
    166171         END DO 
    … …  
    183188         DO jj = 1, jpjm1 
    184189            DO ji = 1, fs_jpim1   ! vector opt. 
    185                ztu(ji,jj) = zeeu(ji,jj) * ( zlt(ji+1,jj  ) - zlt(ji,jj) ) 
    186                zsu(ji,jj) = zeeu(ji,jj) * ( zls(ji+1,jj  ) - zls(ji,jj) ) 
    187                ztv(ji,jj) = zeev(ji,jj) * ( zlt(ji  ,jj+1) - zlt(ji,jj) ) 
    188                zsv(ji,jj) = zeev(ji,jj) * ( zls(ji  ,jj+1) - zls(ji,jj) ) 
     190               ztu(ji,jj,jk) = zeeu(ji,jj) * ( zlt(ji+1,jj  ) - zlt(ji,jj) ) 
     191               zsu(ji,jj,jk) = zeeu(ji,jj) * ( zls(ji+1,jj  ) - zls(ji,jj) ) 
     192               ztv(ji,jj,jk) = zeev(ji,jj) * ( zlt(ji  ,jj+1) - zlt(ji,jj) ) 
     193               zsv(ji,jj,jk) = zeev(ji,jj) * ( zls(ji  ,jj+1) - zls(ji,jj) ) 
    189194            END DO 
    190195         END DO 
    … …  
    194199            DO ji = fs_2, fs_jpim1   ! vector opt. 
    195200               ! horizontal diffusive trends 
    196                zta = zbtr(ji,jj) * (  ztu(ji,jj) - ztu(ji-1,jj) + ztv(ji,jj) - ztv(ji,jj-1)  ) 
    197                zsa = zbtr(ji,jj) * (  zsu(ji,jj) - zsu(ji-1,jj) + zsv(ji,jj) - zsv(ji,jj-1)  ) 
     201               zta = zbtr(ji,jj) * (  ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk)  ) 
     202               zsa = zbtr(ji,jj) * (  zsu(ji,jj,jk) - zsu(ji-1,jj,jk) + zsv(ji,jj,jk) - zsv(ji,jj-1,jk)  ) 
    198203               ! add it to the general tracer trends 
    199204               ta(ji,jj,jk) = ta(ji,jj,jk) + zta 
    … …  
    210215      !                                                ! =============== 
    211216 
    212 #if defined key_diaptr 
    213217      ! "zonal" mean lateral diffusive heat and salt transport 
    214       == >> forced bug NOT implemented, ztv, zsv not 3D arrays 
    215       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     218      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    216219# if defined key_s_coord || defined key_partial_steps 
    217          pht_ldf(:) = prt_vj( ztv(:,:) ) 
    218          pst_ldf(:) = prt_vj( zsv(:,:) ) 
     220         pht_ldf(:) = ptr_vj( ztv(:,:,:) ) 
     221         pst_ldf(:) = ptr_vj( zsv(:,:,:) ) 
    219222# else 
    220          DO jj = 2, jpjm1 
    221             DO ji = fs_2, fs_jpim1   ! vector opt. 
    222               ztv(ji,jj) = ztv(ji,jj) * fse3v(ji,jj,jk) 
    223               zsv(ji,jj) = zsv(ji,jj) * fse3v(ji,jj,jk) 
    224             END DO 
    225          END DO 
    226          pht_ldf(:) = prt_vj( ztv(:,:) ) 
    227          pst_ldf(:) = prt_vj( zsv(:,:) ) 
     223         DO jk = 1, jpkm1 
     224            DO jj = 2, jpjm1 
     225               DO ji = fs_2, fs_jpim1   ! vector opt. 
     226                  ztv(ji,jj,jk) = ztv(ji,jj,jk) * fse3v(ji,jj,jk) 
     227                  zsv(ji,jj,jk) = zsv(ji,jj,jk) * fse3v(ji,jj,jk) 
     228               END DO 
     229            END DO 
     230         END DO 
     231         pht_ldf(:) = ptr_vj( ztv(:,:,:) ) 
     232         pst_ldf(:) = ptr_vj( zsv(:,:,:) ) 
    228233# endif 
    229234      ENDIF 
    230 #endif 
    231235 
    232236   END SUBROUTINE tra_ldf_bilap 

  • trunk/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

    r74 r132  
    2020   USE ldfslp          ! iso-neutral slopes available 
    2121   USE lbclnk          ! ocean lateral boundary condition (or mpp link) 
    22    USE ptr             ! poleward transport diagnostics  
     22   USE diaptr          ! poleward transport diagnostics  
    2323 
    2424   IMPLICIT NONE 
    … …  
    6767      !!   8.5  !  02-08  (G. Madec)  F90: Free form and module 
    6868      !!---------------------------------------------------------------------- 
     69      !! * Modules used 
     70      USE oce           , wk1 => ua,  &  ! use ua as workspace 
     71         &                wk2 => va      ! use va as workspace 
     72 
    6973      !! * Arguments 
    7074      INTEGER, INTENT( in ) ::   kt           ! ocean time-step index 
    … …  
    7377      INTEGER ::   ji, jj, jk                 ! dummy loop indices 
    7478      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
    75          wk1, wk2,            &  ! work array used for rotated biharmonic 
    76          wk3, wk4                ! operator on tracers and/or momentum 
     79         wk3, wk4                ! work array used for rotated biharmonic 
     80         !                       ! operator on tracers and/or momentum 
    7781      !!---------------------------------------------------------------------- 
    7882 
    … …  
    189193         zcof0, zcof1, zcof2,          & 
    190194         zcof3, zcof4 
    191       REAL(wp), DIMENSION(jpi,jpj) ::   & 
    192          zftu, zftv, zfsu, zfsv,       &  ! workspace 
     195      REAL(wp), DIMENSION(jpi,jpj) ::  & 
     196         zftu, zfsu,                   &  ! workspace 
    193197         zdkt, zdk1t,                  & 
    194198         zdks, zdk1s 
    195       REAL(wp), DIMENSION(jpi,jpk) ::   & 
     199      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
     200         zftv, zfsv                       ! workspace (only v components for ptr) 
     201      REAL(wp), DIMENSION(jpi,jpk) ::  & 
    196202         zftw, zfsw,                   &  ! workspace 
    197203         zdit, zdjt, zdj1t,            & 
    … …  
    240246                             +zdk1t(ji+1,jj) + zdkt (ji,jj) )  ) 
    241247 
    242                zftv(ji,jj)= vmask(ji,jj,jk) *   & 
     248               zftv(ji,jj,jk)= vmask(ji,jj,jk) *   & 
    243249                  (  zabe2 *( pt(ji,jj+1,jk) - pt(ji,jj,jk) )   & 
    244250                   + zcof2 *( zdkt (ji,jj+1) + zdk1t(ji,jj)     & 
    … …  
    250256                             +zdk1s(ji+1,jj) + zdks (ji,jj) )  ) 
    251257 
    252                zfsv(ji,jj)= vmask(ji,jj,jk) *   & 
     258               zfsv(ji,jj,jk)= vmask(ji,jj,jk) *   & 
    253259                  (  zabe2 *( ps(ji,jj+1,jk) - ps(ji,jj,jk) )   & 
    254260                   + zcof2 *( zdks (ji,jj+1) + zdk1s(ji,jj)     & 
    … …  
    263269         DO jj = 2 , jpjm1 
    264270            DO ji = 2 , jpim1 
    265                ztah = zftu(ji,jj) - zftu(ji-1,jj) + zftv(ji,jj) - zftv(ji,jj-1) 
    266                zsah = zfsu(ji,jj) - zfsu(ji-1,jj) + zfsv(ji,jj) - zfsv(ji,jj-1) 
     271               ztah = zftu(ji,jj) - zftu(ji-1,jj) + zftv(ji,jj,jk) - zftv(ji,jj-1,jk) 
     272               zsah = zfsu(ji,jj) - zfsu(ji-1,jj) + zfsv(ji,jj,jk) - zfsv(ji,jj-1,jk) 
    267273               plt(ji,jj,jk) = ztah 
    268274               pls(ji,jj,jk) = zsah 
    … …  
    273279      !                                                ! =============== 
    274280 
    275       IF( kaht == 2 .AND. MOD( kt, nf_ptr ) == 0 ) THEN 
    276 #if defined key_diaptr 
    277          ! "zonal" mean diffusive heat and salt transport 
    278 !         pht_ldf(:) = prt_vj( zftv(:,:) ) 
    279 !         pst_ldf(:) = prt_vj( zfsv(:,:) ) 
    280         write(numout,cform_err) 
    281         nstop = nstop + 1 
    282 #endif 
     281      ! "zonal" mean diffusive heat and salt transport 
     282      IF( ln_diaptr .AND. ( kaht == 2 ) .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
     283         pht_ldf(:) = ptr_vj( zftv(:,:,:) ) 
     284         pst_ldf(:) = ptr_vj( zfsv(:,:,:) ) 
    283285      ENDIF 
    284  
    285286 
    286287      !                             ! ************ !   ! =============== 

  • trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90

    r106 r132  
    1919   USE in_out_manager  ! I/O manager 
    2020   USE ldfslp          ! iso-neutral slopes 
    21    USE ptr             ! poleward transport diagnostics 
     21   USE diaptr          ! poleward transport diagnostics 
    2222 
    2323   IMPLICIT NONE 
    … …  
    239239      ENDIF 
    240240 
    241 #if defined key_diaptr 
    242241      !!bug  no separation of diff iso and eiv 
    243       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
    244       ! "zonal" mean lateral diffusive heat and salt transports 
    245          pht_ldf(:) = prt_vj( zftv(:,:,:) ) 
    246          pst_ldf(:) = prt_vj( zfsv(:,:,:) ) 
     242      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
     243         ! "zonal" mean lateral diffusive heat and salt transports 
     244         pht_ldf(:) = ptr_vj( zftv(:,:,:) ) 
     245         pst_ldf(:) = ptr_vj( zfsv(:,:,:) ) 
    247246         ! "zonal" mean lateral eddy induced velocity heat and salt transports 
    248          pht_eiv(:) = prt_vj( zftv(:,:,:) ) 
    249          pst_eiv(:) = prt_vj( zfsv(:,:,:) ) 
     247         pht_eiv(:) = ptr_vj( zftv(:,:,:) ) 
     248         pst_eiv(:) = ptr_vj( zfsv(:,:,:) ) 
    250249      ENDIF 
    251 #endif 
    252250 
    253251   END SUBROUTINE tra_ldf_iso 

  • trunk/NEMO/OPA_SRC/TRA/traldf_iso_zps.F90

    r106 r132  
    1919   USE in_out_manager  ! I/O manager 
    2020   USE ldfslp          ! iso-neutral slopes 
    21    USE ptr             ! poleward transport diagnostics 
     21   USE diaptr          ! poleward transport diagnostics 
    2222 
    2323 
    … …  
    276276 
    277277 
    278 #if defined key_diaptr 
    279278      !!bug  no separation of diff iso and eiv 
    280       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     279      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    281280         ! "zonal" mean lateral diffusive heat and salt transports 
    282          pht_ldf(:) = prt_vj( zftv(:,:,:) ) 
    283          pst_ldf(:) = prt_vj( zfsv(:,:,:) ) 
     281         pht_ldf(:) = ptr_vj( zftv(:,:,:) ) 
     282         pst_ldf(:) = ptr_vj( zfsv(:,:,:) ) 
    284283         ! "zonal" mean lateral eddy induced velocity heat and salt transports 
    285          pht_eiv(:) = prt_vj( zftv(:,:,:) ) 
    286          pst_eiv(:) = prt_vj( zfsv(:,:,:) ) 
     284         pht_eiv(:) = ptr_vj( zftv(:,:,:) ) 
     285         pst_eiv(:) = ptr_vj( zfsv(:,:,:) ) 
    287286      ENDIF 
    288 #endif 
    289287 
    290288   END SUBROUTINE tra_ldf_iso_zps 

  • trunk/NEMO/OPA_SRC/TRA/traldf_lap.F90

    r106 r132  
    1515   USE trdtra_oce      ! ocean active tracer trends 
    1616   USE in_out_manager  ! I/O manager 
    17    USE ptr             ! poleward transport diagnostics 
     17   USE diaptr          ! poleward transport diagnostics 
    1818 
    1919 
    … …  
    149149      ENDIF 
    150150 
    151 #if defined key_diaptr 
    152151      ! "zonal" mean lateral diffusive heat and salt transport  
    153       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     152      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    154153# if defined key_s_coord || defined key_partial_steps 
    155          pht_ldf(:) = prt_vj( ztv(:,:,:) ) 
    156          pst_ldf(:) = prt_vj( zsv(:,:,:) ) 
     154         pht_ldf(:) = ptr_vj( ztv(:,:,:) ) 
     155         pst_ldf(:) = ptr_vj( zsv(:,:,:) ) 
    157156# else 
    158157         DO jk = 1, jpkm1 
    … …  
    164163            END DO 
    165164         END DO 
    166          pht_ldf(:) = prt_vj( ztv(:,:,:) ) 
    167          pst_ldf(:) = prt_vj( zsv(:,:,:) ) 
     165         pht_ldf(:) = ptr_vj( ztv(:,:,:) ) 
     166         pst_ldf(:) = ptr_vj( zsv(:,:,:) ) 
    168167# endif 
    169168      ENDIF 
    170 #endif 
    171169 
    172170   END SUBROUTINE tra_ldf_lap 

  • trunk/NEMO/OPA_SRC/opa.F90

    r109 r132  
    4343   USE flxfwb          !  
    4444 
    45    USE ptr             ! poleward transports           (dia_ptr_init routine) 
     45   USE diaptr          ! poleward transports           (dia_ptr_init routine) 
    4646 
    4747   USE step            ! OPA time-stepping                  (stp     routine) 
    … …  
    199199      CALL flx_fwb_init                     ! FreshWater Budget correction 
    200200 
    201 #if defined key_diaptr 
    202       CALL dia_ptr_init                     ! Poleward TRansports 
    203 #endif 
     201      IF( ln_diaptr ) CALL dia_ptr_init     ! Poleward TRansports initialization 
    204202 
    205203      !                                     ! =============== ! 

  • trunk/NEMO/OPA_SRC/step.F90

    r108 r132  
    109109   USE trdvor          ! vorticity budget                 (trd_vor routine) 
    110110   USE diahdy          ! dynamic height                   (dia_hdy routine) 
    111    USE ptr             ! poleward transports         (dia_ptr_init routine) 
    112    !                   !                                  (dia_ptr routine) 
     111   USE diaptr          ! poleward transports              (dia_ptr routine) 
    113112   USE diahth          ! thermocline depth                (dia_hth routine) 
    114113   USE diafwb          ! freshwater budget                (dia_fwb routine) 
    … …  
    210209            CALL dia_wri_state( 'output.init' ) 
    211210         ENDIF 
    212          IF( lk_diaptr ) CALL dia_ptr_init            ! Poleward TRansport initialization 
    213  
    214211      ENDIF 
    215212 
    … …  
    447444         IF( lk_diahdy  )   CALL dia_hdy( kstp )                 ! dynamical heigh diagnostics 
    448445         IF( lk_diafwb  )   CALL dia_fwb( kstp )                 ! Fresh water budget diagnostics 
    449 #if defined key_diaptr 
    450          IF( kstp == nit000 .OR. MOD( kstp, nf_ptr ) == 0 )   & 
     446 
     447         IF( ln_diaptr .AND. ( kstp == nit000 .OR. MOD( kstp, nf_ptr ) == 0 ) )  & 
    451448                            CALL dia_ptr( kstp )                 ! Poleward TRansports diagnostics 
    452 #endif 
    453449 
    454450         !                                         ! save and outputs 
Note: See TracChangeset for help on using the changeset viewer.