Changeset 1457 for trunk/NEMO/TOP_SRC/PISCES/p4zflx.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/p4zflx.F90 (modified) (7 diffs)

trunk/NEMO/TOP_SRC/PISCES/p4zflx.F90

@@ -23 +23 @@
    USE prtctl_trc
    USE p4zche
+   USE iom
 
    USE lib_mpp
@@ -32 +33 @@
 
    REAL(wp) :: &  ! pre-industrial atmospheric [co2] (ppm)  
-     atcox  = 0.20946 ,   &  !:
-     atcco2 = 278.           !:
-
-   REAL(wp) :: &
-     tco2flx                 !: Total flux of carbon per year
+      atcox  = 0.20946 ,    &  !:
+      atcco2 = 278.            !:
+ 
+   REAL(wp) :: &          
+      xconv  = 0.01/3600.,  &  ! coefficients for conversion       
+      tco2flx                 !: Total flux of carbon per year
 
    !!* Substitution
@@ -58 +60 @@
       INTEGER, INTENT(in) :: kt
       INTEGER  ::   ji, jj, jrorr, nspyr
-      REAL(wp) ::   zttc, zws, zkgwan
+      REAL(wp) ::   ztc, ztc2, ztc3, zws, zkgwan
       REAL(wp) ::   zfld, zflu, zfld16, zflu16, zfact
-      REAL(wp) ::   zph, zah2, zbot, zdic, zalk, zschmitto2, zalka, zschmittco2
+      REAL(wp) ::   zph, zah2, zbot, zdic, zalk, zsch_o2, zalka, zsch_co2
       REAL(wp), DIMENSION(jpi,jpj) ::   zkgco2, zkgo2, zh2co3, ztco2flx
+#if defined key_trc_diaadd && defined key_iomput
+      REAL(wp), DIMENSION(jpi,jpj) ::  zcflx, zoflx, zkg, zdelc
+#endif
       CHARACTER (len=25) :: charout
 
@@ -120 +125 @@
 !CDIR NOVERRCHK
          DO ji = 1, jpi
-
-            zttc = MIN( 35., tn(ji,jj,1) )
-
-      ! Compute the schmidt Number both O2 and CO2
-      ! ------------------------------------------
-
-            zschmittco2 = 2073.1 - 125.62 * zttc + 3.6276 * zttc**2 - 0.043126 * zttc**3
-            zschmitto2  = 1953.4 - 128.0 * zttc + 3.9918 * zttc**2 - 0.050091 * zttc**3
-
+            ztc  = MIN( 35., tn(ji,jj,1) )
+            ztc2 = ztc * ztc
+            ztc3 = ztc * ztc2 
+            ! Compute the schmidt Number both O2 and CO2
+            zsch_co2 = 2073.1 - 125.62 * ztc + 3.6276 * ztc2 - 0.043126 * ztc3
+            zsch_o2  = 1953.4 - 128.0  * ztc + 3.9918 * ztc2 - 0.050091 * ztc3
             !  wind speed 
-            zws  = wndm(ji,jj)
-
-      ! Compute the piston velocity for O2 and CO2
-      ! ------------------------------------------
-
-            zkgwan = (  0.3 * zws * zws    &
-               &   + 2.5 * ( 0.5246 + zttc * ( 0.016256+zttc*0.00049946 ) ) ) &
+            zws  = wndm(ji,jj) * wndm(ji,jj)
+            ! Compute the piston velocity for O2 and CO2
+            zkgwan = 0.3 * zws  + 2.5 * ( 0.5246 + 0.016256 * ztc + 0.00049946  * ztc2 )
 # if defined key_off_degrad
-               &         * facvol(ji,jj,1)      &
-# endif
-               &   / (100. * 3600.)* ( 1.- fr_i(ji,jj) ) * tmask(ji,jj,1)  
-
-            ! COMPUTE GAS EXCHANGE FOR CO2
-            zkgco2(ji,jj) = zkgwan * SQRT( 660./ zschmittco2 )
-            zkgo2(ji,jj)  = zkgwan * SQRT( 660./ zschmitto2 )
-
-         END DO
-      END DO
-
-      ztco2flx(:,:) = 0.
+            zkgwan = zkgwan * xconv * ( 1.- fr_i(ji,jj) ) * tmask(ji,jj,1) * facvol(ji,jj,1)
+#else
+            zkgwan = zkgwan * xconv * ( 1.- fr_i(ji,jj) ) * tmask(ji,jj,1)
+#endif 
+            ! compute gas exchange for CO2 and O2
+            zkgco2(ji,jj) = zkgwan * SQRT( 660./ zsch_co2 )
+            zkgo2 (ji,jj) = zkgwan * SQRT( 660./ zsch_o2 )
+         END DO
+      END DO
+
+#if ! defined key_iomput
       DO jj = 1, jpj
          DO ji = 1, jpi
@@ -167 +164 @@
             tra(ji,jj,1,jpoxy) = tra(ji,jj,1,jpoxy) + ( zfld16 - zflu16 ) / fse3t(ji,jj,1)
 
-# if defined key_trc_diaadd
+# if defined key_trc_diaadd 
             ! Save diagnostics
             trc2d(ji,jj,jp_pcs0_2d    ) = ( zfld - zflu ) * 1000.  * tmask(ji,jj,1)
@@ -177 +174 @@
       END DO
       !
+#else
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            ! Compute CO2 flux for the sea and air
+            zfld = atcco2 * tmask(ji,jj,1) * chemc(ji,jj,1) * zkgco2(ji,jj)
+            zflu = zh2co3(ji,jj) * tmask(ji,jj,1) * zkgco2(ji,jj)
+            tra(ji,jj,1,jpdic) = tra(ji,jj,1,jpdic) + ( zfld - zflu ) / fse3t(ji,jj,1)
+
+            ! compute flux of carbon
+            ztco2flx(ji,jj) = ( zfld - zflu ) * rfact &
+               &             * e1t(ji,jj) * e2t(ji,jj) * tmask(ji,jj,1) * 1000.
+
+            ! Compute O2 flux 
+            zfld16 = atcox * chemc(ji,jj,2) *tmask(ji,jj,1) * zkgo2(ji,jj)
+            zflu16 = trn(ji,jj,1,jpoxy) * tmask(ji,jj,1) * zkgo2(ji,jj)
+            tra(ji,jj,1,jpoxy) = tra(ji,jj,1,jpoxy) + ( zfld16 - zflu16 ) / fse3t(ji,jj,1)
+# if defined key_trc_diaadd 
+            ! Save diagnostics
+            zcflx(ji,jj) = ( zfld - zflu ) * 1000.  * tmask(ji,jj,1)
+            zoflx(ji,jj) = ( zfld16 - zflu16 ) * 1000. * tmask(ji,jj,1)
+            zkg  (ji,jj) = zkgco2(ji,jj) * tmask(ji,jj,1)
+            zdelc(ji,jj) = atcco2 - zh2co3(ji,jj) / ( chemc(ji,jj,1) + rtrn ) * tmask(ji,jj,1)
+# endif
+         END DO
+      END DO
+#endif
 
       ! Total Flux of Carbon
@@ -203 +226 @@
       ENDIF
 
+# if defined key_trc_diaadd && defined key_iomput
+      CALL iom_put( "Cflx", zcflx )
+      CALL iom_put( "Oflx", zoflx )
+      CALL iom_put( "Kg"  , zkg   )
+      CALL iom_put( "DelC", zdelc )
+#endif
 
    END SUBROUTINE p4z_flx