Roger and Mark:

I have remade plots and saved data from Ray's TIEGCM run
/ROBLE/RGR95/TE7ASP1 (and TE7ASP2, TE7ASP3). These are the
same history volumes used in figure 9 of the paper according 
to my files of last August 15th.

There are 2 files on anonymous ftp on the machine
hao.hao.ucar.edu in the directory /pub/foster, as follows:

-rw-r--r--   1 foster   staff     437004 Oct 28 11:13 te7asp.dat
-rw-r--r--   1 foster   staff     940131 Oct 28 11:12 te7asp.ps 

te7asp.ps is a postscript file with 12 frames, and te7asp.dat
is an ascii file with the corresponding data. Postscript frames 
1, 4, 7, and 10 correspond to figure 9 of last August). Note that 
the postscript plots are south polar projections (perimeter lat -40), 
whereas the data are at the global geographic grid (73 lons x 36 
lats). The 12 frames are as follows:

1. tn at zp +2, ut 0	(ps plot also includes wind vectors)
2. un at zp +2, ut 0
3. vn at zp +2, ut 0

4. tn at zp +2, ut 6	(ps plot also includes wind vectors)
5. un at zp +2, ut 6
6. vn at zp +2, ut 6

7. tn at zp +2, ut 12	(ps plot also includes wind vectors)
8. un at zp +2, ut 12
9. vn at zp +2, ut 12

10. tn at zp +2, ut 18	(ps plot also includes wind vectors)
11. un at zp +2, ut 18
12. vn at zp +2, ut 18

Each "frame" in the data file is a 4-line header, followed by
x-axis label, x-axis coordinates, y-axis label, y-axis coordinates,
and finally the 2-d data itself (73x36) (lons rotate 1st). Below 
is a simple IDL procedure that will read the data file, 1 frame at 
a time. It returns a "frame structure" that can be passed to a 
plotting procedure.

BTW, I think the caption for fig 9 might say "at an average height
of approximately 346 km" rather than "at a height of 350 km" since
the plots are actually at a constant pressure surface. The average
heights are shown at the top of each polar plot.

Please let me know if you need anything more.

--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
-----------------------------------------------------------------------


function rdascii,lu,flnm,frame
;
; Read data from ascii file flnm, and return a frame structure:
;
; Four line header, to be used later as follows: 
; line 1: top label
; line 2: x-axis label
; line 3: y-axis label
; line 4: bottom label
;
on_ioerror,bad
hdr = strarr(4)
line = ''
line = string(line,format="(a80)")
nx=1 & ny=1
for i=0,3 do begin
  readf,lu,line
  hdr(i) = strtrim(line,2)
endfor
;
; nx, xx:
;
readf,lu,line
nx = 0
reads,line,nx,format="(3x,i4)"
xx = fltarr(nx)
readf,lu,xx
;
; ny, yy:
;
readf,lu,line
ny = 0
reads,line,ny,format="(3x,i4)"
yy = fltarr(ny)
readf,lu,yy
;
; Now the data:
;
zz = fltarr(nx,ny)
readf,lu,zz
;
frame = {		$
  hdr:hdr,		$
  nx:nx,		$
  xx:xx,		$
  xlab:hdr(1),		$
  ny:ny,		$
  yy:yy,		$
  ylab:hdr(2),		$
  zz:zz,		$
  toplab:hdr(0),	$
  botlab:hdr(3)		}

;
return,0
bad:
  print,'>>> rdascii: i/o error reading from file ',flnm
  return,1
end