2 use machine_post
,only:kind_phys
12 real(kind=kind_phys),
parameter:: con_pi =4.0d0*atan(1.0d0)
13 real(kind=kind_phys),
parameter:: con_sqrt2 =1.414214e+0_kind_phys
14 real(kind=kind_phys),
parameter:: con_sqrt3 =1.732051e+0_kind_phys
16 real(kind=kind_phys),
parameter:: con_rerth =6.3712e+6_kind_phys
17 real(kind=kind_phys),
parameter:: con_g =9.80665e+0_kind_phys
18 real(kind=kind_phys),
parameter:: con_omega =7.2921e-5_kind_phys
19 real(kind=kind_phys),
parameter:: con_rd =2.8705e+2_kind_phys
20 real(kind=kind_phys),
parameter:: con_rv =4.6150e+2_kind_phys
21 real(kind=kind_phys),
parameter:: con_cp =1.0046e+3_kind_phys
22 real(kind=kind_phys),
parameter:: con_cv =7.1760e+2_kind_phys
23 real(kind=kind_phys),
parameter:: con_cvap =1.8460e+3_kind_phys
24 real(kind=kind_phys),
parameter:: con_cliq =4.1855e+3_kind_phys
25 real(kind=kind_phys),
parameter:: con_csol =2.1060e+3_kind_phys
26 real(kind=kind_phys),
parameter:: con_hvap =2.5000e+6_kind_phys
27 real(kind=kind_phys),
parameter:: con_hfus =3.3358e+5_kind_phys
28 real(kind=kind_phys),
parameter:: con_psat =6.1078e+2_kind_phys
29 real(kind=kind_phys),
parameter:: con_sbc =5.6730e-8_kind_phys
30 real(kind=kind_phys),
parameter:: con_solr =1.3533e+3_kind_phys
31 real(kind=kind_phys),
parameter:: con_t0c =2.7315e+2_kind_phys
32 real(kind=kind_phys),
parameter:: con_ttp =2.7316e+2_kind_phys
34 real(kind=kind_phys),
parameter:: con_rocp =con_rd/con_cp
35 real(kind=kind_phys),
parameter:: con_cpor =con_cp/con_rd
36 real(kind=kind_phys),
parameter:: con_rog =con_rd/con_g
37 real(kind=kind_phys),
parameter:: con_fvirt =con_rv/con_rd-1.
38 real(kind=kind_phys),
parameter:: con_eps =con_rd/con_rv
39 real(kind=kind_phys),
parameter:: con_epsm1 =con_rd/con_rv-1.
40 real(kind=kind_phys),
parameter:: con_dldt =con_cvap-con_cliq
41 real(kind=kind_phys),
parameter:: con_xpona =-con_dldt/con_rv
42 real(kind=kind_phys),
parameter:: con_xponb =-con_dldt/con_rv+con_hvap/(con_rv*con_ttp)