UPP/physcons_8f_source.html

module physcons_post

  use machine_post,only:kind_phys

!- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

       implicit none

!

!  Physical constants as set in NMC handbook from Smithsonian tables.

!  Physical constants are given to 5 places.

!  1990/04/30: g and rd are made consistent with NWS usage.

!  2001/10/22: g made consistent with SI usage.

!  Math constants

!  real(kind=kind_phys),parameter:: con_pi      =3.141593e+0          !< pi

  real(kind=kind_phys),parameter:: con_pi     =4.0d0*atan(1.0d0)

  real(kind=kind_phys),parameter:: con_sqrt2   =1.414214e+0_kind_phys

  real(kind=kind_phys),parameter:: con_sqrt3   =1.732051e+0_kind_phys

  real(kind=kind_phys),parameter:: con_rerth   =6.3712e+6_kind_phys

  real(kind=kind_phys),parameter:: con_g       =9.80665e+0_kind_phys

  real(kind=kind_phys),parameter:: con_omega   =7.2921e-5_kind_phys

  real(kind=kind_phys),parameter:: con_rd      =2.8705e+2_kind_phys

  real(kind=kind_phys),parameter:: con_rv      =4.6150e+2_kind_phys

  real(kind=kind_phys),parameter:: con_cp      =1.0046e+3_kind_phys

  real(kind=kind_phys),parameter:: con_cv      =7.1760e+2_kind_phys

  real(kind=kind_phys),parameter:: con_cvap    =1.8460e+3_kind_phys

  real(kind=kind_phys),parameter:: con_cliq    =4.1855e+3_kind_phys

  real(kind=kind_phys),parameter:: con_csol    =2.1060e+3_kind_phys

  real(kind=kind_phys),parameter:: con_hvap    =2.5000e+6_kind_phys

  real(kind=kind_phys),parameter:: con_hfus    =3.3358e+5_kind_phys

  real(kind=kind_phys),parameter:: con_psat    =6.1078e+2_kind_phys

  real(kind=kind_phys),parameter:: con_sbc     =5.6730e-8_kind_phys

  real(kind=kind_phys),parameter:: con_solr    =1.3533e+3_kind_phys

  real(kind=kind_phys),parameter:: con_t0c     =2.7315e+2_kind_phys

  real(kind=kind_phys),parameter:: con_ttp     =2.7316e+2_kind_phys


  real(kind=kind_phys), parameter :: con_rocp = con_rd / con_cp

  real(kind=kind_phys), parameter :: con_cpor = con_cp / con_rd

  real(kind=kind_phys), parameter :: con_rog = con_rd / con_g

  real(kind=kind_phys), parameter :: con_fvirt = con_rv / con_rd - 1.

  real(kind=kind_phys), parameter :: con_eps = con_rd / con_rv

  real(kind=kind_phys), parameter :: con_epsm1 = con_rd / con_rv - 1.

  real(kind=kind_phys), parameter :: con_dldt = con_cvap - con_cliq

  real(kind=kind_phys), parameter :: con_xpona = con_dldt / con_rv

  real(kind=kind_phys), parameter :: con_xponb = -con_dldt / con_rv + con_hvap / (con_rv * con_ttp)

end module