| 1 | !------------------------------------------------------------------------ |
|---|
| 2 | ! Based on shr_const_mod.F90 from CAM |
|---|
| 3 | ! Ported to WRF by William.Gustafson@pnl.gov, Nov. 2009 |
|---|
| 4 | !------------------------------------------------------------------------ |
|---|
| 5 | |
|---|
| 6 | !=============================================================================== |
|---|
| 7 | ! SVN $Id: shr_const_mod.F90 6749 2007-10-04 20:58:20Z jwolfe $ |
|---|
| 8 | ! SVN $URL: https://svn-ccsm-models.cgd.ucar.edu/csm_share/branch_tags/cesm1_0_rel_tags/cesm1_0_rel03_share3_100802/shr/shr_const_mod.F90 $ |
|---|
| 9 | !=============================================================================== |
|---|
| 10 | |
|---|
| 11 | MODULE shr_const_mod |
|---|
| 12 | |
|---|
| 13 | use shr_kind_mod |
|---|
| 14 | |
|---|
| 15 | integer(SHR_KIND_IN),parameter,private :: R8 = SHR_KIND_R8 ! rename for local readability only |
|---|
| 16 | |
|---|
| 17 | !---------------------------------------------------------------------------- |
|---|
| 18 | ! physical constants (all data public) |
|---|
| 19 | !---------------------------------------------------------------------------- |
|---|
| 20 | public |
|---|
| 21 | |
|---|
| 22 | real(R8),parameter :: SHR_CONST_PI = 3.14159265358979323846_R8 ! pi |
|---|
| 23 | real(R8),parameter :: SHR_CONST_CDAY = 86400.0_R8 ! sec in calendar day ~ sec |
|---|
| 24 | real(R8),parameter :: SHR_CONST_SDAY = 86164.0_R8 ! sec in siderial day ~ sec |
|---|
| 25 | real(R8),parameter :: SHR_CONST_OMEGA = 2.0_R8*SHR_CONST_PI/SHR_CONST_SDAY ! earth rot ~ rad/sec |
|---|
| 26 | real(R8),parameter :: SHR_CONST_REARTH = 6.37122e6_R8 ! radius of earth ~ m |
|---|
| 27 | real(R8),parameter :: SHR_CONST_G = 9.80616_R8 ! acceleration of gravity ~ m/s^2 |
|---|
| 28 | |
|---|
| 29 | real(R8),parameter :: SHR_CONST_STEBOL = 5.67e-8_R8 ! Stefan-Boltzmann constant ~ W/m^2/K^4 |
|---|
| 30 | real(R8),parameter :: SHR_CONST_BOLTZ = 1.38065e-23_R8 ! Boltzmann's constant ~ J/K/molecule |
|---|
| 31 | real(R8),parameter :: SHR_CONST_AVOGAD = 6.02214e26_R8 ! Avogadro's number ~ molecules/kmole |
|---|
| 32 | real(R8),parameter :: SHR_CONST_RGAS = SHR_CONST_AVOGAD*SHR_CONST_BOLTZ ! Universal gas constant ~ J/K/kmole |
|---|
| 33 | real(R8),parameter :: SHR_CONST_MWDAIR = 28.966_R8 ! molecular weight dry air ~ kg/kmole |
|---|
| 34 | real(R8),parameter :: SHR_CONST_MWWV = 18.016_R8 ! molecular weight water vapor |
|---|
| 35 | real(R8),parameter :: SHR_CONST_RDAIR = SHR_CONST_RGAS/SHR_CONST_MWDAIR ! Dry air gas constant ~ J/K/kg |
|---|
| 36 | real(R8),parameter :: SHR_CONST_RWV = SHR_CONST_RGAS/SHR_CONST_MWWV ! Water vapor gas constant ~ J/K/kg |
|---|
| 37 | real(R8),parameter :: SHR_CONST_ZVIR = (SHR_CONST_RWV/SHR_CONST_RDAIR)-1.0_R8 ! RWV/RDAIR - 1.0 |
|---|
| 38 | real(R8),parameter :: SHR_CONST_KARMAN = 0.4_R8 ! Von Karman constant |
|---|
| 39 | real(R8),parameter :: SHR_CONST_PSTD = 101325.0_R8 ! standard pressure ~ pascals |
|---|
| 40 | real(R8),parameter :: SHR_CONST_PDB = 0.0112372_R8 ! ratio of 13C/12C in Pee Dee Belemnite (C isotope standard) |
|---|
| 41 | |
|---|
| 42 | real(R8),parameter :: SHR_CONST_TKTRIP = 273.16_R8 ! triple point of fresh water ~ K |
|---|
| 43 | real(R8),parameter :: SHR_CONST_TKFRZ = 273.15_R8 ! freezing T of fresh water ~ K |
|---|
| 44 | real(R8),parameter :: SHR_CONST_TKFRZSW = SHR_CONST_TKFRZ - 1.8_R8 ! freezing T of salt water ~ K |
|---|
| 45 | |
|---|
| 46 | real(R8),parameter :: SHR_CONST_RHODAIR = & ! density of dry air at STP ~ kg/m^3 |
|---|
| 47 | SHR_CONST_PSTD/(SHR_CONST_RDAIR*SHR_CONST_TKFRZ) |
|---|
| 48 | real(R8),parameter :: SHR_CONST_RHOFW = 1.000e3_R8 ! density of fresh water ~ kg/m^3 |
|---|
| 49 | real(R8),parameter :: SHR_CONST_RHOSW = 1.026e3_R8 ! density of sea water ~ kg/m^3 |
|---|
| 50 | real(R8),parameter :: SHR_CONST_RHOICE = 0.917e3_R8 ! density of ice ~ kg/m^3 |
|---|
| 51 | real(R8),parameter :: SHR_CONST_CPDAIR = 1.00464e3_R8 ! specific heat of dry air ~ J/kg/K |
|---|
| 52 | real(R8),parameter :: SHR_CONST_CPWV = 1.810e3_R8 ! specific heat of water vap ~ J/kg/K |
|---|
| 53 | real(R8),parameter :: SHR_CONST_CPVIR = (SHR_CONST_CPWV/SHR_CONST_CPDAIR)-1.0_R8 ! CPWV/CPDAIR - 1.0 |
|---|
| 54 | real(R8),parameter :: SHR_CONST_CPFW = 4.188e3_R8 ! specific heat of fresh h2o ~ J/kg/K |
|---|
| 55 | real(R8),parameter :: SHR_CONST_CPSW = 3.996e3_R8 ! specific heat of sea h2o ~ J/kg/K |
|---|
| 56 | real(R8),parameter :: SHR_CONST_CPICE = 2.11727e3_R8 ! specific heat of fresh ice ~ J/kg/K |
|---|
| 57 | real(R8),parameter :: SHR_CONST_LATICE = 3.337e5_R8 ! latent heat of fusion ~ J/kg |
|---|
| 58 | real(R8),parameter :: SHR_CONST_LATVAP = 2.501e6_R8 ! latent heat of evaporation ~ J/kg |
|---|
| 59 | real(R8),parameter :: SHR_CONST_LATSUB = & ! latent heat of sublimation ~ J/kg |
|---|
| 60 | SHR_CONST_LATICE + SHR_CONST_LATVAP |
|---|
| 61 | real(R8),parameter :: SHR_CONST_OCN_REF_SAL = 34.7_R8 ! ocn ref salinity (psu) |
|---|
| 62 | real(R8),parameter :: SHR_CONST_ICE_REF_SAL = 4.0_R8 ! ice ref salinity (psu) |
|---|
| 63 | |
|---|
| 64 | real(R8),parameter :: SHR_CONST_SPVAL = 1.0e30_R8 ! special missing value |
|---|
| 65 | |
|---|
| 66 | END MODULE shr_const_mod |
|---|