Shengpan:

I have made cgm and ascii data files with O and OH
emission at latitude slices as before, as you requested.

I have put 4 files, 2 each from the equinox and solstice,
runs, in /pub/foster on anon ftp at ftp.hao.ucar.edu, as follows:

rsh hao ls -l \~ftp/pub/foster
total 3728
-rw-r--r--   1 foster   ftp       455040 Apr  5 11:45 shengpan3.eq.cgm
-rw-r--r--   1 foster   ftp       490236 Apr  5 11:45 shengpan3.eq.dat
-rw-r--r--   1 foster   ftp       430560 Apr  5 11:55 shengpan3.sol.cgm
-rw-r--r--   1 foster   ftp       490236 Apr  5 11:55 shengpan3.sol.dat

(As you know, the cgm files are binary graphics metacode, and the
.dat files are the ascii "equivalents"). There are 12 "frames" in each 
file, e.g.:

Latitude slices at ut =  0.00  mtime =  32  0  0
Frame    1: O1       glat=   2.50 min,max= 0.6385E+05 0.5557E+12 iht= 0
Frame    2: O1       glat=   2.50 min,max= 0.3700E+08 0.5552E+12 iht= 1
Frame    3: O1       glat=  17.50 min,max= 0.1970E+06 0.4008E+12 iht= 0
Frame    4: O1       glat=  17.50 min,max= 0.3554E+08 0.3998E+12 iht= 1
Frame    5: O1       glat=  42.50 min,max= 0.1539E+06 0.3781E+12 iht= 0
Frame    6: O1       glat=  42.50 min,max= 0.4556E+08 0.3767E+12 iht= 1
Frame    7: EOH83    glat=   2.50 min,max= 0.2435E-21 0.1210E+04 iht= 0
Frame    8: EOH83    glat=   2.50 min,max= 0.7679E-01 0.5770E+03 iht= 1
Frame    9: EOH83    glat=  17.50 min,max= 0.2241E-21 0.1168E+04 iht= 0
Frame   10: EOH83    glat=  17.50 min,max= 0.9198E-01 0.4656E+03 iht= 1
Frame   11: EOH83    glat=  42.50 min,max= 0.3044E-21 0.1126E+04 iht= 0
Frame   12: EOH83    glat=  42.50 min,max= 0.9515E-01 0.4199E+03 iht= 1

There is one pair of frames per field per latitude, the first with
pressure on y-axis, the second w/ a narrower height scale on y-axis.

BTW, the OH 8-3 emission is a simple calculation as follows:

      parameter (f8=0.29)
      do k=1,id
        pk6n2 = 5.70e-34*(300./tn(k))**2.62
        pk6o2 = 5.96e-34*(300./tn(k))**2.37
        fout(k) = f8*xo(k)*xo2(k)*(pk6n2*xn2(k)+pk6o2*xo2(k)) /
     +            (260.+2.e-11*xo2(k))
      enddo

We also have an "ohrad" code from Sam Makhlouf for OH density at 9 
vibrational levels, and brightness at wavelength bands:

  0  1-0 2-1 3-2 4-3 5-4 6-5 7-6 8-7 9-8
  0   0  2-0 3-1 4-2 5-3 6-4 7-5 8-6 9-7
  0   0   0  3-0 4-1 5-2 6-3 7-4 8-5 9-6
  0   0   0   0  4-0 5-1 6-2 7-3 8-4 9-5
  0   0   0   0   0  5-0 6-1 7-2 8-3 9-4
  0   0   0   0   0   0  6-0 7-1 8-2 9-3

Let me know if you are interested in any of these as well.

--Ben