pro ciret2_itercoef_retrieval, dbz_layer, vdq_layer, temp_layer, dbz, vdq, ext_coeff, iwc, lambda, N0

;;; derive the power law coeffecients from the minnis tau regression - file in c:\mace\data\sgpmmcrC1ciret2MaceAvg.c1\iter_tau_analysis.dat

bm_regress=2.2515733      ; very little information in regression
am_regress=0.014726900  ; these are the mean values from the tau iteration compared to minnis tau. units are cgs

;aa_regress=0.41798814+(0.019587721*dbz_layer)+(0.0022381106*temp_layer)-(0.0037566726*vdq_layer)    ; mean is 0.23294600
;if aa_regress lt 0.01 then aa_regress=0.2329
;if aa_regress gt 0.7 then aa_regress=0.2329
aa_regress=0.2329

;ba_regress=1.5542373+( -0.016941364*dbz_layer)+(-0.00095888698*temp_layer)+(0.0023735435*vdq_layer) ; mean is 1.887
;if ba_regress gt 2.0 then ba_regress=1.887
;if ba_regress lt 1.0 then ba_regress=1.887

ba_regress=1.887

if float(temp_layer) ge 253. then pressure=50000.
if float(temp_layer) lt 253. then pressure=40000.
if float(temp_layer) lt 243. then pressure=30000.
if float(temp_layer) lt 233. then pressure=25000.
if float(temp_layer) lt 223. then pressure=20000.

rho_air=1.e-3*(pressure/(287.04*float(temp_layer)))

az=(0.19)*((6./(3.14159*0.92))^2)*((am_regress)^2)
bz=(2.*bm_regress)-6.
ab=(0.04394+0.06049)/2. & bb=(0.970+0.831)/2.

dyn_visc=(1.7e-5)*(10.)     ; the 1000/100 converts from kg/m*s to g/cm*s
kin_visc=dyn_visc/rho_air
g=981.


bv=(bb*(bm_regress+2.-bm_regress))-1.
av=ab*kin_visc*((2.*am_regress*g/(rho_air*(kin_visc^2)*aa_regress))^bb)

Z=(10.^((double(dbz)/10.d)))*1.d-12
Vdq=(double(vdq))


       ext_coeff=ext_coeff_expo_ZV_bvbzbm(Z, Vdq, double(bm_regress), double(bv),$
        double(bz), double(am_regress), double(av), double(az), double(rho_air))

       iwc=iwc_expo_ZV_bvbzbm(Z, Vdq, bm_regress, bv, bz, am_regress, av, az, ba_regress, rho_air)

       lambda=lambda_expo_ZV_bvbzbm(Z, Vdq,bm_regress, bv, bz, am_regress, av, az, ba_regress, rho_air)

       N0=iwc/(am_regress*bm_regress*(gamma(bm_regress))*(lambda^(-(bm_regress+1))))

return
end