Changeset 3481 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90

Timestamp:

2012-10-02T14:38:42+02:00 (12 years ago)

Author:

cetlod

Message:

NEMOGCM_dev_r3438_LOCEAN15_PISLOB: output a text file for carbon budget

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90 (modified) (3 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90

@@ -57 +57 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: satmco2   !: atmospheric pco2 
 
-   REAL(wp) ::  t_oce_co2_flx               !: Total ocean carbon flux 
-   REAL(wp) ::  t_atm_co2_flx               !: global mean of atmospheric pco2
-   REAL(wp) ::  area                        !: ocean surface
    REAL(wp) ::  xconv  = 0.01_wp / 3600._wp !: coefficients for conversion 
 
@@ -200 +197 @@
 
       t_oce_co2_flx = t_oce_co2_flx + glob_sum( oce_co2(:,:) )            ! Cumulative Total Flux of Carbon
-      IF( kt == nitend ) THEN
-         t_atm_co2_flx = glob_sum( satmco2(:,:) * patm(:,:) * e1e2t(:,:) )            ! Total atmospheric pCO2
-         !
-         t_oce_co2_flx = (-1.) * t_oce_co2_flx  * 12. / 1.e15             ! Conversion in PgC ; negative for out of the ocean
-         t_atm_co2_flx = t_atm_co2_flx  / area                            ! global mean of atmospheric pCO2
-         !
-         IF( lwp) THEN
-            WRITE(numout,*)
-            WRITE(numout,*) ' Global mean of atmospheric pCO2 (ppm) at it= ', kt, ' date= ', ndastp
-            WRITE(numout,*) '------------------------------------------------------- :  ',t_atm_co2_flx
-            WRITE(numout,*)
-            WRITE(numout,*) ' Cumulative total Flux of Carbon out of the ocean (PgC) :'
-            WRITE(numout,*) '-------------------------------------------------------  ',t_oce_co2_flx
-         ENDIF
-         !
-      ENDIF
+      t_atm_co2_flx = glob_sum( satmco2(:,:) * patm(:,:) * e1e2t(:,:) )   ! Total atmospheric pCO2
 
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
@@ -299 +281 @@
       ENDIF
       !
-      area = glob_sum( e1e2t(:,:) )        ! interior global domain surface
-      !
       oce_co2(:,:)  = 0._wp                ! Initialization of Flux of Carbon
       t_atm_co2_flx = 0._wp