Downloading SDO data
Here I give my recommendations for downloading SDO and HMI
data. In
most cases users will want to download cutout data rather than full
disk data, so I discuss this first.
Downloading cutouts
From the JSOC (recommended)
Please check out Procedure
for obtaining SDO cutouts for an EIS data-set for how to use the
JSOC to get cutouts.
The method applies generally as long as you have a pointing and
field-of-view size.
The JSOC website for getting the cutouts is the 'Export data'
webpage.
From Lockheed
LMSAL have a web-based
interface for accessing cutouts, but I prefer to use the IDL
interface to the cutout service. A
typical
call is:
IDL> ssw_cutout_service,'8-feb-2011 23:39','9-feb-2011
00:19',cgiquery,queryout,xcen=-350,ycen=-760,ref_time='9-feb-2011
00:06',fovx=120,fovy=140,wave='193',email=[your e-mail
address],max_frames=41
FOVX and FOVX give the field-of-view in arcsec (not pixels).
The e-mail
address is required so that you receive an e-mail when the query has
been processed. Instead of XCEN and YCEN, REF_HELIO can be specified
which gives in the coordinates in longitude and latitude. For example,
REF_HELIO='S33E29'. I think XCEN and YCEN get converted to longitude
and latitude and rounded to 1 degree, as the cutouts are not precisely
centered on the input coordinates.
If you want multiple wavelengths, make sure to increase
MAX_FRAMES
proportionally to the number of wavelengths.
When the query has been processed, the e-mail links to a
webpage
which gives an IDL command that is used to download the cutouts.
Downloading full disk images
VSO
A typical call is:
IDL> list=vso_search('2010-06-28 00:00','2010-06-28
00:10',inst='aia',wave=171,pixels=4096)
where WAVE= is for specifying the wavelength of the AIA
channel (94,
131, 171, 193, 211, 304, 335), and PIXELS= forces the routine to only
find full resolution data. If you want the AIA synoptic data (1024x1024
images), then use pixels=1024. To take only images every 720 seconds
(for example), use sample=720.
To download the data, do:
IDL> chck=vso_get(list,/rice)
Note that the /rice keyword downloads the data in rice-compressed
format. This routine seems to work out where your computer is and then
choose which file server is best to use. If you have problems with the
default server, then try specifying SITE=. Options include 'nso',
'sdac', 'sao', 'nmsu' and various European sites.
JSOC
My method for getting HMI data from JSOC is go to straight to
the exportdata
page
and in the 'record set' box type something like:
hmi.M_45s[2010.07.25_13:00/5m@45s]
which will look for 45s magnetograms from 25-Jul-2010 from
13:00 for
a
5min period.
You can substitute 'M_45s' for:
- 'hmi.M_720s' (12min magnetograms)
- 'hmi.Ic_45s' (45s continuum images)
- 'hmi.Ic_720s' (12min continuum).
- 'aia.lev1_euv_12s'
- 'aia.lev1_uv_24s'
Where to store data
The SDO team do not give any recommendation on how to store
the SDO
data, unfortunately. I set the environment variable $SDO_TOP to point
to the top level data directory (e.g., /data/sdo), and then I create
sub-directories such as:
/data/sdo/aia/level1/2012/08/05
where all AIA data for 5-Aug-2012 are stored. Thus when I use
the
vso_get() routine I will put the data in here.
The cutout data are generally much smaller in size and just
store
these in the working directory of the project that I'm working on.
What processing level do the data have?
AIA data downloaded from the JSOC and VSO are in level-1 format and so
should be
processed with aia_prep to give a uniform format for all filters.
The LMSAL cutout service data are in level-1.5 format, and so do not
need
to be processed with aia_prep.
Page maintained by Dr Peter R Young.
|