c
c------------------------------------------------------------------
c Begin file /home/sting/foster/timegcm/mksrc/pltlb.f
c------------------------------------------------------------------
c
      subroutine pltlb(zz,nlat,nfield,gcmlat,zplb,htlb,fname,ieqnx,
     +  tndown,ntnd,tndzp,iframe)
c
c Line plot lower boundaries (lat on x-axis, field on y-axis):
c
      dimension zz(nlat,nfield),gcmlat(nlat),tndown(nlat,ntnd),
     +  tndzp(ntnd),viewport(4)
      data spval/1.e36/
      character*8 fname(nfield)
      character*56 toplab
      data viewport /.11,.89,.21,.86/
c
      write(6,"(' ')")
      do ip=1,nfield
c
c Tn, ht, and un lower boundaries are from zatmos at zp = -13,
c others are from cira (tng and tndown), or acd2d (species), and
c are at htlb = 60km)
c
        if (ieqnx.gt.0) then
	  if (fname(ip).eq.'TN      '.or.fname(ip).eq.'HT      '.or.
     +        fname(ip).eq.'UN      ') then
            write(toplab,"('LOWER BOUNDARY (ZP= ',f5.1,
     +        ') FOR FIELD ',a8,' (EQUINOX)',1x)") zplb,fname(ip)
          else
            write(toplab,"('LOWER BOUNDARY (',f6.2,' KM) FOR FIELD ',a8,
     +        ' (EQUINOX)',1x)") htlb,fname(ip)
          endif
        else
	  if (fname(ip).eq.'TN      '.or.fname(ip).eq.'HT      '.or.
     +        fname(ip).eq.'UN      ') then
            write(toplab,"('LOWER BOUNDARY (ZP= ',f5.1,
     +        ') FOR FIELD ',a8,' (SOLSTICE)')") zplb,fname(ip)
          else
            write(toplab,"('LOWER BOUNDARY (',f6.2,' KM) FOR FIELD ',a8,
     +        ' (SOLSTICE)')") htlb,fname(ip)
          endif
        endif
        call displa(1,1,1)
        call anotat('LATITUDE (GCM GRID)',fname(ip),0,1,-1,L)
        call ezxy(gcmlat,zz(1,ip),nlat,toplab)
        iframe = iframe+1
        write(6,"('pltlb frame ',i3,': ip=',i3,' line plot of ',
     +    'lower boundary field ',a8)") iframe,ip,fname(ip)
      enddo
c
c tndown:
      call cpseti('SET',0)
      call cpseti('MAP',0)
      call cpsetr('XC1',gcmlat(1))
      call cpsetr('XCM',gcmlat(nlat))
      call cpsetr('YC1',tndzp(1))
      call cpsetr('YCN',tndzp(ntnd))
      xmid = 0.5*(viewport(1)+viewport(2))
      call cpsetr('ILX',xmid)
      call cpsetr('ILY',-.20)
      call set(viewport(1),viewport(2),viewport(3),viewport(4),
     +         gcmlat(1),gcmlat(nlat),tndzp(1),tndzp(ntnd),1)
c
c Contour tn:
      call contour(tndown,nlat,nlat,ntnd)
      call labrect(gcmlat,nlat,tndzp,ntnd,'LATITUDE (GCM GRID)',
     +  '-LN(P/P0)',0.)
      if (ieqnx.gt.0) then
        write(toplab,"('CIRA TN (DEG K) (EQUINOX)',31x)")
      else
        write(toplab,"('CIRA TN (DEG K) (SOLSTICE)',30x)")
      endif
      call wrlab(toplab(1:lenstr(toplab)),xmid,0.9,0.)
      call frame
      iframe = iframe+1
      write(6,"('pltlb frame ',i3,' contour tndown')") iframe
c
      return
      end