#!/bin/csh


set iRin = 0
while ("$iRin" < 1)
@ iRin++
set rin = `head -$iRin input_rin.list | tail -1`


set incl = 0
while ("$incl" < 4)
@ incl++
set ncl = `head -$incl input_ncl.list  | tail -1`



set idencl = 0
while ("$idencl" < 3)
@ idencl++
set dencl = `head -$idencl input_dencl.list | tail -1`



set idend = 0
while ("$idend" < 3)
@ idend++
set dend = `head -$idend input_dend.list | tail -1`



echo  R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.inp


touch mc.inp


echo idim3 [1=3D,sonst axial]: ibug iblack [1: Qabs=1,Qsca=0]   >>  mc.inp
echo 0 1 0   >>  mc.inp
echo structure: 7 2phase_torus   >>  mc.inp
echo 14   >>  mc.inp
echo number parameter:   >>  mc.inp
echo 9      >>  mc.inp
echo $dencl  >>  mc.inp 
echo $ncl    >>  mc.inp
echo 0.      >>  mc.inp
echo -1      >>  mc.inp
echo $dend   >>  mc.inp
echo 1   >>  mc.inp
echo 0       >>  mc.inp
echo 0  >>  mc.inp
echo 2.   >>  mc.inp
echo iQuelle 2=any LQuelle Tstar distpc nzyk1000   >>  mc.inp
echo 2 3.85d44 10d3 50e6 5000  >>  mc.inp
echo dGW Rinner Router kugelsym   >>  mc.inp
echo -999  $rin -999  1   >>  mc.inp
echo nx ny nz 5 135 135   >>  mc.inp
echo 51 51 51   >>  mc.inp
echo Dust size distribution: lac lec lasi lesi   >>  mc.inp
echo 6 10 7 12   >>  mc.inp
echo Dust type jmrn 0=aCSi 1=fluffy   >>  mc.inp
echo 1   >>  mc.inp
echo anisotropic scattering i_iso 0=gfa 1=HG 2=iso   >>  mc.inp
echo 0   >>  mc.inp
echo Mit jpah jvsg jgas   >>  mc.inp
echo 0 0 0   >>  mc.inp
echo PAH: zcpah hydpah abuPAH 0.1 abuVSG 0   >>  mc.inp
echo 100 0.25 0.125d-20 0.1d-15   >>  mc.inp
echo klump yes=1 nklump 1000 - Zahl der Klumpen   >>  mc.inp
echo 0 1000   >>  mc.inp
echo oldrandom 0=new,1=old Method 0=BW   >>  mc.inp
echo 0 0   >>  mc.inp


mv mc.inp  input/mc.inp

./a.mc >> Mes
./a.spec



mv input/mc.inp Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.inp

mv Mes  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.mes

mv  output/L.Ttau  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.LTtau

mv  output/theta_1.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th1
mv  output/theta_2.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th2
mv  output/theta_3.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th3
mv  output/theta_4.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th4
mv  output/theta_5.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th5
mv  output/theta_6.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th6
mv  output/theta_7.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th7
mv  output/theta_8.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th8
mv  output/theta_9.out  Data/R${rin}_ncl${ncl}_siz1_dencl${dencl}_dend${dend}_nxc0.th9

end
end
end
end
end


