PRO run_remmean print, 'WARNING: This program MUST be run in the directory in which you' print, 'want your files PRODUCED in. If you are not in this directory,' print, 'break the program, exit IDL and mave to the correct directory.' print, ' ' print, ' ' print, 'Enter extension (.new) :' ext = '' read, ext if strmid(ext,0,1) ne '.' then ext='.'+ext print, ' ' print, 'Enter filename of file which contains all of the mlat and mlon :' mfile = '' read, mfile close,1 openr,1,mfile line = '' readf,1,line nstat = fix(strmid(line,0,3)) readf,1,line stats = strarr(nstat) alat = fltarr(nstat) for i=0,nstat-1 do begin readf,1,line stats(i) = string(byte(strmid(line,4,3))+byte(32)) alat(i) = float(strmid(line,36,6)) endfor close,1 print, 'Enter Quiet day #1, Quiet day #2, starting hour :' quiets = '' read, quiets print, 'Enter Start Active day, End Active day :' starts = '' read, starts print, 'Enter cut off latitude for curve / mean (-1 for manual entry) : ' cutlat = 0.0 read, cutlat list = findfile('-l '+'*'+ext) ndir = 0 nold = 0 nnew = 0 dn = 0 flag = 0 start = 0 jflag = 0 for i=0,strlen(list(0))-1 do $ if (strmid(list(0),i,1) eq '/') or $ (strmid(list(0),i,1) eq ' ') then jflag=i for i=0,n_elements(list)-1 do $ list(i)=strmid(list(0),jflag+1,strlen(list(0))-jflag-1) print, list for i=0,n_elements(list)-1 do begin if (cutlat lt 0.0) then begin print, ' ' print, 'File : ',list(i) print, 'Remove Curve (c) or mean (m) ?' mc = '' read, mc endif else begin if abs(alat(i)) le cutlat then mc = 'c' else mc = 'm' endelse close, 1 openw, 1,'temp.input' printf,1, quiets printf,1, 'i' printf,1, starts printf,1, list(i) close,1 if strmid(mc,0,1) eq 'm' then $ spawn,'/galileo/d/ganglu/mags/bin/remmean '+ $ '< temp.input' $ else $ spawn,'/galileo/d/ganglu/mags/bin/remcurv '+ $ '< temp.input' spawn, '/bin/rm temp.input' endfor return end