Is gasesous transmittance at 620nm used elsewhere?

[BaptisteVandecrux]

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**********