program mkmb
      implicit none
      integer,parameter :: nzp=45
      real,parameter :: zpbot=-17., zptop=5.0, dzp=0.5, p0=5.e-4
      real :: zp,pmb
      integer :: k
!
      do k=1,nzp
        zp = zpbot+(k-1)*dzp
        pmb = p0 * exp(-zp) * 1.e-3
        write(6,"('k=',i3,' zp=',f8.2,' pmb=',e12.4)")
     |    k,zp,pmb
      enddo
      end