NAME:
	U_FEM_PLOTTER
 PURPOSE:
	overplot the Ulysses track on an SXT image
 CATEGORY:
	Yohkoh
 CALLING SEQUENCE:
 	.run u_fem_plotter
 INPUTS:
	u_data from ulysses_fem
       xx, yy, rr the limb fit for the image	
	
	NOTE: since this is a main program, you have to conjure up
	those inputs externally. I get xx, yy, rr by trial and error
	from the image and trace_limb, but one could also use any
	other limb-fitting procedure

 OPTIONAL (KEYWORD) INPUT PARAMETERS:
 OUTPUTS:
 COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
 MODIFICATION HISTORY:
	HSH, written 4-Jan-94

   Name: ucon_check

   Purpose: check for ucon routines which are duplicated in other places

 CALLS: ***
	CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], SPLIT_FILES, concat_dir [4]
	dir_list [1], dir_list [2], file_append [1], file_append [2], mk_mapfile [1]
	mk_mapfile [2], path_sw [1], path_sw [2]
   History: (SLF) 18-May-93

 NAME:
	UDN
 PURPOSE:
	Calculate the uncertainties of the SXT data.
 INPUTS:
	index - SXT reformatted index
	data - SXT data
 OPTIONAL INPUTS:
	exp_norm - data is already normalized by exposure time.
	        If supplied, data is not exposure normalized.
	wave - Photon wavelength (Ang).  If not supplied, defaults to 10 aug.
 OUTPUTS:
	uncertain - Uncertainties of the data
 CALLS: ***
	dn2ph, gt_expdur [1], gt_expdur [2], ph2dn
 PROCEDURE:
	Use "dn2ph" & "ph2dn" to convert SXT CCD output data (DN) to/from
       photon counts.
 HISTORY:
	11 Apr, 1994  written by T.Shimizu (Univ. of Tokyo)

 NAME:
	ULYSSES_FEM
 PURPOSE:
	reads in Ulysses ephemeris data (radius, heliolatitude and
	  longitude on 1st day of each month, 1993-1997
	the data were gotten from a public-domain service called
	  in Decnet by RUN NSSDCA::ANON_DIR:[ACTIVE.HELIO]HELICOOR,
 	  which also has ephemeris data for other interplanetary
	  spacecraft if anybody is interested
 CATEGORY:
	Yohkoh 
 CALLING SEQUENCE:
	ulysses_fem, u_data, /plot
 INPUTS:
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
	plot makes a plot of projected central distance 
       stop causes a stop at the end, so that you can see intermediate
         variables
 OUTPUTS:
	structure containing fields as labeled. I'll admit it, I am
         baffled by the coordinates, but I think .synlong has in it
         the heliographic longitude with W limb + 90, E limb 270. 
 CALLS: ***
	RDFILE [1], RDFILE [2], SUN [1], SUN [2], anytim2ints [1], anytim2ints [2]
 CALLED BY:
	HELIO_FOOT, TIM2UPOS
 COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
 MODIFICATION HISTORY:
	HSH, written 4-Jan-94
	HSH, regularized time fields, 11-Feb-94
       HSH, /stop and .synlong, 23-Jul-95

 NAME:
      uncdn
 PURPOSE:
      Given temperature, emission measure, and DN for a set of
      images, get the uncertainty in the DN, by getting the
      number of photons.
 CALLING SEQUENCE:
      unc = uncdn(index,data,te,comp_unc=comp_unc)
 INPUT:
      index = index structure
      data = uncompressed DN, not normalized, a cube or a scalar, or
             an array of (n, n_elements(index))
      te = Log(T), may be the same size as data cube, or a scalar
           If te is not the same size as data, te(0) is used...
 OUTPUT:
      unc = the uncertainty in the DN
 KEYWORDS:
      comp_unc = compression uncertainty in data, same size as data,
                 or a scalar, if not set 0 is used.
 CALLS: ***
	gt_expdur [1], gt_expdur [2], gt_filtb, sxt_flux [1], sxt_flux [2], sxt_flux [3]
 HISTORY:
      25-aug-94, jmm
      27-mar-95, jmm, got rid of em, this is not needed

NAME:
	uniq_fid
PURPOSE:
	Return a list of the unique File IDs in a directory
INPUT:
	dir	- a list of directories to check
OUTPUT:
	Returns a string array of all uniq fileIDs
METHOD:
	Get a directory listing for each valid prefix and
	sort on the file IDs for all files
 CALLS: ***
	PREFIXCK, SPLIT_FILES, UNIQ [1], UNIQ [2], UNIQ [3], file_list [1], file_list [2]
 CALLED BY:
	fndrids, mk_fidlist, mo_fndirck, newfiles [1], newfiles [2], newfiles [3]
HISTORY:
	Written 18-Apr-92 by M.Morrison

 NAME: 
     UNIT_DEM
 PURPOSE:
     Calculate DEM distribution of simple magnetic loop cooling
       by conduction, with idealized radiative losses at loop top.
       This is based on a calculation by Gordon Emslie.
     Hoped to be useful for intercomparing thermal responses of 
       miscellaneous detectors.
 CALLING SEQENCE:
     dem = unit_dem(tmax, alpha = alpha, res = res)
     dem = unit_dem(41)
 OUTPUTS:
     differential emission measure in cm^-3, for each value of
       temperature 0 < T <= TMAX.
 INPUTS:
     tmax is the temperature of loop top in million K;
 KEYWORDS:
     alpha is the parametrized radiative loss function, 
       emissivity = n^2 * T^alpha ergs/cm^2 sec, with default
       value of alpha set at 0.0 to resemble bremsstrahlung 
       (appropriate for high-temperature loops).
     res = divisor of tmax for higher-resolution sampling (e.g.,
       res = 10 means 1/10 MK per step).
 HISTORY:
     HSH, written 5-dec-92
     HSH, set default alpha to 0.0 because of convergence worries

NAME:
 unpack_grause
PURPOSE:
 unpacks hxi_index.grause, the grid use array
CALLING SEQUENCE:
 iuse = unpack_grause(hxi_index)
INPUT:
 hxi_index = the index structure for an HXI image, or an 8-element
             byte array
OUTPUT:
 iuse = a 64 element intarr, set to zero if that mod. collimator
        was not used in the fit, set to 1 if it was.
 CALLS: ***
	bits [1], bits [2]
 CALLED BY:
	hxtpix, makimg0
HISTORY:
 26-nov-1996, jmm

	NAME:
		unsharp_mask
	PURPOSE:
		Improve the contrast in the detailed features of
		an input image via unsharp masking method.
	CALLING SEQUENCE:
       	new_img = unsharp_mask( index, data, [sfd=sfd,
				deltcoef=deltcoef,smooth=smooth, 
				lowdelt=lowdelt, updelt=updelt, 
				lowint=lowint, upint = upint, 
				scale=scale])
	INPUT:
		data	Input image data (can be compressed or decomp).
			If it is decompressed it can not be byte-type.
		index	index or roadmap of input image
			(general defaults on Optional input parameters
			are based on half res. input image when 
			NO input index parameter is present.)
	Optional INPUT:
		sfd	sfd data-type flag for sfd input data.
		deltcoef Multiplier for difference image (img-smooth)
			or masking factor (1, 2, 3, or 4)(default=3)
		smooth	Size "box-car" smoother (default=5)
       	lowdelt Lower delta cutoff, a percentage of minimum, 
			decimal 0.xx. (default is 0)
       	updelt  Upper delta cutoff, a percentage of maximum, 
			decimal 0.xx. (default is 0.5)
       	lowint  Saturation cutoff, don't add back when signal 
			is BELOW a given percentage of the saturation 
			value, decimal 0.xx.
       	upint   Saturation cutoff, don't add back when signal 
			is ABOVE a given percentage of the saturation 
			value, decimal 0.xx.
		scale 	scale flag to ret. img as fltarr with
			a scaling of: bytscl(alog10(img))
	Example Calls:
	 a) For standard "sda" data:
	    IDL> img = unsharp_mask(index(12),data(*,*,12),/scale)
	    IDL> tv, img
		where, data and index are standard parameters from "sda" 
		data files.  The returned parameter "img" is byte-type
		and is already scaled for display (log10).
	 b) For "sfd" data:
	    IDL> img = unsharp_mask(index(9),data(*,*,9),/sfd,/scale)
	    IDL> tv, img
		where, data is compressed data from an "sfd" file.
	 c) For half-res "sfd" data:
	    IDL> img = unsharp_mask(data(*,*,19),/sfd,/scale)
	    IDL> tv, img
		when index is NOT specified half-res data is assumed.
		

 CALLS: ***
	gt_res, sfd_decomp, sxt_decomp [1], sxt_decomp [2]
	HISTORY:
		This routine was written by KTS and GAL 13-Oct-92.
		28-Oct-92, Modified keyword names to be the same
		as the "Original-unsharp masking rountine" called
		"enhancer.pro" which was written by
		A. McAllister and K. Strong, Apr. 1992.

  Unwraps SXT low 8 bit visible light images or other
  low-bit wrapped images.

