Changeset 1830 for branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zflx.F90

Timestamp:

2010-04-12T15:03:51+02:00 (14 years ago)

Author:

cetlod

Message:

Computation of additional diagnostics for PISCES model ( under CPP key key_diaar5 )

• needed for AR5 outputs (vertical inventories, passive tracers at surface,... )
• new output file with suffix dbio_T

File:

• branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zflx.F90 (modified) (1 diff)

branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zflx.F90

@@ -204 +204 @@
           CALL mpp_sum( t_oce_co2_flx )   ! sum over the global domain
         ENDIF
+        ! Conversion in GtC/yr ; negative for outgoing from ocean
+        t_oce_co2_flx = (-1.) * t_oce_co2_flx  * 12. / 1.e15
+        !
         WRITE(numout,*) ' Atmospheric pCO2    :'
         WRITE(numout,*) '-------------------- : ',kt,'  ',t_atm_co2_flx
         WRITE(numout,*) '(ppm)'
-        WRITE(numout,*) 'Total Flux of Carbon :'
-        WRITE(numout,*) '-------------------- : ',t_oce_co2_flx * 12. / 1e15
-        WRITE(numout,*) '(GtC/an)'
+        WRITE(numout,*) 'Total Flux of Carbon out of the ocean :'
+        WRITE(numout,*) '-------------------- : ',t_oce_co2_flx
+        WRITE(numout,*) '(GtC/yr)'
         t_atm_co2_flx = 0.
         t_oce_co2_flx = 0.
+# if defined key_iomput
+        CALL iom_put( "tatpco2" , t_atm_co2_flx  )
+        CALL iom_put( "tco2flx" , t_oce_co2_flx  )
+#endif
       ENDIF
 #endif