Changeset 3481 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.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/p4zsink.F90 (modified) (4 diffs)

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')
       !