pro ebert_curry_radiative_param_longwave, wl, iwp, re, tau, omega, g ; ebert curry ir parameterization constants alphai=[0.0048,0.0054,0.0005,0.0013,0.0030,0.0038,0.0043, $ 0.0003,-0.0013,-0.0020,0.0,0.0028,0.0044,0.0077,0.0074,0.0043, $ -0.0004,0.0072,0.0110,0.0099] gammai=[0.690,0.626,1.279,1.193,1.021,0.949,$ 0.911,1.290,1.452,1.448,1.388,1.099,0.920,0.366,0.603, $ 1.014,1.604,0.926,0.340,0.001] lambdai=[4.1,5.0,6.1,7.0,8.1,9.1,10.0,11.1,12.2,14.1, $ 16.1,18.2,20.0,25.0,31.0,37.0,50.0,75.0,100.0,167.0] ; get into the proper band index=where(abs(wl-lambdai) eq min(abs(wl-lambdai))) ki=alphai[index]+(gammai[index]/re) tau=((iwp)*ki) omega=0.5 g=0.5 return end