PRO unwrap, image, newimage, maxnum=maxnum, noshow=noshow


 INPUT PARAMETERS:
	image = low-8 byte data

 OUTPUT PARAMETERS:
	newimage = integer, unwrapped version

 OPTIONAL INPUT KEYWORDS:
	maxnum = data value at which data wraps.  Default = 256.
	noshow = if set, do not display intermediate images.
		Default = display.

 CALLS: ***
	SPLINE
 CALLED BY:
	DTMAP
 Restrictions:
	Has trouble with sunspots.  Output can be patched with
	PAINT_PIXELS to finish the job.  Once a single clean
	image is made, a data cube can be cleaned with
	CUBE_FILL.
	If the result looks poor, try rotating the input image.

 Written October 1992  Barry LaBonte
 Fixed null array bug  April 11, 1994  BJL
 Added /MAXNUM keyword April 12, 1994  BJL

  pro up16,data,x,y,noscale=noscale,interp= interp

	program to rebin (x16) and display easily
	lwa,  1/30/92
	Modified 5/10/92 to include interp keyword

  pro up2,data,x,y,noscale=noscale,interp= interp

	program to rebin (x2) and display easily
	lwa,  1/30/92
	Modified 5/10/92 to include interp keyword

  pro up4,data,x,y,noscale=noscale,interp= interp

	program to rebin (x4) and display easily
	lwa,  1/30/92
	Modified 5/10/92 to include interp keyword

  pro up8,data,x,y,noscale=noscale,interp= interp

	program to rebin (x8) and display easily
	lwa,  1/30/92
	Modified 5/10/92 to include interp keyword

