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

rm oppotest2.res
rm oppotest2.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 > oppotest2.res
'tmodel020.posdat19'
6., 0.			      	       latsun,fiisun
17.    				       n_latobs
6.,0.  			       lat_obs,alpha_obs	
6.,.0100
6.,.1000
6.,.2500
6.,.5000
6.,.7500
6.,1.000
6.,1.500
6.,2.000
6.,2.500
6.,3.000
6.,3.500
6.,4.000
6.,4.500
6.,5.000
6.,5.500
6.,6.000
5000,10,10				nphot,mx,my
'junk'					parfile (use 'junk')
'junk'					outfile ('junk' -> no)
'junk'					epsfile ('junk' -> no)
'oppotest2.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 oppotest2.res ~/CODE_TO_DONES/EXAMPLES
mv oppotest2.www ~/CODE_TO_DONES/EXAMPLES