Changeset 6806 for branches/NERC

Timestamp:

2016-07-14T13:16:48+02:00 (8 years ago)

Author:

jpalmier

Message:

JPALM -- 14-07-2016 -- Bonnes fetes aux Frenchies

-- convert CO2 flux to kg/m2/s

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (2 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -305 +305 @@
       REAL(wp) ::    f_fco2w, f_BetaD, f_rhosw, f_opres, f_insitut, f_pco2atm, f_fco2atm
       REAL(wp) ::    f_schmidtco2, f_kwco2, f_K0, f_co2starair, f_dpco2, f_kwo2
+      !! jpalm 14-07-2016: convert CO2flux diag in kg/m2/s
+      REAL, PARAMETER :: weight_CO2_mol = 0.0440095
+      REAL, PARAMETER :: secs_in_day = 86400.0
+      REAL, PARAMETER :: CO2flux_conv = weight_CO2_mol / secs_in_day
+
       !!
       INTEGER  ::    iters
@@ -4356 +4361 @@
                      ENDIF
                      IF( med_diag%CO2FLUX%dgsave ) THEN
+                         !! convert from  mol/m2/day to kg/m2/s
+                         f_co2flux2d(:,:) = CO2flux_conv * f_co2flux2d(:,:) 
                          CALL lbc_lnk(f_co2flux2d(:,:),'T',1. )
                          CALL iom_put( "CO2FLUX"  , f_co2flux2d )