Changeset 7550 for branches/2015

Timestamp:

2017-01-12T10:28:57+01:00 (7 years ago)

Author:

davestorkey

Message:

Add calls for total tracer trends diagnostics (one of the CMIP6 diagnostics) in tra_nxt.

File:

• branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90 (modified) (1 diff)

branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

@@ -135 +135 @@
             CALL trd_tra( kt, 'TRA', jp_sal, jptra_zdfp, ztrds )
          ENDIF
+         ! total trend for the non-time-filtered variables. 
+         ztrdt(:,:,:) = 0._wp
+         ztrds(:,:,:) = 0._wp
+         DO jk = 1, jpkm1
+            zfact = 1.0 / rdttra(jk)
+            ztrdt(:,:,jk) = ( tsa(:,:,jk,jp_tem) - tsn(:,:,jk,jp_tem) ) * zfact 
+            ztrds(:,:,jk) = ( tsa(:,:,jk,jp_sal) - tsn(:,:,jk,jp_sal) ) * zfact 
+         END DO
+         CALL trd_tra( kt, 'TRA', jp_tem, jptra_tot, ztrdt )
+         CALL trd_tra( kt, 'TRA', jp_sal, jptra_tot, ztrds )
       ENDIF