Herb:

I have made a new file for day 15 for you to look at:

http://download.hao.ucar.edu/pub/foster/jackman/protons_d15_vmr.nc

There is also an ncdump file in the same directory.
This file has species densities in vmr (units are "cm3-mr").
(Ne electron density is still in cm-3)

For now, I have saved the earlier files in a cm3/ subdirectory,
but will remove these when we're ready to make all new files.

Note that I have added a new dimension zplev, a new coordinate
variable (zplev(zplev)), and a new 4-d geopotential heights 
variable ZP:

	zplev = 45 ;
	float zplev(zplev) ;
		zplev:long_name = "LOG PRESSURE" ;
		zplev:units = "LN(P0/P)" ;
	float ZP(time, zplev, lat, lon) ;
		ZP:long_name = "GEOPOTENTIAL HEIGHT" ;
		ZP:units = "KM" ;
		ZP:zp_interp = "no" ;
		ZP:type = "SCALAR" ;

The zplev coordinate is the full log-pressure scale that is
used in the model (-17 at the bottom to +5 at the top, by 0.5). 
The new 4-d geopotential height variable ZP is provided globally 
at these pressure coords.

ZP is the geopotential from the model that was used by the 
post-processor to interpolate the other fields to the constant 
height surfaces 40 -> 100 by 2 km. If you plot ZP, you'll see 
that 40 to 100 km corresponds to zplev ln(p0/p) from about -15 
up to about -6.

FYI, the conversion from zplev in ln(p0/p) to millibars is:

  mb = p0*exp(-zplev(k))*1.e-3

where p0 = standard pressure = 5.0e-4 (dynes/cm2, I think).

I also added N2 to the file (calculated as 1-o2(mmr)-o(mmr) before 
conversion to vmr).

Ray is planning to redo the model runs, extended to day 25, 
i.e., 2 days beyond the proton data (he will also include H2
on the histories). If this new file looks ok to you, I will wait 
till Ray has extended the runs, then will make new files for all 
days for both runs.

--Ben