Changeset 3481 for branches/2012

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

Location:

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES

Files:

• P4Z/p4zflx.F90 (modified) (3 diffs)
• P4Z/p4zprod.F90 (modified) (2 diffs)
• P4Z/p4zsink.F90 (modified) (4 diffs)
• P4Z/p4zsms.F90 (modified) (4 diffs)
• sms_pisces.F90 (modified) (1 diff)

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

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

@@ -52 +52 @@
    REAL(wp) :: texcret                !: 1 - excret 
    REAL(wp) :: texcret2               !: 1 - excret2        
-   REAL(wp) :: tpp                    !: Total primary production
 
 
@@ -409 +408 @@
      ! Total primary production per year
      tpp = tpp + glob_sum( ( zprorca(:,:,:) + zprorcad(:,:,:) ) * cvol(:,:,:) )
-
-     IF( kt == nitend .AND. jnt == nrdttrc ) THEN
-        IF (lwp) THEN
-           WRITE(numout,*) 'Total PP (Gtc) :'
-           WRITE(numout,*) '-------------------- : ',tpp * 12. / 1.E12
-           WRITE(numout,*) 
-        ENDIF
-      ENDIF
 
      IF( ln_diatrc ) THEN

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

@@ -238 +238 @@
       END DO
 
-      !
-      IF( ln_diatrc ) THEN
+     ! Total primary production per year
+     t_oce_co2_exp = t_oce_co2_exp + glob_sum( ( sinking(:,:,iksed+1) + sinking2(:,:,iksed+1) ) * e1e2t(:,:) * tmask(:,:,1) )
+     !
+     IF( ln_diatrc ) THEN
          zrfact2 = 1.e3 * rfact2r
          ik1  = iksed + 1
@@ -273 +275 @@
       !!                  ***  ROUTINE p4z_sink_init  ***
       !!----------------------------------------------------------------------
+
+      t_oce_co2_exp = 0._wp
+      !
    END SUBROUTINE p4z_sink_init
 
@@ -474 +479 @@
       END DO
 
+     ! Total primary production per year
+     t_oce_co2_exp = t_oce_co2_exp + glob_sum( ( sinking(:,:,:) ) * cvol(:,:,:) )
+     !
       IF( ln_diatrc ) THEN
          !
@@ -655 +663 @@
       END DO
       !
+      t_oce_co2_exp = 0._wp
+      !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_sink_init')
       !

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

@@ -42 +42 @@
    INTEGER ::  numsil  !: logical unit for Si budget
    INTEGER ::  numfer  !: logical unit for Fer budget
+   INTEGER ::  numco2  !: logical unit for co2 budget
 
    !!----------------------------------------------------------------------
@@ -147 +148 @@
       END IF
       !
+       !
+       IF( kt == nitend ) THEN
+         !
+         t_atm_co2_flx  = t_atm_co2_flx / glob_sum( e1e2t(:,:) )
+         t_oce_co2_flx  = t_oce_co2_flx         * 12. / 1.e15 * (-1 )
+         tpp            = tpp           * 1000. * 12. / 1.E15
+         t_oce_co2_exp  = t_oce_co2_exp * 1000. * 12. / 1.E15
+         !
+         IF(lwp) WRITE(numco2,9000) ndastp, t_atm_co2_flx, t_oce_co2_flx, tpp, t_oce_co2_exp
+         !
+      ENDIF
+
       IF( nn_timing == 1 )  CALL timing_stop('p4z_sms')
+      !
+ 9000  FORMAT(i8,f10.5,e18.10,f10.5,f10.5)
       !
    END SUBROUTINE p4z_sms
@@ -368 +383 @@
             CALL ctl_opn( numfer, 'iron.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea )
          ENDIF
+         IF(lwp) CALL ctl_opn( numco2, 'carbon.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea )
       ENDIF
 
@@ -401 +417 @@
          ENDIF
        ENDIF
+       !
  9500  FORMAT(i10,e18.10)     
        !

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/sms_pisces.F90

@@ -46 +46 @@
    REAL(wp) ::   xkmort            !: ???
    REAL(wp) ::   ferat3            !: ???
+
+   !!*  diagnostic parameters 
+   REAL(wp) ::  tpp                !: total primary production
+   REAL(wp) ::  t_oce_co2_exp      !: total carbon export
+   REAL(wp) ::  t_oce_co2_flx      !: Total ocean carbon flux
+   REAL(wp) ::  t_atm_co2_flx      !: global mean of atmospheric pco2
 
    !!* Damping