Changeset 15651

Timestamp:

2022-01-18T13:03:56+01:00 (2 years ago)

Author:

hadjt

Message:

TRA/tradwl.F90

Code to report the trends fixed. This now works without affecting the stratification (unlike the original CO6 code), but hasn't been tested in terms of the trends.

The choice of which code to report the trends too (jptra_qsr) may not be correct.

File:

• NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/TRA/tradwl.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/TRA/tradwl.F90

@@ -74 +74 @@
       !JT USE oce, ONLY :   ztrds => va   ! use va as 3D workspace   
 
+      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrdt, ztrds   ! 3D workspace
+
       !!
       INTEGER, INTENT(in) ::   kt     ! ocean time-step
@@ -109 +111 @@
       !JT    ztrds(:,:,:) = 0.e0
       !JT ENDIF
+
+
+      IF( l_trdtra )   THEN                  !* Save ta and sa trends
+         ALLOCATE( ztrdt(jpi,jpj,jpk) , ztrds(jpi,jpj,jpk) )
+         ztrdt(:,:,:) = tsa(:,:,:,jp_tem)
+         ztrds(:,:,:) = tsa(:,:,:,jp_sal)
+      ENDIF
 !--------------------------------------------------------------------
 !  Set transmissivity
@@ -204 +213 @@
       !JT    !CEODCALL trd_mod( ztrdt, ztrds, jptra_trd_qsr, 'TRA', kt )
       !JT ENDIF
-      !                       ! print mean trends (used for debugging)
+      !                       
+      IF( l_trdtra ) THEN      ! qsr tracers trends saved for diagnostics
+
+         !JT I think I should use jptra_qsr?? 
+
+         !ztrdt(:,:,:) = tsa(:,:,:,jp_tem) - ztrdt(:,:,:)
+         !ztrds(:,:,:) = tsa(:,:,:,jp_sal) - ztrds(:,:,:)
+         DO jk = 1, jpkm1
+            ztrdt(:,:,jk) = tsa(:,:,jk,jp_tem) - ztrdt(:,:,jk)
+            ztrds(:,:,jk) = tsa(:,:,jk,jp_sal) - ztrds(:,:,jk)
+         END DO
+
+         CALL trd_tra( kt, 'TRA', jp_tem, jptra_qsr, ztrdt )
+         CALL trd_tra( kt, 'TRA', jp_sal, jptra_qsr, ztrds )
+         DEALLOCATE( ztrdt , ztrds )
+      ENDIF
+
       IF(ln_ctl)   CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' qsr  - Ta: ', mask1=tmask, clinfo3='tra-ta' )
       !