####################################################################################################
# Same as main script "cloud_albedo_sai" (see header of that file), except the number of streams is 
# reduced from 20 to 2. This simulation is used to test how well a 2-stream model can represent the 
# diffusion-brightening effect.
#
# Note: some simulations for these SAI tests are run with a purely scattering albedo. This is NOT 
# the case for this simulation. Make sure the code has been compiled with the original volcanic 
# aerosol model before running this script.
#
####################################################################################################

outDir=../sbdart_output/                     # name of output directory
outFilename=cloud_albedo_sai_2stream.txt     # name of output file

for aod in 0 0.5 1 1.5 2 ; do
for tau in 0 1 3 5 10 20 50 ; do
for sza in 0 20 40 60 80 ; do

echo "
 &INPUT
    idatm=6          `# Profile: 6 is the US62 Standard atmosphere`
    isalb=4          `# Surface type: 4 is lambertian sea water`
    wlinf=0.2        `# Start wavelength: 0.2 microns`
    wlsup=5.0        `# End wavelength: 5.0 microns`
    wlinc=.005       `# Spectral resolution: 0.005 microns`
    iout=11          `# Output: 11 is broadband fluxes at every level`
    JAER=3           `# Strat. aerosol model: 3 is fresh volcanic`
    zaer=20,-21      `# Strat. aerosol layer: 20-21 km`
    TAERST=$aod, 1   `# Strat. aerosol optical depth: $aod (evenly dist. throughout layer)`
    zcloud=1,-2      `# Cloud layer: 1-2 km`
    Tcloud=$tau, 1   `# Cloud optical depth: $tau (evenly dist. throughout layer)`
    NSTR=2           `# Number of streams: 2, 20 is default in radiance calcs`
    NPHI=19          `# Number of view azimuths: 19`
    PHI=0,180        `# View azimuths: every 10 deg from 0 to 180`
    NZEN=19          `# Number of view zeniths: 19`
    UZEN=0,180       `# View zeniths: every 10 deg from 0 to 180, zero is up`
    SZA=$sza         `# Solar zenith: $sza, 0 is overhead sun`
    SAZA=180         `# Solar azimuth: 180 means forward scattering is at PHI=0`
 /" > INPUT

./sbdart >> ${outDir}${outFilename}

done #sza
done #tau
done #aod