; Read LBH band system parameters: skip='' nlbh = 106 openr, 1, 'lbh.dat' for i=1,3 do readf, 1, skip lbharr=fltarr(6,106) readf,1,lbharr vp=lbharr[0,*] vpd=lbharr[1,*] wlbh=lbharr[2,*]/10. flbh=lbharr[3,*] sguvi=lbharr[4,*] c_o2=lbharr[5,*] close,1 ; Using losintabs to calculate LBH brightness: openr,1,'rt.out' head='' readf,1,head readf,1,jmax,szamod,utmod for ih=0,4 do readf,1,head rta=fltarr(16,80) readf,1,rta close,1 lbhbands=fltarr(nlbh,80) zs=630. ;km zlbh=rta[0,*] ;km re=6370.*1.e5 ;cm rs=re+zs*1.e5 ;cm rlbh=re+zlbh*1.e5 ;cm den=rta[5,*] ;cm-3 for i=0,nlbh-1 do begin eta=rta[15,*]*flbh[i] ;cm-3 s-1 sigma=c_o2(i) ;cm2 losintabs,rs,rlbh,eta,den,sigma,br lbhbands[i,*]=br*1.e-6 ;Rayleighs endfor inds = where(wlbh ge 141.0 and wlbh le 152.8, cnts) indl = where(wlbh ge 167.1 and wlbh le 181.2, cntl) lbhshort = total(lbhbands[inds,*],1) lbhlong = total(lbhbands[indl,*],1) end