#oppotest1.job
#test opposition effect at B=26 degrees
#size distribution 0.20 to 5 meters

rm oppotest1.res
rm oppotest1.www

# go to the directory containing position data
# from size-distribution models

cd ~/CODE_TO_DONES/SIZEDIST_DATA

~/CODE_TO_DONES/mcnewd_oppo.exe<<EOF > oppotest1.res
'tmodel020.posdat19'
26., 0.			      	       latsun,fiisun
17.    				       n_latobs
26.,0.  			       lat_obs,alpha_obs	
26.,.0100
26.,.1000
26.,.2500
26.,.5000
26.,.7500
26.,1.000
26.,1.500
26.,2.000
26.,2.500
26.,3.000
26.,3.500
26.,4.000
26.,4.500
26.,5.000
26.,5.500
26.,6.000
5000,10,10				nphot,mx,my
'junk'					parfile (use 'junk')
'junk'					outfile ('junk' -> no)
'junk'					epsfile ('junk' -> no)
'oppotest1.www'				wwwfile 
4					method 4=HG-spherical-point
-.35,-999.,-999.		        	ppar1,ppar2,ppar3
-1					icheck
-10					idum
1.,1.					sredu,zredu
0.					sizemin
0.,0.					s20,f20
0					iglobe	0-sun 1-saturn 2-south 3-north
128.					distance
100,200					nlat,nfii
3					sat_model 
EOF


#move results to working directory

rm junk
mv oppotest1.res ~/CODE_TO_DONES/EXAMPLES
mv oppotest1.www ~/CODE_TO_DONES/EXAMPLES