NAME:
	update_dpe
PURPOSE:
	Update the DPE value and exposure duration in the index.
	It should really only be used with GET_DC_IMAGE with
	the /INTERPOLATE option.
SAMPLE CALLING SEQUENCE:
	update_dpe, index, new_dpe
INPUT/OUTPUT:
	index	- The SXT index structure. It can be an array.
INPUT:
	new_dpe	- The new DPE value.  It must be a scalar
METHOD:
	The filter A position is modified if necessary to use
	or not use the 10% mask, the mbe value is updated, and
	the measured exposure value is updated.
 CALLS: ***
	HIS_EXIST, gt_dpe, gt_filta
 CALLED BY:
	get_dc_image [2], get_dc_image [3]
HISTORY:
	Written 15-Dec-94 by M.Morrison
	 7-Feb-95 (MDM) - Corrected error in how DPEs which require 
			  the ND filter are updated.

 NAME:
	UPDATE_FROM_PATH
 PURPOSE:
	Updates the YS: directories on Yohkoh from ISASS0:: using
	!path as the pointer as to which directories are to be
	updated.
 CALLING SEQUENCE:  
	.run UPDATE_FROM_PATH
 INPUTS:
	Start time for the "update" in the standard DEC, VMS time format
	The log file may be inhibited is required.
	The execution of the copy must explicitly requested - by default 
	it is inhibited.
 OUTPUTS:

 RESTRICTIONS:
	Designed to work under VMS, updating from ISASS0
	You must have write access to all the directories on the path
	A new branch of the tree on ISASS0 will not be seen.
 PROCEDURE:
	Uses !path to define the directories to be updated. This is
	into a string array using the "," as seperators.
	The equivalent directory of each segment of the path is checked
	for any changes found since the entered "update" date.
	Any change are copied to the equivalent directory.
	By default, all changes are logged to the file UPDATE_LOG.LIS
	Currently updates from the GEN, SXT, HXT and WBS trees.
 MODIFICATION HISTORY:
	RDB	12-March-92
	RDB	19-March-91,	included HXT and WBS
	RDB	03-June-92,	check correct path by number of elements
	RDB	16-June-92,	log files copied to "control" file
				input date made more bullet proof

 NAME:
       UPDATE_GXTOBS
 PURPOSE:
       Uses FTP on isaxc to update ydb files from isass0
 CATEGORY:
 CALLING SEQUENCE:
       .run update_gxtobs
 INPUTS:
 SIDE EFFECTS:
 RESTRICTIONS:
       Works on VMS only; needs FTP
 PROCEDURE:
 MODIFICATION HISTORY:
       RDB  Nov-92     Written
       CDP  Jan 93     Change handling of 'update' variable
                       include ctrl/y handler and $exit in the output file 
       CDP  Feb 93     Trap no files found response from isass0 
                       directory listing of files
       RDB  21-May-93  Changed files copied
	RDB  11-Jun-93	Modified so can copy files stored on directories 
			that are different to file prefix (e.g. gbl, gbe)
			Added gbe and gbl.

 NAME:
       UPDATE_MSSL
 PURPOSE:
       Uses exchange on mssl? update ydb files from sxt1
 CATEGORY:
 CALLING SEQUENCE:
       .run update_mssl_ydb
 INPUTS:
 SIDE EFFECTS:
 RESTRICTIONS:
       Works on VMS only
 PROCEDURE:
 MODIFICATION HISTORY:
       RDB  15-Jul-93	Written
	RDB  19-Jul-93	Added repeated attempts to do directory listing
	RDB  17-Nov-93	Split types into common and less common...
			Purge of copied file types as you go 
	RDB  12-Feb-94	Get copy node from logical COPY_NODE

 Name: UPDATE_SXT_CO_INDEX
 PURPOSE: edit the index.html file to insert the title of the newest nugget.
 CALLING SEQUENCE:  update_sxt_co_index, outfil=outfil
 OPTIONAL INPUT KEYWORD: outfil- specify the output, without which
              index.html is modified.
 SIDE EFFECT:  unless nocopy is set, the newest nugget is copied
              to SXTweekly.html.
 CALLS: ***
	ADDTIME [1], ADDTIME [2], ARR2STR [1], Arr2Str [2], FIND_FILE, FIND_NEWEST_NUGGET
	PLAN2HTML, anytim2ex [1], anytim2ex [2], anytim2weekid [1], anytim2weekid [2]
	int2secarr [1], int2secarr [2], rd_tfile [1], rd_tfile [2]
 HISTORY:  written, 17-Mar-01, NVN
           debugged, 20-Apr-2001, NVN

 NAME:
      UPDATE_TERM
 PURPOSE:
     Update the terminator database structures for use in TERM_SCORE2
 CATEGORY:
      Yohkoh
 CALLING SEQUENCE:
      update_term[,st_date,en_date]
 INPUTS:
      Does not require any 
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
     st_date - starting date if different from date specified in current
               'upto' file names.
     en_date - end date.  If not specified uses current time (!stime).
 OUTPUTS:
 CALLS: ***
	GET_POINTS, anytim2ex [1], anytim2ex [2], restgen [1], restgen [2], savegen [1]
	savegen [2], str_concat [1], str_concat [2]
 COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
 MODIFICATION HISTORY:
      DA, written 26-Jun-1995

 NAME:
	UPDATE_YOHLIB
 PURPOSE:
	Updates the YLIB directories using control file created by 
	the UPDATE_FROM_PATH routine
 CALLING SEQUENCE:  
	.run UPDATE_YOHLIB
 INPUTS:
	file YOHLIB_UPDATE_COM.DAT
 OUTPUTS:

 OPTIONAL OUTPUTS:

 RESTRICTIONS:
	fiel must be on current directory
 PROCEDURE:

 MODIFICATION HISTORY:
	RDB	16-June-92

 NAME:
	ute_max
 PURPOSE:
	To delete pixels in the arrays
	te0,ute0,em0,uem0 which have uncertainties
	greater than the value utemx., this'll also
       delete pixels with te0 gt 2.00
 CALLING SEQUENCE:
	ute_max,te0,ute0,em0,uem0,utemx
 INPUT:
	te0= a fltarr of (ipix,jpix,nexp), Log(T)
	ute0= Uncertainty in te0
	em0= Log(E.M./10^47)
	uem0= Uncertainty in em0
	utemx= the maximum value of ute0 that we'll allow.
 OUTPUT:
	The same arrays, with the offending pixles, i.e., those
	with ute0>utemx set to 0.0
 CALLED BY:
	sxt_fl_analxpro
 HISTORY:
	Written May '92 by J McTiernan
       Added te0<2.0 requirement Feb '93
       Changed units from t6 to te 12-jul-94

