tau=sza_calc_tau(nu,za) Calculate tau at given freq (Hz) and zenith angle (deg) At the moment 0.035 at 26GHz rising as nu^1.5 which is a guess at an appropriate values for Owen's Valley floor.
0001 function tau=sza_calc_tau(nu,za,taurecipe) 0002 % tau=sza_calc_tau(nu,za) 0003 % 0004 % Calculate tau at given freq (Hz) and zenith angle (deg) 0005 % At the moment 0.035 at 26GHz rising as nu^1.5 which is a guess at an 0006 % appropriate values for Owen's Valley floor. 0007 0008 if(~exist('taurecipe')) 0009 taurecipe='initial_valley_floor'; 0010 end 0011 0012 za=za*pi/180; 0013 0014 switch taurecipe 0015 case 'initial_valley_floor' 0016 tau=0.035*(nu/26e9).^1.5; 0017 end 0018 0019 tau=tau.*sec(za); 0020 0021 return 0022 0023 0024