Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

 

Python script to plot the latest Poker Flat data

Matlab script to plot the latest Poker Flat data

IDL script to plot the latest Poker Flat data

Python script to run globalIsprint example

Matlab script to run globalIsprint example

IDL script to run globalIsprint example

 

Step 1: Web exercises

Simple Local Data Access

...

  1. Start at any Madrigal site above (SRI, CEDAR, Millstone, or EISCAT).

  2. Use Full Access Data and then Browse for Individual Madrigal Experiments. Search for all instruments that were running on 2007-03-27. Use the default All Madrigal Sites option. Go back and repeat the search so only Incoherent Scatter Radars are shown.

  3. Select the Poker Flat World Day experiment that was running on 2007-03-27, and do the following with that experiment:

    1. Determine how many data files there are for that experiment, and how they differ.

    2. For one of the files, choose "View description from the catalog and/or header records", and read the summary information stored there. Ask questions if anything is unclear.

    3. Download one of the files is ascii format by using the Download file link, and sticking with the default Simple column-formated ascii option. Remember that this option does not filter the data, and no derived parameters will be included. Open the downloaded file with a text editor to make sure its easy to understand and parse.

    4. Madrigal administrators can add plots and links to each Madrigal experiment. Click on the Geometry Plot link to see the beams that were used. Examine other plots from that experiment.

    5. For the default file for Alternating Code (AC16-30), choose "Print file as ascii (isprint)". This link allows you print both measured and derived parameters. For this file, choose time parameters (year, month, day, hour, min, sec), geographic parameters (elm, azm, gdlat, glon, gdalt), geophysical parameter (kp), and I. S. Radar Basic Parameters (ne, dne, ti, dti, te, dte). Which of these parameters are in the file, and which are derived?

    6. Repeat the above with headers off and missing data replaced with the string NaN.

    7. Save the result in a file using the Save text to file button.

    8. Determine what the Autocorrelation parameter NUMTXAEU means.

    9. Next, we'll try to filter the data. There are some standards filter at the top of the web page, such as elevation or altitude. Just under them are free-form filters that allow filtering using any parameter. First, apply a filter to eliminate elevations under 75 degrees. Look at the resulting data to be sure all data with elevation less than 75 degrees is eliminated.

    10. Add a filter so that only data where te/dte is greater than 10.0 is printed (that is, the error in te is less than one-tenth the measured value).

Run Models
  1. Use the Madrigal page Run Models->Calculate any Madrigal parameter for a given time and range of lat, lon, and alt to calculate the shadow height (SDWHT) and magnetic field vector (BN,BE,BD) 1000 km directly above PFISR (lat 65.130, lon -147.471) at 2007-03-27 12:00:00 UT.

...

  1. Make sure you have python installed. If not, install it from python.org.

  2. To run the plotting test script in the next exercise, run the next two steps. The non-plotting parts of the Madrigal python API will run without those third-party modules.

  3. Install the package numpy from numpy.scipy.org if you don't already have it. This package does numerical arrays and linear algebra.

  4. Install the ploting package matplotlib from matplotlib.sourceforge.net if you don't already have it.

  5. To install the remote python madrigal api from www.openmadrigal.org, you will need to uncompress the file. Then you will need to run python setupMadrigalWeb.py install. If you are installing on Windows, you will need to open a command prompt, and then change directory to where you uncompressed the file. From windows the python command may be something like c:\python25\python.exe setupMadrigalWeb.py install. All the command line scripts (globalIsprint.py, madrigalPColor.py, and madrigalScatter.py) will be put into your PATH and so can be called directly from the command line.

...

  • Download either the python, Matlab, or IDL version of the script to run the globalIsprint example (see top of page).

  • Go to the SRI or the CEDAR Madrigal site, and choose Full data access then Global Madrigal search.

  • Choose the Poker Flat radar and the date range of March 10-20, 2007.

  • Select Show advanced filters.

  • Choose Alternating code for kind of data.

  • Add three filters:

 

Mnemonic

Lower limit

Upper limit

Kp

4.0

leave blank

gdalt

240

260

ne

2e11

leave blank

 

  • Hit the Select parameters button.

  • Choose the parameters year,month,day,hour,min,sec,elm,azm,gdalt,gdlat,glon,kp,ne,dne and hit Done.

  • Choose Continue at the bottom.

  • Select the language you are using, and hit OK.

  • Verify the script runs successfully as is. If there is a problem, compare to the scripts at the top.

  • Modify the script as suggested below. You can either directly edit your script, or regenerate it using the web interface.

  • Instead of using Kp as a filter, filter on Ap3 > 15.

  • Add Te and Ti as parameters

  • Filter also on Te > 1000