pro uterrplot,time,plus,minus,skip=inskip,horizontal=horizontal

NAME:
     UTERRPLOT
PURPOSE:
     Plot Error bars on a UTPLOT
CATEGORY:
CALLING SEQUENCE:
     uterrplot,time,plus,minus
INPUTS:
   IF HORIZONTAL=0:
     time = time vector 
     plus = the data value at the top of the error range
     minus = the data value at the bottom of the error range
   IF HORIZONTAL=1:
     time = time at earliest point in error range 
     plus = time at latest point in error range
     minus = the data value 
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
     /horizontal = Plot a horizontal error bar rather than a vertical error
                   bar. 
     skip = integer specifying how frequently an error bar should be
            plotted. Default is 1 which means plot every error bar.  If skip
            is set to 2, then every other error bar si plotted, etc.
OUTPUTS:
 CALLS: ***
	OUTPLOT [1], OUTPLOT [2], OUTPLOT [3]
 CALLED BY:
	HXT_LCURVE
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
     You must call utplot first to plot the data points.  This routine only
     adds error bars to a data plot which is assumed to already be plotted.
PROCEDURE:
     First call utplot to plot your data, then call uterrplot to add error
     bars. 
MODIFICATION HISTORY:
     T. Metcalf 2000-April-05

   Name: uvf_interval

   Purpose: return UV Flood interval

   Calling Sequence:
      uvfi=uvf_interval(times)		; return interval in seconds
      uvfi=uvf_interval(times,/string) ; return as string from common table
      uvfi=uvf_interval(fid=fid)	; Specify file id's instead of time


   Input Parameters:
      times - scaler or vector of desired dates
   
   RETURNS:
      Output is UV flood interval in seconds (if /string is not specified)

   OPTIONAL INPUT KEYWORDS:
      fid	- Times may be specified by file id
      strings - if set, return as string (default is integer seconds)

   PROCEDURE:
      Calls gtab_comm to get the common table output and "greps" it.
 CALLS: ***
	gtab_comm, str2cols [1], str2cols [2], wc_where [1], wc_where [2]
   Restrictions: assumes fixed format for string (output from gtab_comm)

   History:
     12-Jul-1994 (SLF) Written.
     12-jul-1994 (JRL) Added the fid=fid keyword option