Changeset 7712


Ignore:
Timestamp:
2017-02-22T12:26:36+01:00 (4 years ago)
Author:
dford
Message:

Update comment relating to pCO2 to fCO2 conversion.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r4650_general_vert_coord_obsoper_surf_bgc/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

    r7625 r7712  
    16811681         pco2_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabm_o3pc) 
    16821682         pco2(:,:) = pco2_3d(:,:,1) 
    1683          ! Now, convert pCO2 to fCO2, based on SST in K 
    1684          ! This follows the method used by SOCAT to convert pCO2 to fCO2, see: 
    1685          ! http://www.socat.info/Uniform_format_surface_fCO2_database.pdf 
     1683         ! Now, convert pCO2 to fCO2, based on SST in K. This follows the standard methodology of: 
     1684         ! Pierrot et al. (2009), Recommendations for autonomous underway pCO2 measuring systems 
     1685         ! and data reduction routines, Deep-Sea Research II, 56: 512-522. 
     1686         ! and 
     1687         ! Weiss (1974), Carbon dioxide in water and seawater: the solubility of a non-ideal gas, 
     1688         ! Marine Chemistry, 2: 203-215. 
     1689         ! In the implementation below, atmospheric pressure has been assumed to be 1 atm and so 
     1690         ! not explicitly included - atmospheric pressure is not necessarily available so this is 
     1691         ! the best assumption. 
     1692         ! Further, the (1-xCO2)^2 term has been neglected. This is common practice 
     1693         ! (see e.g. Zeebe and Wolf-Gladrow (2001), CO2 in Seawater: Equilibrium, Kinetics, Isotopes) 
     1694         ! because xCO2 in atm is ~0, and so this term will only affect the result to the 3rd decimal 
     1695         ! place for typical values, and xCO2 would need to be approximated from pCO2 anyway. 
    16861696         fco2(:,:) = pco2(:,:) * EXP((-1636.75                                                                               + & 
    16871697            &                         12.0408      * (tsn(:,:,1,jp_tem)+rt0)                                                 - & 
Note: See TracChangeset for help on using the changeset viewer.