Changeset 9450 for branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.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/p4zpoc.F90 (modified) (9 diffs)

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

@@ -67 +67 @@
       REAL(wp), POINTER, DIMENSION(:,:,:)   :: zremipoc, zremigoc, zorem3, ztremint
       REAL(wp), POINTER, DIMENSION(:,:,:,:) :: alphag
+      REAL(wp), POINTER, DIMENSION(:,:,:)   :: zfolimi
       !!---------------------------------------------------------------------
       !
@@ -74 +75 @@
       CALL wrk_alloc( jpi, jpj,      totprod,  totthick, totcons )
       CALL wrk_alloc( jpi, jpj, jpk, zremipoc, zremigoc, zorem3, ztremint )
+      CALL wrk_alloc( jpi, jpj, jpk, zfolimi )
       ALLOCATE( alphag(jpi,jpj,jpk,jcpoc) )
 
@@ -98 +100 @@
       orem  (:,:,:)   = 0.
       ztremint(:,:,:) = 0.
+      zfolimi(:,:,:)  = 0.
 
       DO jn = 1, jcpoc
@@ -217 +220 @@
                   tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zorem2
                   tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer2
+                  zfolimi(ji,jj,jk)   = zofer2
                END DO
             END DO
@@ -247 +251 @@
                   tra(ji,jj,jk,jpgop) = tra(ji,jj,jk,jpgop) - zopop2 * (1. + solgoc)
                   tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) - zofer2 * (1. + solgoc)
+                  zfolimi(ji,jj,jk)   = zofer2
                END DO
             END DO
@@ -422 +427 @@
                     tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zorem
                     tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zofer
+                    zfolimi(ji,jj,jk)   = zfolimi(ji,jj,jk) + zofer
                   ENDIF
                END DO
@@ -447 +453 @@
                 tra(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + zopop 
                 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer 
+                zfolimi(ji,jj,jk)   = zfolimi(ji,jj,jk) + zofer
              END DO
            END DO
@@ -457 +464 @@
           CALL iom_put( "REMINP" , zremipoc(:,:,:) * tmask(:,:,:) )  ! Remineralisation rate
           CALL iom_put( "REMING" , zremigoc(:,:,:) * tmask(:,:,:) )  ! Remineralisation rate
+          CALL iom_put( "REMINF" , zfolimi(:,:,:) * tmask(:,:,:)  * 1.e+9 * zrfact2 )  ! Remineralisation rate
         ENDIF
      ENDIF
@@ -468 +476 @@
       CALL wrk_dealloc( jpi, jpj,      totprod,  totthick, totcons )
       CALL wrk_dealloc( jpi, jpj, jpk, zremipoc, zremigoc, zorem3, ztremint )
+      CALL wrk_dealloc( jpi, jpj, jpk, zfolimi )
       DEALLOCATE( alphag )
       !