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