New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7156 for branches/NERC – NEMO

Changeset 7156 for branches/NERC


Ignore:
Timestamp:
2016-10-28T13:04:10+02:00 (8 years ago)
Author:
jpalmier
Message:

JPALM -- 28-10-2016 -- correct CO2 flux conversion in MEDUSA

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r7017 r7156  
    305305      REAL(wp) ::    f_fco2w, f_BetaD, f_rhosw, f_opres, f_insitut, f_pco2atm, f_fco2atm 
    306306      REAL(wp) ::    f_schmidtco2, f_kwco2, f_K0, f_co2starair, f_dpco2, f_kwo2 
    307       !! jpalm 14-07-2016: convert CO2flux diag in kg/m2/s 
    308       REAL, PARAMETER :: weight_CO2_mol = 0.0440095 
    309       REAL, PARAMETER :: secs_in_day = 86400.0 
    310       REAL, PARAMETER :: CO2flux_conv = weight_CO2_mol / secs_in_day 
    311  
     307      !! jpalm 14-07-2016: convert CO2flux diag from mmol-C/m2/s to kg-CO2/m2/s 
     308      REAL, PARAMETER :: CO2flux_conv = 44.0095 * 1.e-6   !! CO2 weight per mol (kg/mmol) 
    312309      !! 
    313310      INTEGER  ::    iters 
     
    14991496#  endif 
    15001497                  !! 
     1498                  !! JPALM -- 28-10-2016 -- save co2flux in mmol/m2/s in 2D field 
     1499                  !!                     before it is converted to mmol/m3/d 
     1500                  IF ( lk_iomput ) THEN 
     1501                      IF( med_diag%CO2FLUX%dgsave ) THEN 
     1502                         f_co2flux2d(ji,jj) = f_co2flux 
     1503                      ENDIF 
     1504                  ENDIF 
     1505                  !! 
    15011506                  !! mmol/m2/s -> mmol/m3/d; correct for sea-ice; divide through by layer thickness 
    15021507                  f_co2flux = (1. - fr_i(ji,jj)) * f_co2flux * 86400. / fthk 
     
    15781583                  !! 
    15791584                  !! store 2D outputs 
     1585                  !! JPALM -- 28-10-2016 -- removed co2flux from the list below 
     1586                  !!                     it is now a little bit above, because a  
     1587                  !!                     conversion is done in between. 
    15801588                  IF ( lk_iomput ) THEN 
    15811589                      IF( med_diag%ATM_PCO2%dgsave ) THEN 
     
    15841592                      IF( med_diag%OCN_PCO2%dgsave ) THEN 
    15851593                         f_pco2w2d(ji,jj) = f_pco2w 
    1586                       ENDIF 
    1587                       IF( med_diag%CO2FLUX%dgsave ) THEN 
    1588                          f_co2flux2d(ji,jj) = f_co2flux 
    15891594                      ENDIF 
    15901595                      IF( med_diag%TCO2%dgsave ) THEN 
     
    43714376                     ENDIF 
    43724377                     IF( med_diag%CO2FLUX%dgsave ) THEN 
    4373                          !! convert from  mol/m2/day to kg/m2/s 
     4378                         !! convert from  mmol-C/m2/s to kg-CO2/m2/s 
    43744379                         f_co2flux2d(:,:) = CO2flux_conv * f_co2flux2d(:,:)  
    43754380                         CALL lbc_lnk(f_co2flux2d(:,:),'T',1. ) 
Note: See TracChangeset for help on using the changeset viewer.