c
      program mklevs
      parameter(ncplzp=52,p0t21=5.412154e-4,p0maccm=4.15872e-4)
      real cplzp(ncplzp),pmbt21(ncplzp),pmbmaccm(ncplzp)
c
      do k=1,ncplzp
        cplzp(k) = -20.5+(k-1)*.5
        pmbt21(k) = p0t21*exp(-cplzp(k))*1.e-3
        pmbmaccm(k) = p0maccm*exp(-cplzp(k))*1.e-3
      enddo
      write(6,"('LEV      ZP        PMB_T21      PMB_MACCM')")
      do k=ncplzp,1,-1
        write(6,"(i3,f8.2,2(1pe15.6))") k,cplzp(k),pmbt21(k),pmbmaccm(k)
      enddo
      stop
      end