Changeset 1830 for branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zopt.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/p4zopt.F90 (modified) (2 diffs)

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

@@ -61 +61 @@
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zekg, zekr, zekb
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   ze1 , ze2 , ze3, ze0
-#if defined key_trc_diaadd && defined key_iomput
-     REAL(wp), DIMENSION(jpi,jpj)      ::   zw2d
-     REAL(wp), DIMENSION(jpi,jpj,jpk)  ::   zw3d
-#endif
       !!---------------------------------------------------------------------
 
@@ -238 +234 @@
 # else
       ! write diagnostics 
-      zw2d(:,:  ) =  heup(:,:  ) * tmask(:,:,1)
-      zw3d(:,:,:) =  etot(:,:,:) * tmask(:,:,:)
-      IF( jnt == nrdttrc ) CALL iom_put( "Heup", zw2d )               
-      IF( jnt == nrdttrc ) CALL iom_put( "PAR" , zw3d )
+      IF( jnt == nrdttrc ) then 
+         CALL iom_put( "Heup", heup(:,:  ) * tmask(:,:,1) )  ! euphotic layer deptht
+         CALL iom_put( "PAR" , etot(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation
+      ENDIF
 # endif
 #endif