Differences between revisions 2 and 18 (spanning 16 versions)
Revision 2 as of 2008-07-18 04:18:12
Size: 6375
Editor: l4-m0
Comment:
Revision 18 as of 2009-01-29 06:19:24
Size: 3188
Editor: l4-m0
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
= Making CR Synoptic Charts = = Making MDI Magnetic Field Synoptic Charts =
== 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.
Line 3: Line 5:
== DESCRIPTION: == == INSTRUCTIONS: To make the magnetic field charts: ==
Line 5: Line 7:
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 (4?) days of data after the official end of the CR.  1. On '''your own machine''' (for speed), cd /synoptic/new_scripts/synop_qual.
Line 7: Line 9:
== INSTRUCTIONS: ==  1. Run new_cr [CR] for the new Carrington Rotation. You can also run new_cr_noscale to highlight mostly one-pixel cosmic ray effects.
  . The magnetograms for the entire CR (plus a few days after) will pop up in DS9 for you to look at.
  . To set the Cosmic Ray bits edit the file CR####.####-#### where the numbers are the CR, and the first and the last day. Indicate in the file the day and image numbers for which the bit must be set. When finished, run the file and save it in the directory CR.BitsSet.Lists NoteOnCosmicRays
  . If any other trouble bits need to be set, run 'udpate' in a '''sonar''' window. It will prompt you for the day and image numbers, and the quality bit(s) you want to set. Quality bit definitions are at http://soi.stanford.edu/production/QUALITY/.
