real(kind=kind_phys), parameter physcons::rhosnow = 100._kind_phys