Changeset 1457 for trunk/NEMO/TOP_SRC/PISCES/p4zlys.F90

Timestamp:

2009-05-23T10:16:38+02:00 (15 years ago)

Author:

cetlod

Message:

distribution of iom_put in TOP routines, see ticket:437

File:

• trunk/NEMO/TOP_SRC/PISCES/p4zlys.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/PISCES/p4zlys.F90

@@ -22 +22 @@
    USE sms_pisces
    USE prtctl_trc
+   USE iom
 
    IMPLICIT NONE
@@ -61 +62 @@
       REAL(wp) ::   zomegaca, zexcess, zexcess0
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zco3
+#if defined key_trc_dia3d && defined key_iomput
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zw3d
+#endif
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
@@ -151 +155 @@
 
 # if defined key_trc_diaadd &&  defined key_trc_dia3d
-      DO jk = 1, jpk
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-               trc3d(ji,jj,jk,jp_pcs0_3d    ) = hi  (ji,jj,jk)          * tmask(ji,jj,jk)
-               trc3d(ji,jj,jk,jp_pcs0_3d + 1) = zco3(ji,jj,jk)          * tmask(ji,jj,jk)
-               trc3d(ji,jj,jk,jp_pcs0_3d + 2) = aksp(ji,jj,jk) / calcon * tmask(ji,jj,jk)
-            ENDDO
-         ENDDO
-      ENDDO
+#  if ! defined key_iomput
+      trc3d(:,:,:,jp_pcs0_3d    ) = hi  (:,:,:)          * tmask(:,:,:)
+      trc3d(:,:,:,jp_pcs0_3d + 1) = zco3(:,:,:)          * tmask(:,:,:)
+      trc3d(:,:,:,jp_pcs0_3d + 2) = aksp(:,:,:) / calcon * tmask(:,:,:)
+#  else
+      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 )
+#  endif
 # endif
       !