real(kind=kind_dyn), parameter physcons::con_hfus_dyn =3.3358e+5_kind_dyn