c
c------------------------------------------------------------------
c Begin file /home/sting/foster/timesdif/pltglb.f
c------------------------------------------------------------------
c
      subroutine pltglb(utc,fld)
c
      include 'gettime.h'
      include 'timesdif.h'
      character*56 toplab,lab
      dimension plt(imx,jmx),viewport(4),fld(imx,kmx,jmx,nplfld)
      data viewport /.12,.90,.2,.70/
c
      xmid = 0.5*(viewport(1)+viewport(2))
      call cpseti('SET',0)
      call cpseti('MAP',1)
      call cpsetr('XC1',gcmlon(1))
      call cpsetr('XCM',gcmlon(imx))
      call cpsetr('YC1',gcmlat(1))
      call cpsetr('YCN',gcmlat(jmx))
      call cpsetr('ILX',xmid)
      call cpsetr('ILY',-.25)
      call cpseti('ILP',0)
      call mkproj('CE',viewport,0.,0.,0.,dum,dum)
c
c Field loop: 
      ixp = 0
      do 200 ip=1,ntimefld
        if (ipltime(ip).le.0) goto 200
        ixp = ixp+1
c
c Selected pressure loop:
        do 100 izp = 1,npls
          ixzp = ixfind(gcmzp,kmx,spls(izp),dzp)
          if (ixzp.le.0) then
            write(6,"('>>> pltglb: bad selected pressure=',f10.3,
     +        ' izp=',i3)") spls(izp),izp
            goto 100
          endif
c
          plt(:,:) = fld(:,ixzp,:,ixp)
          if (icolor.le.0) then
            call contour(plt,imx,imx,jmx,cint(ip),cmin(ip),cmax(ip))
          else
            call conclr(plt,imx,imx,jmx,cint(ip),cmin(ip),cmax(ip))
          endif
          if (icont.gt.0) call maplot
          call labrect(gcmlon,imx,gcmlat,jmx,'LONGITUDE','LATITUDE',.04)
c
c Top label:
          if (wtime(ip).gt.0.) then
            write(toplab,"(a8,' (% DIFFERENCE)  ZP=',f5.1,' UT = ',
     +        f6.2,11x)") timelab_short(ip),spls(izp),utc
          else
            write(toplab,"(a8,' (RAW DIFFERENCE)  ZP=',f5.1,' UT = ',
     +        f6.2,9x)") timelab_short(ip),spls(izp),utc
          endif
          call wrlab(toplab(1:lenstr(toplab)),xmid,viewport(4)+.05,0.)
c
c History vols label:
          write(lab,"(2x,'Difference of TIGCM histories ',a8,
     +      ' minus ',a8,1x)") cpertvol,ccntrvol
          call wrlab(lab(1:lenstr(lab)),xmid,viewport(3)-.1,.012)
          call frame
          iframe = iframe+1
          write(6,"('pltglb diffs: frame ',i4,' field ',a8,
     +      ' zp=',f5.1)") iframe,timelab_short(ip),spls(izp)
c
c End selected pressure loop: izp=1,npls
 100    continue
c
c End field loop: ip=1,ntimefld
 200  continue
c
      return
      end