Changeset 1880 for branches/CMIP5_IPSL/NEMO/OPA_SRC/TRA/traadv_tvd.F90

Timestamp:

2010-05-21T09:39:46+02:00 (14 years ago)

Author:

aclsce

Message:

Put in CMIP5_IPSL branch bug fixes for diaptr code, see ticket:#361

File:

• branches/CMIP5_IPSL/NEMO/OPA_SRC/TRA/traadv_tvd.F90 (modified) (2 diffs)

branches/CMIP5_IPSL/NEMO/OPA_SRC/TRA/traadv_tvd.F90

@@ -179 +179 @@
       END DO
 
+      ! "zonal" mean advective heat and salt transport
+      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN
+         pht_adv(:) = ptr_vj( ztv(:,:,:) )
+         pst_adv(:) = ptr_vj( zsv(:,:,:) )
+      ENDIF
 
       ! Save the intermediate i / j / k advective trends for diagnostics
@@ -366 +371 @@
       ! "zonal" mean advective heat and salt transport
       IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN
-         pht_adv(:) = ptr_vj( ztv(:,:,:) )
-         pst_adv(:) = ptr_vj( zsv(:,:,:) )
+         pht_adv(:) = ptr_vj( ztv(:,:,:) ) + pht_adv(:)
+         pst_adv(:) = ptr_vj( zsv(:,:,:) ) + pst_adv(:)
       ENDIF
       !