#tilttest2.job
#0.50-5 m

#tilt-effect for alpha=6 degrees

rm tilttest2.res
rm tilttest2.ttt

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

cd ~/CODE_TO_DONES/SIZEDIST_DATA

~/CODE_TO_DONES/mcnewd_tilt.exe<<EOF > tilttest2.res
'tmodel050.posdat19'
2.0000, 0.			      	       latsun,fiisun
1              				       n_latobs
0.0000,6.000				       testlat,	testfii (see below*)
1000,10,10				nphot,mx,my
'junk'					parfile (use 'junk')
'junk'					outfile ('junk' -> no)
'junk'					epsfile 
'tilttest2.www'				wwwfile 
8					method 8 callisto-spherical-point
-0.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
50,100					nlat,nfii
5					sat_model -zones 
1,1					alpha_choice,lat_choice (see below*)
11
1.
2.
4.
6.
9.
12.
15.
18.
21.
24.
26.74
'tilttest2.ttt'
EOF

#     intepretation of *marked variables:
#     alpha_choice = 0 -> testfii intepreted as dfii
#     alpha_choice = 1 -> testfii intepreted as alpha     
#     lat_choice   = 0 -> testlat intepreted as they are
#     lat_choice   = 1 -> testlat added to bsun


#move results to working directory

rm junk
rm tilttest2.www
mv tilttest2.res ~/CODE_TO_DONES/EXAMPLES
mv tilttest2.ttt ~/CODE_TO_DONES/EXAMPLES