real(kind=kind_phys), parameter physcons::con_jcal =4.1855E+0_kind_phys