Vince and Ken:

I have made an ascii data file containing data from a time-gcm 
run, interpolated to the Bear Lake location. FYI, the run was
a spring equinox (94080), and the history volumes on the ncar
mss are /ROBLE/RGR95/TSESDR1 through TSESDR7. This is the same
run Gordon Shepard has been looking at. f107d=67., f107a=72.,
ctpoten=30., and hempower=3., so this is a solar minimum case.

The file contains ascii data for 11 "frames", one frame per field.
Each frame begins with a 4-line header followed by the "x" values
(0 to 24 ut, nx=25 ) and the "y" values (35 to 120km, ny=86).
There is a 2d array f(nx,ny) immediately following the y values,
with nx varying first. The fields are given in each header; they
are t,u,v,o2,o3,o,oh,ho2,h,w,n2, in that order. T is in deg K,
u, v, and w are in m/s, and all species are number densities.
If you contour these, it is useful to look at log10 of the number
densities with height on the y-axis.

The file is on the machine hao.hao.ucar.edu, in the directory
/pub/foster, with the file name bearlake.dat. It's only .33Mb,
so I did not bother to compress it. Below is a crude fortran
program to help you read it initially. Please let me know when
you have retrieved the file, and if you need any other data or
any help, etc.

(p.s., Gordon wants to talk with Vince (gshep@ncar.ucar.edu))

--Ben

-----------------------------------------------------------------------
Ben Foster		      	High Altitude Observatory (HAO)
foster@ncar.ucar.edu		phone: 303-497-1595  fax: 303-497-1589  
Nat. Center for Atmos. Res.     P.O. Box 3000 Boulder CO 80307 USA
-----------------------------------------------------------------------


      program rdascii
c
      character*40 fname
      character*80 histvol,field,lab,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
      read(lurd,"(a)",end=900) histvol
      read(lurd,"(a)") field
      read(lurd,"(6x,i4)") iframe
      read(lurd,"(a)") lab
      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,"(' ')")
      write(6,"(a)") histvol(1:lenstr(histvol))
      write(6,"(a)") field(1:lenstr(field))
      write(6,"('Frame ',i4)") iframe
      write(6,"(a)") lab(1:lenstr(lab))
      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
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c