Changeset 1830 for branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zlys.F90

Timestamp:

2010-04-12T15:03:51+02:00 (14 years ago)

Author:

cetlod

Message:

Computation of additional diagnostics for PISCES model ( under CPP key key_diaar5 )

• needed for AR5 outputs (vertical inventories, passive tracers at surface,... )
• new output file with suffix dbio_T

File:

• branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zlys.F90 (modified) (4 diffs)

branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zlys.F90

@@ -67 +67 @@
 #if defined key_trc_dia3d && defined key_iomput
       REAL(wp) ::   zrfact2
-      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zcaldiss, zw3d
+      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zcaldiss
 #endif
       CHARACTER (len=25) :: charout
@@ -94 +94 @@
                   ! SET DUMMY VARIABLE FOR TOTAL BORATE
                   zbot  = borat(ji,jj,jk)
+
+                  ! SET DUMMY VARIABLE FOR TOTAL BORATE
+                  zbot  = borat(ji,jj,jk)
                   zfact = rhop (ji,jj,jk) / 1000. + rtrn
 
@@ -171 +174 @@
 #  else
       zrfact2 = 1.e3 * rfact2r
-      zw3d(:,:,:) = hi  (:,:,:)                    * tmask(:,:,:)
-      CALL iom_put( "PH", zw3d )
-      zw3d(:,:,:) = zco3(:,:,:)                    * tmask(:,:,:)
-      CALL iom_put( "CO3", zw3d )
-      zw3d(:,:,:) = aksp(:,:,:) / calcon           * tmask(:,:,:)
-      CALL iom_put( "CO3sat", zw3d )
-      zw3d(:,:,:) = zcaldiss(:,:,:) * zrfact2 * tmask(:,:,:)
-      CALL iom_put( "Dcal", zw3d )
+      CALL iom_put( "PH"    , hi      (:,:,:)           * tmask(:,:,:) )
+      CALL iom_put( "CO3"   , zco3    (:,:,:)           * tmask(:,:,:) )
+      CALL iom_put( "CO3sat", aksp    (:,:,:) / calcon  * tmask(:,:,:) )
+      CALL iom_put( "DCAL"  , zcaldiss(:,:,:) * zrfact2 * tmask(:,:,:) )
+#if defined key_diaar5
+      CALL iom_put( "PHSFC" , hi     (:,:,1)           * tmask(:,:,1) )
+#endif
 #  endif
 # endif
@@ -232 +234 @@
    END SUBROUTINE p4z_lys
 #endif 
-
    !!======================================================================
 END MODULE  p4zlys