;****************************************
; MC_ASYOBS_READ.saveinfo
; 04/01/01 (06/12/99 HS)
;****************************************

Explanations for the variables stored in the savefile 
produced by mc_asyobs_read.pro


runide		identification for the dynamical simulation
	        (all particle position files have this prefix)

allfile		identification of the photometric simulation
		(name of savfile is allfile+'.save')

pos1,pos2	defines the posdat-files used from dynamical 
	        simulation (runide.pos1 - runide.pos2)

postable	indices of posdat-files 
	        (currently: pos1,pos1+1,  ..pos2)
npos		number of posdat-files

azitable	table of azimuths in photometric simulation
nazi		number of azimuths
		(ccw from sub-observer point, 270 corresponds to Eastern ansa)

bsun,lonsun	Sun elevation and longitude wrt sub-observer point

bobs,lonobs 	Observer elevation and longitude (lonobs=0 always)

nphot		number of photons / single position file / single azimuth

mx,my		number of sub-divisions in search of intersectors
		in the photon path when using spatial partitioning
		(1,1 -> not used)

cpu		total cpu time of the photometric run
cpu_comp        cpu time / 1000 photons/ single position file / single azimuth
	        (this to compare speed of various versions/ cases)

npartori	number of particles in dynamical simulations
npart		number of particles in photometric simulation, includes
		dublicates of particles extending over the boundaries 

npot_mean	mean number of particles for which intersection is actually
		searched for when using sub-divisions (instead of npart)

nbox_mean	mean number of original + image boxes photon has traveled

nscat_mean	mean number of scatterings before escape

tablew		raw values of weights (nazi * npos elements)
tabled	        error estimate of weights (nazi * npos elements)
tablen		number of scatterings contributing to weights 
		(nazi * npos elements)
tablem          raw values of weights for various orders of scatterings
		(nazi * npos elements * 10 elements ; for the last index
		0-8 refers to 1-9 scatterings, 9 refers to the rest)
		
meanw		same as above except that they are means overs
meand		various npos posdat-files
meann
meanm

iperf		normalization factor to convert
		meanw to I/F

ifmin,ifmax	min and max of I/F, assuming albedo=1
amp		(ifmax-ifmin)/ifmin    -"-

ifmax50,ifmin50,amp50,meanw50 	same for albedo 0.5


-----------------------------------------
table_fdirect   fraction of photons passing directly through
		(nazi * npos elements)


method		  phase-function model (1-lambert-element, 8 callisto etc)
smethod		  string variable describing method
ppar1,ppar2,ppar3 parameters of phase-functions


variables related to Saturn-shine
iglobe		  ligh-source model (1-sun, 2-saturn)
sglobe		  string variable describing light-source
		  (SUN, SATURN, SATURN-NORTH, SATURN-SOUTH)
smodel		  string variable describing Saturn-model


flux_sat_azi	  saturn-flux (nazi elements)
nphot_sat_azi	  number of saturn-photons (nazi elements)

abox	          Planetocentric distance
rsaturn		  Planet radius
nlat_sat,nfii_sat number of divisions in Saturn-grid

-----------------------------------------
variables related to geometric correction:

d_per_a		  = observing distance/Saturnocentric distance
	             0  -> infinity
		    >0  -> make geometric correction 

suntable          = Sun longitudes with respect to simulation x-axis
		    in degrees, ccw, nazi-values	 
obstable	  = observer longitudes with respect to simulation x-axis
		    without geometric correction

testfii_tod	  = observer longitudes with respect to simulation x-axis
		    with geometric correction
testlat_tod       = observer latitudes 
		    with geometric correction

alpha		  = phase angle (nazi-values)
tau_n		  = photometric normal optical depth (nazi-values)
dtau_n		  = formal error in tau_n
taumean 	  = mean (tau_n)
tau_dyn		  = dynamical optical depth


----------------------------------------------------------------------
;SOME EXAMPLES

1) To plot brightnes vs azimuth:  plot,azitable,iperf*meanw
This is for albedo=1.

2) For some other value of albedo do like:

alb=0.5
meanw50=meanw*0.

for i=0,8 do begin
meanw50=meanw50+alb^(i+1.)*meanm(*,i)
endfor

plot,azitable,iperf*meanw50

3) Different orders of scattering

plot,azitable,iperf*meanm(*,0)*alb	;single scattered
oplot,azitable,iperf*meanm(*,1)*alb^2	;scattered twice etc.