the new script calculates gaseous transmittance at 620nm
Is that something we wnt to add as output?
c ***********************STEP 2************************
c calculation of ozone gaseous transmittance at 620nm
c snow reflectance at 620nm in absence of ozone absorption:
DARM=4.*pi * 8.58e-9/ (1.e-3 * 0.620)
reska=exp(-sqrt( DARM* al))
rozone=r0*exp(-u1*u2/r0*sqrt( DARM* al))
c atmospheric spectral characteristics at 620nm:
tauaer=aot*(0.62/0.5)**(-anna)
taumol=0.008735/0.62**(4.08)
okota=height/6000.
if (okota.gt.0.) taumol=taumol*exp(-okota)
if (MOLEC.eq.1) go to 99
taumol=0.0053/0.62**(4.0932)
99 continue
c total aerosol optical thickness at 620nm:
tau=tauaer+taumol
c atmospheric reflectance, total transmittance and spherical albedo
refatm= func1(tau)
tatm= func2(tau)
albatm= func3(tau)
c gaseous transmittance at 620nm:
c T620=(ctoa(7)-refatm)*(1.-reska*albatm)/tatm/rozone
TOAOZ=refatm+tatm*rozone/(1.-reska*albatm)
T620=ctoa(7)/TOAOZ
c write(*,*) TOAOZ,T620,ctoa(7)
c estimation of gaseous transmittance T620 at 620nm is finished
c ************END of STEP 2**********
the new script calculates gaseous transmittance at 620nm Is that something we wnt to add as output?