Format of "history" synoptic maps. These data are very strange! They do not represent true pB. They should be used for morphology only! MK????.LE Ascii files containing 34 days worth of data. (lower height) MK????.LW Ascii files containing 34 days worth of data. MK????.ME Ascii files containing 34 days worth of data. (middle height) MK????.MW Ascii files containing 34 days worth of data. Data format is: 1496 65 176.64 203.94 1.0 1.5 5.0 0.0 East File 1496 65 190.28 217.58 2.0 1.5 5.0 0.0 West File CR YR CRB CRE L H I ? (see descrip. below) Then there is a short header for each day, followed by 37 data points of pB*1.e+8 going from NORTH to SOUTH in 5 degree increments, starting at 0 degrees north and going to 180 degrees south PA (37 points) The west limb maps are the opposite direction. MK????.LEib.X Ascii files containing a synoptic map. Data format is self explanatory. These files were created from the binary data below. MK????.LEib Binary files containing 30 days worth of data. (lower height) MK????.LWib Binary files containing 30 days worth of data. MK????.MEib Binary files containing 30 days worth of data. (middle height) MK????.MWib Binary files containing 30 days worth of data. These files are created from the MK????.LE (ascii) files. They have been "processed", interpolated, and changed to "binary" data. I don't like this format or this data. The format is set up to match the published Mark I, II and III "history" synoptic maps in Tech Note form, available from Liz Boyd (303) 497-1535 (liz@ncar.ucar.edu). --Alice Lecinski (303) 497-1523 (mlso_data_requests@ncar.ucar.edu) The first 8*32 bits are 8 floats indicating the following: 1496 65 176.64 203.94 1.0 1.5 5.0 0.0 East File 1496 65 190.28 217.58 2.0 1.5 5.0 0.0 West File CR YR CRB CRE L H I ? CR = Carrington Rotation number YR = Year For East File: CRB = First day of CR data (CMday - 6.82) CRE = Last day of CR data (CMday - 6.82 + 27.3) For West File: CRB = First day of CR data (CMday + 6.82) CRE = Last day of CR data (CMday + 6.82 + 27.3) L = Limb (East(1.0) or West(2.0)) H = Height above sun center I = Increment of the data (Every 5 degrees) ? = unknown Note: 1) the height is incorrect! For MARK III data, the data in the heightL directory is from diode 22 (1.35 Rsun) the data in the heightM directory is from diode 58 (1.73 Rsun) Both have been multiplied by wierd factors to make them look like MARK I and II data. For MARK I and II data, the data is ***assumed*** to be at 1.5 Rsun 2) CRB, and CRE are incorrect! Error is +-0.5 days. The next 30 bytes are an ascii key for the 30 days of data. An octal 30 (decimal 24) indicates real data existed for this day. An octal 31 (decimal 25) indicates NO data existed for this day. Periods with less than two day data gaps were interpolated across. This is why the ascii key is useful. The next 2 bytes are filler material and should be 0's The next 37x30x16bits are the 30day "map". Time is arranged backwards, (longitude going forwards) so that the last day is first, and the first day last. This makes the data easier to display in typical "map" fashion. For the map, the first 37 * 16 bits are pB*1.e+8 data as shorts (16 bits) for the last day of the 30 day "map". For an east limb map, the first point is centered at 90 degrees south latitude, then next at 85 deg. south latitude, and so on. For a west limb map, the first point is centered at 90 degrees NORTH, latitude, then next at 85 deg. NORTH latitude, and so on. The next 37 * 16 bits are pB*1.e+8 data as shorts (16 bits) for the second to last day of the 30 day "map". This point (the second day) is close to where a "real" east limb map would end. Each 30day map was created to start at least 1.5 days before the above CRB. (A typical days data was collected at 19:00 UT on the day in question.) A quick and dirty way (+-7 degrees longitude) to make a map is to assume that the map simply ends here, and each next 37 points are the next 13 degrees of longitude. That is, the second series of 37 points covers longitudes 0:12, the third longitudes 13:24, etc. Here is my idl code that reads a map and then displays it: (This is the code 'seehist.pro'. To produce postscript contour maps, use the idl code 'conthist.pro' ) ; March 21, 1992 Alice Lecinski ; ; This idl procedure can be used to read in and display ; ANY MK____.[EW]ib file. ; if ( n_elements(histfilenm) eq 0 ) then begin ;{ histfilenm=' ' ahdr=fltarr(8) ckey=bytarr(30) cpad=bytarr(2) histmap=fltarr(30,37) bigmap=fltarr(360,181) window,xsize=360,ysize=181 endif ;} mini=0.2 & maxi=14.0 imfilenm=' ' print,'Enter the filename which contains the list of MK????.[EW]ib files to see--> ' read,imfilenm close,11 openr,11,imfilenm print,' ' print,'The images listed in file: ',imfilenm,' will now be displayed.' print,' ' print,'You can use the mouse to move backwards and forwards in your image list.' print,'Place the mouse in the map window, then:' print,' To move forwards, use the left most mouse button.' print,' To move backwards, use the middle mouse button.' print,' ' print,' ' while ( not eof(11) ) do begin ;{ readf,11,histfilenm close,1 openr,1,histfilenm readu,1,ahdr readu,1,ckey readu,1,cpad readu,1,histmap close,1 print,histfilenm doy = 1 ; Have to skip doy=0 since the doy=0 vector is before ; the image actually begins ( or map ends ) for i=0,350,13 do begin ;{ One day is ~13 degrees... lat = 0 bigmap(i:i+12,0:2) = histmap(doy,lat) lat = lat + 1 for j=3,181-4,5 do begin ;{ bigmap(i:i+12,j:j+4) = histmap(doy,lat) lat = lat + 1 endfor ;} bigmap(i:i+12,178:180) = histmap(doy,lat) doy = doy +1 endfor ;} bigmap(351:359,178:180) = histmap(doy,lat) if(ahdr(4) eq 2.) then bigmap=rotate(bigmap,7) ; Flip the west limb map. tv,bytscl(bigmap,min=mini,max=maxi) cursor,x,y,/normal endwhile ;} close,11 end