real(kind=kind_dyn), parameter physcons::con_cp_dyn =1.0046e+3_kind_dyn