#Erstellt aus template.inp verschiedene uvspec input und output files, 
#anschließend wird mit template.py jedes einzelne output file geplottet.

#!/bin/csh

set template=/athome/hube_is/libRadtran_exercises/Lecture_exercises/6_Impact_of_nstr/uvspec_withAerosol_nstrtemplate.inp

set outdir=/athome/hube_is/libRadtran_exercises/Lecture_exercises/6_Impact_of_nstr

foreach nstr ( 04 12 16 20 28 32 )

    echo ... nstr = $nstr

    gawk -vs=$nstr '{gsub("VARIABLE",s); print}' $template > /athome/hube_is/libRadtran_exercises/Lecture_exercises/6_Impact_of_nstr/uvspec_withAerosol_nstr$nstr.inp

    uvspec <uvspec_withAerosol_nstr$nstr.inp >uvspec_withAerosol_nstr$nstr.out


end 

set template=/athome/hube_is/libRadtran_exercises/Lecture_exercises/6_Impact_of_nstr/radfield_polar_nstrtemplate.py

set outdir=/athome/hube_is/libRadtran_exercises/Lecture_exercises/6_Impact_of_nstr

foreach nstr ( 04 12 16 20 28 32 )


    gawk -vs=$nstr '{gsub("VARIABLE",s); print}' $template > /athome/hube_is/libRadtran_exercises/Lecture_exercises/6_Impact_of_nstr/radfield_polar_nstr$nstr.py

    python radfield_polar_nstr$nstr.py


end


exit 0