see restart-trace for sources of restart variable definitions

Master restart file (unit nrg) and 
primary history buffer restart file (also unit nrg for all but abs/emis data(nrg2))
are always written.

Secondary (nrg2) only written if timestep after restart is not a abs/emis calc time,
i.e. if the primary history file write frequency, NHTFRQ(1), 
is not a multiple of the absorptivity/emissivity calculation frequency, IRADAE 

I count 58 variables here

I did not include global integrals printed out, since I can't 
imagine that they are prognostic variables.  However, they are
read in for restarting the model, so we'll need to include
them in some form other than the state vector.

Table entry definition:
====================================================
restart variable   definition
masterlist name    definition
                   progn(ostic)/diag(nostic)/stat(ic)), confidence (sure/hi/med/low)
                   basic restart or expanded (when restart is not a physics update time)
                   comments 
====================================================

from write_restart_dynamics (nrg) calls  dynamics/eul/restart_dynamics.F90

====================================================
u3(t-1)		
u3(t-2)		
U		Zonal wind component
		prog high
		basic
------------
v3(t-1)		
v3(t-2)		
V		Merdional wind component
		prog high
		basic
------------
t3(t-1)		
t3(t-2)		
T		Temperature
		prog high
		basic
------------
q3(t-1,plndlv*(pcnst+pnats))		
q3(t-2,plndlv*(pcnst+pnats))		
Q		
		prog high
		basic
		pcnst is advected constituents, pnats is non-advected
------------
ps(t-1)		
ps(t-2)		
PS		Surface pressure
		prog high
		basic
------------
vort(t-1)	vorticity
vort(t-2)
XXX		no "vorticity" in master
		prog low
		basic
------------
div(t-1)	divergence
div(t-2)
XXX		no "divergence" in master
		prog low
		basic
------------
dpsl		longitudinal pressure gradient
XXX		NONE
		prog low
		basic
------------
dpsm		meridional pressure gradient
XXX		NONE
		prog low
		basic
------------
dps 		pressure gradient
XXX		NONE
		prog low
		basic
------------
phis		surface geopotential
XXX		NONE
		stat low
		basic
------------
omga		vertical velocity
OMEGA		vertical pressure velocity
		prog low
		basic
------------
lammp   	Longitude coordinates of the trajectory mid-points of the
		parcels that correspond to the global grid points contained
		in the latitude slice being forecasted.  On entry lammp
		is an initial guess.
XXX
                prog med
                basic low
------------
phimp   	Latitude coordinates of the trajectory mid-points of the
		parcels that correspond to the global grid points contained
		in the latitude slice being forecasted.  On entry phimp
		is an initial guess.
XXX
                prog med
                basic low
------------
sigmp   	Hybrid value at the trajectory midpoint for each gridpoint
		in a vertical slice from the global grid.  On entry sigmp is
		an initial guess.
XXX
                prog med
                basic low
------------
qfcst		slt moisture forecast
XXX
		prog med
		basic low
------------
====================================================
from write_restart_physics (nrg) calls  physics/cam1/restart_physics.F90

====================================================
pblht		PBL height
PBLH		Height of planetary boundary layer
		prog med
		basic med                
------------
tpert		temperature perturbation (PBL)
TPERT		PBL plume temperature perturbation
		prog med
		basic med
------------
qrs		shortwave radiative heating rate
QRS		Solar heating rate
		prog med
		basic med
------------
qrl		shortwave radiative heating rate
QRL		Longwave heating rate
		prog med
		basic med
------------
qpert		moisture perturbation (PBL)
QPERT		PBL plume moisture perturbation
		prog med
		basic med
------------
cld		Fractional cloud cover
CLOUD ?		cloud fraction     
CLDTOT ?	total cloud cover
                prog low
                basic med
------------
qcwat(1,1,begchunk,n3tmp)
                Grid average condensed water mixing ratio
CWAT            Grid average condensed water mixing ratio
yyCWAT ?	yyCWAT is passive tracer quanitity
		but yy are not Q,T, or L (q,t,l)
		prog high
		basic med
		cld, qcwat, and tcwat are physics things,
		but have dynamics time levels
