function euvac37, f107, f107a common csized, d nbins=37 ; read solar minimum reference spectrum (first time only) sized=size(d) if sized[0] eq 0 then begin print,'Reading euvac_37.dat' str=' ' d=fltarr(4,nbins) close,1 openr,1,'/home/stans/sun/euvac/lowres/euvac_37.dat' for i=1,2 do readf,1,str readf,1,d close,1 endif waves=reform(d[0,*]) wavel=reform(d[1,*]) sfmin=reform(d[2,*]) afac=reform(d[3,*]) ; EUVAC formula pind=0.5*(f107+f107a) sflux=sfmin*(1.+afac*(pind-80.)) sflux=sflux>0.1*sfmin ; return array model=fltarr(3,nbins) model[0,*]=waves model[1,*]=wavel model[2,*]=sflux return, model end