Changeset 9450 for branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90

Timestamp:

2018-03-30T18:08:29+02:00 (6 years ago)

Author:

aumont

Message:

debug PISCES code

File:

• branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90 (modified) (5 diffs)

branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90

@@ -71 +71 @@
       REAL(wp) :: zgrazp, zgrazm, zgrazsd
       REAL(wp) :: zgrazmf, zgrazsf, zgrazpf
-      REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d, zfezoo
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zzligprod
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
@@ -77 +78 @@
       IF( nn_timing == 1 )  CALL timing_start('p4z_micro')
       !
-      CALL wrk_alloc( jpi, jpj, jpk, zgrazing )
+      CALL wrk_alloc( jpi, jpj, jpk, zgrazing, zfezoo )
+      IF (ln_ligand) THEN
+         CALL wrk_alloc( jpi, jpj, jpk, zzligprod )
+         zzligprod(:,:,:) = 0._wp
+      ENDIF
       !
       DO jk = 1, jpkm1
@@ -139 +144 @@
                tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zgrarem - zgrarsig
                !
-               IF( ln_ligand ) tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + (zgrarem - zgrarsig) * ldocz
+               IF( ln_ligand ) THEN
+                  tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + (zgrarem - zgrarsig) * ldocz
+                  zzligprod(ji,jj,jk) = (zgrarem - zgrarsig) * ldocz
+               ENDIF
                !
                tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2ut * zgrarsig
                tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zgrafer
+               zfezoo(ji,jj,jk) = zgrafer
                tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zgrapoc
                prodpoc(ji,jj,jk) = prodpoc(ji,jj,jk) + zgrapoc
@@ -184 +193 @@
               CALL iom_put( "GRAZ1", zw3d )
            ENDIF
+           IF( iom_use( "FEZOO" ) ) THEN
+              zw3d(:,:,:) = zfezoo(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:)   !
+              CALL iom_put( "FEZOO", zw3d )
+           ENDIF
+           IF( iom_use( "LPRODZ" ) .AND. ln_ligand )  THEN
+              zw3d(:,:,:) = zzligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:)
+              CALL iom_put( "LPRODZ"  , zw3d )
+           ENDIF
            CALL wrk_dealloc( jpi, jpj, jpk, zw3d )
          ENDIF
@@ -194 +211 @@
       ENDIF
       !
-      CALL wrk_dealloc( jpi, jpj, jpk, zgrazing )
+      CALL wrk_dealloc( jpi, jpj, jpk, zgrazing, zfezoo )
+      IF (ln_ligand) CALL wrk_dealloc( jpi, jpj, jpk, zzligprod )
       !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_micro')