pro sfc_coords, stt=stt, ent=ent if n_elements(stt) eq 0 then stt='28-aug-96 00:00' if n_elements(ent) eq 0 then ent=strmid(systime(0),8,2)+'-'+ $ strmid(systime(0),4,3)+'-'+strmid(systime(0),22,2)+ $ ' '+strmid(systime(0),11,8) print,stt,ent infil=findfile('/ydb/sfc/*') rd_roadmap,infil,roadmap help,roadmap sel=sel_timrange(roadmap,stt,ent,/bet) help,sel if n_elements(sel) ge 1 then roadmap=roadmap(sel) else begin print,'No terminators in the sfc data base in the specified period' return endelse ;; You may or may not have read in the sfc data from /ydb/sfc using yodat. ;get coords of FFI images in full res pixel values: num = n_elements(roadmap) ;num = number of terminator images in index,data coords=sxt_cen(roadmap) print,'*** Setting duff terminators 3 and 52 to filters 0,0 ***' sfc0 = create_struct('coords',fltarr(2,num), $ 'time',lonarr(num),'day',intarr(num), $ 'filta',lonarr(num),'filtb',lonarr(num)) help,coords sfc0.coords = coords sfc0.time = gt_time(roadmap) sfc0.day = gt_day(roadmap) sfc0.filta = gt_filta(roadmap) sfc0.filtb = gt_filtb(roadmap) ;; sfc0.filta(3)=0 ;; sfc0.filtb(3)=0 ;; image 3 is 1,2 filter ;; sfc0.filta(52)=0 ;; sfc0.filtb(52)=0 ;; image 52 is 1,6 filter ; save as a genx file savegen,sfc0,file='sfc_coords.genx' return end