####################################################################################################
# Same as main script "cloud_albedo_sai" (see header of that file), except iout = 22. This option 
# will output radiances and fluxes at every level, used to understand the angular distribution of 
# the radiances that are contributing to the overall irradiances. Only overhead sun with cloud 
# optical depth of 10 (the "idealized" SAI scenario) is calculated.
#
# 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_withRadiances.txt     # name of output file

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=22          `# Output: 22 is broadband radiances and fluxes at every level`
    JAER=3           `# Strat. aerosol model: 3 is fresh volcanic`
    zaer=20,-21      `# Strat. aerosol layer: 20-21 km`
    TAERST=1, 1      `# Strat. aerosol optical depth: 1 (evenly dist. throughout layer)`
    zcloud=1,-2      `# Cloud layer: 1-2 km`
    Tcloud=10, 1     `# Cloud optical depth: 10 (evenly dist. throughout layer)`
    NSTR=20          `# Number of streams: 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=0            `# Solar zenith: 0 is overhead sun`
    SAZA=180         `# Solar azimuth: 180 means forward scattering is at PHI=0`
 /" > INPUT

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