EIT Software Examples


 
a) Useful routines (Jeff's top 10):

EIT_CATRD (and EIT_XCAT = widget form) -   
  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)

EIT_DISPLAY -   
       This function will display EIT movies (IDL Save Sets) and
       EIT images (FITS files) using eit$\_$movie and eit$\_$image
       respectively.

EIT_DUMP - produce postscript or gif files
  
EIT_FILE2PATH - 
    translate eit file name to "standard path" on local system
      (vms<->unix, handle collapsed or yyyy/mm/dd trees)

EIT_FXPAR - Obtain the value of a parameter in a FITS header

EIT_PREP -
       A preliminary version of a processing routine for analyzing
       EIT images. Presently this version reads in the named FITS
       image and produces a background subtracted, degridded and
       flat-fielded output array or FITS file. In the future this
       routine will also vignette correct. 
 
MK_EIT_L1 -
       Create Level 1 FITS files via eit_prep. Demo program to show usage.
       Works on 2D/3D conbinations.

READ_EIT -  read eit image(s) into 2D or 3D array and 
    map header->structure  

MOVIE_MAKER - 
      This procedure will create EIT movies (IDL Save Sets).
      The user interactively selects a series of FITS images,
      scales them, and saves them.

EIT_FLUX and EIT_TEMP -
   Specify Emission ratio (e.g. 195/171) ==> compute Te
   Specify Te + filter + DEM ==> compute DN 
   Specify Te + filter + DN ==> compute DEM


b) Other "eit specific"

   read_eit.pro		; read EIT files into IDL structure array and data cube
			; (if no data, then only read headers for speed)

   eit_file2path.pro	; translate file names to local (online names)

   lasteit		; display most recent NN images for selected waves

   eit_colors.pro	; load selected color table (via wavelength)

   eit_catrds.pro	; span days for 'catrd_daily' and concatenate output
			; (include concatenation of /raw structures)

   misc...
   eit_fill_cube	; fill images in cube from nearest neighbors
			  (ex: try: lasteit,info,data,10,/w304,/fill)
   eit_data_scale	; "standard" 1 line scaling
   eit_qklcat		; ascii catalog of all online images 

c) "generic"            

   get_infox.pro 	; IDL structure->string array (1 line per structure)
			; permits user flexibility in output tags and order

   str2cols.pro		; string array -> table via fast algorithm
			; (example, catrd_daily string array -> columns)

   utplot,tim2dset,	; tons (literally) tools to act upon the IDL 
			; structures returned by read_eit 

   xstepper.pro		; data cube reviewer  (sxt,spartan,eit,lasco,...)
			; includes zoom, pixmap, etc, options.

------- demo "problem" - read online 171 images taken between 1-may and 5-may

IDL> dailys=eit_catrds('1-may-96','5-may-96',file=file)  ; read daily cats

IDL> fullnames=eit_file2path(file,exist,ecnt)      ; convert names->online 
IDL> found=where(exist)			           ; (only those found locally)
IDL> read_eit, fullnames(found), eit0		   ; read FITS headers
						   ; and map ->IDL struct

IDL> ss=where(eit0.wavelnth eq 171 and $	   ; apply some filter
   eit0.naxis1 eq 1024 and eit0.naxis2 eq 1024)    ; (future widget app.)

IDL> wanted=found(ss)				   ; filtered AND online
IDL> read_eit, fullnames(wanted), eit1, data	   ; read images-> cube

******** some "generic" routines shown acting on the EIT variables above ****

1) get_infox - provide user defined ascii summary of an IDL structure

                           |--user selected list---|

IDL> array=get_infox(eit1,'exptime,wavelnth,filter',/fmt_tim,/more)
     START TIME      EXPTIME  WAVELNTH  FILTER
 1-MAY-96  01:01:48     5.00     00171   Clear
 1-MAY-96  09:54:10     5.00     00171   Clear
 1-MAY-96  18:32:56     5.00     00171   Clear
 2-MAY-96  11:20:06     5.00     00171   Clear
 (--- etc ---)
 5-MAY-96  08:52:40     5.00     00171   Clear
 5-MAY-96  16:18:15     5.00     00171   Clear

2)  str2cols - break a string array into columns (eit catalog for example)

IDL> cattable=str2cols(dailys) & help, cattable
CATTABLE        STRING    = Array(8, 58)
IDL> more, cattable(6,wanted) + cattable(1,wanted) + cattable(0,wanted) 
 EFR19960501.010148  01:01:48 1996-May-01
 EFR19960501.095410  09:54:10 1996-May-01
 EFR19960501.183256  18:32:56 1996-May-01
 EFR19960502.112006  11:20:06 1996-May-02
 (--- etc ---)
 EFR19960505.085240  08:52:40 1996-May-05
 EFR19960505.161815  16:18:15 1996-May-05

3) utplot, outplot, evt_grid, plot_goes, tim2dset, sel_timrange, etc.
   correlate EIT instrument events or other SOHO/YOHKOH parameters
   select structures, plot structure tags and databases against UT


RE: SSW/EIT remote installations 
Cron job was installed on penumbra this week to generate one compressed 
tar file per SSW instrument or branch 
(see $SSW/offline/swminat/tar)

    http://www.space.lockheed.com/eit/ssw_install.html

FORM interface - user selects SSW instruments and a customized 
installation script is generated - user executes on local machine
to install SSW.


Return to the EIT User's Guide Index

Go back to Chapter 1 (EIT and Solar Software)

Return to the EIT Home Page

The SDAC home page

The SOHO home page