Paul:

I have put 3 time-gcm netcdf history files on the hao ftp site

ftp://ftp.hao.ucar.edu/pub/foster/hays

These files contain 3d model results on the time-gcm grid:

 foster   ftp      40610492 Oct  9 09:33 ROBLE.tgcm24d.sdshays.nc
 foster   ftp      40610492 Oct  9 09:33 ROBLE.tgcm24d.seqhays.nc
 foster   ftp      40610492 Oct  9 09:33 ROBLE.tgcm24d.sjshays.nc

I have also put output from the netcdf command "ncdump" for each of 
the files in the same /pub/foster/hays directory (*.ncd). These are 
text files describing the contents of the *.nc files.

sdshays files are for december solstice,
seqhays files are for march equinox,
sjshays files are for june solstice. 

I made some runs of our post-processor, and saved *.cgm and *.ps
files, also on the ftp site:

 foster   ftp      4178880 Oct  9 09:33 sdshays.log.cgm
 foster   ftp      4518720 Oct  9 09:33 seqhays.log.cgm
 foster   ftp      4222080 Oct  9 09:33 sjshays.log.cgm
 foster   ftp      10089549 Oct  9 09:33 sdshays.log.ps
 foster   ftp      10981610 Oct  9 09:33 seqhays.log.ps
 foster   ftp      10209962 Oct  9 09:33 sjshays.log.ps

These are horizontal maps and vertical slices of the fields.
(vertical slices are plotted in log10) The ps files have same plots 
as the cgm files, so you can use either.

The 3-d fields on the histories are:

TN  : neutral temperature (note that bottom boundary of TN is stored
      at the top level)
XO2N: O2 number density
XO3N: O3 number density
XON : O number density
EO200A: O2 (0-0) band emission
XLO : 
O21SIG: O2(1SIG) number density
O2  : O2 mass mixing ratio
O1  : O mass mixing ratio
Z   : geopotential

The O2 emission was calculated as follows:

c Given tn, xo2,xo,xn2 (species in cm-3), dimensioned id, return
c   O2 (0-0) band emission:
c   (from Ian McDade, (mcdade@windic.yorku.ca) Mar 7, 1994 mail
c    communication to Roble)
c
c The nighttime O2 Atmospheric (0-0) band emission rates in the region may be
c similarly modelled using:
c
c              A1 k1 [O]^2 {[N2]+[O2]} [O2]
c Vat=--------------------------------------------------(photons cm-3 sec-1)
c     {A2 + k2a[O2] + k2b[N2] + k2c[O]} {D'[O2] + D"[O]}
c
c where [O], [O2] and [N2] are number densities cm-3
c A1 = 0.079 sec-1
c A2 = 0.083 sec-1
c k1 = 4.7E-33 (300/T)^2  cm6 sec-1
c k2a = 4.0E-17 cm3 sec-1
c k2b = 2.2E-15 cm3 sec-1
c k2c = 8.0E-14 cm3 sec-1
c D' = 6.6 dimensionless
c D" = 19  dimensionless
c
      dimension tn(id),xo2(id),xo(id),xn2(id),fout(id)
      parameter(a1=.079, a2=.083, rk2a=4.e-17, rk2b=2.2e-15,
     +          rk2c=8.e-14, dp=6.6, dpp=19.)
c
      do k=1,id
        rk1 = 4.7e-33*(300./tn(k))**2
        fout(k) = (a1*rk1*xo(k)**2 * (xn2(k)+xo2(k)) * xo2(k)) /
     +    ((a2 + rk2a*xo2(k) + rk2b*xn2(k) + rk2c*xo(k)) *
     +    (dp*xo2(k)+dpp*xo(k)))
      enddo


I'm assuming you have experience reading netcdf files. If you need
any help, let me know. I can also make more runs of the post-
processor -- could make ascii data files of specific slices if you
want. Let me know when you have copied the files so I can free that
disk. Good Luck,

--Ben