Changeset 1459

Timestamp:

2009-05-23T10:53:26+02:00 (15 years ago)

Author:

cetlod

Message:

distribution of iom_put in CFC and C14b models, see ticket:437

Location:

trunk/NEMO/TOP_SRC

Files:

• C14b/trcsms_c14b.F90 (modified) (4 diffs)
• CFC/trcsms_cfc.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/C14b/trcsms_c14b.F90

@@ -19 +19 @@
    USE trdmld_trc_oce
    USE trdmld_trc
+   USE iom
 
    IMPLICIT NONE
@@ -125 +126 @@
        zpv  ,        &      !: piston velocity 
        zdemi, ztra
+#if defined key_trc_dia3d  && defined key_iomput
+      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d
+#endif
       !!----------------------------------------------------------------------
 
@@ -266 +270 @@
             qint_c14(ji,jj) = qint_c14(ji,jj) + qtr_c14(ji,jj) * rdt
 
-# if defined key_trc_diaadd
+# if defined key_trc_diaadd && ! defined key_iomput
             ! Save 2D diagnostics
             trc2d(ji,jj,jp_c14b0_2d    ) = qtr_c14 (ji,jj)
@@ -284 +288 @@
 #endif
                tra(ji,jj,jk,jpc14) = tra(ji,jj,jk,jpc14) - ztra / rdt
-
-# if defined key_trc_dia3d
+#if defined key_trc_dia3d
                ! Save 3D diagnostics
+# if ! defined key_iomput
                trc3d(ji,jj,jk,jp_c14b0_3d ) = ztra    !  radioactive decay
+# else 
+               zw3d(ji,jj,jk) = ztra    !  radioactive decay
 # endif
+#endif
             END DO
          END DO
       END DO
 
+#if defined key_trc_diaadd  && defined key_iomput
+      CALL iom_put( "qtrC14b"  , qtr_c14  )
+      CALL iom_put( "qintC14b" , qint_c14 )
+#endif
+#if defined key_trc_dia3d  && defined key_iomput
+      CALL iom_put( "fdecay" , zw3d )
+#endif
       IF( l_trdtrc ) THEN
          CALL trd_mod_trc( tra(:,:,:,jpc14), jpc14, jptrc_trd_sms, kt )   ! save trends

trunk/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

@@ -20 +20 @@
    USE trdmld_trc_oce
    USE trdmld_trc
+   USE iom
 
    IMPLICIT NONE
@@ -165 +166 @@
                ! cumulation of surface flux at each time step
                qint_cfc(ji,jj,jl) = qint_cfc(ji,jj,jl) + qtr_cfc(ji,jj,jl) * rdt
-
-# if defined key_trc_diaadd
-               ! Save diagnostics , just for CFC11
-               js = 1 
-               trc2d(ji,jj,jp_cfc0_2d    ) = qtr_cfc (ji,jj,js)
-               trc2d(ji,jj,jp_cfc0_2d + 1) = qint_cfc(ji,jj,js)
-# endif              
                !                                               !----------------!
             END DO                                             !  end i-j loop  !
@@ -179 +173 @@
       !                                                     !----------------!
 
+#if defined key_trc_diaadd 
+      ! Save diagnostics , just for CFC11
+# if ! defined key_iomput
+      trc2d(:,:,jp_cfc0_2d    ) = qtr_cfc (:,:,1)
+      trc2d(:,:,jp_cfc0_2d + 1) = qint_cfc(:,:,1)
+# else
+      CALL iom_put( "qtrCFC11"  , qtr_cfc (:,:,1) )
+      CALL iom_put( "qintCFC11" , qint_cfc(:,:,1) )
+# endif
+#endif
+
       IF( l_trdtrc ) THEN
           DO jn = jp_cfc0, jp_cfc1