pro me,wavelen,r,c,radius1,radius2,qabs,qext,qsca

; This program assembles a matrix qabs of emissivities
; as a function of wavelength by using the procedure
; john for all wavelengths

s=size(wavelen)
ilow=fix(radius1/.1d)
ihigh=fix(radius2/.1d)
nradii=ihigh-ilow+1

qabs=dblarr(nradii,s(1))
qext=dblarr(nradii,s(1))
qsca=dblarr(nradii,s(1))
qratio=dblarr(nradii,s(1))

for irad=0,nradii-1 do begin
  radius=0.1d*(ilow+irad)
  for j=0,s(1)-1 do begin
    john,r(j),c(j),radius,wavelen(j),qabs(irad,j),qext(irad,j),qsca(irad,j)
    endfor
endfor
end