real(kind=kind_phys), parameter physcons::rhowater = 1000._kind_phys