pro swea_cnts_to_flux, sat, factor factor = 1./5.8e-3 gf = swea_geom_factor(sat,0d,0d) mygf = fltarr(80) for i=0,7 do mygf[i] = gf[2*i,0] + gf[2*i+1,0] mygf[8:23] = gf[*,1] mygf[24:39] = gf[*,2] mygf[40:55] = gf[*,3] mygf[56:71] = gf[*,4] for i=0,7 do mygf[72+i] = gf[2*i,5] + gf[2*i,5] factor = factor/mygf return end