c
      subroutine conlat(f,xx,yy,nx,ny,vp,cmin,cmax,cint,
     +  scalefac,iht,icolor,ibox_clabs,gkswid,zmin,zmax,ciu)
c
c Contour latitude slice f(nx,ny), where x-axis is longitude, 
c   y-axis is zp (if iht=0), or height (if iht=1). Add axes and axis labels.
c zmin,zmax,ciu are output.
c
      real f(nx,ny),xx(nx),yy(ny),vp(4)
      integer gkswid
c
      call cpseti('SET',0)
      call cpseti('MAP',0)
      call cpsetr('XC1',xx(1))
      call cpsetr('XCM',xx(nx))
      call cpsetr('YC1',yy(1))
      call cpsetr('YCN',yy(ny))
      call set(vp(1),vp(2),vp(3),vp(4),xx(1),xx(nx),yy(1),yy(ny),1)
      call contour(f,nx,nx,ny,cmin,cmax,cint,
     +  scalefac,icolor,ibox_clabs,gkswid)
      call cpgetr('CIU',ciu)	! contour interval used
      if (cmin.ge.cmax) then
        call cpgetr('ZMN',zmin)	! field min
        call cpgetr('ZMX',zmax)	! field max
      else
        zmin = cmin
        zmax = cmax
      endif
c
c Add axis labels:
c
      chsize = 0.
      if (iht.le.0) then
        call labrect(xx,nx,yy,ny,'LONGITUDE (DEG)','ZP (LN(P0/P))',
     +    chsize)
      else
        call labrect(xx,nx,yy,ny,'LONGITUDE (DEG)','HEIGHT (KM)',
     +    chsize)
      endif
      return
      end