Changeset 8453 for branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90

Timestamp:

2017-08-22T18:54:47+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM --22-08-2017-- repair trends for passive tracers

File:

• branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90 (modified) (3 diffs)

branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90

@@ -91 +91 @@
                                    xxi
       USE trc,               ONLY: med_diag, tra
+      USE trd_oce,           ONLY: jptra_sms 
+      USE trdtrc
 
    !!* Substitution
@@ -119 +121 @@
       !! temporary variables
       REAL(wp) :: fq0
+
+      !! trend temporary array:
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrmed
 
       !!==========================================================
@@ -734 +739 @@
       ENDDO
 
+      !! MEDUSA SMS trends:
+      IF( l_trdtrc ) THEN
+          CALL wrk_alloc( jpi, jpj, jpk, ztrmed )
+          DO jn = jp_msa0, jp_msa1
+            ztrmed(:,:,:) = tra(:,:,:,jn)
+            CALL trd_trc( ztrmed, jn, jptra_sms, kt )   ! save trends
+          END DO
+          CALL wrk_dealloc( jpi, jpj, jpk, ztrmed )
+      END IF
+
       DO jj = 2,jpjm1
          DO ji = 2,jpim1