| 1 | module module_gfs_physcons |
|---|
| 2 | use module_gfs_machine,only:kind_phys |
|---|
| 3 | ! Physical constants as set in NMC handbook from Smithsonian tables. |
|---|
| 4 | ! Physical constants are given to 5 places. |
|---|
| 5 | ! 1990/04/30: g and rd are made consistent with NWS usage. |
|---|
| 6 | ! 2001/10/22: g made consistent with SI usage. |
|---|
| 7 | ! Math constants |
|---|
| 8 | real(kind=kind_phys),parameter:: con_pi =3.1415926535897931 ! pi |
|---|
| 9 | real(kind=kind_phys),parameter:: con_sqrt2 =1.414214e+0 ! square root of 2 |
|---|
| 10 | real(kind=kind_phys),parameter:: con_sqrt3 =1.732051e+0 ! square root of 3 |
|---|
| 11 | ! Primary constants |
|---|
| 12 | real(kind=kind_phys),parameter:: con_rerth =6.3712e+6 ! radius of earth (m) |
|---|
| 13 | real(kind=kind_phys),parameter:: con_g =9.80665e+0! gravity (m/s2) |
|---|
| 14 | real(kind=kind_phys),parameter:: con_omega =7.2921e-5 ! ang vel of earth (1/s) |
|---|
| 15 | real(kind=kind_phys),parameter:: con_rd =2.8705e+2 ! gas constant air (J/kg/K) |
|---|
| 16 | real(kind=kind_phys),parameter:: con_rv =4.6150e+2 ! gas constant H2O (J/kg/K) |
|---|
| 17 | real(kind=kind_phys),parameter:: con_cp =1.0046e+3 ! spec heat air @p (J/kg/K) |
|---|
| 18 | real(kind=kind_phys),parameter:: con_cv =7.1760e+2 ! spec heat air @v (J/kg/K) |
|---|
| 19 | real(kind=kind_phys),parameter:: con_cvap =1.8460e+3 ! spec heat H2O gas (J/kg/K) |
|---|
| 20 | real(kind=kind_phys),parameter:: con_cliq =4.1855e+3 ! spec heat H2O liq (J/kg/K) |
|---|
| 21 | real(kind=kind_phys),parameter:: con_csol =2.1060e+3 ! spec heat H2O ice (J/kg/K) |
|---|
| 22 | real(kind=kind_phys),parameter:: con_hvap =2.5000e+6 ! lat heat H2O cond (J/kg) |
|---|
| 23 | real(kind=kind_phys),parameter:: con_hfus =3.3358e+5 ! lat heat H2O fusion (J/kg) |
|---|
| 24 | real(kind=kind_phys),parameter:: con_psat =6.1078e+2 ! pres at H2O 3pt (Pa) |
|---|
| 25 | real(kind=kind_phys),parameter:: con_sbc =5.6730e-8 ! stefan-boltzmann (W/m2/K4) |
|---|
| 26 | real(kind=kind_phys),parameter:: con_solr =1.3533e+3 ! solar constant (W/m2) |
|---|
| 27 | real(kind=kind_phys),parameter:: con_t0c =2.7315e+2 ! temp at 0C (K) |
|---|
| 28 | real(kind=kind_phys),parameter:: con_ttp =2.7316e+2 ! temp at H2O 3pt (K) |
|---|
| 29 | real(kind=kind_phys),parameter:: con_jcal =4.1855E+0 ! JOULES PER CALORIE () |
|---|
| 30 | ! Secondary constants |
|---|
| 31 | real(kind=kind_phys),parameter:: con_rocp =con_rd/con_cp |
|---|
| 32 | real(kind=kind_phys),parameter:: con_cpor =con_cp/con_rd |
|---|
| 33 | real(kind=kind_phys),parameter:: con_rog =con_rd/con_g |
|---|
| 34 | real(kind=kind_phys),parameter:: con_fvirt =con_rv/con_rd-1. |
|---|
| 35 | real(kind=kind_phys),parameter:: con_eps =con_rd/con_rv |
|---|
| 36 | real(kind=kind_phys),parameter:: con_epsm1 =con_rd/con_rv-1. |
|---|
| 37 | real(kind=kind_phys),parameter:: con_dldt =con_cvap-con_cliq |
|---|
| 38 | real(kind=kind_phys),parameter:: con_xpona =-con_dldt/con_rv |
|---|
| 39 | real(kind=kind_phys),parameter:: con_xponb =-con_dldt/con_rv+con_hvap/(con_rv*con_ttp) |
|---|
| 40 | end module module_gfs_physcons |
|---|
