Making MDI Synoptic Charts in Intensity

DESCRIPTION:

The magnetic and intensity Carrington Rotation synoptic charts are generated as soon as the definitive MDI data are processed for the period. They require 3 (or 4?) days of data after the nominal end of the CR.

INSTRUCTIONS: For Making the intensity charts

Option A: For a CR that had '''NO events''' that would affect the intensity (ie, exposure change, focus changes, tuning changes, spacecraft roll or proton storms):

1. On SONAR, cd /synoptic/new_scripts/Ic_synop.

2. Run get_data_cr <CR> for the new Carrington Rotation.

   • This will retrieve the lev1.8 files you need (which may take awhile) and make a mapfile, map.[CR] for you to run.

   • It will create a subdirectory CR<CR> with some ancilliary information in it.

3. Run the mapfile pe -v -A map=map.<CR>

   • (pe -v -A map=map.2067, for example). The mapfile does the following:
   1. makes lev1.8 data
   2. removes limb darkening
   3. makes a median filter using all of the Ic images without limb darkening for the CR
   4. makes lev2 data by applying the filter to the images without limb darkening
   5. remaps the corrected images
   6. creates the synoptic charts
4. Make the synoptic gif images for the fits files.
   • Run /ssw/gen/setup/ssw_idl on your machine.

     • IDL> .run /synoptic/new_scripts/plot_Ic_synop

   • Run plot_Ic_synop for both fits files:

     • IDL> plot_Ic_synop, in='/synoptic/carrot/Ic/2066/synop_Ic_N=3.2066.fits', /gif, /next, /carrot

     • IDL> plot_Ic_synop, in='/synoptic/carrot/Ic/2066/synop_Ic_N=5.2066.fits', /gif, /next, /carrot

5. Check the images to make sure they look OK. You may need to set trouble bits if it looks like a bad image was used. You can check the files online at

   • http://soi.stanford.edu/magnetic/index6.html

6. Copy the mapfile map.[CR] file to CR[CR] and then mv the original to /synoptic/carrot/Ic/[CR]
7. Once the synoptic charts are done, make the individual lev2 GIF images.

   1. On SONAR run /synoptic/new_scripts/Ic_synop/make_gifs_cr <CR>

      • This generates two scripts to run on your own machine

   2. On your own machin run CR<CR>/do_dark.<CR>

   3. On your own machin run CR<CR>/do_int.<CR>

Option B: To make the intensity charts for a CR that '''HAD events''' that would affect the intensity (ie, focus changes, tuning changes, spacecraft rolls or proton storms):

The median filter requires ~30 days of data that have no intensity changes due to events listed above. So for example, if there was a proton storm in the middle of the CR with an obvious change in intensity in the images, you would need 2 median filters - one for before the proton storm, and one for after the proton storm - and each filter would require 30 days of data. Note that each of those ~ 30-day chunks also needs to have no events. This isn't always possible to do because of other events...a roll 5 days after the storm, for instance. If you can't get the 30-day chunks of data to make separate median filters, just process the data as if there were no events.

• Check for relevant events at http://mdisas.nascom.nasa.gov/events/events.html.

• If you have 30-day chunks of good data for each median filter you need:
• On sonar, cd /synoptic/new_scripts/Ic_synop.
• Run get_data_cr [CR] for the new Carrington Rotation.
  • This will retrieve the lev1.8 files you need (which may take awhile) and make a preliminary mapfile for you to edit.
• Edit the mapfile so it creates two median filters (med.2066.1.fits and med.2066.2.fits, for example) and uses them on the two sets of input data.
  • An example of an edited mapfile is at /synoptic/carrot/Ic/2066/map.2066.example
  • As this will use additional data, it may take some time for the data to come on line
• Run the edited mapfile (pe -v -A map=map.2066, for example).
• Make the gif images for the fits files.
  • Run /ssw/gen/setup/ssw_idl on your machine.

    • IDL> .run /synoptic/new_scripts/plot_Ic_synop

  • Run plot_Ic_synop for both fits files:

    • IDL> plot_Ic_synop, in='/synoptic/carrot/Ic/2066/synop_Ic_N=3.2066.fits', /gif, /next, /carrot

    • IDL> plot_Ic_synop, in='/synoptic/carrot/Ic/2066/synop_Ic_N=5.2066.fits', /gif, /next, /carrot

• Check the images to make sure they look OK.
  • You may need to set trouble bits if it looks like a bad image was used.

  • You can check the files online at http://soi.stanford.edu/magnetic/index6.html

• Copy the map.[CR] file to CR[CR] and then mv it to /synoptic/carrot/Ic/[CR]
• Once the synoptic charts are done, make the lev2 GIF images on your machine.
  • See instructions above