Changeset 2072 for branches/devmercator2010/NEMO/TOP_SRC/PISCES/p4zsed.F90

Timestamp:

2010-09-08T16:29:38+02:00 (14 years ago)

Author:

cbricaud

Message:

add change from DEV_r1784_3DF

File:

• branches/devmercator2010/NEMO/TOP_SRC/PISCES/p4zsed.F90 (modified) (2 diffs)

branches/devmercator2010/NEMO/TOP_SRC/PISCES/p4zsed.F90

@@ -96 +96 @@
       REAL(wp) :: zrfact2
 # if defined key_iomput
+     REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zw3d 
      REAL(wp), DIMENSION(jpi,jpj)    ::    zw2d 
 # endif
@@ -331 +332 @@
       trc2d(:,:,jp_pcs0_2d + 12) = znitrpot(:,:,1) * 1.e-7 * zrfact2 * fse3t(:,:,1) * tmask(:,:,1)
 # else
-      ! surface downward net flux of iron
-      zw2d(:,:)   =  ( zirondep(:,:,1) + ironsed(:,:,1) * rfact2 ) * zrfact2 * fse3t(:,:,1) * tmask(:,:,1) 
+      ! write diagnostics
+      zw2d(:,:)   =  ( zirondep(:,:,1) + ironsed(:,:,1) * rfact2 ) &
+      &            * zrfact2 * fse3t(:,:,1) * tmask(:,:,1)     
       IF( jnt == nrdttrc ) CALL iom_put( "Irondep", zw2d )
-      ! nitrogen fixation at surface
-      zw2d(:,:)   =  znitrpot(:,:,1) * 1.e-7 * zrfact2  * fse3t(:,:,1) * tmask(:,:,1)
-      IF( jnt == nrdttrc ) CALL iom_put( "Nfix" , zw2d )
-# endif
+      zw3d(:,:,:) = znitrpot(:,:,:) * 1.e-7 * zrfact2  * fse3t(:,:,:) * tmask(:,:,:)
+      IF( jnt == nrdttrc ) CALL iom_put( "Nfix", zw3d  ) 
+# endif
+
 # endif
       !