module module_gfs_physcons
  use module_gfs_machine,only:kind_phys





  real(kind=kind_phys),parameter:: con_pi      =3.1415926535897931 
  real(kind=kind_phys),parameter:: con_sqrt2   =1.414214e+0 
  real(kind=kind_phys),parameter:: con_sqrt3   =1.732051e+0 

  real(kind=kind_phys),parameter:: con_rerth   =6.3712e+6 
  real(kind=kind_phys),parameter:: con_g       =9.80665e+0
  real(kind=kind_phys),parameter:: con_omega   =7.2921e-5 
  real(kind=kind_phys),parameter:: con_rd      =2.8705e+2 
  real(kind=kind_phys),parameter:: con_rv      =4.6150e+2 
  real(kind=kind_phys),parameter:: con_cp      =1.0046e+3 
  real(kind=kind_phys),parameter:: con_cv      =7.1760e+2 
  real(kind=kind_phys),parameter:: con_cvap    =1.8460e+3 
  real(kind=kind_phys),parameter:: con_cliq    =4.1855e+3 
  real(kind=kind_phys),parameter:: con_csol    =2.1060e+3 
  real(kind=kind_phys),parameter:: con_hvap    =2.5000e+6 
  real(kind=kind_phys),parameter:: con_hfus    =3.3358e+5 
  real(kind=kind_phys),parameter:: con_psat    =6.1078e+2 
  real(kind=kind_phys),parameter:: con_sbc     =5.6730e-8 
  real(kind=kind_phys),parameter:: con_solr    =1.3533e+3 
  real(kind=kind_phys),parameter:: con_t0c     =2.7315e+2 
  real(kind=kind_phys),parameter:: con_ttp     =2.7316e+2 
  real(kind=kind_phys),parameter:: con_jcal    =4.1855E+0 

  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 module_gfs_physcons