[Previous]
[Next]
NAME:
edac_summary
PURPOSE:
To plot the EDAC error rates
SAMPLE CALLING SEQUENCE:
dirs = file_list('/md84/fits/sci5k/lev0/','2000*', /cd)
edac_summary, file_list(dirs, '9605*'), img, daytim
edac_summary, '9*_000000.fits', img, daytim ;one file a day
edac_summary, ff, img, daytim, print_list=[2,3]
edac_summary, xx, img, daytim, log_list='$MDI_CAL_INFO/edac_summary.txt'
edac_summary,sttim='1-feb',entim='2-feb'
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], RFITS2 [1], RFITS2 [2], SXPAR [1]
SXPAR [2], SXPAR [3], UNSIGN, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4]
UTPLOT [5], UTPLOT [6], anytim2ints [1], anytim2ints [2], break_file [4]
fid2ex [1], fid2ex [2], file_list [1], file_list [2], fmt_tim [1], fmt_tim [2]
plottime [1], plottime [2], pprint [1], pprint [2], sel_filetimes [1]
sel_filetimes [2], ut_time [1], ut_time [2]
OPTIONAL KEYWORD INPUT:
print_list - the list of pages to print info about.
0=page 16, 1=page 14, 2=page 10, 3=page 12
HISTORY:
Written Mar-96 by M.Morrison
V2.0 30-May-96 (MDM) - Turned into a procedure
- Only plot sums
13-Jun-96 (MDM) - Changed to accept /ynozero type of plotting
16-Oct-96 (MDM) - Used /data13 instead of /small1
V2.1 21-Oct-96 (MDM) - Reordered the pages to put page 10 on the
bottom
V2.2 6-Jan-97 (MDM) - Unsign the count
- Added unsign
- Added "log_list" option
V2.3 9-Jan-97 (MDM) - Reversed the "log_list" labels for single/double
error counts
V2.4 4-Feb-97 (MDM) - Added STTIM/ENTIM option
V2.5 26-Mar-97 (MDM) - Corrected "log_list" output (not quite right order)
V2.6 7-Jan-98 (RIB) - changed dirs to 200* from 2000* to pickup 20020000 dpc
[Previous]
[Next]
Project:
SOHO - SUMER
Name:
edges
Purpose:
Apply distortion correction to sumer data
Use:
edg_dat = edges(data,slit,RefPix,detector)
Inputs
Data - data to be analyzied - 1 full frame
RefWav - Refernce wavelength.
Keyword Input:
Quiet - supress the general information printed out at start of program
Output
edg_dat - spectrum edges
Restrictions:
This program should NOT be used on data which has been binned or
subject to compression schemes 6 and up.
Written:
Tom Moran, Sept. 22, 1997
[Previous]
[Next]
Project : SOHO -LASCO/EIT
Name : EDITFRAME
Purpose : This tool allows the user to edit a particular image "on the
fly," while the WRUNMOVIE3 procedure is running.
Use : EDITFRAME, in, out, [moviev]
Arguments : in - The input, or original image
out - The output, or final image
moviev - a structure generated by the WRUNMOVIE3 program,
containing information about the movie. If moviev
is not provided, then certain features of the program,
such as referencing images by "Current frame", will
not be supported.
Calls : ***
EDITHELP, EDIT_GETFILENAME, EDIT_LOADIMG, LASCO_READFITS [1], LASCO_READFITS [2]
LZ_GETLASCODIR, OFFSET_BIAS, QL_GETLASCODIR, SWAP [1], SWAP [2], point_filter
Side effects : None
Category : Image Processing
Written : Scott Hawley, NRL Jul 15, 1996
Version :
See Also : WRUNMOVIE3.PRO
[Previous]
[Next]
NAME:
ef2time
PURPOSE: convert input EIT QKL or LZ filename into corresponding time.
CALLING SEQUENCE:
times = ef2time(filenames,format='ccsds')
INPUTS:
filenames = QKL or LZ file names, e.g. 'efz200001.123456'
may be an array
OPTIONAL INPUTS:
None
KEYWORDS
format = choice of output format, default = UTC
allowed values : 'ecs','ccsds','vms','external','stime'
CAVEATS: The file names do not always correspond to the
"CORRECTED OBS_TIME" in the header. The times may be up to
5 minutes off. Do not use this if higher tie precision is needed.
MODIFICATION HISTORY
17-Jan-2001 - J.S.Newmark
CALLS:
[Previous]
[Next]
NAME:
eit_bakeouts
PURPOSE
returns a list of bakeouts
INPUTS
none
OUTPUTS
a string array containing the start and stop dates for each bakeout
KEYWORDS
none
WARNINGS
needs to be regularly updated.
may not be up to date for the latest bakeout
CREATION
18-May-2001. F. Auchere
MODIFICATIONS
Updated : 11-Dec-2001. F. Auchere
Updated : 14-Jan-2002. F. Auchere
Updated : 19-Jul-2002. F. Auchere
Updated : 04-Sep-2002. F. Auchere
Updated : 24-Dec-2002. K. Schenk
Updated : 04-Feb-2003. S. McIntosh
Updated : 04-Feb-2004: F. Auchere
Updated : 13-jul-2004: F. Auchere
Updated : 08-mar-2005: K. Schenk
Updated : 30-mar-2005: K. Schenk
Updated : 29-mar-2006: K. Schenk
Updated : 12-jun-2006: K. Schenk
Updated : 28-Aug-2006: K. Schenk
Updated : 26-Oct-2006: K. Schenk
CALLED BY
eit_findcalgroup, eit_norm_response [1]
[Previous]
[Next]
NAME:
eit_calgroups
PURPOSE
returns a list of valid averaged calibration lamps file names
INPUTS
none
OUTPUTS
a string array containing the list of calibration lamps
KEYWORDS
none
WARNINGS
needs to be regularly updated
may not include the most recent calibration lamps
CREATION
8-Jun-2001. F. Auchere
MODIFICATIONS
Updated: Mon Dec 24 EIT planner
Updated: Thu Jan 10 EIT planner
Updated: Mon Feb 25 EIT planner
Updated: Mon Feb 25 EIT planner
Updated: Wed Apr 3 EIT planner
Updated: Fri Jul 19 EIT planner
Updated: Wed Jul 24 F. Auchere
Updated: Wed Sep 4 F. Auchere
Updated: Wed Sep 25 EIT planner
Updated: Tue Dec 24 EIT planner
Updated: Mon Jan 6 EIT planner
Updated: Wed Jan 8 EIT planner
Updated: Tue Feb 4 EIT planner
Updated: Mon Mar 10 EIT planner
Updated: Tue Mar 11 EIT planner
Updated: Thu Sep 4 EIT planner
Updated: Wed Feb 4 EIT planner
Updated: Wed Feb 11 EIT planner
Updated: Thu Feb 12 EIT planner
Updated: Fri May 14 EIT planner
Updated: Mon May 17 EIT planner
Updated: Sat Jun 12 EIT planner
Updated: Mon Jun 14 EIT planner
Updated: Tue Jul 13 EIT planner
Updated: Mon Jul 19 EIT planner
Updated: 02-sep-2004 F. Auchere : removed cal19990717.fts from the list - bad file
Updated: Tue Oct 12 2004 EIT planner
Updated: Thu Oct 21 2004 EIT planner
Updated: Tue Nov 16 2004 EIT planner
Updated: Fri Jan 7 2005 EIT planner
Updated: Tue Mar 8 2005 EIT planner
Updated: Thu Mar 31 2005 EIT planner
Updated: Sun Apr 3 2005 EIT planner
Updated: Fri Apr 22 2005 EIT planner
Updated: Fri Jun 10 2005 EIT planner
Updated: Tue Jul 19 2005 EIT planner
Updated: Tue Jul 19 2005 EIT planner
Updated: Tue Jul 19 2005 EIT planner
Updated: Wed Aug 3 2005 EIT planner
Updated: Thu Aug 4 2005 EIT planner
Updated: Mon Aug 8 2005 EIT planner
Updated: Mon Aug 8 2005 EIT planner
Updated: Sun Sep 11 2005 EIT planner
Updated: Mon Oct 31 2005 EIT planner
Updated: Sun Nov 27 2005 EIT planner
Updated: Sun Nov 27 2005 EIT planner
Updated: Mon Dec 26 2005 EIT planner
Updated: Mon Dec 26 2005 EIT planner
Updated: Mon Dec 26 2005 EIT planner
Updated: Tue Jan 24 2006 EIT planner
Updated: Wed Mar 22 2006 EIT planner
Updated: Wed Mar 29 2006 EIT planner
Updated: Sun Apr 2 2006 EIT planner
Updated: Sun Jun 4 2006 EIT planner
Updated: Mon Jun 12 2006 EIT planner
Updated: Mon Jun 12 2006 EIT planner
Updated: Wed Jun 14 2006 EIT planner
Updated: Wed Jun 14 2006 EIT planner
Updated: Thu Jun 15 2006 EIT planner
Updated: Thu Jun 15 2006 EIT planner
Updated: Fri Jun 16 2006 EIT planner
Updated: Sun Jul 16 2006 EIT planner
Updated: Sun Jul 16 2006 EIT planner
Updated: Fri Jul 21 2006 EIT planner
Updated: Thu Oct 26 2006 EIT planner
Updated: Tue Oct 31 2006 EIT planner
Updated: Wed Nov 29 2006 EIT planner
CALLED BY
eit_find_last_cal, eit_findcalgroup
[Previous]
[Next]
NAME :
EIT_CATRD
PURPOSE : Read the daily catalogs and return entries as requested output
CATEGORY : Catalogs
EXPLANATION : This function reads the daily QKL catalogs or the single
Level Zero catalog.
SYNTAX : list = eit_catrd(in_date, group=group, interactive=interactive,
date, long=long, raw=raw, yester=yester, lz=lz,
no_next=no_next, short=short, input_file=input_file,
records=records, filename=filename, tai=tai,
start_date=start_date,stop_date=stop_date,object=object,
wave=wave,filter=filter,nx=nx,ny=ny,xbin=xbin,ybin=ybin,
timerange=timerange,sci_obj=sci_obj,infov=infov,hihg_c=high_c)
CALLED BY:
EIT_DISPLAY, EIT_PA, EIT_PA_ev, EIT_POSTAGE, EIT_SUBFIELD, EIT_SUB_MOVIE, EIT_XCAT
MK_EIT_L1, MOVIE_MAKER, SUBFIELD_EV, get_eit_data
EXAMPLES : return string listing of all LZ 304 images from Jun 2, 1996
to June 22, 1996 taken in the clear filter and return
filenames,
IDL> file=''
IDL> list=eit_catrd(start_date='2-jun-96',stop_date=$
'22-jun-96',wave='304',filter='clear',file=file)
CALLED BY : top level
CALLS TO : timegrid, str_copy_tags, cat_open, cat_directory,
cat_to_text, where_arr, anytim2utc, utc2doy, doy2utc,
xmessage, utc2tai
ARGUMENTS :
DATE IF present, a string date in the format "yyyy/mm/dd"
of the daily catalog to be searched.
Default is the date of the most recent catalog.
KEYWORDS :
INPUT_FILE IF present, this is the catalog which is read in.
This may be any of the catalog: "merged", QKL daily or
LZ daily. Default is the "merged" catalog.
LZ IF present read from the merged Level-zero catalog
YESTER IF present, indicates the number of 'days/catalogs' backward
in time you want as the catalog input. Note that this is
really the number of backward catalogs, so days may not
line up as expected, but they should.
New feature: print,catrd_daily("y=3") for three days ago.
NO_NEXT IF present and set THEN return a null string IF the specIFied
date does not have a catalog.
SHORT_TEXT IF present and set and TEXT is the output mode then the
displayed text will fit in 80 characters w/ less info.
RECORDS IF present, then return catalog file record numbers matching
the entries in the returned variable. This is the VMS
record number.
LONG IF present and set, then all catalog entries are considered,
including duplicates and entries are left in unsorted,
reform order.
RAW IF present and set , return the catalog entries as the raw
data structure instead of the default of a strarr of text
entries.
NOTE: the default is to return a string array of entries.
FILENAME IF present, then return the fits file names matching
the entries in the returned variable.
TAI IF present, then return the TAI time matching
the entries in the returned variable.
START_DATE IF present and a LZ listing then sets start date of list
STOP_DATE IF prsent and a LZ listing then sets end date of list
WAVE Return matches from LZ/QKL catalog (171,195,284,304 or 1-4)
FILTER Return matches from LZ/QKL catalog (0-4 or one of:
['AL+1', 'BLK EST', 'CLEAR', 'BLK WST', 'AL+2' ] )
NX Return matches from LZ/QKL catalog for xsize in pixels
NY Return matches from LZ/QKL catalog for ysize in pixels
XBIN Return matches from LZ/QKL catalog for x summing in pixels
YBIN Return matches from LZ/QKL catalog for y summing in pixels
TIMERANGE 2 element vector containing start and stop dates
OBJECT IF present select object, dark, cal lamp, full FOV,
partial FOV
CAT_STRUCT An alternate input, a raw catalog structure previsouly
created - useful to get search tools, if input is
a index structure will return output structure not text
SS Return subscripts of input index structure instead of
selected structure
SCI_OBJ Return matches from LZ/QKL catalog for entered keyword
INFOV Return matches (partial FOV) if entered coordinates,
Solar arcsecs, is in field of view, e.g, [1,-960]=south pole
HIGH_C Set if want high cadence (defined to be greater than 5
images per file) subfield images
OUTPUTS :
Return a string arrar of catalog entries, unless the /RAW keyword
is set then return raw structure of entries.
IF the /LONG keyword is set then return all catalog entries and
do not sort.
IF input is index structure returns output structure
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], CATDISP_EVENT, CAT_DIRECTORY, CAT_OPEN [1]
CAT_OPEN [2], CAT_TO_TEXT, CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
CW_BGROUP, CW_FIELD, DATATYPE [1], DATATYPE [2], DATATYPE [3], DELVARX [1]
DELVARX [2], DELVARX [3], DELVARX [4], DOY2UTC, EIT_FXPAR, EIT_PIXSIZE
GET_DFONT [1], GET_DFONT [2], HELPCAT, IS_GSFCVMS, KILLHELP, SEARCH_CAT, STR_SEP
TAG_EXIST [1], TAG_EXIST [2], UTC2DOY, UTC2SEC, UTC2TAI, XMANAGER, XMESSAGE
concat_dir [4], delvarx [5], get_logenv [1], get_logenv [2], str_copy_tags [1]
str_copy_tags [2], str_copy_tags [3], str_copy_tags [4], timegrid, where_arr [1]
where_arr [2]
COMMON : cat_params, cat_par2, helpshare
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1996 March 31, added /NO_NEXT keyword
1996 May 13, allow YESTER to use numerical values
V2, Jeff Newmark - complete revision, LZ catalog, widget
interface, etc...
1996 Dec 23 - add selection over times as well as dates
J Newmark
1996 Jan 15 - changed date variable to indate due to
conflict with IDLASTRO function date!
1996 Jan 16 - use eit_fxpar to obtain selection parameters
allow cat_struct input to be index structures
created by read_eit as well as raw catalog structures
V3.0
1997 Mar 20 J. Newmark - add in sci_obj field for catalog
1997 Jun 03 J. Newmark - fix bug in wavelength selection
1998 Jan 20 J. Newmark - add infov keyword
V3.1
1998 Apr 01 J. Newmark - automatic selection of LZ vs. QKL
1998 Apr 07 J. Newmark - allow mixed LZ + QKL
1998 Oct 23 J. Newmark - add high_c keyword
1999 Oct 15 R.D.Bentley - Try EIT_LZ_CAT before defaulting to EIT_LZ for catalog
2000 Jan 19 J. Newmark - bug fix in little used path
Version 2.1
2003 February 13 D.M. fecit - Added handling of big-endian OSes.
CONTACT : gurman@gsfc.nasa.gov
[Previous]
[Next]
NAME :
EIT_CATRD
PURPOSE : Read the daily catalogs and return entries as rEQuested output
CATEGORY : Catalogs
EXPLANATION : This function reads the daily QKL catalogs or the single
Level Zero catalog.
SYNTAX : list = eit_catrd(in_date, group=group, interactive=interactive,
date, long=long, raw=raw, yester=yester, lz=lz,
no_next=no_next, short=short, input_file=input_file,
records=records, filename=filename, tai=tai,
start_date=start_date,stop_date=stop_date,object=object,
wave=wave,filter=filter,nx=nx,ny=ny,xbin=xbin,ybin=ybin,
timerange=timerange,sci_obj=sci_obj,infov=infov,hihg_c=high_c)
CALLED BY:
EIT_DISPLAY, EIT_PA, EIT_PA_ev, EIT_POSTAGE, EIT_SUBFIELD, EIT_SUB_MOVIE, EIT_XCAT
MK_EIT_L1, MOVIE_MAKER, SUBFIELD_EV, get_eit_data
EXAMPLES : return string listing of all LZ 304 images from Jun 2, 1996
to June 22, 1996 taken in the clear filter and return
filenames,
IDL> file=''
IDL> list=eit_catrd(start_date='2-jun-96',stop_date=$
'22-jun-96',wave='304',filter='clear',file=file)
CALLED BY : top level
CALLS TO : timegrid, str_copy_tags, cat_open, cat_directory,
cat_to_text, where_arr, anytim2utc, utc2doy, doy2utc,
xmessage, utc2tai
ARGUMENTS :
DATE IF present, a string date in the format "yyyy/mm/dd"
of the daily catalog to be searched.
Default is the date of the most recent catalog.
KEYWORDS :
INPUT_FILE IF present, this is the catalog which is read in.
This may be any of the catalog: "merged", QKL daily or
LZ daily. Default is the "merged" catalog.
LZ IF present read from the merged Level-zero catalog
YESTER IF present, indicates the number of 'days/catalogs' backward
in time you want as the catalog input. Note that this is
really the number of backward catalogs, so days may not
line up as expected, but they should.
New feature: print,catrd_daily("y=3") for three days ago.
NO_NEXT IF present and set THEN return a null string IF the specIFied
date does not have a catalog.
SHORT_TEXT IF present and set and TEXT is the output mode then the
displayed text will fit in 80 characters w/ less info.
RECORDS IF present, then return catalog file record numbers matching
the entries in the returned variable. This is the VMS
record number.
LONG IF present and set, then all catalog entries are considered,
including duplicates and entries are left in unsorted,
reform order.
RAW IF present and set , return the catalog entries as the raw
data structure instead of the default of a strarr of text
entries.
NOTE: the default is to return a string array of entries.
FILENAME IF present, then return the fits file names matching
the entries in the returned variable.
TAI IF present, then return the TAI time matching
the entries in the returned variable.
START_DATE IF present and a LZ listing then sets start date of list
STOP_DATE IF prsent and a LZ listing then sets end date of list
WAVE Return matches from LZ/QKL catalog (171,195,284,304 or 1-4)
FILTER Return matches from LZ/QKL catalog (0-4 or one of:
['AL+1', 'BLK EST', 'CLEAR', 'BLK WST', 'AL+2' ] )
NX Return matches from LZ/QKL catalog for xsize in pixels
NY Return matches from LZ/QKL catalog for ysize in pixels
XBIN Return matches from LZ/QKL catalog for x summing in pixels
YBIN Return matches from LZ/QKL catalog for y summing in pixels
TIMERANGE 2 element vector containing start and stop dates
OBJECT IF present select object, dark, cal lamp, full FOV,
partial FOV
CAT_STRUCT An alternate input, a raw catalog structure previsouly
created - useful to get search tools, if input is
a index structure will return output structure not text
SS Return subscripts of input index structure instead of
selected structure
SCI_OBJ Return matches from LZ/QKL catalog for entered keyword
INFOV Return matches (partial FOV) if entered coordinates,
Solar arcsecs, is in field of view, e.g, [1,-960]=south pole
HIGH_C Set if want high cadence (defined to be greater than 5
images per file) subfield images
OUTPUTS :
Return a string arrar of catalog entries, unless the /RAW keyword
is set then return raw structure of entries.
IF the /LONG keyword is set then return all catalog entries and
do not sort.
IF input is index structure returns output structure
COMMON : cat_params, cat_par2, helpshare
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1996 March 31, added /NO_NEXT keyword
1996 May 13, allow YESTER to use numerical values
V2, Jeff Newmark - complete revision, LZ catalog, widget
interface, etc...
1996 Dec 23 - add selection over times as well as dates
J Newmark
1996 Jan 15 - changed date variable to indate due to
conflict with IDLASTRO function date!
1996 Jan 16 - use eit_fxpar to obtain selection parameters
allow cat_struct input to be index structures
created by read_eit as well as raw catalog structures
V3.0
1997 Mar 20 J. Newmark - add in sci_obj field for catalog
1997 Jun 03 J. Newmark - fix bug in wavelength selection
1998 Jan 20 J. Newmark - add infov keyword
V3.1
1998 Apr 01 J. Newmark - automatic selection of LZ vs. QKL
1998 Apr 07 J. Newmark - allow mixed LZ + QKL
1998 Oct 23 J. Newmark - add high_c keyword
1999 Oct 15 R.D.Bentley - Try EIT_LZ_CAT before defaulting to EIT_LZ for catalog
2000 Jan 19 J. Newmark - bug fix in little used path
CONTACT : newmark@eitv2.nascom.nasa.gov
[Previous]
[Next]
Name: eit_catrd
Purpose: front end to catrd_daily (read time range)
Input Paramters:
time0 - start time
time1 - stop time
Keyword Parameters:
raw - if set, return catalog structures (default=string)
logn - switch (see catrd_daily documentation)
filename (output) - eit filenames
records (output) - (see catrd_daily doc)
Calling Sequence:
cats=eit_catrds(time0 [,time1, /raw, /long, records=records, filename=filename]
CALLS: ***
CATRD_DAILY, EIT_CATRDS, data_chk [1], data_chk [2], str_concat [1]
str_concat [2], timegrid, ut_time [1], ut_time [2]
CALLED BY:
EIT_DISPLAY, EIT_PA, EIT_PA_ev, EIT_POSTAGE, EIT_SUBFIELD, EIT_SUB_MOVIE, EIT_XCAT
MK_EIT_L1, MOVIE_MAKER, SUBFIELD_EV, get_eit_data
History:
24-may-1996 S.L.Freeland
[Previous]
[Next]
Name: eit_colors
Purpose: load EIT color tables
Input Parameters:
table - EIT table number or associated wavelength
Output Parameters:
r,g,b - optional return RGB
Calling Examples:
eit_colors,0 ; load the first EIT table (relative)
eit_colors,41 ; load the first EIT table (absolute)
eit_colors,284,r,g,b ; map wavelenght-> table, return RGB
Method:
call ssw_colors ; contains links to SSW color table files
CALLS: ***
is_member [1], is_member [2], ssw_colors
CALLED BY:
EIT_MKMOVIE, LOAD_EIT_COLOR, do_disp_month, eit_fulldisk, lasteit
mk_lasteit_movie, ssw_track_demo
History:
18-Feb-1996 (S.L.Freeland)
[Previous]
[Next]
EIT_COORD_FITPARAM
Purpose - return the parameters of the Fourier fitting of the 195
pointing data
Inputs: named variables that will contain the values of the amplitude
and phase of the Fourier components as well as time boundaries
Output: the values for these variables
Created: 29 Nov 2001. F. Auchere
Updated: Mon Feb 25 EIT planner
Updated: Wed Apr 3 EIT planner
Updated: Sat Jul 20 EIT planner
Updated: Mon Aug 12 EIT planner
Updated: Tue Aug 13 EIT planner
Updated: Thu Aug 15 EIT planner
Updated: Wed Sep 4 EIT planner
Updated: Wed Sep 11 EIT planner
Updated: Wed Dec 25 EIT planner
Updated: Tue Mar 11 EIT planner
Updated: Wed Sep 3 EIT planner
Updated: Mon Oct 20 EIT planner
Updated: Wed Feb 4 EIT planner
Updated: Thu Feb 5 EIT planner
Updated: Wed Feb 11 EIT planner
Updated: Wed Feb 11 EIT planner
Updated: Thu Feb 12 EIT planner
Updated: Fri May 14 EIT planner
Updated: Mon May 17 EIT planner
Updated: Sat Jun 12 EIT planner
Updated: Sat Jun 12 EIT planner
Updated: Mon Jun 14 EIT planner
Updated: Tue Jul 13 EIT planner
Updated: Mon Jul 19 EIT planner
Updated: Wed Oct 13 EIT planner
Updated: Thu Oct 21 EIT planner
Updated: Wed Nov 17 EIT planner
Updated: Fri Jan 7 EIT planner
Updated: Wed Mar 9 EIT planner
Updated: Fri Apr 1 EIT planner
Updated: Wed Aug 3 EIT planner
Updated: Thu Aug 4 EIT planner
Updated: Mon Aug 8 EIT planner
Updated: Sun Sep 11 EIT planner
Updated: Mon Oct 31 EIT planner
Updated: Sun Nov 27 EIT planner
Updated: Mon Dec 26 EIT planner
Updated: Mon Dec 26 EIT planner
Updated: Tue Dec 27 EIT planner
Updated: Wed Jan 25 EIT planner
Updated: Thu Mar 23 EIT planner
Updated: Wed Mar 29 EIT planner
Updated: Sun Apr 2 EIT planner
Updated: Mon Jun 5 EIT planner
Updated: Mon Jun 12 EIT planner
Updated: Fri Jun 16 EIT planner
Updated: Fri Jun 16 EIT planner
Updated: Sun Jul 16 EIT planner
Updated: Fri Oct 27 EIT planner
Updated: Tue Oct 31 EIT planner
Updated: Thu Nov 2 EIT planner
Updated: Wed Nov 29 EIT planner
CALLED BY
eit_point
[Previous]
[Next]
Project : SOHO - LASCO/EIT
Name : EIT_DAILY
Purpose : Get a list of EIT images on a selected date and then create
a postscript file of a wavelength sequence of 4.
Explanation : creates a database query to get a list of files for a single day and
then brings up a list widget to select the desired files and they
will then be scaled and output to a postscript file
Use : EIT_DAILY
Inputs : none: program will ask a date to be typed in
Keywords : bw is set to create a black and white output
Outputs : None: program creates a postscript file
Calls : ***
ANYTIM2CAL, ANYTIM2UTC [1], ANYTIM2UTC [2], EIT_CT, EIT_DEGRIDN, EIT_LISTER
EIT_SCALE, EIT_TRANSFER, FXPAR [1], FXPAR [2], GREP, HEADFITS [1], HEADFITS [2]
HEADFITS [3], L_RDFTS, REMCHAR [1], REMCHAR [2], REMCHAR [3], STR_INDEX [1]
STR_INDEX [2]
Restrictions: Uses data in NRL archive
Side effects: Creates postscript file in current working directory
Category : Image Display.
Written : Ken Dere September 1996.
Modified : Ken Dere, Dec. 1996 to include black and white output
Version :
[Previous]
[Next]
EIT_DARK
Purpose - returns a fourier series fit to the measured offset values
Inputs: date
Output: the offset value
Created: Mar 2001. F. Auchere
Modified: 14 Dec 2001. F. Auchere - New fit with Fourier series
Updated: Fri Dec 14 EIT planner
Updated: Mon Dec 17 EIT planner
Modified: 12 Jan 2002. F. Auchere - Handles different offset in binned images
CALLS:
CALLED BY
EIT_DISPLAY, EIT_PREP, MOVIE_MAKER, eit_dump, eit_fulldisk, eit_readcal, lasteit
[Previous]
[Next]
EIT_DARK_FITPARAM
PURPOSE: Returns the coefficients of the fourier
series fit to the measured offset values
Inputs: date
Output: the offset value
Created: 20-dec-2001. F. Auchere
Updated: Mon Feb 25 EIT planner
Updated: Wed Apr 3 EIT planner
Updated: Sat Jul 20 EIT planner
Updated: Mon Aug 12 EIT planner
Updated: Tue Aug 13 EIT planner
Updated: Thu Aug 15 EIT planner
Updated: Wed Sep 4 EIT planner
Updated: Wed Sep 11 EIT planner
Updated: Wed Dec 25 EIT planner
Updated: Tue Mar 11 EIT planner
Updated: Thu Sep 4 EIT planner
Updated: Mon Oct 20 EIT planner
Updated: Wed Feb 4 EIT planner
Updated: Wed Feb 11 EIT planner
Updated: Fri May 14 EIT planner
Updated: Mon May 17 EIT planner
Updated: Sat Jun 12 EIT planner
Updated: Tue Jul 13 EIT planner
Updated: Wed Oct 13 EIT planner
Updated: Thu Oct 21 EIT planner
Updated: Wed Nov 17 EIT planner
Updated: Sat Jan 8 EIT planner
Updated: Wed Mar 9 EIT planner
Updated: Wed Aug 3 EIT planner
Updated: Wed Dec 28 EIT planner
Updated: Mon Jun 5 EIT planner
Updated: Fri Jun 16 EIT planner
Updated: Sun Jul 16 EIT planner
Updated: Fri Oct 27 EIT planner
Updated: Thu Nov 2 EIT planner
CALLED BY
eit_dark
[Previous]
[Next]
Name: eit_data_scale
Purpose: scale EIT data to some standard (and provide single point maint)
Input Parameters:
data - data to scale
wave - (optional) - wavelength dependent scale
History:
17-Feb-1996 (S.L.Freeland)
Restrictions:
For now, only one algorithm, used for all wavelengths
[Previous]
[Next]
Name:
eit_degrid
Purpose:
Degrid an EIT image, either full FOV or POV and/or binned.
Input Parameters:
image = image w/ missing blocks raised to the detector offset, and
detector offset subtracted.
header = FITS header or EIT index structure for image
Note that image can be full-res or a binned image.
Output:
clean_image = degridded image
Keywords:
final = set if using LZ data
image_no = set if image is from a 3-D LZ FITS file
verbose = set if messages discriminating new read or pre-read wanted
Calling Sequence:
clean_image = eit_degrid(image, fits_header)
clean_image = eit_degrid(image, index)
clean_image = eit_degrid(image, fits_header, image_no=image_no)
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
DATATYPE [1], DATATYPE [2], DATATYPE [3], EIT_FXPAR, IS_GSFCVMS, concat_dir [4]
CALLED BY:
EIT_IMAGE, EIT_PREP, MAKE_IMAGE_EIT img hdr
Restrictions:
Can't handle a vector of images {yet}.
If you pipe an LZ file to this, you'd better have the fits header.
If you're doing this with single-image files, it will work just
fine if have a FITS header structure.
History:
23-Mar-96 - (DMF) - Written
8-Apr-96 - (JSN) - Fix for incorrect entry of pixel size.
22-Jul-96 - (JSN) - Added capability of export SSW
20-Aug-96 - (JSN) - Added image_no keyword for 3D LZ data.
28-Aug-96 - (JSN) - changed to always XDR format
29-Aug-96 - (DMF) - XDR too slow! Backed off to block I/O files in
OVMS interal representation
9-Oct-96 - (BNH) - Now accepts FITS headers or index structures.
Cleaned a few things up and dispatched FIND_KEYWORD
for good. Pulled all the history of this routine
into the "History" section. Still need to sort
out how to make index structures to cope
with the new-improved "LZ" file structure.
16-Oct-96 - (BNH) - Add a common block for the currently-used grid.
1997/04/11 - JBG - Added VERBOSE keyword.
1997/11/10 - JSN - Use newly calculated grid (after July 1997)
1998/02/04 - JSN - New AL +1 grids
3-Jun-1998 - Zarro (SAC/GSFC)
- Added calls to temporary for improved
memory management.
1999/07/14 - JSN - Streamline code
28-Sep-1999 - N.Rich - Adapt for LASCO fits header: object, wave, DATE-OBS
18-May-2000 - JSN - online LASCO compatibility
17-May-2001 - JSN - only original grids
[Previous]
[Next]
NAME:
eit_degrid_fft
PURPOSE:
Degrid an EIT image via an FFT transform. EIT_DEGRID_FFT assumes
missing blocks have been raised to the detector offset, and the
detector offset subtracted.
CALLING SEQUENCE:
clean_image = eit_degrid_fft(image, header [, /final])
INPUTS:
image = Dark-subtracted EIT image
header = FITS header
OPTIONAL KEYWORDS:
final = ????
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], find_fits_keyword
RESTRICTIONS:
The input image must have the missing blocks raised to detector
offset, and the detector offset subtracted. (Fecit's words)
At this time (9-Apr-1996) only supports full-frame full-resolution
images.
TODO:
Need to stamp the FITS header so we know this image was degridded
with <mumble> version of EIT_DEGRID_FFT.
MODIfICATION HISTORY:
23-Mar-1996 - (DMF) - Written
9-Apr-1996 - (BNH) - Cleaned up and documented
[Previous]
[Next]
NAME:
eit_degrid_smooth
PURPOSE:
Smoothing algorithm for degridding EIT images. Adaptation of
A. Maucherat's early routines.
CALLING SEQUENCE:
clean_image = eit_degrid_smooth(dirty_image [,residue=residue,
grid=grid])
INPUTS:
dirty_image = dark-subtracted image.
OUTPUT:
clean_image = alog10(degridded image)
OPTIONAL OUTPUTS:
residue = residue of degridding algorithm
grid = Calculated grid
CALLS: ***
REVERSE
RESTRICTIONS:
Can't handle images smaller than 42x42 pixels.
TO DO:
* This routine doesn't conserve counts. Maybe it never will.
* I wonder how well this code deals with small images.
* Need to keep a history of this in the FITS header. EIT_DEGRID
also suffers from this.
MODifICATION HISTORY:
15-Mar-1996 - (AM) - Written
9-Apr-1996 - (BNH) - Cleaned up to make slightly more readable,
turned into a function call, Added some
comments and streamlined some code.
10-Apr-1996 - (BNH) - Added support for sub_images.
- Image is converted to a FLOAT on entry.
10-Oct-1996 - (BNH) - Corrected an entry in the median array.
Thanks to RAS for noticing this.
[Previous]
[Next]
NAME:
eit_despeckle
PURPOSE
a simple despecklimg routine
INPUTS
img - image to despeckle
OUTPUTS
returns the despeckled image
KEYWORDS
none
PROCEDURE
median filter with threshold
CREATION
18-May-2001. F. Auchere
MODIFICATIONS
[Previous]
[Next]
NAME:
EIT_DISPLAY
PURPOSE:
This function will display EIT movies (IDL Save Sets) and
EIT images (FITS files) using eit_movie and eit_image
respectively.
CATEGORY:
WIDGET interface
CALLING SEQUENCE:
EIT_DISPLAY
INPUTS: None, presently this routine reads predefined directories
namely the logical REF_DIR (FITS images) or MOVIE_DIR
(3-dim IDL Save Sets).
KEYWORD PARAMETERS: None
OUTPUTS:
This routine displays images. It does not return any outputs.
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3]
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], CW_BGROUP, CW_ZOOM, DATATYPE [1]
DATATYPE [2], DATATYPE [3], DEVICELIB [1], DEVICELIB [2], DISPHD, DOY2UTC
EITDISP_EVENT, EIT_CATRD [1], EIT_CATRD [2], EIT_FXPAR, EIT_IMAGE, EIT_MOVIE
FIND_FILE, GET_DFONT [1], GET_DFONT [2], HEADFITS [1], HEADFITS [2], HEADFITS [3]
HELPOUT, IMAGELIB, IS_GSFCVMS, KILLHELP, KILLZOOM, LABEL_IMAGE, LOADCT, MKZOOM, OBSLOG
PUT, RDPIX, READFITS [1], READFITS [2], READFITS [3], REVERSE, SAVE_PSDIR, SHOWEIT
STR_SEP, UTC2DOY, XIN, XLOADCT [1], XLOADCT [2], XLOADCT [3], XMANAGER, XMESSAGE
break_file [4], concat_dir [4], eit_catrd [3], eit_dark, eit_dump, eit_file2path
is_fits [1], is_fits [2], is_ss [1], is_ss [2], qmenu
COMMON BLOCKS:
filestuff,dirspec,flist,index,list7,text6 - this common block is
used to pass back and forth to the event handler the
directory, filelist, file index, file widget and directory widget.
helpshare,helpbase - a common block for the help widget
SIDE EFFECTS:
unsure if common block must be re-initialized after each run.
RESTRICTIONS:
Serious memory hog, especially if display movies
PROCEDURE:
A widget is displayed. The user can choose either IDL Save
Sets (assumed therefore to be movies) or FITS files (assumed
to be images). The user then selects a file to display using
the appropriate mechanism. The user can zoom an image. This
calss the widget cw_zoom.
MODIFICATION HISTORY:
Written by: J. Newmark Date. Feb. 1996
1996 march 14 - E. Einfalt - changed title banner
1996 june 12 - J. Newmark - use catalog for image
listings, add RDPIX function
1996 june 20 - J. Newmark - enable use of "efz" files
1996 july 8 - J. Newmark - add print button
1996 july 19 - J. Newmark - added calls to eit_file2path
to get directory names and pathways, allow
for collapsed data trees
1996 aug 20 - J. Newmark - added ability to display
LZ 3-D data
1996 nov 26 - D.M. fecit - added /flat to all EIT_IMAGE
calls
1996 dec 3 - J. Newmark - added button for MOVIE_ARCHIVE
1997 mar 12 - D.M. fecit - changed exposure time strings
to floating point
2003 may 2 - C.A. Young - Added year droplist to select
the year subdirectory in new movie archive
2003 may 5 - E Einfalt - Added VMS compatibility to yearly
subdirectory code
2003 may 5 - D.M. fecit - Modified annual movie directory
dropdown code to accommodate years other than
2003.
2003 jul 29 - D.M. fecit - Modified default YSIZE of
file list in main widget to allow view of more
file names.
[Previous]
[Next]
Name:
eit_dump
Purpose: produce postscript or gif files
Input Parameters:
filenam - FITS file of image
Output Parameters:
file - name format is efr'date'_wave.'ps or gif'
Keyword Parameters:
post - set for postscript output - (default)
gif - set for gif output -image only no labels
jpeg - set for JPEG output - image only, no labels
encap - set for encapsulated postscript output
printer - set for printing to value
outdir - set for output directory of saved file
window - set to dump current window
surround - set to replace missing blocks with surrounding ones
file - set for output file name
annotate -set to string if wish annotated PS output for window dump
image_no -set for image_no in 3d LZ file
Calling Sequence:
eit_dump,'efr19960702.123456',/post
Method:
read FITS file, load proper color table, subtract background,
replace missing blocks, degrid, bytscl(alog10(image)), save
in selected format
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
EIT_FXPAR, EIT_IMAGE, HEADFITS [1], HEADFITS [2], HEADFITS [3], LOADCT, WRITE_GIF
concat_dir [4], eit_dark, eit_file2path, is_fits [1], is_fits [2]
CALLED BY:
EIT_DISPLAY
History:
2-Jul-1996 - J. Newmark - created
22-Jul-1996 - J. Newmark - allow for capture of IDL window
23-Jan-1997 - J. Newmark - allow file keyword for output
18-Feb-1997 - J. Newmark - call eit_image for processing, add
annotate keyword
19-May-1997 - J. Newmark - add image_no keyword
1997 May 22 - D.M. fecit - added JPEG keyword
[Previous]
[Next]
NAME:
eit_eff_area
PURPOSE:
Return the EIT effective area as a function of wavelength
CALLING SEQUENCE:
Area = eit_eff_area(index,waveout=waveout)
Area = eit_eff_area(roadmap,wave) ;User supplies wavelengths
Area = eit_eff_area(index,wave)
Area = eit_eff_area(Sector,wave) ;Sector is an integer array
Area = eit_eff_area(Sector,wave,filter=filt) ;Filt is an integer array
Area_structure = eit_eff_area()
Area_structure = eit_eff_area(filein=filein)
INPUTS:
Item - EIT index structure or roadmap. May be a scalar
or an array. Sector and Filters are extracted.
or - an integer scalar or array with the EIT sector number.
OPTIONAL INPUTS:
wave - Wavelengths at which the values of the area are
desired. If not specified, the wavelengths are
generated internally and returned in waveout.
OPTIONAL INPUT KEYWORDS:
filters - To override the values of filters in Item
verbose - Set to provide additional information
filein - To read a specific EIT effective area file. By
default reads the following file:
concat_dir('SSW_EIT_RESPONSE','era*.genx')
Use filein='your_file_name' to override the default.
text_area - If file is read, return the text information (string) array
OUTPUTS:
Returns the EIT effective area.
Output size will be (N_elements(wave), N_elements(item))
OPTIONAL OUTPUT KEYWORDS:
waveout - If wave is not provided, then wavelengths are
defined by eit_eff_area and returned in waveout.
OPTIONAL INPUT/OUTPUT KEYWORD:
area_str - The contents of the EIT effective area file is
returned in this structure variable. If this
variable is defined, the routine will not bother
to re-read the file on subsequent calls.
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2]
CONCAT_DIR [3], DSPLINE, FILE_EXIST [2], break_file [4], concat_dir [4]
file_exist [1], file_exist [3], gt_EITfilter, gt_EITsector, restgen [1]
restgen [2]
CALLED BY:
plot_eit_area, plot_eit_entrance, plot_eit_filters, plot_eit_mirror
RESTRICTIONS:
If filters is defined, it must be of length 1 or length n_elements(Item)
NOTE: It is assumed that area.filter(0) = Entrance filter
area.filter(1) = CCD fixed filter
area.filter(2) = Filter wheel C+Al filter
area.filter(3) = Filter wheel thin Al filter
This is order is established in the routine mk_eit_area.pro
METHOD:
The EIT effective area file is read. If wave is specified, the
values are interpolated in Log(area).
HISTORY:
9-Jun-95, J. R. Lemen (LPARL), Written
[Previous]
[Next]
Name: eit_file2path
Purpose: translate eit file name to "standard path" on local system
Input Parameters:
files - one or more eit file names (with or without path)
Output Parameters:
exist - boolean vector - <files> online?
exist_count - count where(<files>) online
Keyword Parameters:
topeit - top level directory (default=EIT_QUICKLOOK)
curdir - if set, set topeit to current directory
lz - if set, set topeit to EIT_LZ and use YYYY/MM format
collapse - if set, assme all files in single directory <topeit>
default is GSFC standard topeit/YYYY/MM/DD/filename
gavroc - set if running on gavroche, special directory structure
otherwise will check
Calling Sequence:
eitpath=eit_file2path(eitfilenames [,/collapse, /curdir, topeit='xxx'])
eitpath=eit_file2path(filenames,exist,count) ; exist = online?
Calling Examples:
print,eit_file2path('efr19960521.043112') ; "GSFC-like" tree
....ate/data/processed/eit/quicklook/1996/05/21/efr19960521.043112
print,eit_file2path('efr19960521.043112',/collapse) ; collapsed tree
...ate/data/processed/eit/quicklook/efr19960521
print,eit_file2path('efr19960521.043112',/curdir) ; local version
/usr/users/freeland/dev/eitpath/efr19960521.043112
CALLED BY:
EIT_DISPLAY, EIT_IMAGE, EIT_MKMOVIE, EIT_POSTAGE, EIT_PREP, EIT_SUB_MOVIE, MK_EIT_L1
MOVIE_MAKER, eit_dump, eit_fulldisk, eit_fulldiskdb, eit_readcal, get_eit_data
lasteit, read_eit
History:
21-May-1996 (S.L.Freeland)
22-May-1996 (S.L.Freeland) - added 2nd param (exist)
01-Aug-1996 (J. Newmark) - added lz keyword for level-zero data
16-Aug-1996 (J. Newmark) - added special paths for GAVROC
21-Jan-1997 (J. Newmark) - use is_gsfcvms function
04-Mar-1997 (J. Nemwark) - allow input of short catalog listing
28-Mar-1997 (S. Freeland) - check 'EIT_DATA_STYLE' to allow
site specification ('collapsed')
12-May-1998 (J. Newmark) - allow combination of QKL and LZ
fix archive location on IS_GSFCVMS
22-Jul-1998 (J. Newmark) - fix multiple yr/month in LZ data
on VMS
27-Jan-1998 (J. Newmark) - bug in LZ/QKL combination
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3]
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], FILE_EXIST [2], IS_GSFCVMS
break_file [4], concat_dir [4], curdir [1], curdir [2], data_chk [1], data_chk [2]
file_exist [1], file_exist [3], get_logenv [1], get_logenv [2], prstr [1]
prstr [2], strmids [1], strmids [2]
Restrictions:
assume filenames in given call have same length prefix
[Previous]
[Next]
Name: eit_files
Purpose: return online eit file names
Input Parameters:
t0, t1 - Optional start time and stop time
Keyword Parameters:
last - number of days to look at (Most recent N days)
collapse - SWITCH, if set, all files in one directory
full_disk - if $EIT_FULLDISK is defined, use that tree
matchfiles - instead of t0/t1, pass in list of filenames/times
and get best match (at same or different wavelength)
delta - only if MATCHFILES input, deltaT(sec) between desired/return
Calling Sequence:
files=eit_files
CALLS: ***
BOX_MESSAGE, CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], FILE2TIME
FILE_EXIST [2], Month_id [1], Month_id [2], STR_SEP, TIME2FILE, TRIM, UNIQ [1]
UNIQ [2], UNIQ [3], concat_dir [4], curdir [1], curdir [2], data_chk [1]
data_chk [2], file_exist [1], file_exist [3], file_list [1], file_list [2]
get_logenv [1], get_logenv [2], sel_timrange [1], sel_timrange [2], ssw_strsplit
str_replace [1], str_replace [2], tim2dset [1], tim2dset [2], time_window
timegrid, ut_time [1], ut_time [2], wc_where [1], wc_where [2]
CALLED BY:
EIT_COPY, EIT_GBO, EIT__DEFINE, eit_fulldisk, eit_genx_cat, eit_mirror
eit_proton_summary, lasteit, ssw_track_demo, trace_movie_context
Restrictions:
assume environmental EIT_QUICKLOOK points to collapsed tree or top of tree
History:
15-jan-1996 (S.L.Freeland) for LPARL access
14-feb-1996 (S.L.Freeland) allow uncollapsed tree , LAST keyword
3-Jun-1996 (S.L.Freeland) check EIT_QUICKLOOKnn
7-aug-1996 (J. Newmark) add lz (level-zero) keyword and associated
lz queries, allow one or more QUICKLOOK archives
25-mar-1997 (S.L.Freeland) - add COLLAPSE keyword and function
recurse w/collapse when required
add T0 and T1 parameters
6-apr-1997 (S.L.Freeland) - trim files to true t0 and t1 if supplied
1-may-1997 (S.L.Freeland) - fix problem if auto-collapse search tried
(recursive logic)
3-jun-1997 (S.L.Freeland) - add /FULLDISK and WAVE keywords
(required $EIT_FULLDISK defined)
11-jun-1997 (S.L.Freeland) - fix of "minor" 3-jun bug (/full)
26-Jun-1997 (S.L.Freeland) - add MATCHFILES(in) and DELTA(out)
13-oct-1997 (S.L.Freeland) - fixed problem with LAST
8-Apr-1998 (S.L.Freeland) - use current directory if topdir not set
10-sep-1998 (S.L.Freeland) - permit only t0
15-feb-1999 (S.L.Freeland) - use qkl if lz not found
23-apr-1999 (S.L.Freeland) - fix problem for small time ranges
which span a month (subdirectory) boundry
Allow time vector as t0 (use matchfiles logic)
2-Jan-1999 slf - quick y2k fix (still need to look at year spanning)
15-Aug-2001 Zarro (EITI/GSFC) - modified to return blank string when
no files found
1-Sept-2001, Zarro (EITI/GSFC) - added /QUIET
14-May-2003 - S.L.Freeland - replace 'strsplit' with ssw_strsplit
[Previous]
[Next]
Name: eit_fill_cube
Purpose: fill missing data with data from neighbors
Input Parameters:
cube - data cube (output if only one parameter passed)
Output Parameters:
ocube - optional - output (filled cube) - default replaces input
Keyword Parameters:
missing - if set, pixel value considered 'missing' (default=0)
times - vector of image times (not yet implemented)
track_progress - switch, if set, print correction status info
window - limit window of neighbors to check (+/- limit)
(for large cubes, dont wander further than required
Output Parameters:
orig_miss - percent missing data in original cube [fltarr(nimages)]
final_miss - percent missing data in final cube
Calling Sequence:
eit_fill_cube,cube ; overwrite input with filled version
eit_fill_cube,cube,ocube ; filled version in ocube (orig saved)
Calling Examples:
eit_fill_cube,cube, orig=orig, final=final
eit_fill_cube,cube, mising='aa'x ; AA(hex) is "missing" (ex: fill data)
eit_fill_cube,cube, window=2 ; check neghbors within +2/-2
CALLED BY:
eit_fulldisk, eit_proton_summary, lasteit
Restrictions:
Assumes images are kindof like their neighbors
(normalization, wavelenth, structure, registration, time, whatever)
History:
17-feb-1996 S.L.Freeland
18-Feb-1996 S.L.Freeland - add WINDOW keyword and function
Nothing about this is EIT specific (fills any data cube) so it might->GEN
CALLS: ***
data_chk [1], data_chk [2]
Side Effects:
Input array is filled ('clobbered') if only one parameter (for space)
[Previous]
[Next]
NAME:
eit_fillmissingblocks
PURPOSE
a fast removal of missing blocks designed for calibration lamps
INPUTS
img = input image
OUTPUTS
returns the image with missing blocks replaced
KEYWORDS
bin = binning of the input image. must be 16 (2x2 binned image) or 32 (full res)
PROCEDURE
replace missing blocks with the nearest valib blocks
CREATION
18-May-2001. F. Auchere
MODIFICATIONS
12-Aug 2002. F. Auchere bug fixed
28-jul-2003. F. Auchère changed replacement by nearby missing block by
a bilinear interpolation
CALLED BY
eit_find_last_cal, eit_getcal
[Previous]
[Next]
NAME:
eit_find_last_cal
PURPOSE
finds the latest group of calibration lamps
INPUTS
date = date of observation
OUTPUTS
returns the latest averaged calibration lamp
KEYWORDS
none
CREATION
12-Dec-2001. F. Auchere
MODIFICATIONS
12-Aug-2002. F. Auchere : added fill in of missing blocks.
11-Jul-2003. J. Newmark : Windows compatibility
CALLS:
CALLED BY
eit_getcal
[Previous]
[Next]
NAME:
eit_findcalgroup
PURPOSE
finds an averages groups of calibration lamps taken within 24 hours
INPUTS
date = date of observation
dir = search direction. dir>=0 : forward, dir<0 : backwards
OUTPUTS
returns an averaged calibration lamp
KEYWORDS
start_date = date of the first calibration lamp
end_date = date of the last calibration lamp
CREATION
18-May-2001. F. Auchere
MODIFICATIONS
2-Jun-2001. F. Auchere : Completely re-written in order to use the new cal lamps data base
17-Sep-2001. F. Auchere : Bug fixed
11-Jul-2003. J. Newmark : Windows compatibility
CALLS:
CALLED BY
eit_getcal
[Previous]
[Next]
NAME:
eit_flat
PURPOSE:
CALLING SEQUENCE:
flat_fielded_img = eit_flat(img, fits_header, image_no=image_no)
INPUTS:
img = Raw image, detector offset subtracted and missing pixels raised
to baseline level
fits_header = index structure/fits header for img
OPTIONAL INPUTS:
None
KEYWORDS
image_no = image number in LZ file
verbose = set if messages discriminating new read or pre-read wanted
no_copy = don't make dual copies of input image
(only use if output image replaces input image)
CALLED BY:
EIT_IMAGE, EIT_PREP, MAKE_IMAGE_EIT img hdr
EXAMPLE:
flat_img = eit_flat(img, fits_header) ; 2-D version
flat_img = eit_flat(img, fits_header, image_no) ; 3-D version
MODIFICATION HISTORY
The epoch - (JSN) - Written
17-Oct-96 - (BNH) - Now sports a document header
- Also sports a common block
1997/04/11 - D.M. - Added VERBOSE keyword.
1998/04/13 - D.M. - Added /VAX_FLOAT keyword in OpenVMS OPENR statement
for IDL 5.1 compatibility.
3-Jun-1998 - Zarro (SAC/GSFC)
- Added calls to temporary for improved
memory management.
28-Sep-1999 - N.Rich - adapt for LASCO fits header: object, wave, DATE-OBS
18-May-2000 - J.S.Newmark - online LASCO compatibility
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
DATATYPE [1], DATATYPE [2], DATATYPE [3], EIT_FXPAR, IS_GSFCVMS, concat_dir [4]
COMMON BLOCKS:
eit_flat_blk = filename and flat field for last image flat-fielded (to
speed up operation)
[Previous]
[Next]
NAME:
eit_flux
PURPOSE:
Specify Te + filter + DEM ==> compute DN
Specify Te + filter + DN ==> compute DEM
CALLING SEQUENCE:
dn_s = eit_flux(171,1.e6,filter='CLEAR',edensity=1e9)
or
dem = eit_flux(195,1.5e6,filter='AL+1',edensity=1.5e9,dn=1000,
expotime=5)
INPUTS:
wav_band: EIT bandpass, one of 171,195,284,304
te: Temperature, may be a vector
OUTPUTS:
This function returns DN/s/pixel.
If DN keyword is set, will return Column Emission Measure cm^-5
this value represents an isothermal emission measure
approximation. We include a factor 0f 0.7 to account
for the Potasch approximation that lines really form
over +/- 0.15 dex. If one wanted to compare this to a
differential emission measure (DEM) then you should
multiply this by 1./0.3
OPTIONAL INPUT KEYWORDS:
DN: if set returns Emission Measure, if vector should correspond
to elements of temperature
filter: default is CLEAR position can also specify AL+1 or AL+2
wavestep: wavelength bin for spectrum calculation, default=0.2Ang
Specify one of below:
Edensity: electron density in emitting region (cm^-3)
Pressure: electron pressure (cm^-3 K)
Note- these can be single valued or arrays corresponding to
the specified temperatures
these are used in the population calculations
Can specify both of below:
inspec: an input spectrum, useful if have spectra already calculated
waves: wavelengths corresponding to inspec
dem_name = differential emission measure file - temp vs dem, can
also be a float value=alog10(dem) which correspond to temps.
set to "unit" if just need unit emission measure - UNIT=1e26
wmin = minimum wavelength of interest (default = 165A)
wmax = maximum wavelength of interest (default = 350A)
expotime = exposure time if data set not normalized
OPTIONAL OUTPUT KEYWORDS:
CALLS: ***
DATATYPE [1], DATATYPE [2], DATATYPE [3], EIT_PIXSIZE, PICKFILE, READ_DEM, SPLINE
eit_parms, ssw_packages
CALLED BY:
eit_temp
RESTRICTIONS:
**** This routine requires the CHIANTI package. Please run the
following command before using this routine:
ssw_packages,/chianti,/prepend
Presently this reads the Instrument Data from the EIT Sensitivity
files on the SSW as provided by J. Newmark. This section is
modularized so as to facilitate updates or even new instruments.
This program runs fairly slow since it is calculating the spectrum
from 165-350 A. A better way to run this is to calculate the spectrum
using ISOTHERMAL once for the density and temps (or pressure) you wish
and pass it into this program.
MODIFICATION HISTORY:
1996 Oct 30, J.S. Newmark, Written.
1998 Jan 26 J.S. Newmark, V. 2.0 delete wavestep in
calculation - should not be included, auto add chianti
to path
1999 May 11 J.S. Newmark, V 3.0 New calibration
1999 Jun 08 J.S. Newmark, Clarify return value of Emission measure
1999 Nov 17 J.S. Newmark, V 3.1 New Calibration
2000 Sep 26 J.S. Newmark, change chianti_spec -> isothermal,/cont
[Previous]
[Next]
Name: eit_fulldisk
Purpose: read, return and optionally display Full Disk EIT images
Input Parameters:
t0 - start time of interest
t1 - stop time of interest
Output Parameters:
index - EIT instrument parameters (structures from read_eit)
data - EIT image or image cube
Keyword Parameters:
/w171, /w195, /w284, /304 - desired wavelength (default=171)
/wall - show all four wavel.
/norebin - return full size (default reduced to 512x512)
rebin - size to rebin to (assumed sqare!)
/notv - just return the data
/display - if set, display data (via xstepper(
nmissb - maxmimum number of missing blocks to accept
quality - (empirical data check) minimum quality "accepted"
(see data_quality.pro, or, if set as flag, default=90%)
sizes - sizes to accept (ex: 256,512,1024) - default=all full fov
hours - if set, limit cadence to this number of hours
(see grid_data.pro and timegrid.pro for generic details)
count - number of files found matching criteria
Calling Sequence:
eit_fulldisk, t0, t1, index, data [, /w171, /w195, /w284, /304]
Calling Examples:
eit_fulldisk,t0,t1,index,data,/w304 ; 304 full between t0&t1
eit_fulldisk,t0,t1,index,data,hour=4,/w195,nmiss=0 ; 195, perfect only,
; 4-hour cadence
eit_fulldisk,t0,t1,index,data,hour=4,/w195, $ ; same as previous
nmissb=0, sizes=[512,1024] ; but ignore 256x256 etc
eit_fulldisk,t0,t1,index,data,/w171, $ ; good to perfect 171s
nmissb=2,sizes=1024 ; only 1024x1024
eit_fulldisk,t0,t1,index,/w284 ; INDEX only (no data read)
History:
29-may-1997 - S.L.Freeland - from 'lasteit.pro'
(and asap, lasteit will just call this routine)
2-Apr-1998 - S.L.Freeland - added COUNT output keyword
add /LZ search if nothing in qkl
Restrictions:
Full disk only, Files in $EIT_QUICKLOOK
Calls: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], DELVARX [1], DELVARX [2]
DELVARX [3], DELVARX [4], WDEF [1], WDEF [2], anytim [1], anytim [2], anytim [3]
anytim [4], anytim [5], break_file [4], data_chk [1], data_chk [2], data_quality
delvarx [5], eit_colors, eit_dark, eit_file2path, eit_files, eit_fill_cube
get_eit_info, grid_data [1], grid_data [2], read_eit, reltime [1], reltime [2]
timegrid, uniqo [1], uniqo [2], ut_time [1], ut_time [2], xstepper [1], xstepper [2]
Common Blocks:
eit_fulldisk_blk - speed up subsequent wavelenths for same time interva;
[Previous]
[Next]
Name: eit_fulldiskdb
Purpose: generate EIT full disk data base (calls eit_fulldisk)
Input Parameters:
t0, t1 - start and stop time
Keyword Parameters:
copy - switch - if set, COPY the files to $EIT_FULLDB/YYYYMM
links - switch - if set, symbolic links only
nmissb - missing block limit (see eit_fulldisk.pro)
hours - limit cadence in hours (see eit_fulldisk.pro, grid_data.pro)
sizes - limit sizes (256, 512, 1024, etc. see eit_fulldisk.pro)
waves - optional wave array (default=[171,195,284,304]
eit_prep - if set, prep the output (overrides COPY and LINKS)
History:
29-may-1997 - S.L.Freeland, written
2-April-1998 - S.L.Freeland - check COUNT return by eit_fulldisk
Restrictions:
many - assumed running on a machine with all EIT data online
CALLS: ***
BOX_MESSAGE, CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], EIT_FDSPAWN
FILE_EXIST [2], TIME2FILE, UNIQ [1], UNIQ [2], UNIQ [3], concat_dir [4]
data_chk [1], data_chk [2], eit_file2path, file_exist [1], file_exist [3]
get_logenv [1], get_logenv [2], gt_day [1], gt_day [2], gt_tagval [1]
gt_tagval [2], prstr [1], prstr [2], str_replace [1], str_replace [2], strjustify
timegrid
Side Effects:
may create required subdirectories ($EIT_FULLDB/YYYYMM)
[Previous]
[Next]
Project : SOHO - EIT
Name :
EIT_FXPAR
Purpose :
Obtain the value of a parameter in a FITS header.
Explanation : This is a shell around the IDLASTRO routine FXPAR.
The procedure can parse for keywords embedded into
comments. Usageis identical to FXPAR.
Use :
Result = EIT_FXPAR( HDR, NAME [, ABORT ] [,IMAGE_NO=IMAGE_NO] )
e.g. IDL> hdr=headfits('efr19960715.123746')
IDL> print,eit_fxpar(hdr,'EXPTIME')
7.000
IDL> print,eit_fxpar(hdr,'N_MISSING_BLOCKS')
4
IDL> hdr=headfits('efz19960623.132832')
IDL> print,eit_fxpar(hdr,'FLTR',image_no=23)
Al + 1
Inputs :
HDR = FITS header string array (e.g. as returned by FXREAD). Each
element should have a length of 80 characters
Optionaly, this can be an IDL structure CREATED from a
FITs header as by read_eit
NAME = String name of the parameter to return. If NAME is of the
form 'keyword*' then an array is returned containing values
of keywordN where N is an integer. The value of keywordN
will be placed in RESULT(N-1). The data type of RESULT will
be the type of the first valid match of keywordN found.
If keyword is in COMMENT field will return keyval.
Opt. Inputs :
ABORT = String to output if can not return value
IMAGE_NO = Image number for files containing multiple subimages,
set to 'ALL' to return values for all subimages
Outputs :
The returned value of the function is the value(s) associated with the
requested keyword in the header array.
Keywords :
COUNT = Optional keyword to return a value equal to the number of
parameters found by FXPAR.
COMMENTS= Array of comments associated with the returned values.
INITIALIZE = Used by read_eit only, parse comments for 3D index
structures, else pull out of tag value
CAT_WAVES = return all wavelngths in catalog structure, else
returns first wavelength
Calls : ***
ANYTIM2UTC [1], ANYTIM2UTC [2], DATATYPE [1], DATATYPE [2], DATATYPE [3]
FXPAR [1], FXPAR [2], TAG_EXIST [1], TAG_EXIST [2], UTC2TAI, gt_tagval [1]
gt_tagval [2]
CALLED BY:
EIT_CATRD [1], EIT_DISPLAY, EIT_IMAGE, EIT_POSTAGE, EIT_PREP, EIT_SUB_UTIL_FILE
MK_EIT_L1, MOVIE_MAKER, NRL2EIT, eit_degrid, eit_dump, eit_flat, eit_getlimb
eit_gt_corner, eit_gt_filter, eit_gt_wave, eit_lzsort, eit_norm_response [1]
eit_ratio, eit_readcal, eit_ring_flux, read_eit
Common :
None.
Restrictions:
None.
Side effects:
Category :
Data Handling, I/O, FITS, Generic.
Written : 1996 July 15, J. Newmark
Modified : 1996 Aug 23 J. Newmark allow input if IDL structure (or
array of structures) created from FITS headers.
1996 Oct 18. J. Newmark - allow use of Single image keywords
for multiple image LZ files, return scalar for 3D keywords
if 1 element array.
; 1996 Nov 5, S.L. Freeland - force 'have_sub' to scalar value
1996 Dec 27 J. Newmark - fixed index structure input.
Version : 1.2
1997 Jan 14 J. Newmark - added capability to handle
raw catalog structure input, add /initialize keyword
cat_waves keyword
Version : 2.0
1997 Jan 24 J. Newmark - changed many if-thens to CASE
1997 Mar 20 J. Newmark - add in sci_obj field for catalog
1997 Apr 14 J. Newmark - add support for second ASCII table
1997 Aug 12 J. Newmark - return INT/FLT where expected, 3D
1998 Mar 26 J. Newmark - add corrected date_obs field
1998 Oct 23 J. Newmark - add n_repeat field for raw catalog
[Previous]
[Next]
Name: eit_genx_cat
Purpose: generate and/or read local eit qkl catalog
Input Parameters:
NONE
Output Parameters:
index - IDL 'eit_struct.pro' style structures for last few days
text - text version of index (via get_infox)
Keyword Parameters:
cat_name - (output) - name of current catalog
name_only - (input) - switch, just return catalog name via CAT_NAME
full_fov - (input) - if set, read FULL FOV catalog
generate - (input) - switch, GENERATE catalog (default is just read)
(generally , via CRON job)
last - (input) - only with GENERATE, number of days to include
(most recent NN=LAST days, default=2)
Calling Sequence:
eit_genx_cat, eitstructs, eitinfo ; most recent few days
eit_genx_cat, eitstructs, eitinfo, /full_fov ; same, only full FOV (speed)
eit_genx_cat,/generate ; Make it, via priviledged cron job
CALLED BY:
lasteit
History:
3-mar-1997 - generate/read EIT_STRUCT (via read_eit.pro) based catalog
(replace old "temporary" ascii catalog)
11-mar-1997 - make a smaller full fov also
23-apr-1997 - add filter for zero block files (and send email warning)
Calls: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], STR2ARR [1], STR2ARR [2]
concat_dir [4], eit_files, file_size [1], file_size [2], get_infox, get_logenv [1]
get_logenv [2], mail [1], mail [2], pr_status [1], pr_status [2], read_eit
restgen [1], restgen [2], savegen [1], savegen [2], str_replace [1]
str_replace [2]
[Previous]
[Next]
NAME:
eit_getcal
PURPOSE
computes a calibration lamp for a given date
INPUTS
date_obs = date of observation
OUTPUTS
returns a computed calibration lamp interpolated or extrapolated form the closed calibration lamps
KEYWORDS
none
WARNINGS
still in test
may no return correct values for the most recent bakeout periods
PROCEDURE
interpolates or extrapolates from the closest calibration lamps
CREATION
18-May-2001. F. Auchere
MODIFICATIONS
12-Dec 2001. F. Auchere call to eit_find_last_cal added to avoid
blow up when no recent cal lamp is found
12-Jan 2003. F. Auchere changed a b and c variables type to float
06-Sep 2004. F. Auchere correct calls to eit_fillmissingblocks
CALLS:
CALLED BY
eit_norm_response [1]
[Previous]
[Next]
Name: eit_getlimb
Purpose: get the limb parameters from EIT FITS or fit limb if not present
Input Parameters:
eitfiles - EIT FITS file list
Output Parameters:
xx - sun center, full res EIT pixels
yy - sun center, full res EIT pixels
rr - S radius, full res EIT pixels
Keyword Parameters:
fit - switch, if set, call fit_limb (default if no FITS params)
find - switch, if set, call find_limb
guess - switch, if set, return guess (for fast prototyping)
update - switch, if set, update FITS file (add point info to header)
re_fit - switch, if set redo fitting even if in header
Calling Sequence:
eit_getlimb, eitfiles, xx,yy,rr [/fit, /find, /guess, /update]
Method:
get parameters from FITS file, fit_limb or find_limb , or guess
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3]
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], EIT_FXPAR, EIT_PIXSIZE
FXADDPAR [1], FXADDPAR [2], FXPAR [1], FXPAR [2], FXWRITE [1], FXWRITE [2]
HEADFITS [1], HEADFITS [2], HEADFITS [3], READFITS [1], READFITS [2], READFITS [3]
break_file [4], concat_dir [4], ssw_strsplit
CALLED BY:
sxt2eit
History:
1-Feb-1996 - S.L.Freeland
15-Feb-1996 - S.L.Freeland - minor adjustments
3-Sep-1996 - J. Newmark -updated for whole mission paramters
14-May-2003 - S.L.Freeland - replace strsplit with ssw_strsplit
[Previous]
[Next]
Name: eit_getobs
Purpose: interactive selection of EIT obsserving blocks (32x32)
Input Parameters:
NONE
Output Parameters:
obs - command blocks selected [** currently 0-1023 from lower left **]
Keyword Parameters:
magnification - image (size/1024.) - default derived from current window
Calling Sequence:
eit_getobs, obs [magnification=magnification]
Method:
read cursor postion from current X-window containing an EIT or
EIT-aligned image - permit mouse selection of EIT command regions
"Press LEFT button to SELECT or UNSELECT a region"
"Press/Drag CENTER button to SELECT or UNSELECT adjacent regions"
"Press RIGHT button to EXIT"
History:
30-jan-1996 (S.L.Freeland)
CALLS: ***
MORE [1], MORE [2], WDEF [1], WDEF [2], grid, strjustify, tbeep [1], tbeep [2]
tbeep [3]
Restrictions:
assumes EIT or (EIT co-aligned image) is in current X windows display
[Previous]
[Next]
Name:
eit_gt_corner
Purpose:
Return the defining corners of an EIT image.
Input Parameters:
header = FITS header or EIT index structure for image
Output:
corners = [P1_X, P2_X, P1_Y, P2_Y]
Calling Sequence:
corners = eit_gt_corner(index)
corners = eit_gt_corner(fits_header)
CALLS: ***
EIT_FXPAR, data_chk [1], data_chk [2]
Restrictions:
The one defining reason for this routine is to work around the
COMMENT header definitions. If P[12]_[XY] ever were to become
real FITS header keywords instead of comments, it is my hope that
we would only have to fix this one spot to deal with it.
History:
31-October-96 - (BNH) - Written
[Previous]
[Next]
Name:
eit_gt_filter
Purpose:
Return the filter of a 2- or 3-D index structure or
FITS header.
Input Parameters:
header = FITS header or EIT index structure for image
Input Keywords:
image_no = Image number from within above header
Output:
Filter of chosen image
Calling Sequence:
filter = eit_gt_filter(fits_header)
filter = eit_gt_filter(fits_header,image_no=7)
filter = eit_gt_filter(index [,image_no=4])
Running Commentary:
Another workaround for fits keywords disguised as comments.
CALLS: ***
EIT_FXPAR, data_chk [1], data_chk [2], gt_tagval [1], gt_tagval [2]
History:
3-November-96 - (BNH) - Written
[Previous]
[Next]
Name:
eit_gt_wave
Purpose:
Return the wavelength of a 2- or 3-D index structure or
FITS header.
Input Parameters:
header = FITS header or EIT index structure for image
Input Keywords:
image_no = Image number from within above header
Output:
wavelength of chosen image
Calling Sequence:
wave = eit_gt_wavelnth(fits_header)
wave = eit_gt_wavelnth(fits_header,image_no=7)
wave = eit_gt_wavelnth(index [,image_no=4])
Running Commentary:
Another workaround for fits keywords disguised as comments.
CALLS: ***
EIT_FXPAR, EIT_GT_WAVELNTH, data_chk [1], data_chk [2], gt_tagval [1]
gt_tagval [2]
History:
3-November-96 - (BNH) - Written
[Previous]
[Next]
NAME:
EIT_IMAGE
PURPOSE:
Read a SOHO Extreme ultraviolet Imaging Telescope (EIT) image
from a FITS file, and return the byte-scaled image.
CALLING SEQUENCE:
scaled_image = EIT_IMAGE(filename, [original_image =
original_image,] [scale_factor = scale_factor,]
[top = top], [flat=flat],
[image_time = image_time,] [header = header],
[show_image = show_image], [replacement_value =
replacement_value], [fill_block = fill_block],
[surround = surround], [adjust = adjust],
[dark_current = dark_current], [v_offset = v_offset],
[no_degrid = no_degrid], [half = half],
[quarter = quarter], [no_rescale = no_rescale],
[xpos = xpos, ypos = ypos], [image_no=image_no])
INPUTS:
FILENAME = String containing the name of the FITS file to be read.
Optionally, this can instead be the string returned
by EIT_CATRD for a file.
OUTPUTS:
Result = byte-scaled image array (usually 1024 x 1024)
;
KEYWORDS:
show_image - if present and non-zero, display the image on the
currently selected graphics device
original_image - original, floating-point version of contents of file
scale_factor = scaling for logarithmic scaling on the current X
device
top - maximum value to which to scale (useful when scaling a series
of images to the same maximum)
image_time = string concatenating DATE-OBS and TIME-OBS from the
FITS header
header = FITS header (string array)
replacement_value - DN with which to replace missing data block pixels
fill_block - intarr(32, 32) with which to replace missing data blocks
surround - if set, replace missing block data with the average of
surrounding blocks
adjust - multiplicative factor by which to adjust for strange exposure
times
dark_current - dark current value to be subtracted when image is
scaled (Note: if replacement_value is defined and
dark_current is not, the dark current is set to the
replacement value.)
v_offset - if show is set, the vertical offset for display (useful
for displays incapable of displaying the full image)
no_degrid - if set, do not remove the grid, otherwise remove by
the method of F. Clette
half - if set, rebin image to half size
quarter - if set, rebin to quarter size
xpos = x position of image window
ypos = y position of image window
image_no = selected image from a multiple image LZ file
flat - if set flat field the image
no_rescale - set if image is not to be rescaled into 0 - 255 range
CALLED BY:
EIT_DISPLAY, MOVIE_MAKER, eit_dump
EXAMPLE:
Find the most current EIT image on the EIT image directory,
and display it on the current device:
s = findfile("reformat_dir:ef*.*")
n = n_elements(s)
a = EIT_IMAGE(s(n - 1), /show_image)
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], EIT_FXPAR, EIT_RESCALE, HEADFITS [1], HEADFITS [2]
HEADFITS [3], LOADCT, RAISE_MISSING_BLOCKS, READFITS [1], READFITS [2]
READFITS [3], eit_degrid, eit_file2path, eit_flat, eit_norm_response [1]
eit_norm_response [2], is_fits [1], is_fits [2]
RESTRICTIONS:
None.
NOTES:
None.
PROCEDURES USED:
READFITS (Astronomy library)
RAISE_MISSING_BLOCKS
EIT_DEGRID
EIT_FLAT
MODIFICATION HISTORY:
D.M. fecit, from SXT_IMAGE, 25 October, 1993.
Changed to ignore v_offset unless n_y eq 1024. D.M. fecit, 1996/02/25.
Call EIT_DEGRID instead of reinventing the wheel. D.M. fecit,
1996/03/30.
1996-Aug-20 JSN added support for 3-D LZ data
Changed filename argument handling to accept "brief" catalog lines.
D.M. fecit, 1996/09/01.
1996-Sep-25 J Newmark added flat fielding
1996/10/20 Changed scaling to use only southernmost 15/16
of a full-frame image (to avoid scaling 284 A,
clear filter pos. images to the light leak).
D.M. fecit.
1996/12/04 J. Newmark - removed last traces of find_keyword
1996/12/23 Added QUARTER keyword. D.M. fecit.
1997/02/05 Added NO_RESCALE keyword. D.M. fecit.
1998/03/27 J Newmark change code to use NSLICE keyword in
READFITS routine for 3D images
1998 Jun 3 - D. Zarro (SAC/GSFC)
- sprinkled liberally with calls to temporary
2000 Aug 31 J. Newmark - include flat fielding
2001 Jun 4 J. Newmark - only use original grid
2003 July 17 - A. Young - compensate for spacecraft roll
[Previous]
[Next]
NAME:
eit_lampexpt
PURPOSE
corrects the calibration lamp exposure time for the response time of the lamp's circuit.
INPUTS
t = exposure time as recorded in the header
OUTPUTS
returns the corrected exposure time
KEYWORDS
none
CREATION
29-Jul-2002. F. Auchere
MODIFICATIONS
CALLED BY
eit_readcal
[Previous]
[Next]
Name: eit_lzsort
Purpose: get index of subimages from a level zero FITS file which
match selected wavelength
Input Parameters:
filename - EIT level zero FITS file
wave - specified wavelength of subimages
Output Parameters:
Function return indices of subimages matching specified wavelength
make_array - set if want to return data
img - if make_array specified returns subimage data
Keyword Parameters: none
Calling Sequence:
indices=eit_lzsort(filename,wave [,img=img, /make_array])
Method:
Check if FITS file, call eit_fxpar to read in subimage info
returns data cube if asked for via:
IDL> img=readfits(file)
IDL> img=img(index)
CALLS: ***
EIT_FXPAR, HEADFITS [1], HEADFITS [2], HEADFITS [3], READFITS [1], READFITS [2]
READFITS [3], is_fits [1], is_fits [2]
CALLED BY:
MOVIE_MAKER
History:
18-July-1996 - J. Newmark - created
[Previous]
[Next]
Name: eit_mirror
Purpose; make and submit EIT quicklook mirror job
6-feb-1996 (S.L.Freeland)
CALLS:
[Previous]
[Next]
NAME:
EIT_MKMOVIE
PURPOSE:
Create a GIF or MPEG movie from an Input FITS file list or byte-scaled
data cube
CATEGORY:
Animation
CALLING SEQUENCE:
eit_mkmovie,input [,outsize=outsize][,/mpeg][,/gif][,wave=wave][movie_name=name]
INPUTS:
list I = listing of LZ or QKL data files or a bytescaled datacube as
may come from an existing IDL save set.
OPTIONAL INPUT KEYWORD PARAMETERS:
wave I = wavelength
outsize I = 2 dimensional array for output size of image
movie_name I = output movie name (default derived from input data)
gif I = select if want gif animation
mpeg I = select if want mpeg animation
java I = select if want javascript animation
outidr I = set for specified output directory
OUTPUTS:
Creates either a GIF or MPEG file
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
DATATYPE [1], DATATYPE [2], DATATYPE [3], EIT_PREP, FILE_EXIST [2], LOADCT, RM_FILE
WDISPLAY, XLOADCT [1], XLOADCT [2], XLOADCT [3], XMESSAGE, concat_dir [4]
curdir [1], curdir [2], data2files, eit_colors, eit_file2path, file_exist [1]
file_exist [3], get_logenv [1], get_logenv [2], image2movie, is_fits [1]
is_fits [2], read_eit
COMMON BLOCKS: none.
SIDE EFFECTS: temporarily creates GIF files for each image
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
Written by: J. Newmark Date. September 1997
1999 February 11 J. Newmark Add Java, ooutdir, clean up
[Previous]
[Next]
NAME:
eit_norm_response
PURPOSE
compute a normalized response for each eit band
INPUTS
date_obs = date of observation
inwave = wavelength of observation
fits_header = header to obtain pointing/size information
OUTPUTS
returns normalization factor to be applied to measured DN level
on given date such that all EIT data is normalized to 2-Feb-1996.
KEYWORDS
no_index = set if do not want to use recommended Index correlation
for response scaling
verbose = set for testing
WARNINGS
this routine may not always have the most recent normalizations
especially just following a bakeout.
PROCEDURE
Uses EUV flat field determined by scaling from Calibration lamps
CREATION
17-May-2001. J. Newmark and F. Auchere
MODIFICATIONS
10-Dec-2001. J. Newmark inclusion of trend corrections
28-Feb-2002. J. Newmark inclusion of wavelength dependence
of flat field
11-Jul-2003. J. Newmark : Windows compatibility
02-Sep-2004. F. Auchère : added /nan keyword to total() when commputing
the degrad parameter
30-May-2006. J. Newmark : Set extrapolation after last calibrated day equal
to last calibrated day
CALLS:
CALLED BY
EIT_IMAGE, EIT_PREP, MAKE_IMAGE_EIT img hdr
[Previous]
[Next]
NAME:
eit_norm_response
PURPOSE
compute a normalized response for each eit band
INPUTS
date_obs = date of observation
inwave = wavelength of observation
fits_header = header to obtain pointing/size information
OUTPUTS
returns normalization factor to be applied to measured DN level
on given date such that all EIT data is normalized to 2-Feb-1996.
KEYWORDS
no_index = set if do not want to use recommended Index correlation
for response scaling
verbose = set for testing
WARNINGS
this routine may not always have the most recent normalizations
especially just following a bakeout.
PROCEDURE
Uses EUV flat field determined by scaling from Calibration lamps
CREATION
17-May-2001. J. Newmark and F. Auchere
MODIFICATIONS
10-Dec-2001. J. Newmark inclusion of trend corrections
28-Feb-2002. J. Newmark inclusion of wavelength dependence
of flat field
11-Jul-2003. J. Newmark : Windows compatibility
02-Sep-2004. F. Auchère : added /nan keyword to total() when commputing
the degrad parameter
CALLED BY
EIT_IMAGE, EIT_PREP, MAKE_IMAGE_EIT img hdr
[Previous]
[Next]
NAME : EIT_PA
PURPOSE : Interface used to select specific files from EIT catalogs
CATEGORY : EIT subfield
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_PA_ev
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS :
KEYWORDS : none
OUTPUTS : none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], EIT_CATRD [1], EIT_CATRD [2], EIT_PA_UTIL
EIT_PA_ev, XMANAGER, eit_catrd [3]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_PA_ev
PURPOSE : Event handler for EIT_PA
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_PA
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], EIT_CATRD [1], EIT_CATRD [2], EIT_PA, EIT_PA_UTIL
XMANAGER, eit_catrd [3]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_PA_PROCESS
PURPOSE : Parse info from file header, display, and set up coord system
CATEGORY :
EXPLANATION :
SYNTAX :
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS :
Cat_text text array containing catalog output of the selected files
KEYWORDS :
OUTPUTS : none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], LOADCT, READFITS [1], READFITS [2], READFITS [3]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
2001/06/13 modified for EMS movie files
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_PA_UTIL
PURPOSE : Routines to handle various user actions from EIT_PA_EV
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_PA, EIT_PA_ev
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (EIT)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_PA_VIEW
PURPOSE : No frills display of EIT image
CATEGORY :
EXPLANATION :
SYNTAX :
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
LOADCT, READFITS [1], READFITS [2], READFITS [3]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME:
eit_parms
PURPOSE:
Generally used as a sub-routine in EIT_FLUX, EIT_TEMP. This
routine reads the instrument calibration parameters.
CALLING SEQUENCE:
Area = eit_parms(waves,band,filter,text,header,[units=units],
INPUTS:
waves: wavelength array of intest
band: bandpass (e.g. 171,195,284,304)
filter: Filter wheel selection, default = 'CLEAR=clr',
can select 'al1', 'al2'
OUTPUTS:
area: The effective area per wavelength, default units =
cm^2 dn sr Ang photons^-1
OPTIONAL INPUT KEYWORDS:
OPTIONAL OUTPUT KEYWORDS:
text: returns text description of calibration
header: returns file header information
units: returns output units
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], INTERPOL, MINMAX [1], MINMAX [2]
READFITS [1], READFITS [2], READFITS [3], concat_dir [4]
CALLED BY:
eit_flux
RESTRICTIONS:
The SOHO EIT calibration is presently being upgraded. This file
may change. Updates will be posted on EIT homepage.
MODIFICATION HISTORY:
1999 May 11, J.S. Newmark, Written.
1999 November 12, J.S. Newmark, New Calibration -
Dere et al. Sol.Phys. 2000
Modifications from this include far wings and
second order He II in 304 bandpass. These modifications
have almost negligible effect on final output.
Graphically see calibration, Dere = solid, extra = dashed:
http://umbra.nascom.nasa.gov/eit/images/calib.ps
[Previous]
[Next]
NAME:
EIT_PIXSIZE
return pixel scale, i.e. arcsec/pixel
"/pix=(0.206265/focal length)*pixel size microns
=(0.206265/1.652)*21=2.62
However, a problem arises with comparison of features with other
SOHO instruments. Part of the problem is caused by the physical
size difference of the solar disk in the photoshphere and transition
region and corona.
Reference: Auchere, DeForest, Artzner, Ap. J. Letters accepted, 2000.
CALLED BY
EIT_CATRD [1], EIT_PREP, NRL2EIT, eit_flux, eit_getlimb, fixinate_eit, read_eit
sxt2eit
[Previous]
[Next]
EIT_POINT
Purpose - return EIT pointing = Sun center, for a specific date and
wavelength
Inputs: date, wavelength (uses average if not defined)
Option Inputs: bin - if set uses binning (default = 1)
Output: XCEN, YCEN in coordinate system starting at (0,0)
Usage:
xy = eit_point('23-dec-1996',195)
print,xy
505.130 514.410
Created: 17 Feb 1998. J. Newmark
Modified: 23 Jul 1998 J. Newmark - completely redone
Modified: 17 Aug 1998 J. Newmark - new fit
Modified: 04 Dec 1998 J. Newmark - added pointing after recovery
- added bin keyword
Modified: 03 Feb 1999 J. Newmark - new fit
Modified: 06 Aug 1999 J. Newmark - new fit
Modified: 20 Oct 1999 J. Newmark - new fit
Modified: 02 Feb 2000 J. Newmark - new fit
Modified: 19 Sep 2000 J. Newmark - new fit
Modified: 28 Nov 2000 J. Newmark - new fit
Modified: 29 Nov 2001 F. Auchere - new fitting function : fourier series and handles arrays
CALLS:
CALLED BY
EIT_PREP, read_eit
[Previous]
[Next]
NAME :
EIT_POSTAGE
PURPOSE : display small images of whole day on one screen
CATEGORY : image summary
EXPLANATION :
SYNTAX : eit_postage
EXAMPLES : none
CALLED BY :
CALLS TO : put, eit_catrd, label_image, eit_file2path, readfits
ARGUMENTS :
INDATE - string of the day to display
KEYWORDS :
LZ - if set then display level zero data, not quicklook data
GIF - if set then output a GIF file
USE_CT - if set creates a color table for each wavelength
compare - if set, use half width of [my] screen
OUTPUTS : IF GIF set then a file name eit_YYYY-DD-MM.gif
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], CONGRID [1], CONGRID [2], CONGRID [3]
DEVICELIB [1], DEVICELIB [2], EIT_CATRD [1], EIT_CATRD [2], EIT_FXPAR, IMAGELIB
LABEL_IMAGE, LOADCT, PUT, READFITS [1], READFITS [2], READFITS [3], WRITE_GIF
eit_catrd [3], eit_file2path
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS:
If USE_CT keyword is set will mess up the color tables. User
should reset them by running loadct,0 after using this routine
PROJECT : SOHO - EIT
HISTORY : V1, Jeff Newmark, 24-jan-1997
1997 Feb 21 - J. Newmark use multiple color tables
1997 Mar 17 - E. Einfalt add titles
1997 Feb 10 - Einfalt added compare
CONTACT : newmark@eitv2.nascom.nasa.gov
[Previous]
[Next]
NAME:
EIT_PREP
PURPOSE:
Process an EIT image.
Processing routine for preparing Level-1 EIT data.
Presently this version reads in the named FITS
image and produces a background subtracted, degridded,
flat-fielded, and degradation corrected output array or FITS file.
In the future this routine will also vignette correct.
The input may also be the output index structure and data cube as
produced by read_eit.
CATEGORY:
FITS processing
CALLING SEQUENCE:
EIT_PREP,filename,[outheader,outimage] [data=data] [,/outfits]
[,/verbose] [,/cosmic] [,/no_calibrate] [,/no_index]
[,/save_zero] [,/nrl] [n_block = n_block]
[,/surround] [,/nodata] [,outdir = outdir]
[,image_no = image_no]
INPUTS:
filename - FITS file name or the output of eit_catrd
may also be index structure (e.g. output of read_eit)
OPTIONAL INPUT KEYWORD PARAMETERS:
outfits - set if want to produce an output FITS file
(or set string filename)
outdir - Destination directory for prepped FITS files
(if outdir is set, /outfits is assumed)
nrl - set if input file was processed at NRL
surround - set to replace missing block data with the average of
surrounding blocks.
image_no - Image number(s) (or 'ALL') to process in multiple image LZ file.
fill_block - intarr(32,32) which is used to replace missing data blocks
save_zero- set to retain images that have only 0 counts.
verbose - set for a few more messages
data - data array(cube) corresponding to index structure
nodata - set if don't want output array (auto set if lt 3 params)
cosmic - set if want to remove cosmic rays - may remove small
real features
no_calibrate - set if want "raw" returned images, i.e. only background subtracted
no_index - set if don't want to use recommended Index correlation
for final response scaling
no_roll - do not correct for SOHO roll (mainly for testing)
OPTIONAL OUTPUTS:
outheader- FITS header or output index structure
outimage - Processed EIT image. Default data type is Float
OPTIONAL OUTPUT KEYWORD PARAMETERS:
n_block - Number of blocks repaired
CALLS: ***
ANYTIM2CAL, ANYTIM2UTC [1], ANYTIM2UTC [2], CONCAT_DIR [1], CONCAT_DIR [2]
CONCAT_DIR [3], COSMICR, DATATYPE [1], DATATYPE [2], DATATYPE [3], DELVARX [1]
DELVARX [2], DELVARX [3], DELVARX [4], EIT_FXPAR, EIT_PIXSIZE, FILE_STAT [1]
FILE_STAT [2], FXADDPAR [1], FXADDPAR [2], FXPAR [1], FXPAR [2], HAVE_TAG
HEADFITS [1], HEADFITS [2], HEADFITS [3], IS_STRING, NRL2EIT, RAISE_MISSING_BLOCKS
READFITS [1], READFITS [2], READFITS [3], REP_TAG_VALUE, ROT_FITS_HEAD
WRITEFITS [1], WRITEFITS [2], concat_dir [4], delvarx [5], eit_dark, eit_degrid
eit_file2path, eit_flat, eit_norm_response [1], eit_norm_response [2], eit_point
file_stat [3], get_logenv [1], get_logenv [2], is_fits [1], is_fits [2]
str_copy_tags [1], str_copy_tags [2], str_copy_tags [3], str_copy_tags [4]
struct2fitshead
CALLED BY:
EIT_MKMOVIE, EIT_SUB_MOVIE, MAKE_IMAGE_EIT img hdr, MKMOVIE, MK_EIT_L1, MK_EIT_MAP
MK_GIF, MK_IMG, eit_proton_summary, eit_ratio, eit_ring_flux, lasteit
ssw_fov_context, ssw_track_demo, todays_targets, trace_movie_context
COMMON BLOCKS: none.
SIDE EFFECTS:
If /outfits is set, then a FITS file is written.
RESTRICTIONS:
1) Processes only one 3-Dim file at a time, either a single image
vector, or the entire file.
2) Does not process a mixture of 2D and 3D files
3) If filename is a vector, there is no checks to make sure
that all the images are the same size. User must be careful
to check this.
PROCEDURE:
The FITS file is read in. Any missing blocks are replaced
with the present dark noise value. The dark noise is
subtracted off the entire image. The image is degridded.
The image is flat fielded and corrected for degradation.
The image is optionally returned.
A new FITS file with an updated header is written out
with the format eit_l1_yyyymmdd_hhmmss.fits
Note: the De-gridding and Flat-fielding are multiplicative
operations with an AVERAGE value over the full field
of 1.0
MODIFICATION HISTORY:
Written by: J. Newmark Date. July 1996
1996 Sep 24 - J. Newmark - add flat fielding
1996 Oct 8 - J. Newmark - allow NRL format files
1996 Oct 16, J. R. Lemen, Changed calling arguments.
1996 Dec 09 - J. Newmark - add usage of image_no throughout.
1997 Feb 10 - J. Newmark - allow input if index structure and data
cube
1997 Mar 09 - J. Newmark - write out FITS files from indx/data
structures
1997 Mar 10 - J. Newmark - fix output FITS file names.
1997 Apr 09 - J. Newmark - account for new exptime
1997 Apr 12 - J. Newmark - account for correction of Al+1 to ClEAR
1997 Aug 14 - J. Newmark - implemented Handy requests of outdir, use
of nodata
1997 Aug 18 - J. Newmark - updated conversion of struct->fitshead
fix up code for /surround
clean up code for output file name
1997 Oct 06 - J. Newmark - added preliminary response normalization
1997 Oct 30 - J. Newmark - update header for LZ Partial Fov images
1997 Nov 10 - J. Newmark - add new grids after July 1997
Version 2.0
1998 Mar 27 - J. Newmark - support ALL or vector images in single 3D
file, speed up 3D processing,
clean up history, pathway checks,
update documentation, fix n_block bug
added adjust keyword
1998 Jun 3 - D. Zarro (SAC/GSFC)
- sprinkled liberally with calls to temporary
1998 Jul 22 - J. Newmark - add phot_norm keyword
1998 Jul 23 - J. Newmark - add cosmic keyword
1998 Dec 08 - J. Newmark - change adjust keyword
2000 Jan 05 - J. Newmark - new response normalization, deleted
phot_norm keyword
2000 Aug 04 - J. Newmark - NOPROMPT keyword for batch processing
large data sets, cleaned up memory managment
2001 Mar 15 - F. Auchere - added time varying dark current
Version 3.0
2001 May 18 - J. Newmark/F. Auchere - new response normalization
2001 Jun 21 - N. Rich - Add NRL check at (near) line 332
2001 Aug 15 - Zarro (EITI/GSFC) - add optional FITS out file name in
OUTFITS
Version 4.0
2001 Dec 10 - J. Newmark - changed default options,
new calibration, new default options lead
to default output type as float
2002 Jan 12 - F. Auchere - modified call to eit_dark to handle
different offset in binned images
2003 Apr 17 - D. Biesecker - modified call to is_fits to correctly
handle collapsed directory path names
2003 Jul 17 - C.A. Young - added 180 degree roll correction
2003 Sep 13 - C.A. Young - added fix for multifile list
2004 Feb 4 - D.M. Zarro - Update header for rotation correction
[Previous]
[Next]
Name: eit_proton_summary
Purpose: "standard" processing pipeline/summary plot for proton events
Calls: ***
ARR2STR [1], Arr2Str [2], BOX_MESSAGE, CONCAT_DIR [1], CONCAT_DIR [2]
CONCAT_DIR [3], CT2RGB, DELVARX [1], DELVARX [2], DELVARX [3], DELVARX [4], EIT_PREP
EVENT_MOVIE2, FILE2TIME, FILE_EXIST [2], LOADCT, PLOT_GOES, TIME2FILE, UTPLOT [1]
UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], WDEF [1], WDEF [2]
anytim [1], anytim [2], anytim [3], anytim [4], anytim [5], concat_dir [4], confac
cosmic_stat, cube_edit, data_chk [1], data_chk [2], delvarx [5], deriv_arr [1]
deriv_arr [2], eit_files, eit_fill_cube, evt_grid, file_append [1]
file_append [2], file_exist [1], file_exist [3], get_logenv [1], get_logenv [2]
grid_data [1], grid_data [2], gt_tagval [1], gt_tagval [2], html_doc
linecolors [1], linecolors [2], mkthumb, mreadfits_fixup, plot_goesp, rd_gxd [1]
rd_gxd [2], rd_gxd [3], read_eit, safe_log10, savegen [1], savegen [2]
sel_timrange [1], sel_timrange [2], special_movie, ssw_deltat, ssw_unspike_cube
str_replace [1], str_replace [2], strtab2html, thumbnail_table_html, time_window
zbuff2file [1], zbuff2file [2]
[Previous]
[Next]
Project : SOHO
Name : EIT_QUALITY
Purpose : flag problematic EIT file
Category : synoptic gbo
Syntax : IDL> flag=eit_quality(file)
Inputs : FILE = EIT FITS file
Outputs : 1 if good quality science data, 0 if problem with EIT files (camera err, engineering, etc)
CALLS: ***
ARR2STR [1], Arr2Str [2], IS_STRING, LOC_FILE [1], LOC_FILE [2], LOC_FILE [3]
read_eit
History : Written 20 Dec 2005, D. Zarro (EITI/GSFC)
Contact : dzarro@solar.stanford.edu
[Previous]
[Next]
Name:
eit_ratio
Purpose: create a ratio of two images and optionally display it
Input Parameters:
file1 - FITS file of image 1 = numerator
file2 - FITS file of image 2 = denominator
Output Parameters:
eit_ratio returns ratio array
Keyword Parameters:
show - set if want output displayed
xpixs, ypixs- set if want a subset pixel range of image,
format is xpixs=[xstart,xstop]
magnify - set to integer multiple for output REBIN size
image[1,2]_no - set for image # of 3D LZ file
temp_bar - set to display a temperature bar on bottom of image
Calling Sequence:
eit_ratio,'efr19960924.123456','efr19960924.654321',/show,$
magnify=3,xpixs=[100,300],ypixs=[100,300]
CALLS: ***
COLOR_BAR, DEVICELIB [1], DEVICELIB [2], EIT_FXPAR, EIT_PREP, FAKE_MISSING_BLOCKS
IMAGELIB, LOADCT, eit_temp
RESTRICTIONS:
the temperature calibration is preliminary and not very
accurate.
Method:
read FITS file, load proper color table, subtract background,
replace missing blocks, degrid, flat field, take ratio, subset if
desired, rebin if desired,optionally display
History:
24-Sep-1996 - J. Newmark - created
08-Dec-1996 - J. Newmark - took out all processing steps and
directly called routine EIT_PREP
1997/02/06 Added missing block replacement (unit ratio) and
MIN_RATIO, MAX_RATIO keywords for display. D.M. fecit
1997/02/10 J. Newmark - added "rough" temperature bar.
1997/09/19 D.M. fecit - added PIXMAP keyword
1997/10/01 J. Newmark - new call to eit_temp
[Previous]
[Next]
NAME:
eit_readcal
PURPOSE
reads a calibration lamp
INPUTS
file = file name
OPTIONAL INPUTS
hdr = a variable to recieve the FITS header
OUTPUTS
returns a calibration lamp
KEYWORDS
nonorm : if present, then does not normalize to exposure time.
CREATION
18-May-2001. F. Auchere
MODIFICATIONS
29-Jul-2002. F. Auchere Takes into account the calibration lamp exposure time correction
30-Jul-2002. F. Auchere Added the optional "hdr" parameter to return the header of the file
23-Aug-2002. F. Auchere Added /nonorm keyword
CALLS:
[Previous]
[Next]
NAME:
eit_ring_flux
PURPOSE:
Calculate the total flux (DN) per second with a given annulus in
an EIT full field image.
CALLING SEQUENCE:
flux = eit_ring_flux(fits_file, inner_rad, outer_rad)
INPUTS:
infile = FITS file of full field EIT image
inner_rad = inner edge for annulus in solar radii
outer_rad = outer edge for annulus in solar radii
OUTPUTS:
flux = total flux per second in given annulus
OPTIONAL INPUTS:
eastonly = set if want semi-annulus on east limb
westonly = set if want semi-annulus on west limb
EXAMPLES:
Return flux within annulus 1.13-1.17 Rsun:
flux = eit_ring_flux(infile,1.13,1.17)
Return all flux on disk:
flux = eit_ring_flux(infile,0,1.0)
Return all flux off disk:
flux = eit_ring_flux(infile,1.001,2.0)
MODIFICATION HISTORY
1997 Jan 08 - (JSN) - Written
1997 Feb 13 - J. Newmark - add east and west keywords
1997 Mar 17 - J. Newmark - account for change in exptime parameter
1997 Apr 10 - J. Newmark - use normalize exptime, normalize filter
in eit_prep
CALLS: ***
DIST_CIRCLE, EIT_FXPAR, EIT_PREP
COMMON BLOCKS: none
NOTES: this procedure will fill in missing blocks with the /surround
option, i.e. with average of surrounding blocks.
If image taken using AL+1 Filter will normalize to CLEAR.
[Previous]
[Next]
Name: eit_struct
Purpose: return eit_structure
Input Parameters:
number (optional) - number structures returned - default is one
Keyword Parameters:
catalog - if set, return catalog structure (subset)
ncomment - if set, number of COMMENT records to include
Calling Sequence:
str=eit_struct( [number , ncomment=xx ] )
CALLS: ***
REM_TAG [1], REM_TAG [2], REP_TAG_VALUE
CALLED BY:
read_eit
History:
27-apr-1996 S.L.Freeland (map FITs header)
24-may-1996 S.L.Freeland (add .VERSION, make .MJD long)
15-aug-1996 S.L.Freeland (Version 2, naxis3 - add /CATALOG switch)
27-oct-1996 S.L.Freeland (typo (BCALE->BSCALE))
11-dec-1996 S.L.Freeland (add /CATALOG keyword for testing)
14-mar-1997 J.S.Newmark (add history to obtain version number)
11-apr-1997 S.L.Freeland (add NCOMMENT keyword , a few new tags)
14-apr-1997 J.S.Newmark (add sci-obj field to catalog)
05-may-1997 J.S. Newmark (add s/c pointing info tags)
26-mar-1998 J.S. Newmark (add handling for corrected date_obs)
02-jun-1998 J.S. Newmark (add cmp_no keyword)
[Previous]
[Next]
NAME : EIT_SUB_BLKPIX
PURPOSE : Provide locs pixel (P1,P1/P2,P2) of current block, for port B
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_DRAW
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME:
EIT_SUB_CONVOFF
PURPOSE:
Convert between EIT grid subfield offset (0-1023) and offset
coordinate block (column,row). And vice verse.
Or when either is input along with the subfield pattern then
[xcen, ycen, xsize, ysize] is returned as arcsec.
CATEGORY:
EIT SUBFIELD
CALLING SEQUENCE:
result = EIT_SUB_CONVOFF(offset=offset,[sub_ids=sub_ids, /text])
result = EIT_SUB_CONVOFF(coor=[column,row],[sub_ids=sub_ids])
CALLED BY:
Any
CALLS TO:
none
INPUTS:
none
OPTIONAL INPUTS:
OFFSET When present this scalar is the individual EIT grid block
ranging from 0 to 1023, that is the offset.
When both OFFSET and COOR are present OFFSET take precedence.
COOR When present this 2-D vector is the column and the row of
the EIT offset.
When both OFFSET and COOR are present OFFSET take precedence.
SUB_IDS An array of size 1024 representing which of the EIT blocks
are selected. (1 represents select, 0 is not selected)
TEXT When present with OFFSET then the (column and row are
returned as a string in the form "column, row"
OUTPUTS:
When OFFSET is input, returns the COOR vector (intarr)
if TEXT is set then return string "col, row"
When COOR is input, returns the OFFSET (integer)
OPTIONAL OUTPUTS:
none
COMMON BLOCKS:
none
SIDE EFFECTS:
none
RESTRICTIONS:
none
PROCEDURE:
When OFFSET is input, modulo by the number of blocks in a column to
get the column and divide by the number of blocks in a row to get
the row.
When COOR is input, multiply row by number of blocks in a row and
add to the column value to get the OFFSET.
MODIFICATION HISTORY:
1995 February 17, Elaine Einfalt (HSTX)
[Previous]
[Next]
NAME : EIT_SUB_DISP
PURPOSE : Display Fresh image and draw block pattern.
CATEGORY : EIT Subfield
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUBFIELD, EIT_SUB_DRAW, EIT_SUB_MDI, EIT_SUB_PRINT, EIT_SUB_RESIZE
EIT_SUB_UTIL_AREA, SUBFIELD_EV
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_TOOL
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_DRAW
PURPOSE :
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_DRAW_EV
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], ARCMIN2HEL, CV_COORD, EIT_SUB_BLKPIX, EIT_SUB_DISP
EIT_SUB_ROTATE, EIT_SUB_UTIL_AREA, EIT_SUB_UTIL_BOX
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
2003/07/21 _ Einfalt - Spacecraft antenna problem result
in spending part of the year rolled 180
Images are still displayed as taken,
but Arcsec, Helio, and R/theta outputs
need to be adjusted.
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_DRAW_EV
PURPOSE : Handle EIT subfield tool draw window events
CATEGORY : EIT Subfield
EXPLANATION : When the draw_base was define in EIT_SUBFIELD the keyword
EVENT_PRO was set to this routine name, so this will be the
event handle for all draw window activities.
Those activities are either movements within the window,
including mouse button down and up, or darw window resizing.
Resizing events were defined in the draw_base definition,
by setting the TLB_SIZE_EVENT keyword.
SYNTAX : Indirectly through the XMANAGER by the EVENT_PRO keyword.
EXAMPLES : none
CALLED BY : XMANAGER whenever there is mouse activity in/at the window
CALLS TO : none
ARGUMENTS :
EVENT The widget sturcture as passed by the XMANAGER
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_DRAW, EIT_SUB_RESIZE
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, 1995 Dec - Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_FIXOFF
PURPOSE :
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_MOVE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
RESPOND_WIDG, STR_SEP
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1999-march Einfalt - - passing group to more routines for V5 modal
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_MAN_INPUT
PURPOSE : User enters coordinates to place a mark.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUBFIELD, SUBFIELD_EV
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], EIT_SOHO_PLOT_COLORS, EIT_SUB_IO_TESTER
EIT_SUB_MAN_EV, EIT_SUB_UTIL_RO, HEL2ARCMIN [1], HEL2ARCMIN [2], STR_SEP, XMANAGER
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, 1997 Dec, Elaine Einfalt (HSTX)
2003/07/21 _ Einfalt - Spacecraft antenna problem result
in spending part of the year rolled 180
Images are still displayed as taken,
but Arcsec and Helio "X" need to be adjusted.
2004/06/14 - einfalt - changed a call from anytim to anytim2utc
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_MDI
PURPOSE : Toggle the displaying of the MDI FOV
CATEGORY : eit subfield
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUBFIELD, SUBFIELD_EV
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_DISP
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME:
EIT_SUB_MOVIE
PURPOSE:
Create a movie of a subfield image from a sequence of full
field iages, tracking the region of interest with solar
rotation.
CATEGORY:
FITS processing
CALLING SEQUENCE:
EIT_SUB_MOVIE outheader,outmovie,start time,stop time [,/arcsec]
[,ref_helio=ref_helio] [,ref_date=ref_date] [,wave=wave]
[,outsize=outsize] [,object=object] [,list=list]
INPUTS:
tstart I = start time
tstop I = stop time
OPTIONAL INPUT KEYWORD PARAMETERS:
ref_helio I = heliocentric Center of FOV (lon,lat)
arcsec I = set if ref_helio is in arcsec coords (lon-lat) instead of
helio lon-lat
ref_date I = reference time/date of center
outsize I = number of pixels for output image (x,y)
wave I = wavelength of images (default=195)
object I = object type (e.g. 'full fov') default = all images
list I = input list instead of time range
phot_norm I = set for using phot_norm keyword in eit_prep
OUTPUTS:
header O = output header (structure)
movie O = output processed subfield 3D data array
CALLS: ***
EIT_CATRD [1], EIT_CATRD [2], EIT_PREP, FAKE_MISSING_BLOCKS, eit_catrd [3]
eit_file2path, read_eit, ssw_track_fov
COMMON BLOCKS: none.
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
Read in catalog data. Read in selected images and headers.
Force solar pointing over short period of images to be constant,
this eliminates the "bad" limb fits. Extract selected FOV, tracking
solar rotation.
MODIFICATION HISTORY:
Written by: J. Newmark Date. September 1997
24-Nov-1997 J. Newmark - Account for CRPIX1,2 offset in FITS header
07-Apr-1998 J. Newmark - update
30-jul-1999 J. Newmark - add phot_norm keyword
[Previous]
[Next]
NAME : EIT_SUB_PRINT
PURPOSE : Print the current subfield display
CATEGORY : EIT subfield
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUBFIELD, SUBFIELD_EV
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : 0 = b/w print, 1=soho_color3 (cheap color), 2= $$$ color
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_DISP, PSPLOT [1], PSPLOT [2], RESPOND_WIDG
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, 1996/10/15 Elaine Einfalt (HSTX)
1999-march Einfalt - - passing group to more routines for V5 modal
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_RESIZE
PURPOSE : Resize window to user spec. and display image at that size.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_DRAW_EV
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_DISP
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_ROTATE
PURPOSE : Rotate subfield click and mark points.
CATEGORY :
EXPLANATION : Single middle click to draw window marks a point for rotation
(stored in common rotation_click). Double click middle button
unstores the click...(well, loads an out of range value to
be rejected by this routine).
SYNTAX :
CALLED BY:
EIT_SUB_DRAW
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], ARCMIN2HEL, EIT_SUB_UTIL_ROT, RESPOND_WIDG
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1999-march Einfalt - - passing group to more routines for V5 modal
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_SC_PIXEL
PURPOSE : Read EIT S/C pixel file return latest X and Y S/C pix, w/date
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_COOR
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : OLD_POINT if set, use pre recentering (sun more north) coords
OUTPUTS : none
CALLS: ***
STR_SEP, rd_tfile [1], rd_tfile [2]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1998 Jan - added old_point keyword for doing coords
on pre-center images.
2004-Jun-02 - einfalt - added getenv call to get stat_dir
modified to read last line smarter
and other mac related conversions
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_TOOL
PURPOSE : Select one or more EIT subfields.
CATEGORY : EIT Subfields
EXPLANATION :
SYNTAX : result = EIT_SUB_TOOL([file])
CALLED BY:
EIT_SUBFIELD, EIT_SUB_DISP, SUBFIELD_EV
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS :
FILE When present, the data file that will be displayed under the grid.
If not present, routine uses last receive full sun image.
OLD_POINT if set, use pre recentering (sun more north) coords.
OUTPUTS :
Returns the selected subfields, with the subfield offset substracted.
CALLS: ***
EIT_SUB_UTIL_BLK, EIT_SUB_UTIL_FILE, EIT_SUB_UTIL_LINE, EIT_SUB_UTIL_MOVE
EIT_SUB_UTIL_OFF, XREGISTERED
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS:
When un-highlighting a subfield block, then the part of the image
that was the little "nubby" is lost as well.
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1998 Jan - eitnfalt- added old_point keyword
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_AREA
PURPOSE :
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_DRAW
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_DISP, EIT_SUB_UTIL_BOX
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_BLK
PURPOSE : Process just these new blocks, adding/subtracting as needed.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_TOOL
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
BLOCKS The EIT subfield indices to be processed.
SUB_IDS BYRARR(1024) with indices of currently highlighted
blocks set to 1b, and all the rest set to 0b.
OFFSET The current offset value, may be shifted to fit.
OUTPUTS : none
CALLS: ***
EIT_SUB_UTIL_HIGH
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_BOX
PURPOSE :
CATEGORY :
EXPLANATION :
SYNTAX :
result = eit_sub_util_box(/init,size=size)
result = eit_sub_util_box(x,y)
result = eit_sub_util_box(/get_coord)
CALLED BY:
EIT_SUB_DRAW, EIT_SUB_UTIL_AREA
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS :
SIZE The size (squared) the area boxing will be incremented
OUTPUTS : none
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_COOR
PURPOSE : Set of the arcsec data coord. system for EIT images.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_FILE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : OLD_POINT if set, use pre recentering (sun more north) coords.
OUTPUTS : none
CALLS: ***
EIT_SUB_SC_PIXEL
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1998 Jan - eitnfalt- added old_point keyword
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_FILE
PURPOSE : Request to get a new EIT file and display it.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_TOOL
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : OLD_POINT if set, use pre recentering (sun more north) coords.
OUTPUTS : none
CALLS: ***
CONGRID [1], CONGRID [2], CONGRID [3], EIT_FXPAR, EIT_SUB_UTIL_COOR
EIT_SUB_UTIL_HELIO, EIT_SUB_UTIL_HIGH, EIT_SUB_UTIL_LINE, EIT_SUB_UTIL_MDI
EIT_SUB_UTIL_RO, FXPAR [1], FXPAR [2], LOADCT, READFITS [1], READFITS [2]
READFITS [3], REVERSE, linecolors [1], linecolors [2]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1996 Sept 25 Einfalt - added over_file ability (and c1)
1998 Jan - eitnfalt- added old_point keyword
2003-july-21 Einfalt - Added substruc.antenna_roll value
2004-may Einfalt - Added Mac/SSW compatablity
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_HELIO
PURPOSE : Draw solar heliographic coords. over top of EIT image.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_FILE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], EIT_SUB_SOLAR_GRID
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
2003/07/21 _ Einfalt - Spacecraft antenna problem result
in spending part of the year rolled 180
Images are still displayed as taken,
but Helio tricked into upside down b angle.
2004/06/14 - einfalt- change call from anytim to anytim2utc,
in order to minimize unnecessary code calls.
2004/06/15 - einfalt - changed keyword from FILE= to DATE=
[Previous]
[Next]
NAME : EIT_SUB_UTIL_HIGH
PURPOSE :
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_BLK, EIT_SUB_UTIL_FILE, EIT_SUB_UTIL_MOVE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS :
COLOR 0=black line
1=white line
BLOCKS The subfields to be processed (highlighted or unhighlighted)
SUB_IDS The array of all current field states, 0=un, 1=highlighted
OUTPUTS : none
CALLS: ***
EIT_SUB_UTIL_SIDE
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_LINE
PURPOSE : Assign vector values to be grid and subfield highlight lines.
CATEGORY : EIT Subfield
EXPLANATION :
Whenever the draw window size has changed (interactive user) then the
grid lines and block line sizes need to be redetermined.
SYNTAX :
CALLED BY:
EIT_SUB_TOOL, EIT_SUB_UTIL_FILE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS :
PIX_PER The number of device pixels per a single EIT block size.
GRID When present, draw the EIT grid over.
OUTPUTS : none
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_MDI
PURPOSE : DRAW THE MDI FOV
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_FILE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SOHO_PLOT_COLORS
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_MOVE
PURPOSE : Redraw EIT subfield pattern at the input offset location.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_TOOL
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
SUB_IDS BYRARR(1024) with indices of currently highlighted
blocks set to 1b, and all the rest set to 0b.
OFFSET A new offset value, may be shifted to fit.
OUTPUTS : none
CALLS: ***
EIT_SUB_FIXOFF, EIT_SUB_GETOFF, EIT_SUB_UTIL_HIGH, XREGISTERED
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
1999-march Einfalt - - passing group to more routines for V5 modal
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_NUB
PURPOSE :
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_OFF, EIT_SUB_UTIL_SIDE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
OUTPUTS : none
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_OFF
PURPOSE : Indicate the EIT subfield offset subfield with "nubbies".
CATEGORY : EIT Subfield
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_TOOL
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_GETOFF, EIT_SUB_UTIL_NUB
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_RO
PURPOSE : Plot camera readout area
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_MAN_INPUT, EIT_SUB_UTIL_FILE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS :
CAM_RO Intarr of camera readout area in blocks (0-31)
96/5/22 (p1row,p1col,p2row,p2col) a.k.a (y1,x1), (y2,x2)
(p1row,p1col,p2row,p2col) a.k.a (x1,y1), (x2,y2)
PIX_PER The number of device pixels per a single EIT block size.
OUTPUTS : none
CALLS: ***
EIT_SOHO_PLOT_COLORS
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_ROT
PURPOSE : Display rotation point of single middle click in subfield tool
CATEGORY :
EXPLANATION :
eit_sub_util_rot,ref_helio=[-20,15],rotat=24,date='10-oct-97',
arcsec=arcsec
print,arcsec = -416 428
SYNTAX :
CALLED BY:
EIT_SUB_ROTATE
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGUMENTS : none
KEYWORDS :
ref_helio - heliographic coords (latitude, longitude) -
rotat - number of hours to rotate
date - date of image
OUTPUTS :
arcsec - position of rotated point in arcsecs (latitude, longitude)
CALLS: ***
DIFF_ROT [1], DIFF_ROT [2], HEL2ARCMIN [1], HEL2ARCMIN [2]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Jeff Newmark
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME : EIT_SUB_UTIL_SIDE
PURPOSE : Draw 4 sides around highlighted EIT subfields.
CATEGORY :
EXPLANATION :
SYNTAX :
CALLED BY:
EIT_SUB_UTIL_HIGH
EXAMPLES : none
CALLED BY :
CALLS TO : none
ARGURMENTS : none
KEYWORDS : none
OUTPUTS : none
CALLS: ***
EIT_SUB_UTIL_NUB
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : V1, Elaine Einfalt (HSTX)
CONTACT : eit@xanado.nascom.nasa.gov
[Previous]
[Next]
NAME:
EIT_SUBFIELD
PURPOSE:
Create a widget for the selection of EIT program subfields
CATEGORY:
EIT Planning
CALLING SEQUENCE:
result = EIT_SUBFIELD
CALLED BY:
EIT_PROGRRAM_EV
CALLS TO:
SUBFIELD_EV (through the XMANAGER)
ARGUMENTS:
NCATDAYS If present, the number of catalog to search to
obtain full field images. (default = 5)
KEYWORDS:
SUBFIELDS The most recently accessed EIT subfield array
OVER_FILE String or string array (?) of a fits file name
with full path included.
CALLED BY:
SUBFIELD_EV
Note: if over_file is a LASCO image include
keyword LASCO
LASCO If over_file is set, then LASCO may be set to
1 for C1, 2 for C2, or 3 for C3 to provide flips
OPTIONAL INPUTS:
none
OUTPUTS:
Returns an INTARR of the subfield numbers selected by user
OPTIONAL OUTPUTS:
none
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3]
CW_PDMENU, DATATYPE [1], DATATYPE [2], DATATYPE [3], EIT_CATRD [1], EIT_CATRD [2]
EIT_SUB_DISP, EIT_SUB_MAN_INPUT, EIT_SUB_MDI, EIT_SUB_PRINT, EIT_SUB_TOOL, GET_UTC
RESPOND_WIDG, REVERSE, STR_SEP, SUBFIELD_EV, UTC2TAI, XLOADCT [1], XLOADCT [2]
XLOADCT [3], XMANAGER, break_file [4], eit_catrd [3], linecolors [1]
linecolors [2]
COMMON BLOCKS:
none
SIDE EFFECTS:
none
RESTRICTIONS:
none
PROCEDURE:
Create the EIT subfield widget where the draw window displays the
EIT subfield grid for user to select desired subfields.
MODIFICATION HISTORY:
1994 December 2, Elaine Einfalt (HSTX)
1996 July 3, Einfalt - search catalog for full field images only
1996 Sept 25 Einfalt - added over_file keyword (full path)
and LASCO keywords
1999-march Einfalt - - passing group to more routines for V5 modal
2003-july Einfalt - Added substruc.antenna_roll to structure
2004-may Einfalt - Added Mac/SSW compatablity
2004-june einfalt - removed call to UTIL_TRUNC_FILE
[Previous]
[Next]
NAME:
eit_temp
PURPOSE:
Specify 195/171 or 284/195 ratio or individual DN values ==> compute TEMP,
emission measure
CALLED BY:
eit_ratio
NOTE: THIS IS MAINLY A PLACEHOLDER, HOWEVER IT MAY BE USEFUL.
CALLING SEQUENCE:
temp=eit_temp(dn171=dn171,dn195=dn195,expotime=expotime,$
inratio=inratio,dn284=dn284]
INPUTS:
MUST SPECIFY EITHER
dn171,dn195,dn284 = observed flux in image, best if corrected by
eit_prep, must have background subtracted
OR
inratio = calculated 195/171 or 284/195 ratio, e.g. from eit_ratio
pre_calc = MUCH FASTER - use pre-calculated spectra if online.
show = set if want to see temperature dependence upon ratio
expotime = total exposure time (including shutter close time)
if specifying non-normalized DN values
r195_to_171 = (default) return temperature based upon 195/171 ratio
r284_to_195 = return temperature based upon 284/195 ratio
OUTPUTS:
This function returns temperature and emission measure (cm^-5).
IF a dn195 level is not input the output emission measure is
normalized to 1 DN a observed in 195.
temperature = computed temperature for given ratio
dem = column emission measure (cm^-5)
value returned can be compared to what one
would get from a true DEM calculation, i.e.
we return the Isothermal Emission measure *0.7/0.3
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], INTERPOL, concat_dir [4], eit_flux
interp2d [1], interp2d [2], interp2d [3], ssw_packages
RESTRICTIONS:
Default: Pressure for level populations = (Ne*T) 1e15 K cm^-3
(not very sensitive to this)
VERY Slow if calculating spectra from scratch.
Pre_calc: this uses a set of pre-calculated spectra for a
pressure of 1e15 cm^-3 K and spanning emperatures
from 30,000-9.0e6 K. These are then interpolated to
obtain the final answer.
MODIFICATION HISTORY:
1997 Feb 05, J.S. Newmark, Written.
1997 Sep 29, Version 2.0 Substantial re-write to work with CHIANTI
J. Newmark
1997 Oct 06, Version 2.1 J. Newmark, use fact that temperature is
slowly varying and dem is linearly proportional to
signal, use histograms - speed up for large inputs
1998 Feb 13, Add error checking, put floor on arrays of 1 DN, max
ratio = 5
1999 May 07, J.S. Newmark - new calibration
1999 Jun 08, Version 3.0, J. Newmark
BUG FIX for Emission measure calculation
1999 Nov 17, J.S. Newmark - new calibration
[Previous]
[Next]
Name:
EIT_TIME_CORRECTION
Purpose:
To interpolate or extrapolate the most recent OBE - LOBT time difference
from the values in TIME_DIFFERENCE_DB.
Input Parameters:
None
Output:
None
Keywords:
CORRECTION_STRING - Used to return an ASCII string with the
time difference.
INTERP - Turn on interpolation
Calling Sequence:
dt = EIT_TIME_CORRECTION(obe_time, CORRECTION_STRING = CORRECTION_STRING)
CALLS: ***
ANYTIM2UTC [1], ANYTIM2UTC [2], CONVERT_TIME, INTERP_TIME_CORR, WORD_SWAP
Restrictions:
If large jumps in the difference occur between realtime contacts, this
routine could return inaccurate values.
History:
1997 April 17 - D.M. fecit.
1997 August 27 - Added CORRECTION_STRING keyword D.M. fecit.
1999 Feb 7 - Added binary search and indexing to speed finding
the right record, created unix version for Solaris - DW
1999 Feb 8 - Added INTERP keyword - IMPORTANT due to time offset jumps
during out of contact periods, interpolation can produce
totally bogus answers - DW
@(#)eit_time_correction.pro 1.2 02/09/99 NRL IDL LIBRARY
[Previous]
[Next]
NAME:
eit_vis2euv
PURPOSE
applies the visible to euv relationship to a vector or array
INPUTS
img = vector or array to convert
OUTPUTS
the converted vector or array
KEYWORDS
degrad : degradation index defined as the average of a normalized calibration lamp
reverse : if present, then conputes the euv to visible conversion
CREATION
18-May-2001. F. Auchere
MODIFICATIONS
04-Sep-2002. F. Auchere Added the reverse keyword to make the euv to visible conversion
CALLED BY
eit_norm_response [1]
[Previous]
[Next]
NAME:
eit_witness
PURPOSE:
Return the product of the reflectivities of the
primary & secondary witness mirrors
CALLING SEQUENCE:
eit_witness, 171, wave, reflect
eit_witness, 195, wave, reflect
eit_witness, 284, wave, reflect
eit_witness, 304, wave, reflect
INPUTS:
Sector = 171, 195, 284 or 304
OUTPUTS:
wave = Vector of wavelengths
reflect = Vector of reflectivities
CALLED BY:
plot_eit_area, plot_eit_mirror
HISTORY:
5-Jun-95, J. R. Lemen and R. C. Catura, Written
[Previous]
[Next]
NAME :
EIT_XCAT
PURPOSE : Interactive Widget interface to LZ catalog
CATEGORY : Catalogs
EXPLANATION : This function calls EIT_CATRD with interactive + LZ switch
SYNTAX : list = eit_xcat(filename=filename,tai=tai)
EXAMPLES :
IDL> file=''
IDL> list=eit_xcat(file=file)
CALLED BY : top level
CALLS TO : eit_catrd
ARGUMENTS : none
KEYWORDS :
FILENAME IF present, then return the fits file names matching
the entries in the returned variable.
TAI IF present, then return the TAI time matching
the entries in the returned variable.
CAT_STRUCT -Alternate catalog input, raw catalog previously
created, as by read_eit
OUTPUTS :
Return a string array of catalog entries,or a raw IDL structure
of entries. Optional outputs include list of filenames or tai times
CALLS: ***
EIT_CATRD [1], EIT_CATRD [2], eit_catrd [3]
COMMON : none
RESTRICTIONS: none
SIDE EFFECTS: none
PROJECT : SOHO - EIT
HISTORY : J. Newmark 22-Aug-96
CONTACT : newmark@eitv2.nascom.nasa.gov
[Previous]
[Next]
Name: eitoversxt
Purpose: draw EIT contours on SXT image; optionally overlay EIT 32x32 grid
Calling Sequence:
History: proto-version, 23-jan-1996
CALLS: ***
CONGRID [1], CONGRID [2], CONGRID [3], FILE_EXIST [2], FIND_LIMB, FXPAR [1]
FXPAR [2], PB0R, READFITS [1], READFITS [2], READFITS [3], WDEF [1], WDEF [2]
coal_image, data_chk [1], data_chk [2], file_exist [1], file_exist [3]
fmt_tim [1], fmt_tim [2], gt_day [1], gt_day [2], ocontour [1], ocontour [2]
ocontour [3], ocontour [4], ocontour [5]
Restrictions:
only full sun EIT and SXT , simple FITS for now
TODO - rotate SXT image to input time
[Previous]
[Next]
NAME:
ELTHEORY
PURPOSE:
This procedure computes various properties of the Thomson scattered
light by an electron in the solar corona.
CATEGORY:
CME
CALLING SEQUENCE:
ELTHEORY, Rin, T, R, B, Bt, Br, Pol
INPUTS:
Rin: Impact Distance (in solar radii)
T: Angle from plane of sky (in degrees)
OUTPUTS:
R: Distance of electron from sun center (in solar radii)
B: Total Brightness for one electron at R,T
Bt: Tangential Brightness for one electron at R,T
Br: Radial Brightness for one electron at R,T
Pol: Polarization for one electron at R,T
KEYWORDS:
LIMB: The limb darkening coefficient. If not set then use
0.63 (550 nm).
CENTER: If set, then the output is returned as the central solar
disk brightness. THe default is mean solar brightness.
CALLED BY:
CALC_CME_MASS, VDHCOEFF
COMMON BLOCKS:
None
SIDE EFFECTS:
None
RESTRICTIONS:
None
PROCEDURE:
This procedure uses the Thomson electron scattering equations
to compute the total brightness, the radial and tangential
components and the polarization associated with one electron.
Limb darkening of 0.63 is included.
The results are returned in units of mean solar brightness.
MODIFICATION HISTORY:
Created IDL version 16 May 1996, RA Howard, NRL
RAH 11/25/98 Added common block for VDH coefficients
RAH 09/01/99 Added keyword for limb darkening
RAH 10/02/99 Changed output to be mean solar brightness
RAH 10/27/99 Added keyword for central/mean solar brightness
@(#)eltheory.pro 1.5 11/29/05 LASCO IDL LIBRARY
[Previous]
[Next]
Name: embed
Purpose: build ORs from multiple pfis - crude version
Input Parameters:
data = image array from sda file
rmap = structure (roadmap or sxt index record corresponding to
the image array
Output Parameters
function returns image
CALLS:
[Previous]
[Next]
@(#)emieit.pro 1.1 09/22/96 LASCO IDL LIBRARY
CALLS:
[Previous]
[Next]
@(#)emilook.pro 1.1 09/22/96 LASCO IDL LIBRARY
CALLS:
[Previous]
[Next]
@(#)emiview.pro 1.1 09/22/96 LASCO IDL LIBRARY
CALLS:
[Previous]
[Next]
PROJECT:
SOHO - SUMER
NAME:
ENCOFIT()
PURPOSE:
Determines the fitted encoder value for a given step
CATEGORY:
EXPLANATION:
SYNTAX:
Result = encofit()
CALLED BY:
CONVERS5
EXAMPLES:
INPUTS:
step - step
stepboarder - step border for polynoms switching
p1 - Coeff
p2 - Coeff
p3 - Coeff
p4 - Coeff
OPTIONAL INPUTS:
None.
OUTPUTS:
None.
OPTIONAL OUTPUTS:
None.
KEYWORDS:
None.
CALLS: ***
POLY
COMMON:
None.
RESTRICTIONS:
None.
SIDE EFFECTS:
None.
HISTORY:
Version 1, September 12, 1996, Dietmar Germerott, MAPE Lindau. Written
CONTACT:
Dietmar Germerott, MAPE Lindau (germerott@linax1.mpae.gwdg.de)
[Previous]
[Next]
NAME:
EPHEMERIS
PURPOSE:
To do basic calculation of solar ephemeris. The Date
for which the calculations should be performed can be
entered by an input variable or they will be prompted
by the program. The ephemeris are displayed on the screen
if none of the various output keywords is specified, else
the output of the procedure is given by keywords (see there).
CATEGORY:
PICO
CALLING SEQUENCE:
EPHEMERIS
INPUTS:
None
OPTIONAL INPUTS:
date: The date (in UT) for which the ephemeris are to
be calculated. date is a string matching to
the IDL time format; it has the same form
as the result of SYSTIME(). See there for
the correct format.
KEYWORD PARAMETERS:
There are quite a lot of Keywords to specify
with which you can retrieve the calculated
values. If no output keyword is specified,
the result of the calculations will be displayed
in the Log Window.
INPUT KEYWORDS:
===============
YEAR: gives the year (four digits requested, f.ex.
1994 instead of 94)
MONTH: gives the month (1 for JAN, 12 for DEC etc.)
DAY: gives the day.
HOUR: gives the hour in UT (from 0 to 23)
MINUTE: gives the respective minute
NOW: If specified, the calculations are done for
the systemtime given by SYSTIME().
HEADER: Date and time are taken out of an image header
MODIFY_HEADER: If set, no output is given on the screen
but the header on input is modified in a way
that all the calculated results are written
in. A header must be given.
OUTPUT KEYWORDS:
================
AZ_SUN: The Azimut of the sun, measured from the
south about the west. Azimut and height are
calculated for Pic Du Midi.
B0: The latitude of the solar disk center
DEL_SUN: The declination of the sun in degrees
HANGLE: The hourangle of the sun.
H_SUN: The height above the horizon. Height and
Azimut are calculated for Pic Du Midi.
INCLINATION: Returns the inclination of the solar
rotational axis.
JULDAT: Named variable which returns the Julian day
for that instant exact to 1/4 of a day. Up to
now it is not known how to perform a better
precision.
L0: The longitude of the solar disk center
LONGSUN: Returns the longitude of the sun in the
ecliptical plane.
NORTHANGLE: The angle between the solar rotational
axis and the direction to the zenith.
OBLIQUITY: Returns the obliquity of the ecliptic
(angle between ecliptical plane and equator)
RA_SUN: The right ascension of the sun (in degrees)
SIDERIAL_TIME: The siderial time for the date
SOLDIST: A named variable returning the distance of
the earth to the sun in Astronomical units.
OUTPUTS:
The output of the results is mainly performed by named
keywords. However, if no keywords are specified, the
output is sent to the Log window.
OPTIONAL OUTPUTS:
None
CALLED BY:
SUNGRID
EXAMPLE:
To get the basic ephemeris data about the sun for this
instant enter at the prompt
EPHEMERIS,SYSTIME()
however, if you need the B0 and L0 values for
any further treatment, call Ephemeris by
EPHEMERIS, date, B0=clatitude, L0=clongitude
CALLS: ***
SIGN
COMMON BLOCKS:
None
SIDE EFFECTS:
Unknown
RESTRICTIONS:
None
PROCEDURE:
This procedure has been transposed from a FORTRAN
routine written by Jean Michel Niot, spring 1994.
The numerical formulae are taken from Astronomical
Algorithms, Jean Meeus.
MODIFICATION HISTORY:
written by Jean Michel Niot, spring 1994
adapted to IDL by Alexander Epple, 26-OCT-1994, Pic Du Midi
[Previous]
[Next]
NAME:
EPHEMERIS
PURPOSE:
To do basic calculation of solar ephemeris. The Date
for which the calculations should be performed can be
entered by an input variable or they will be prompted
by the program. The ephemeris are displayed on the screen
if none of the various output keywords is specified, else
the output of the procedure is given by keywords (see there).
CATEGORY:
PICO
CALLING SEQUENCE:
EPHEMERIS
INPUTS:
None
OPTIONAL INPUTS:
date: The date (in UT) for which the ephemeris are to
be calculated. date is a string matching to
the IDL time format; it has the same form
as the result of SYSTIME(). See there for
the correct format.
KEYWORD PARAMETERS:
There are quite a lot of Keywords to specify
with which you can retrieve the calculated
values. If no output keyword is specified,
the result of the calculations will be displayed
in the Log Window.
INPUT KEYWORDS:
===============
YEAR: gives the year (four digits requested, f.ex.
1994 instead of 94)
MONTH: gives the month (1 for JAN, 12 for DEC etc.)
DAY: gives the day.
HOUR: gives the hour in UT (from 0 to 23)
MINUTE: gives the respective minute
NOW: If specified, the calculations are done for
the systemtime given by SYSTIME().
HEADER: Date and time are taken out of an image header
MODIFY_HEADER: If set, no output is given on the screen
but the header on input is modified in a way
that all the calculated results are written
in. A header must be given.
OUTPUT KEYWORDS:
================
AZ_SUN: The Azimut of the sun, measured from the
south about the west. Azimut and height are
calculated for Pic Du Midi.
B0: The latitude of the solar disk center
DEL_SUN: The declination of the sun in degrees
HANGLE: The hourangle of the sun.
H_SUN: The height above the horizon. Height and
Azimut are calculated for Pic Du Midi.
INCLINATION: Returns the inclination of the solar
rotational axis.
JULDAT: Named variable which returns the Julian day
for that instant exact to 1/4 of a day. Up to
now it is not known how to perform a better
precision.
L0: The longitude of the solar disk center
LONGSUN: Returns the longitude of the sun in the
ecliptical plane.
NORTHANGLE: The angle between the solar rotational
axis and the direction to the zenith.
OBLIQUITY: Returns the obliquity of the ecliptic
(angle between ecliptical plane and equator)
RA_SUN: The right ascension of the sun (in degrees)
SIDERIAL_TIME: The siderial time for the date
SOLDIST: A named variable returning the distance of
the earth to the sun in Astronomical units.
OUTPUTS:
The output of the results is mainly performed by named
keywords. However, if no keywords are specified, the
output is sent to the Log window.
OPTIONAL OUTPUTS:
None
CALLED BY:
SUNGRID
EXAMPLE:
To get the basic ephemeris data about the sun for this
instant enter at the prompt
EPHEMERIS,SYSTIME()
however, if you need the B0 and L0 values for
any further treatment, call Ephemeris by
EPHEMERIS, date, B0=clatitude, L0=clongitude
CALLS: ***
EPHEMERIX, SIGN
COMMON BLOCKS:
None
SIDE EFFECTS:
Unknown
RESTRICTIONS:
None
PROCEDURE:
This procedure has been transposed from a FORTRAN
routine written by Jean Michel Niot, spring 1994.
The numerical formulae are taken from Astronomical
Algorithms, Jean Meeus.
MODIFICATION HISTORY:
written by Jean Michel Niot, spring 1994
adapted to IDL by Alexander Epple, 26-OCT-1994, Pic Du Midi
[Previous]
[Next]
Project : SOHO - CDS
Name : ERASE_NIS_WIN
Purpose : Erase the NIS line-windows on the technical planning plot.
Explanation : Allows updating of the display in the technical planning
routine MK_RASTER to reflect the deletion of a line window
from the current list.
Use : IDL> erase_nis_win ,[num]
Inputs : None
Opt. Inputs : num - the particular line window to erase, otherwise all.
Outputs : None
Opt. Outputs: None
Keywords : None
Calls : ***
TP_BSCALE
CALLED BY:
MK_RASTER, TP_DISP_IEF
Restrictions: None
Side effects: None
Category : Technical planning
Prev. Hist. :
Written : C D Pike, RAL, 14-May-1993
Modified : Replace calls to NINT by ROUND. CDP, 17-Jun-95
Version : Version 2, 17-Jun-95
[Previous]
[Next]
NAME:
ERROR_MESSAGE
PURPOSE:
The purpose of this function is to have a device-independent
error messaging function. The error message is reported
to the user by using DIALOG_MESSAGE if widgets are
supported and MESSAGE otherwise.
AUTHOR:
FANNING SOFTWARE CONSULTING
David Fanning, Ph.D.
1645 Sheely Drive
Fort Collins, CO 80526 USA
Phone: 970-221-0438
E-mail: davidf@dfanning.com
Coyote's Guide to IDL Programming: http://www.dfanning.com/
CATEGORY:
Utility.
CALLING SEQUENCE:
ok = Error_Message(the_Error_Message)
INPUTS:
the_Error_Message: This is a string argument containing the error
message you want reported. If undefined, this variable is set
to the string in the !Error_State.Msg system variable.
KEYWORDS:
NONAME: If this keyword is set the name of the calling routine
is not printed along with the message.
TRACEBACK: Setting this keyword results in an error traceback
being printed to standard output with the PRINT command.
In addition, any keyword appropriate for the MESSAGE or DIALOG_MESSAGE
routines can also be used.
OUTPUTS:
Currently the only output from the function is the string "OK".
CALLS: ***
STR_SEP
CALLED BY:
FSC_DROPLIST, FSC_WINDOW, PICKCOLORNAME, TVREAD [3]
RESTRICTIONS:
The "Warning" Dialog_Message dialog is used by default. Use keywords
/ERROR or /INFORMATION to select other dialog behaviors.
EXAMPLE:
To handle an undefined variable error:
IF N_Elements(variable) EQ 0 THEN $
ok = Error_Message('Variable is undefined', /Traceback)
MODIFICATION HISTORY:
Written by: David Fanning, 27 April 1999.
Added the calling routine's name in the message and NoName keyword. 31 Jan 2000. DWF.
Added _Extra keyword. 10 February 2000. DWF.
Forgot to add _Extra everywhere. Fixed for MAIN errors. 8 AUG 2000. DWF.
Adding call routine's name to Traceback Report. 8 AUG 2000. DWF.
Switched default value for Dialog_Message to "Error" from "Warning". 7 OCT 2000. DWF.
[Previous]
[Next]
Name:
ex2dow
Purpose:
To convert from the standard 7 element time array
(HH,MM,SS,MSEC,DD,MM,YY) to a day of week number
for that date. (0=sunday, 1=monday,6=saturday).
Input:
timarr - Standard "ex" time array
(HH,MM,SS,MSEC,DD,MM,YY)
year is assumed to be of the form 91,
not 1991.
Output:
day of week
0 = sunday
1 = monday
str_dow = optional output, the day of week in
a string.
CALLS: ***
Ex2Int [1], Ex2Int [2]
CALLED BY:
anytim2weeks [1], anytim2weeks [2], do_reg_backup, ex2week [1], ex2week [2]
fits_disp_month [1], fits_disp_month [2], mk_fem [1], mk_fem [2]
mk_fem_file [1], mk_fem_file [2], week_loop [1], week_loop [2], week_loop [3]
History:
Written 12-Oct-91 by M.Morrison
[Previous]
[Next]
NAME:
ex2fid
PURPOSE:
Given a 7-element time array, generate the fileID (YYMMDD.HHMM)
for that time.
CALLING SEQUENCE:
fid = ex2fid(tarr)
fid = ex2fid(tarr, /sec)
INPUT:
tarr
OPTIONAL INPUT:
sec - If present, return (YYMMDD.HHMMSS)
CALLED BY:
BCS_BROWSE24, DISPLOI_MON5K, FIND_CAL [1], FIND_CAL [2], HXTPIXON, HXT_QLOOK
OBS_PLOT, bbso_name_fix, cfl_summary [1], cfl_summary [2], check_compile [1]
check_compile [2], check_dumps [1], check_process [1], check_process [2]
check_sci5k_trunc, daily_forecast [2], disp_sci5k, disp_therm_rs232
do_reg_backup, event_movie [1], event_movie [2], fits2time [1], fits2time [2]
fits2time [3], get1hk_info [1], get1hk_info [2], get_hk_info [1]
get_hk_info [2], get_mk3 [1], get_mk3 [2], get_sfm, go_batch [1], go_batch [2]
goes3sec_copy, gtab_file, hxtimg_accum, ip_que_dmpver, jitter_gif_xyimg
kp_name_fix, laststat [1], laststat [2], mk_gbl, mk_gev [1], mk_gev [2], mk_gsn_obs
mk_gsn_obs_s1, mk_gx, mk_hst_summary, mk_lasteit_movie, mk_mdi_fits, mk_mdi_iap
mk_pnt, mk_sfc [1], mk_sfc [2], mk_trace_i1, mk_week_file [1], mk_week_file [2]
mk_ydbtape [1], mk_ydbtape [2], new_disp_sci5k [1], new_disp_sci5k [2]
nob_img_copy, nob_name_fix, norik_name_fix, pref_super, rd_gbl_raw [1]
rd_gbl_raw [2], rd_hirsaiso, rd_sci5k, rd_so_at_ftr, rd_therm_rs232
rd_week_file [1], rd_week_file [2], rd_week_file [3], rd_week_file [4]
redo_disploi, redo_mon_sci5k, reslot, sda2fits, selsis_copy [1], selsis_copy [2]
sft2sfc, soon2fits, ssw_check_contrib, sxt2file, sxt2mpeg, timeline, topsdb [1]
topsdb [2], web_seq, weekid [2], wfile, xdisp_fits, xdisp_trace [1], xdisp_trace2
xdisp_trace3, xhkplot, xset_chain [1], xset_chain [2], xwrite_hist, ydb_exist [2]
HISTORY:
Written 11-Dec-91 by M.Morrison
14-Oct-92 (MDM) - Modified to accept an array of times
[Previous]
[Next]
Name:
Ex2Int
Purpose:
Convert conventional date and time into days since 1979
and milliseconds of day.
Calling Sequence:
Ex2Int, time, msod, ds79
Inputs:
time= 7 element integer array containing, in order,
hr min sec msec day mon yr
Output:
msod= 4-byte integer: milliseconds of the day
ds79= 2-byte integer: number of days since 1-Jan-1979
CALLS: ***
JDCNV [1], JDCNV [2], JDCNV [3]
CALLED BY:
ADDTIME [1], ADDTIME [2], ALIGN_AR, AdjUsrTim, BCS_24HR_PLOT [1]
BCS_24HR_PLOT [3], BCS_ADDGOES, BCS_BROWSE24, BCS_CREATE_CAT, BCS_LOAD_CAT
EUTPLOT, LIST_BDA, MANY_DAYS, OP_ADD_SPECIAL, OUTPLOT [1], OUTPLOT [2], OUTPLOT [3]
OUT_GRSPHL1, OUT_GRSPHL2, OUT_HXSPH, RD_AR, READ_SPFDGRS, READ_SPFDHXS, SEL_AR
UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], XMOVIE_SFM, anytim2dd79
anytim2ints [1], anytim2ints [2], bcs_bpcfs, contact_sum, dp_clock_ref, ex2dow [1]
ex2dow [2], ex2week [1], ex2week [2], get_dn_rate [1], get_dn_rate [2], gettime
go_get_sirius, hxi_interp, input_time [1], input_time [2], int2sec [1]
int2sec [2], interp_img, make_32 [1], make_32 [2], mk_fem [1], mk_fem [2]
mk_fem_aoslos [1], mk_fem_aoslos [2], mk_fem_file [1], mk_fem_file [2]
mk_orbit [1], mk_orbit [2], op_edit [1], op_edit [2], op_same_time [1]
op_same_time [2], op_terminator [1], op_terminator [2], op_time_sort, pr_evn [2]
pr_fem, pr_gbe, pr_gev, pr_nel, pr_trav_summ, rd_week_file [1], rd_week_file [2]
rd_week_file [3], rd_week_file [4], rdbda_smm, sxt_interp [1], sxt_interp [2]
tim2intstr [1], tim2intstr [2], valid_pass4 [1], valid_pass4 [2], weekid [2]
ydb_exist [2]
Side Effects:
Restrictions:
None
History:
version 1.0, was adopted from Ex2Int.FOR (SMM software),
written by GAL, 15-Feb-1991
8-oct-91, Updated, JRL: Make output vectors 1-d if the input
time is 2-d, or a scalar if the input
time is 1-d.
17-mar-92, Modified, JRL: Made for loop index long type.
Modified to use all vector operations, ras, 93/6/7
Modified to correctly deal with years GE 2000, jmm, 7/28/94
[Previous]
[Next]
Name:
ex2week
Purpose:
To convert from any standard times to a week number.
Method:
The definition of week starts on a Sunday, and the
first week of the year is = 1. For the first
week, there can be less than 7 days in the week
The number of weeks in a year is 53 because
of the way the first and last weeks are counted.
Input:
timarr - Standard "ex" time array
(HH,MM,SS,MSEC,DD,MM,YY)
year is assumed to be of the form 91,
not 1991.
CALLS: ***
Ex2Int [1], Ex2Int [2], Int2Ex [1], Int2Ex [2], UNIQ [1], UNIQ [2], UNIQ [3]
anytim2ints [1], anytim2ints [2], ex2dow [1], ex2dow [2]
CALLED BY:
GetObsID, anytim2weekid [1], anytim2weekid [2], anytim2weekinfo
anytim2weeks [1], anytim2weeks [2], fmt_rasm, goes3sec_copy, mk_fem [1]
mk_fem [2], mk_fem_file [1], mk_fem_file [2], mk_pnt, mk_sd2, mk_sdc [2], mk_sdc [3]
mk_sdc [4], mk_sdl, mk_sl [1], mk_sl [2], mk_sot, mk_ssl, mk_week_file [1]
mk_week_file [2], nts_copy [1], nts_copy [2], week2ex [1], week2ex [2]
week_loop [1], week_loop [2], week_loop [3], wfile
OPTIONAL KEYWORD INPUT:
wid - If set, then return the string variable for
the week ID for each input time.
get_uniq- If set, then return all of the uniq week
numbers.
Output:
returns week number
year - the year number
HISTORY:
Written Oct-91 by M.Morrison
17-May-93 (MDM) - Modified to allow any time formats in
- Added /WID and /GET_UNIQ option
- Corrected the handling of vector inputs
[Previous]
[Next]
NAME:
EXP_CORR
PURPOSE:
This procedure modifies the header to give the correct exposure time and
adds comments to give the old time and the correction factor.
CATEGORY:
LASCO EXPFAC
CALLING SEQUENCE:
EXP_CORR,Hdr
INPUTS:
Hdr: LASCO image header
OUTPUTS:
The image header is modified
CALLS: ***
GET_EXP_FACTOR [1], GET_EXP_FACTOR [2], OFFSET_BIAS
CALLED BY:
IMG2DNPERSEC
RESTRICTIONS:
The LASCO image structure is required.
PROCEDURE:
MODIFICATION HISTORY:
Written by: RA Howard, NRL, 2/6/98
@(#)exp_corr.pro 1.1 04/17/98 LASCO IDL LIBRARY
[Previous]
[Next]
@(#)exp_fit.pro 1.1 09/22/96 LASCO IDL LIBRARY
[Previous]
[Next]
@(#)exp_fit_prog.pro 1.1 09/22/96 LASCO IDL LIBRARY
[Previous]
[Next]
Name: exp_num
Purpose: return sxt exposure number from roadmap or sxtindex
structure
CALLED BY
embed
[Previous]
[Next]
FUNCTION: exp_scale, min_exp, max_exp
Calculates the exposures for the light transfer data set from minimum and
maximum exposure times. The exposures are closer together at the extreams.
CALLS:
[Previous]
[Next]
FUNCTION: exp_scale, min_exp, max_exp
Calculates the exposures for the light transfer data set from minimum and
maximum exposure times. The exposures are closer together at the extreams.
CALLS:
[Previous]
[Next]
Project : SOHO - CDS
Name : EXPINT_DEMO
Purpose : Demonstrate effect of packetization module patch
Explanation :
Use : IDL> expint_demo
Inputs : None
Opt. Inputs : None
Outputs : None
Opt. Outputs: None
Keywords : None
Calls : ***
CIRCLE_SYM, CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], concat_dir [4]
Common : None
Restrictions: None
Side effects: None
Category :
Prev. Hist. : None
Written : C D Pike, RAL, 4-Jun-97
Modified :
Version : Version 1, 4-Jun-97
[Previous]
[Next]
Project : SOHO - CDS
Name : EXPORT_PLAN
Purpose : Exports a CDS science plan from the database
Explanation : This routine exports a CDS science plan covering a given period
of time from the database to a text file in a specific export
format. This text file can then be used by the routine
IMPORT_PLAN to import the plan into another database. This
allows the plan to be moved between different sites, or from a
user's personal database to the official CDS database.
Use : EXPORT_PLAN, START_TIME, END_TIME
Inputs : START_TIME = The start time of the period, in any CDS time
format.
END_TIME = The end time of the period.
Opt. Inputs : None.
Outputs : The complete plan definition, including any study definitions,
is written to a file. The filename is derived from the start
date, includes a version number, and has the extension ".plan";
for example, "19970410_123456_001.plan".
Opt. Outputs: None.
Keywords : PATH = The directory to write the file to. If not
defined, then the file will be written to the
current directory.
FILENAME = Returns the name of the output file.
ERRMSG = If defined and passed, then any error messages will
be returned to the user in this parameter rather
than depending on the MESSAGE routine in IDL. If
no errors are encountered, then a null string is
returned. In order to use this feature, ERRMSG
must be defined first, e.g.
ERRMSG = ''
EXPORT_PLAN, ERRMSG=ERRMSG, ...
IF ERRMSG NE '' THEN ...
Calls : ***
ANYTIM2CAL, ANYTIM2TAI, ANYTIM2UTC [1], ANYTIM2UTC [2], BREAK_FILE [1]
BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
DATATYPE [1], DATATYPE [2], DATATYPE [3], EXPORT_STUDY, GET_DETAIL, GET_UTC
LIST_DETAIL, LIST_PLAN, TRIM, break_file [4], concat_dir [4]
CALLED BY:
XPORT
Common : None.
Restrictions: The user must have write privilege in the current default
directory.
Side effects: There are certain assumptions which are made in order for this
software to be able to work correctly:
1. It is assumed that the GSET database is exactly the same on
both the export and import sides. There's no checking of
the validity of the GSET_ID number.
2. There's no checking of the durations of the various
studies. In particular, any implied assumptions about
whether or not GSETs need to be loaded are carried over
without checking.
3. It is assumed that any engineering studies (STUDY_ID=0) are
defined identically on both the export and import sides.
4. Campaign numbers are assumed to have exactly the same
meaning on both the export and import sides.
Category : Planning, Databases
Prev. Hist. : None.
Written : William Thompson, GSFC, 10-Apr-1997
Modified : Version 1, William Thompson, GSFC, 10-Apr-1997
: Version 2, Zarro, GSFC, 22-Apr-1997 -- fixed ERRMSG
Version 3, William Thompson, GSFC, 5 May 1997
Changed filename convention to include a version
number, and to put an "_" character between the date
and the time.
Version 4, William Thompson, GSFC, 17 May 1999
Fixed bug involving TRACKING and GET_GSET
Version : Version 3, 5 May 1997
[Previous]
[Next]
Project : SOHO - CDS
Name : EXPORT_STUDY
Purpose : Exports a study definition from the database
Explanation : This routine exports a study definition from the database to a
text file in a specific export format. This text file can then
be used by the routine IMPORT_STUDY to import the software into
another database. This allows study definitions to be moved
between different sites, or from a user's personal database to
the official CDS database.
Use : EXPORT_STUDY, STUDY_ID, STUDYVAR
EXPORT_STUDY, TITLE_ID, STUDYVAR, /USE_TITLE
Inputs : STUDY_ID
Opt. Inputs : STUDY_ID = Study ID number, or
TITLE_ID = Title ID number
STUDYVAR = Study variation index
Outputs : The complete study definition, including raster description,
linelists, and data window lists, is written to a file. The
filename is derived from the study name (OBS_PROG) and the
study variation number, and has the extension ".def"; for
example, "synop_1.def".
Opt. Outputs: None.
Keywords : USE_TITLE = If set, then the input parameter is the title ID
rather than the study ID. Not required if OBS_PROG
is passed.
PATH = The directory to write the file to. If not
defined, then the file will be written to the
current directory.
FILENAME = Returns the name of the output file.
UNIT = If passed, then contains the unit number of an
already opened file. This option is used by the
routine export_plan. When used, the header is not
written to the file, and the file is left open when
the routine is done.
ERRMSG = If defined and passed, then any error messages will
be returned to the user in this parameter rather
than depending on the MESSAGE routine in IDL. If
no errors are encountered, then a null string is
returned. In order to use this feature, ERRMSG
must be defined first, e.g.
ERRMSG = ''
EXPORT_STUDY, ERRMSG=ERRMSG, ...
IF ERRMSG NE '' THEN ...
Calls : ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], GET_DATAWIN, GET_LINELIST
GET_RASTER, GET_STUDY, GET_UTC, TRIM, concat_dir [4]
CALLED BY:
EXPORT_PLAN, MK_STUDY, XSTUDY
Common : None.
Restrictions: The user must have write privilege in the current default
directory.
Side effects: None.
Category : Planning, Databases
Prev. Hist. : None.
Written : William Thompson, GSFC, 16 May 1995
Modified : Version 1, William Thompson, GSFC, 23 May 1995
Version 2, William Thompson, GSFC, 20 June 1995
Added item WIN_FLAG.
Version 3, William Thompson, GSFC, 27 October 1995
Added item IEF_ID.
Version 4, William Thompson, GSFC, 2 October 1996
Added keyword PATH
Version 5, William Thompson, GSFC, 10 April 1997
Added keyword UNIT
Version 6, William Thompson, GSFC, 18 April 1997
Fixed bug introduced in version 5
Version : Version 6, 18 April 1997
[Previous]
[Next]
Project : SOHO - CDS
Name : EXPSLICE()
Purpose : Returns slice of a Spectral Window at specific exposure
Explanation :
Use : exposure = EXPSLICE(QLDS,window,exposure_no)
Inputs : QLDS: A CDS QL data structure.
window: The number or name of a detector window.
exposure_no: The number of the desired exposure.
Opt. Inputs : VALID: Supply to return a value of 0 instead of
generating an error message.
Outputs : Returns a 1- (GIS) or 2- (NIS) dimensional array
corresponding to the data from a specific exposure
in a specified detector data window for a specific
exposure. (A slice in the SOLAR_Y - WAVELNTH plane).
Opt. Outputs: VALID: Returns TRUE (NOT 0) if the slice went ok
and 0 if some error was detected.
Keywords : NOCHECK :
Set to skip QLMGR parameter check.
Calls : ***
DATATYPE [1], DATATYPE [2], DATATYPE [3], QLMGR, QL_DIMENSIONS, TRIM, WINDOWNO
Common : None.
Restrictions: QLDS must be a valid CDS QL data structure
WINDOW must be a scalar integer or string,
with the number or name of an existing detector
window.
EXPOSURE_NO must be <1, ... , nexposures(QLDS)-1 >
Side effects:
Category : Data_Handling, QuickLook, Interface.
Prev. Hist. : Requested by Paal Brekke 19 October 1993
Made by modifying WAVESLICE() - bugs here
should be checked out in WAVESLICE(), and
vice versa.
Written : Stein Vidar Hagfors Haugan, 19 October 1993
Modified : SVHH, 21 November 1993 - Changed name:
rstslice -> expslice
SVHH, 1-June-1994, Fixed minor bug in calling windowno.
SVHH, 20-Feb-1995, New data storage formats causes
some changes. Needs to be looked at. What
does EXPSLICE mean with data of dimensions
(SOLAR_X,SOLAR_Y,DEL_TIME) ?
For now I'm assuming (WAVELNTH,SOLAR_X,SOLAR_Y).
SVHH, 21-Aug-1995, Version 1.2
Fixed a mis-interpretation when both solar_x
AND dimensions.del_time are singular.
Version : 1.2, 21-Aug-1995
[Previous]
[Next]
PROJECT:
SOHO - SUMER
NAME:
EXPSTA()
PURPOSE:
returns start of exposure in CDS_EXT_TIME Structure
CATEGORY:
EXPLANATION:
SYNTAX:
Result = expsta()
EXAMPLES:
INPUTS:
header - BYTARR(92) Image Header
OPTIONAL INPUTS:
None.
OUTPUTS:
None.
OPTIONAL OUTPUTS:
None.
KEYWORDS:
None.
CALLS: ***
OBT2TAI, TAI2UTC
COMMON:
None.
RESTRICTIONS:
None.
SIDE EFFECTS:
None.
HISTORY:
Version 1, October 7, 1996, Dietmar Germerott, MAPE Lindau. Written
Version 2, April 15, 1998, GERMEROTT, MPAE LINDAU
modified to use the time functions obt2tai and tai2utc
CONTACT:
Dietmar Germerott, MAPE Lindau (germerott@linax1.mpae.gwdg.de)
[Previous]
[Next]
PROJECT:
SOHO - SUMER
NAME:
EXPTIM()
PURPOSE:
returns exposure time in seconds
CATEGORY:
EXPLANATION:
SYNTAX:
Result = exptim()
EXAMPLES:
INPUTS:
header - BYTARR(92) Image Header
OPTIONAL INPUTS:
None.
OUTPUTS:
None.
OPTIONAL OUTPUTS:
None.
KEYWORDS:
None.
COMMON:
None.
RESTRICTIONS:
None.
SIDE EFFECTS:
None.
HISTORY:
Version 1, October 7, 1996, Dietmar Germerott, MAPE Lindau. Written
CONTACT:
Dietmar Germerott, MAPE Lindau (germerott@linax1.mpae.gwdg.de)
[Previous]
[Next]
Project : SOHO - LASCO
Name : EXSUNC
Purpose : Extract sun's x,y and radius from image header
Category :
Explanation :
Syntax :
CALLED BY:
WR_BEAM, WR_BEAM1
Examples :
Inputs : None
Opt. Inputs : None
Outputs : None
Opt. Outputs: None
Keywords : None
CALLS: ***
chk_h4s
Common :
Restrictions:
Side effects: Not known
History : Version 1, 03-Oct-1995, B Podlipnik. Written
Contact : BP, borut@lasco1.mpae.gwdg.de
[Previous]
[Next]
NAME:
ext_subimg
PURPOSE:
Given an SXT full frame image, extract the portion of the PFI image
described in the second parameter.
CALLING SEQUENCE:
subimg = ext_subimg(index1, data1, index2)
INPUT:
index1 - The FFI index structure
data1 - The FFI data array (NX x NY)
index2 - The PFI index structure
OUTPUT:
returns the PFI portion of the FFI image
CALLS: ***
gt_pfi_ffi, gt_res
OPTIONAL KEYWORD INPUT:
use_avg - If set and the overlap of the PFI and the FFI is not
complete, then find the average of the signal level in the
edge portion of the FFI image, and set all non-overlapping
portions in the PFI image to that value. Used for dark
current subtraction method.
qdebug - If set, print some debugging messages
OPTIONAL KEYWORD OUTPUT:
line_map - A byte array of length equal to the number of lines in the
PFI. A value of 1 is set for every line that has overlap
with the FFI.
HISTORY:
Written 7-Apr-93 by M.Morrison (taking DARK_SUB work as the beginning)
[Previous]
[Next]
NAME: extract_ima2.pro
PURPOSE: extract a image from another
CATEGORY: General tools high level routine
CALLING SEQUENCE: extract_ima2, iin, hin, x1, x2, y1, y2, iout, hout
INPUTS: iin image array
hin header of image array
x1,x2 x interval of rectangle
y1,y2 y interval of rectangle
OPTIONAL INPUT PARAMETERS: None
KEYWORD PARAMETERS: None
OUTPUTS:iout extracted image array
hout header for extracted image
OPTIONAL OUTPUT PARAMETERS: None
COMMON BLOCKS: None
SIDE EFFECTS: opens a catalog of darks
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY: defined by ALL 6/24/93
corrected by M.B 11/03/93 for header's filling
SCCS variables for IDL use
@(#)extract_ima2.pro 1.0 25/6/93 :LAS
[Previous]
[Next]
Project : SOHO - CDS
Name : EZMGN
Purpose : To fit multiple Gaussians to NIS data and compare to CHIANTI
Explanation : A front end to D Zarro's multiple Gauss fitting routines
that also allows a visual comparison with a CHIANTI synthetic
spectrum. (ie Easy Multiple Gauss, NIS) cf EZFIT
A maximum of 20 profiles can be fitted simulstaneously.
Use : IDL> ezmg, w, sp
Inputs : w - wavelength array (as returned eg by pix2wave)
sp - the spectral data
Opt. Inputs : chfile - if supplied the CHIANTI spectrum is loaded from
there rather than the default.
Outputs : Results listed to screen and file "ezmg.result" if wanted.
Opt. Outputs: None
Keywords : WAVELENGTH - range of wavelengths to plot
AR - use CHIANTI active region data (default = QS)
WIDTH - fix width of line profiles to this value
BACKFIX - fix background to cursor value
HARDCOPY - do hardcopy of the screen plots
ILIMIT - set limit for intensity of identified lines
in CHIANTI spectrum (default = 1.0e10)
HELP - give instructions on use of cursor
NO_ID - do not put line IDs on hardcopy output.
Calls : FITSPEC
Common : None
Restrictions: The cursor interaction is querky as follows
(LEFT or MIDDLE buttons only)
1) when the spectrum is displayed you can either set a fixed
or first guess background level with the first click
(whether it is a fixed parameter or not is determined by
the /BACKFIX keyword) ..... or ....
2) if the first 2 clicks are below the X-axis they will be
taken to redefine the range of data to be fitted (if this
option is taken the NEXT click will determine the
background).
3) thereafter click at the peak (wavelength and intensity
important) to define each line you want fitted.
4) when done, click to the left of the Y axis to indicate
you are finished.
(RIGHT button only)
5) If pressed at any time within the DATA plot will give a
possible ID of the line.
6) if pressed to left of Y-axis will ABORT the procedure (ie
no fitting is attempted)
NB if you just want a hardcopy of the plots with no fitting
then enter ezmgn with the /HARD keyword and immediately click
the left button to the left of the Y-axis.
Side effects: None
Category : Spectral analysis
Prev. Hist. : None
Written : C D Pike, RAL, 15-Oct-96
Modified : Include errors in hardcopy. CDP, 4-Feb-97
Fix bug in total intensity calculation. CDP, 28-Feb-97
Add help and clean plot at end. CDP, 5-Mar-97
Allow data to be negative and make line identification
on hardcopy optional. CDP, 13-Jun-97
Version 6, 08-Jul-1997, William Thompson, GSFC
Changed format statements to allow larger dynamic range
Ignore sign of width
Explicitly trap > 20 profiles. CDP, 1-Sep-97
Version : Version 7, 01-Sep-1997