; ======================================================================================== pro rev_diffs,rev_dat,nrev_dat,rev_base,nrev_base,nremove,rev_remove,nadd,rev_add ; ======================================================================================== ; This procedure finds orbit discrepancies between a summary file and data list. ; --------------------------------------------------------------------------------------- ; --------------------------------------------------------------------------------------- ; Orbits that are in the summary files but not in the current data set. To be removed ; --------------------------------------------------------------------------------------- rev_remove=lonarr(nrev_base) nremove=0 for i=0,nrev_base-1 do begin x=where(rev_dat eq rev_base[i],nx) if(nx eq 0) then begin rev_remove[nremove]=rev_base[i] nremove+=1 endif endfor if(nremove gt 0) then rev_remove=rev_remove[0:nremove-1] else rev_remove=[0] ; --------------------------------------------------------------------- ; Orbits that are in the data set but missing from the summary files. ; --------------------------------------------------------------------- rev_add=lonarr(nrev_dat) nadd=0 for i=0,nrev_dat-1 do begin x=where(rev_base eq rev_dat[i],nx) if(nx eq 0) then begin rev_add[nadd]=rev_dat[i] nadd+=1 endif endfor if(nadd gt 0) then rev_add=rev_add[0:nadd-1] else rev_add=[0] return end