module diffusion_h implicit none real*8, parameter :: Pdiff = 1. ! pressure below which diffusion is computed real*8, parameter :: tdiffmin = 5d0 real*8, parameter :: dzres = 2d0 ! grid resolution (km) for diffusion end module diffusion_h