Changeset 7617 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.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/p4zlys.F90 (modified) (4 diffs)

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

@@ -67 +67 @@
       CHARACTER (len=25) :: charout
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zco3, zcaldiss, zhinit, zhi, zco3sat
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zwork
       !!---------------------------------------------------------------------
       !
@@ -72 +73 @@
       !
       CALL wrk_alloc( jpi, jpj, jpk, zco3, zcaldiss, zhinit, zhi, zco3sat )
+      CALL wrk_alloc( jpi, jpj, jpk, zwork )
       !
       zco3    (:,:,:) = 0.
       zcaldiss(:,:,:) = 0.
+      zwork   (:,:,:) = 0.
       zhinit(:,:,:)   = hi(:,:,:) * 1000. / ( rhop(:,:,:) + rtrn )
       !     -------------------------------------------
@@ -133 +136 @@
          !
       IF( lk_iomput .AND. knt == nrdttrc ) THEN
-         IF( iom_use( "PH"     ) ) CALL iom_put( "PH"    , -1. * LOG10( hi(:,:,:) )          * tmask(:,:,:) )
+         IF( iom_use( "PH"     ) ) CALL iom_put( "PH"    , -1. * LOG10( hi(:,:,:) + rtrn )          * tmask(:,:,:) )
          IF( iom_use( "CO3"    ) ) CALL iom_put( "CO3"   , zco3(:,:,:) * 1.e+3               * tmask(:,:,:) )
          IF( iom_use( "CO3sat" ) ) CALL iom_put( "CO3sat", zco3sat(:,:,:) * 1.e+3            * tmask(:,:,:) )
          IF( iom_use( "DCAL"   ) ) CALL iom_put( "DCAL"  , zcaldiss(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) )
-
+         IF( iom_use( "AOU"    ) ) THEN 
+             zwork(:,:,:) = MAX(0., (chemo2(:,:,:) - trn(:,:,:,jpoxy) ) ) * 1.E6 * tmask(:,:,:)
+             CALL iom_put( "AOU"   , MAX(0., zwork(:,:,:) ) )
+         ENDIF
       ELSE
          IF( ln_diatrc ) THEN
@@ -154 +160 @@
       !
       CALL wrk_dealloc( jpi, jpj, jpk, zco3, zcaldiss, zhinit, zhi, zco3sat )
+      CALL wrk_dealloc( jpi, jpj, jpk, zwork )
       !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_lys')