####################################################################################################
# Same as "cloud_albedo_sai_noAbs" (see header of that file), except the cloud optical depths are 
# scaled, via trial-and-error such that the albedo at cloud top is equal to the value from the 
# simulation without SAI. Only overhead sun with cloud optical depth of 10 (the "idealized" SAI 
# scenario) is calculated.
#
# Note: Make sure the code has been recompiled with the overwritten single scattering albedo for 
# the volcanic aerosol model before running this script.
#
####################################################################################################

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


# cloud optical depth was scaled via trial and error, to 3 significant figures, until the best 
# match was found to the albedo at cloud top with the non-SAI simulation.
#
#TARGET ALBEDOS from NON-SAI SIMULATION:
#Cloud albedo without SAI, tau=0: 0.03607167293876154
#Cloud albedo without SAI, tau=1: 0.07627014581734458
#Cloud albedo without SAI, tau=3: 0.16640373197625105
#Cloud albedo without SAI, tau=5: 0.25237170830519917
#Cloud albedo without SAI, tau=10: 0.4148089305125846
#Cloud albedo without SAI, tau=20: 0.5846045665122436
#Cloud albedo without SAI, tau=50: 0.7498691527641479
#
#FINAL ALBEDOS from SCALED COD SAI SIMULATION:
#Cloud albedo with COD-scaled SAI, tau=0: 0.03896264394719596
#Cloud albedo with COD-scaled SAI, tau=0.544: 0.07625821378991209
#Cloud albedo with COD-scaled SAI, tau=2.07: 0.16647064145480867
#Cloud albedo with COD-scaled SAI, tau=3.76: 0.25214691270404543
#Cloud albedo with COD-scaled SAI, tau=8.14: 0.41478915405012406
#Cloud albedo with COD-scaled SAI, tau=16.76: 0.584573637036428
#Cloud albedo with COD-scaled SAI, tau=41.16: 0.7498855112514805

for tau in 0 0.544 2.07 3.76 8.14 16.76 41.16 ; 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: .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=1, 1      `# Strat. aerosol optical depth: 1 (evenly dist. throughout layer)`
    zcloud=1,-2      `# Cloud layer: 1-2 km`
    Tcloud=$tau, 1   `# Cloud optical depth: $tau (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}

done #tau