Line 9: Line 14:
=== To make the magnetic field charts: ===  1. Once quality bits are set, run 'make_magsynop.pl [CR] [FIRST] [LAST] on '''n02''' in the directory /synoptic/new_scripts
  . Here CR is the Carrington Rotation and FIRST and LAST are the beginning and end times of the CR. There is an alias called 'carrot' (run on sonar) that will get those times for you:
   . alias carrot "echo \!* > /tmp/cr;echo \!*;~kehcheng/work/synop/carrot < /tmp/cr"
    . sonar.Stanford.EDU: carrot 2067
     . 2067
     . T_FIRST = 2008.02.21_17:34:00_TAI 982690440.315253 5529.73
     . T_ROT = 2008.03.06_09:44:30_TAI 983871870.041679 5543.41
     . T_LAST = 2008.03.20_01:40:33_TAI 985052432.670730 5557.07
     . T_EARTH = 2008.03.06_09:31:46_TAI 983871105.862488
  . So to make the synoptic charts for CR 2067, on n02 run make_magsynop.pl 2067 5529.73 5557.07.
  . make_magsynop.pl puts the new files in /synoptic/carrot/M/[CR]
   . You can ignore the error messages about not being able to read the overview.fits files.
   . Missing magnetograms also generate a complaint
  This takes roughly 15 minutes to run.
Line 11: Line 29:
 1. On your machine, cd /synoptic/new_scripts/synop_qual.  1. Make the gif images for the fits files that were created by make_magsynop.pl.
  . Run /ssw/gen/setup/ssw_idl on your machine.
   . IDL> .run /synoptic/new_scripts/plot_synop
   . IDL> .run /synoptic/new_scripts/synop_gif_loop
  . To make each of the 0, 15 E/W, 30 E/W, 45 E/W, and 60 E/W gif files for both line of sight and radial, do
   . IDL> synop_gif_loop,####,####
  . this calls, for example plot_synop, in='/synoptic/carrot/M/####/synop_Ml_30W.####.fits', /gif, /next, /carrot
   . IDL> exit
Line 13: Line 38:
 2. Run new_cr [CR] for the new Carrington Rotation.
      The magnetograms for the entire CR (plus a few
      days after) will pop up in DS9 for you to look
      at. If any trouble bits need to be set, run
      'udpate' in a sonar window. It will prompt you
      for the day and image numbers, and the quality
      bit(s) you want to set. Quality bit definitions
      are at http://soi.stanford.edu/production/QUALITY/.
 1. Check the images to make sure they look OK. You may need to set additional trouble bits if it looks like a bad image was used, and rerun make_magsynop.pl and remake the gif images. You can check the files online at
  . http://soi.stanford.edu/magnetic/index6.html
Line 22: Line 41:
 3. Once quality bits are set, run 'make_magsynop.pl [CR] [FIRST] [LAST]
      on n02, where CR is the Carrington Rotation and FIRST and LAST
      are the beginning and end times of the CR. I have an alias called
      'carrot' (run on sonar) that will get those times for you:

      alias carrot "echo \!* > /tmp/cr;echo \!*;~kehcheng/work/synop/carrot < /tmp/cr"

      sonar.Stanford.EDU: carrot 2067
      2067
      T_FIRST = 2008.02.21_17:34:00_TAI 982690440.315253 5529.73
      T_ROT = 2008.03.06_09:44:30_TAI 983871870.041679 5543.41
      T_LAST = 2008.03.20_01:40:33_TAI 985052432.670730 5557.07
      T_EARTH = 2008.03.06_09:31:46_TAI 983871105.862488
   
      So to make the synoptic charts for CR 2067, on n02 run
      make_magsynop.pl 2067 5529.73 5557.07.

      make_magsynop.pl puts the new files in /synoptic/carrot/M/[CR]

 4. Make the gif images for the fits files that were created by make_magsynop.pl.
        Run /ssw/gen/setup/ssw_idl on your machine.
        IDL> .run /synoptic/new_scripts/plot_synop
        For each of the 0, 30 E/W and 60 E/W fits files, and for each l and r, run plot_synop:
        IDL> plot_synop, in='/synoptic/carrot/M/2066/synop_Ml_30W.2066.fits', /gif, /next, /carrot
        IDL> exit

 5. Check the images to make sure they look OK. You may need to set additional
      trouble bits if it looks like a bad image was used, and rerun make_magsynop.pl and
      remake the gif images. You can check the files online at
      http://soi.stanford.edu/magnetic/index6.html

 6. Update /synoptic/carrot/latest to reflect the latest CR that's available when you have updated the charts.

=== To make the intensity charts ===

For a CR that had no events that would affect the
intensity (ie, focus changes, tuning changes, spacecraft rolls or proton storms):

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

  2. Run get_data[CR] for the new Carrington Rotation. This will retrieve the lev1.8
      files you need (which may take awhile) and make a mapfile for you to run.

  3. Run the mapfile (pe -v -A map=map.2067, for example). The mapfile does the following:
        i. makes lev1.8 data
       ii. removes limb darkening
      iii. makes a median filter using all of the Ic images without limb darkening for the CR
       iv. makes lev2 data by applying the filter to the images without limb darkening
        v. remaps the corrected images
       vi. creates the synoptic charts

  4. 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
        IDL> exit

  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. Once the synoptic charts are done, make the lev2 GIF images on your machine.



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:

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

  2. Run get_data[CR] for the new Carrington Rotation. This will retrieve the lev1.8
      files you need (which may take awhile) and make a mapfile for you to run.

  3. 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.

  4. Run the mapfile (pe -v -A map=map.2066, for example).

  5. 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
        IDL> exit

  6. 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

  7. Once the synoptic charts are done, make the lev2 GIF imageson your machine.
 1. Update /synoptic/carrot/latest to reflect the latest CR that's available when you have updated the charts.

Making MDI Magnetic Field Synoptic Charts

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: To make the magnetic field charts:

  1. On your own machine (for speed), cd /synoptic/new_scripts/synop_qual.

  2. Run new_cr [CR] for the new Carrington Rotation. You can also run new_cr_noscale to highlight mostly one-pixel cosmic ray effects.
    • The magnetograms for the entire CR (plus a few days after) will pop up in DS9 for you to look at.
    • To set the Cosmic Ray bits edit the file CR####.####-#### where the numbers are the CR, and the first and the last day. Indicate in the file the day and image numbers for which the bit must be set. When finished, run the file and save it in the directory CR.BitsSet.Lists NoteOnCosmicRays

    • If any other trouble bits need to be set, run 'udpate' in a sonar window. It will prompt you for the day and image numbers, and the quality bit(s) you want to set. Quality bit definitions are at http://soi.stanford.edu/production/QUALITY/.

  3. Once quality bits are set, run 'make_magsynop.pl [CR] [FIRST] [LAST] on n02 in the directory /synoptic/new_scripts

    • Here CR is the Carrington Rotation and FIRST and LAST are the beginning and end times of the CR. There is an alias called 'carrot' (run on sonar) that will get those times for you:
      • alias carrot "echo \!* > /tmp/cr;echo \!*;~kehcheng/work/synop/carrot < /tmp/cr"

        • sonar.Stanford.EDU: carrot 2067
          • 2067
          • T_FIRST = 2008.02.21_17:34:00_TAI 982690440.315253 5529.73
          • T_ROT = 2008.03.06_09:44:30_TAI 983871870.041679 5543.41
          • T_LAST = 2008.03.20_01:40:33_TAI 985052432.670730 5557.07
          • T_EARTH = 2008.03.06_09:31:46_TAI 983871105.862488
    • So to make the synoptic charts for CR 2067, on n02 run make_magsynop.pl 2067 5529.73 5557.07.
    • make_magsynop.pl puts the new files in /synoptic/carrot/M/[CR]
      • You can ignore the error messages about not being able to read the overview.fits files.
      • Missing magnetograms also generate a complaint
      This takes roughly 15 minutes to run.
  4. Make the gif images for the fits files that were created by make_magsynop.pl.
    • Run /ssw/gen/setup/ssw_idl on your machine.
      • IDL> .run /synoptic/new_scripts/plot_synop

      • IDL> .run /synoptic/new_scripts/synop_gif_loop

    • To make each of the 0, 15 E/W, 30 E/W, 45 E/W, and 60 E/W gif files for both line of sight and radial, do
      • IDL> synop_gif_loop,####,####

    • this calls, for example plot_synop, in='/synoptic/carrot/M/####/synop_Ml_30W.####.fits', /gif, /next, /carrot
      • IDL> exit

  5. Check the images to make sure they look OK. You may need to set additional trouble bits if it looks like a bad image was used, and rerun make_magsynop.pl and remake the gif images. You can check the files online at
  6. Update /synoptic/carrot/latest to reflect the latest CR that's available when you have updated the charts.

WSO_Wiki: SynopticChartReductions (last edited 2013-05-01 09:59:02 by localhost)