------------
tcwat 		?
tcwato		cloud water old temperature
tcwatn		cloud water new temperature
XXX		see qcwat(?)
		prog high
		basic med
------------
lcwat		?
lcwato		cloud liquid water old q
lcwatn		cloud liq. water new q


		prog low
		basic low
------------
Surface flux info from Comsrf (COUP_CSM)
------------
fsnt		Net column abs solar flux at model top
FSNT		Net downward solar flux at top of model
                stat low
                basic med
------------
fsns		surface absorbed solar flux
FSNS            Net downward solar flux at surface
	 	prog low
                basic low
------------
flnt		Net outgoing lw flux at model top
FLNT		Net upward longwave flux at top of model
		prog low
		basic low
------------
flns		Srf longwave cooling (up-down) flux
FLNS		Net upward longwave flux at surface
		prog low
		basic low
------------
These (non-"ice") are components of srfflx_state_2d, 
written to nrg
------------
a[sl]di[rf]{ice}
                albedo: short/long wave, diffuse/direct, {for ice}
XXX		no 'albedo' anywhere!
                prog low
                basic low
------------
tsice		
TSICE		Sea-ice surface Temperature
		stat low
		basic med
------------
lwup		longwave up radiative flux
XXX             
		prog med
		basic low
------------
landfrac	fraction of sfc area covered by land
LANDFRAC  	Fraction of grid-square covered by land
                prog med
                basic med
------------
landm		land/ocean/sea ice flag
LANDM 		Land-Ocean transition mask:(Output as instantaneous field by default)
                prog med
                basic med
------------
sgh		land/ocean/sea ice flag
       yes, that's what's in code; looks like accidental copy of landm descr.
SGH 		Standard deviation of orography
		static low
		basic med
------------
ts		sfc temp
TS 		Surface temperature
                prog med
		basic med
------------
tssub		cam surface/subsurface temperatures
XXX
		prog low
		basic med
------------
sicthk		cam sea-ice thickness (m)
XXX
		prog low
		basic med
------------
snowhland	snow depth (liquid water) ovr lnd
SNOWHLND	Water equivalent snow depth over land
		prog low
		basic med
------------
snowhice	snow depth (liquid water) ovr ice
SNOWHICE	Water equivalent snow depth over sea-ice
		prog low
		basic med
------------
flwds		Down longwave flux at surface
???		Many upward fluxes
------------
sols		Direct solar rad on surface (< 0.7)
SOLS		Downward visible direct solar to surface
		prog low
		basic low
------------
soll		Direct solar rad on surface (>= 0.7)
SOLL		Downward near IR direct solar to surface.
		prog low
		basic low
------------
solsd		Diffuse solar rad on surface (< 0.7)
SOLSD		Downward visible diffuse solar to surface
		prog low
		basic low
------------
solld		Diffuse solar rad on surface (>= 0.7)
SOLLD		Downward near IR diffuse solar to surface
		prog low
		basic low
------------
trefmxav	diagnostic: tref max over the day
TREFMXAV	Daily maximum reference height
                temperature. The output frequency and
                type of the file determines if it's an
                average or minimum/maximum over
                another time frequency. If the output
                file is a monthly average this will be the
                monthly average of the daily minimum temperature.

                diag high
                basic med
------------
trefmnav	diagnostic: tref min over the day
TREFMNAV	Daily minimum reference height
                temperature. The output frequency and
                type of the file determines if it's an
                average or minimum/maximum over
                another time frequency. If the output
                file is a monthly average this will be the
                monthly average of the daily minimum temperature.

                diag high
                basic med
------------
icefrac		fraction of sfc area covered by seaice
ICEFRAC		Fraction of grid-square covered by sea-ice
                prog low
                basic med
------------
abstot_3d	Non-adjacent layer absorptivites
XXX		
		prog low
		expanded high
------------
absnxt_3d	Nearest layer absorptivities
XXX		
		prog low
		expanded high
------------
emstot_3d	Total emissivity
XXX		
		prog low
		expanded high
------------
