pro check_4_velmap,year,norev=norev,all=all ; ; checks to see if SYNOP atlas images have the OV velocity map ; Up to last 100 checked unless /ALL given ; headcat,cat,y=year n = where(strpos(cat.obs_prog,'SYNOP') ge 0, c) if c gt 0 then begin if keyword_set(norev) then begin f = cat(n).filename endif else begin f = reverse(cat(n).filename) endelse if not keyword_set(all) then begin nf = n_elements(f) f = f(0:(99 < (nf-1))) nf = n_elements(f) endif nf = n_elements(f) for i = 0,nf-1 do begin break_file,f(i),disk,dir,file,ext giff = '/archive/private/data/cds/atlas/'+file+'.jpg' if file_exist(giff) then begin show_jpg,giff,arr print,'Max for: ',giff,' is ',max(arr(450:500,200:250)) if max(arr(450:500,200:250)) lt 100 then begin add_vel_map,file print,'Amending file: ',file endif endif else begin add_vel_map,file print,'Adding new file: ',file endelse wdelete,!d.window endfor endif end