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 74 for trunk/NEMO/OPA_SRC/TRA/traadv_cen2.F90 – NEMO

Ignore:
Timestamp:
2004-04-22T14:38:42+02:00 (20 years ago)
Author:
opalod
Message:

CT : BUGFIX048 : Bug correction of Poleward Transport diagnostics

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/TRA/traadv_cen2.F90

    r32 r74  
    1313   USE oce             ! ocean dynamics and active tracers 
    1414   USE dom_oce         ! ocean space and time domain 
    15    USE trdtra_oce     ! ocean active tracer trends 
     15   USE trdtra_oce      ! ocean active tracer trends 
    1616   USE flxrnf          ! 
    1717   USE trabbl          ! advective term in the BBL 
    18    USE ocfzpt           ! 
     18   USE ocfzpt          ! 
    1919   USE lib_mpp 
    20    USE lbclnk 
     20   USE lbclnk          ! ocean lateral boundary condition (or mpp link) 
    2121   USE in_out_manager  ! I/O manager 
     22   USE ptr             ! poleward transport diagnostics 
    2223 
    2324   IMPLICIT NONE 
     
    299300      !                                                ! =============== 
    300301 
    301       IF( l_ctl .AND. lwp ) THEN 
     302      IF(l_ctl) THEN 
    302303         zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    303304         zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     
    310311      IF( MOD( kt, nf_ptr ) == 0 ) THEN 
    311312# if defined key_s_coord || defined key_partial_steps 
    312          pht_adv(:,:) = prt_vj( zwy(:,:,:) ) 
    313          pst_adv(:,:) = prt_vj( zwz(:,:,:) ) 
     313         pht_adv(:) = prt_vj( zwy(:,:,:) ) 
     314         pst_adv(:) = prt_vj( zwz(:,:,:) ) 
    314315# else 
    315316         DO jk = 1, jpkm1 
     
    321322            END DO 
    322323         END DO 
    323          pht_adv(:,:) = prt_vj( zwy(:,:,:) ) 
    324          pst_adv(:,:) = prt_vj( zwz(:,:,:) ) 
     324         pht_adv(:) = prt_vj( zwy(:,:,:) ) 
     325         pst_adv(:) = prt_vj( zwz(:,:,:) ) 
    325326# endif 
    326327      ENDIF 
     
    394395      END DO 
    395396 
    396       IF( l_ctl .AND. lwp ) THEN 
     397      IF(l_ctl) THEN 
    397398         zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    398399         zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
Note: See TracChangeset for help on using the changeset viewer.