pro ebert_curry_radiative_param_shortwave, wl, iwp, re, tau, omega, g ; ebert curry vis parameterization constants ai=[3.448e-3,3.448e-3,3.448e-3,3.448e-3,3.448e-3] bi=[2.431,2.431,2.431,2.431,2.431] ci=[0.00001,0.00011,0.01240,0.03779,0.46658] di=[0.0,1.405e-5,6.867e-4,1.284e-3,2.050e-5] ei=[0.7661,0.7730,0.7865,0.8172,0.9595] fi=[5.851e-4,5.665e-4,7.204e-4,7.463e-4,1.076e-4] wi=[0.4686,0.3577,0.1149,0.0372,0.0216] ec_wl_s=[0.25,0.7,1.3,1.9,2.5] ec_wl_l=[0.7,1.3,1.9,2.5,4.4] ; get into the proper band i=1 while ec_wl_s[i] lt wl and ec_wl_l[i] lt wl do begin i=i+1 endwhile tau=((iwp)*(ai[i]+(bi[i]/re))) omega=(1.-(ci[i]+(di[i]/re))) if max(omega) ge 1. then omega[*]=0.99 if min(omega) le 0.0 then omega[*]=0.01 g=(ei[i]+(fi[i]*re)) if max(g) ge 1. then g[*]=0.99 if min(g) le -1. then g[*]=-0.99 return end