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