function ratime, value hour = long(abs(value)/15.0) min = long((abs(value)/15.0 - hour)*60.0 + 0.2) return, string(hour, min, format="(i2.2,'h',i2.2,'m')") end function decdeg, value deg = long(value) min = long(abs((value - deg)*60.0)+0.1) return, string(deg, min, format="(i3.2,'!9%!X',i2.2)")+string("222B) end pro plotso set_plot, 'z' loadct, 39 !P.background = 255 SPAWN, 'wc names', wc_string num_points = long(STRMID(wc_string(0),0,8)) num_cols = long(STRMID(wc_string(0),8,8))/num_points data = STRARR(num_cols,num_points) OPENR, unit, 'names', /GET_LUN READF, unit, data CLOSE, unit FREE_LUN, unit cor = read_data('catalog') box00 = read_data('obox211094600.10.0') box01 = read_data('obox211094600.11.0') box02 = read_data('obox211145240.18.1') box03 = read_data('obox211145240.19.0') latvals = 62.5333333333333 + findgen(11) * 0.216666666666667 lonvals = -22.2437770122798 + findgen(11) * 0.5 latlbs = strarr(11); lonlbs = strarr(11); for i=0,10 do begin latlbs[i] = decdeg(latvals[i]) if (lonvals[i] le 0.0) then $ lonlbs[i] = ratime(lonvals[i]) $ else lonlbs[i] = ratime(lonvals[i]-360.0) end map_set, /gnomic, 63.5333333333333, -20, /isotropic, $ limit=[62.5333333333333, -22.2437770122798, 64.5333333333333, -17.7562229877202], $ title='Sky Map Centered at ASM Error Box Intersection', color=0 map_grid, label=1, latalign=0.5, latlab=-18.0062229877202, lonalign=0.5, $ lonlab=62.6416666666667, lats=latvals, lons=lonvals, charsize=0.7, $ lonnames=lonlbs, latnames=latlbs, color=0 for i=0,num_points-1 do begin xyouts, -cor(0,i)+0.5, cor(1,i), data(0,i), color=0 endfor oplot, -cor(0,*), cor(1,*), psym=4, symsize=0.6, color=0 oplot, -box00(0,*), box00(1,*), color=220 oplot, -box01(0,*), box01(1,*), color=220 oplot, -box02(0,*), box02(1,*), color=220 oplot, -box03(0,*), box03(1,*), color=220 filename = 'map.gif' write_gif, filename, tvrd() end