Mike:

I have rerun the processor and saved ascii data and
ps files. This time there are 18 frames, as follows:

Model time (day:hr:min)  16: 0: 0 on volume /ROBLE/RGR97/TSMAHB9

Map projections at ut =  0.0  mtime =  16  0  0
UN+VN    low,high,scale =     0.00    0.00    0.00 vmag_len=  0.0000
wrascii: opened file tgcmproc.dat with lu  22
Frame    1: TN     global map: Ht=  80.0: min,max= 1.74E+02,  2.21E+02
Frame    2: TN     global map: Ht= 100.0: min,max= 1.53E+02,  1.92E+02
Frame    3: TN     global map: Ht= 120.0: min,max= 3.26E+02,  4.73E+02
Frame    4: NO     global map: Ht=  80.0: min,max= 1.51E+04,  3.62E+08
Frame    5: NO     global map: Ht= 100.0: min,max= 8.61E+05,  8.47E+07
Frame    6: NO     global map: Ht= 120.0: min,max= 3.10E+05,  3.95E+07
Frame    7: O3     global map: Ht=  80.0: min,max= 1.56E+07,  1.22E+09
Frame    8: O3     global map: Ht= 100.0: min,max= 3.29E+05,  9.48E+07
Frame    9: O3     global map: Ht= 120.0: min,max= 9.94E+00,  1.20E+03

Longitude slices at ut =  0.0  mtime =  16  0  0
Frame   10: TN       lon slice at glon= -180.00 (lat vs ht)
Frame   11: TN       lon slice at glon=    0.00 (lat vs ht)
Frame   12: TN       lon slice at zonal means (lat vs ht)
Frame   13: NO       lon slice at glon= -180.00 (lat vs ht)
Frame   14: NO       lon slice at glon=    0.00 (lat vs ht)
Frame   15: NO       lon slice at zonal means (lat vs ht)
Frame   16: O3       lon slice at glon= -180.00 (lat vs ht)
Frame   17: O3       lon slice at glon=    0.00 (lat vs ht)
Frame   18: O3       lon slice at zonal means (lat vs ht)

There are 2 files on hao.ucar.edu in /pub/foster:

-rw-r--r--   1 foster   ftp       559008 Jul 29 09:16 mahrsi.dat
-rw-r--r--   1 foster   ftp      1575199 Jul 29 09:16 mahrsi.ps

Below is a small f77 program that will read the ascii
data file. There may be some very minor differences between
these plots and the ones I sent earlier, but the numbers in
the ascii file should correspond exactly to these latest plots.

c
      program rdascii
c
c Read ascii data file generated by tgcmproc:
c
      character*80 fname
      character*80 hdr(4),labx,laby
      character*1 ans
      dimension xx(100),yy(100),f(100,100)
      data lurd/9/
c
      do 200 iii=1,1000
      write(6,"(' ')")
      write(6,"('Enter file name (zzz to stop): ',$)")
      read(5,*) fname
      if (fname(1:3).eq.'zzz') stop 'bye'
      open(lurd,file=fname,status='OLD',err=901)
      write(6,"('Ask about printing for every frame in the',
     +  ' file? (y/n): ',$)")
      read(5,*) ans
      iprompt = 0
      if (ans.eq.'y') iprompt = 1
      nfr = 0
c
c Read data from each frame:
c
 100  continue
      do i=1,4
        read(lurd,"(a)",end=900) hdr(i)
      enddo
      read(lurd,"(a)") labx
      read(labx,"(3x,i4)") nx
      if (nx.gt.1) read(lurd,"(6e13.5)") (xx(i),i=1,nx)
      read(lurd,"(a)") laby
      read(laby,"(3x,i4)") ny
      if (ny.gt.1) read(lurd,"(6e13.5)") (yy(i),i=1,ny)
c
      if (nx.gt.1.and.ny.gt.1) then
        read(lurd,"(6e13.5)") ((f(i,ii),i=1,nx),ii=1,ny)
      elseif (nx.eq.1.and.ny.gt.1) then
        read(lurd,"(6e13.5)") (f(1,i),i=1,ny)
      elseif (ny.eq.1.and.nx.gt.1) then
        read(lurd,"(6e13.5)") (f(i,1),i=1,nx)
      endif
      nfr = nfr+1
c
c Conditionally echo to stdout:
c
      write(6,"(' ')")
c     write(6,"(a)") histvol(1:lenstr(histvol))
c     write(6,"(a)") field(1:lenstr(field))
c     write(6,"('Frame ',i4)") iframe
c     write(6,"(a)") lab(1:lenstr(lab))
      do i=1,4
        write(6,"(a)") hdr(i)
      enddo
      write(6,"(a)") labx(1:lenstr(labx))
      if (nx.gt.1) then
        if (nfr.eq.1.or.iprompt.gt.0) then
          write(6,"('Print x-axis values? (y/n/q): ',$)")
          read(5,*) ans
        else
          ans = 'n'
        endif
        if (ans.eq.'y') write(6,"(6e13.5)") (xx(i),i=1,nx)
        if (ans.eq.'q') goto 902
      endif
      write(6,"(a)") laby(1:lenstr(laby))
      if (ny.gt.1) then
        if (nfr.eq.1.or.iprompt.gt.0) then
          write(6,"('Print y-axis values? (y/n/q): ',$)")
          read(5,*) ans
        else
          ans = 'n'
        endif
        if (ans.eq.'y') write(6,"(6e13.5)") (yy(i),i=1,ny)
        if (ans.eq.'q') goto 902
      endif
      if (nfr.eq.1.or.iprompt.gt.0) then
        write(6,"('Print data values? (y/n/q): ',$)")
        read(5,*) ans
      else
        ans = 'n'
      endif
      if (ans.eq.'y') then
        if (nx.gt.1.and.ny.gt.1) then
          write(6,"(6e13.5)") ((f(i,ii),i=1,nx),ii=1,ny)
        elseif (nx.eq.1.and.ny.gt.1) then
          write(6,"(6e13.5)") (f(1,i),i=1,ny)
        elseif (ny.eq.1.and.nx.gt.1) then
          write(6,"(6e13.5)") (f(i,1),i=1,nx)
        endif
      endif
      if (ans.eq.'q') goto 902
      goto 100
c
      stop
 900  continue
      write(6,"('EOF encountered on file ',a,' (lu ',i3,')')")
     +  fname,lurd
 902  close(lurd)
      goto 200
 901  write(6,"('Error opening file ',a)") fname
 200  continue
      stop
      end
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
      function lenstr(str)
      character*(*) str
c
c Return index to last non-blank char in str
c
      length = len(str)
      do i=length,1,-1
        if (str(i:i).ne.' ') then
          lenstr = i
          return
        endif
      enddo
      lenstr = 0
      return
      end
c
      subroutine clearstr(str)
c 
c Set given string to all blanks 
c 
      character*(*) str
      length = len(str)
      do i=1,length
        str(i:i) = ' '
      enddo
      return
      end