Changeset 7617 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90

Timestamp:

2017-01-27T18:17:01+01:00 (7 years ago)

Author:

aumont

Message:

update diagnostics + changes in quota code

File:

• branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90 (modified) (3 diffs)

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90

@@ -59 +59 @@
       REAL(wp) ::   zlgwp, zlgwpr, zlgwr, zlablgw, zrfepa, zfepr
       REAL(wp) ::   zstep, zstep2, zrfact2
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zligrem, zligpr, zrligprod
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
@@ -89 +90 @@
                zlgwpr = prlgw * zstep * etot(ji,jj,jk) * trn(ji,jj,jk,jplgw) * (1. - fr_i(ji,jj))
                tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zlgwp - zlgwr - zlgwpr
+               zligrem(ji,jj,jk)   = zlgwr
+               zligpr(ji,jj,jk)    = zlgwpr
+               zrligprod(ji,jj,jk) = zlgwp
             END DO
          END DO
@@ -111 +115 @@
          END DO
       END DO
+
+      !  Output of some diagnostics variables
+      !     ---------------------------------
+      IF( ln_diatrc .AND. lk_iomput ) THEN
+         IF( jnt == nrdttrc ) THEN
+            CALL iom_put("LIGREM"   , zligrem(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) )
+            CALL iom_put("LIGPR"    , zligpr(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) )
+            CALL iom_put("LPRODR"   , zrligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) )
+         ENDIF
+      ENDIF
+
+      CALL wrk_dealloc( jpi, jpj, jpk, zligrem, zligpr, zrligprod )
 
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)