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

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