PROGRAM: NAME
	sxt_satscreen

 PURPOSE
	To identify images for which the only saturation is dark spike
	saturation, i.e., no solar x-ray saturation.  This allows these
	images to be processed into movie images without compositing.

 PROCEDURE
	The program takes advantage of the fact that dark spike saturation
	is only a single column wide in the x-direction.  Images with no
	x-ray saturation features above row 5 (to eliminate "bottom row"
	problems) are accepted as x-ray saturation free.  The rare cases
	of genuine single-column saturated solar features will not be caught
	by this code -- but this is not much of a problem as they
	will still be in the movie.
	*** TAKES ABOUT 20 SEC FOR (512,512,36) DATA CUBE ***
	*** WASTED A LOT OF TIME TRYING TO MAKE IT FASTER  :-(  ***

	Test is implemented on gt_res or keyword RES because saturation level
	for FR images is 235 DN.  Use either index structure or set
	/fr keyword flat.  Default is saturation=255 DN.

 CALLING SEQUENCE
	outlist = sxt_satscreen(index,data)
	outlist = sxt_satscreen(xx,data,/fr,/timer
	
 INPUT
	index = sxt index structure OR any non-structure input to ignore.
	****OK to use on mixed-resolution data cube if index is passed in.***
	data = raw, sxt_compressed image or data cube.

 OPTIONAL INPUT KEYWORD
	/fr specifies that all input data are full resolution
	/timer will print out 'SXT_SATSCREEN took:   ....." message.
 OUTPUT
	outlist = logic vector, unsaturated images are 1's.

 HISTORY
	16-Oct-96 [LWA], written.
 	19-Oct-96 [LWA], added timer keyword.
	23-Oct-96 [LWA], changed handle FR case, too.
 CALLS:

	PROGRAM: NAME
		DUSTOFF
	PURPOSE
		A program to interpolate out the 2 big dust specks
		in a half resolution FFI.
	CALLING SEQUENCE
		out = dustoff(img)
	HISTORY
		Written by LWA on 6-Aug-93
 CALLS:

	PROGRAM: NAME
		SXT_DIFFUSER2
		The difference from SXT_DIFFUSER is a variable smoothing
		width set by the value of iterate..
	PURPOSE
		To return a flat-fielded diffuser image for the purpose
		of monitoring x-ray damage to the CCD.
	CALLING SEQUENCE
		out=sxt_diffuser2(index,lo8_data,comp_data, $
			[iterate=iterate,nn=nn,spike=spike, $
			 dust=dust,ff_smooth=ff_smooth]
		out=sxt_diffuser2(index,lo8_data,comp_data)
		out=sxt_diffuser2(index,lo8_data,com_data,interate=6, $
                               nn=25,spike=300,/dust)
	INPUT
		index = index of either image, used for dark subtraction.
		lo8_data = diffuser image taken in lo8 mode.
		comp_data = compressed mode diffuser image.
	OPTIONAL KEYWORD INPUT
		interate = number of iterations in smoothing
			default = 6 iterations
		nn = initial smoothing window
			default = 17
		spike = amplitude of spike to be removed, default=75.
			set spike = 1000 to keep from removing spikes
		dust = set this to prevent removal of 2 big dust spots
		ff_smooth = do final smoothing with FFT.
	METHOD
		1.  Form reconstructed 12 bit image.
		2.  Subtract background, remove dust and fix edges.
		3.  Remove single pixel spikes.
		4.  Iteratively prepare a smoothed image.
		5.  Subtract smooth image from original image.
		6.  Further smoothing with Fourier Transform.
	HISTORY
		Written by LWA on 6-Aug-93
		Modified to use SPIKES and BIG_SMOOTH_DFSR.   4/29/94  LWA
 CALLS:

 NAME:
     NAR2COORD 
 PURPOSE:
     generate SXT pointing information from the NOAA AR number 
 CATEGORY:
 CALLING SEQUENCE:
 INPUTS:
     time, the string time to rotate coordinates to
     name, the string name of the AR sought (defaults to the first)
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
     angle allows input of angle coordinates 
     /gif writes a .gif file to /home/www/sxt_co/obs_coord.gif
       instead of to X
 OUTPUTS:
     prints to screen and makes an image. The three verticals on
       the image correspond to the NAR time, the LASTSFD time, and
       the prediction time
 CALLS: ***
	CONV_A2H [1], CONV_A2H [2], LOADCT, RD_NAR, WDEF [1], WDEF [2], WRITE_GIF, anytim [1]
	anytim [2], anytim [3], anytim [4], anytim [5], conv_h2p [1], conv_h2p [2]
	fmt_tim [1], fmt_tim [2], lastsfd [1], lastsfd [2]
 COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
 MODIFICATION HISTORY:
     4-Apr-00, written (HSH)

NAME       : ncongrid (pro)
 PURPOSE :
	shrink/expand multiple images in norikura format
 CATEGORY :
	idl/lib/nkr
 CALLING SEQUENCE :
	ncongrid,files,xfact=xfact,yfact=yfact,nx=nx,ny=ny,	$
		outfiles=outfiles,fact=fact
 INPUTS :
	files    : NKR type filenames
 OPTIONAL INPUT PARAMETERS : 
 KEYWORD PARAMETERS :
 OUTPUTS :
	none
 CALLS: ***
	CONGRID [1], CONGRID [2], CONGRID [3], NCONGRID, NEXTFILE, NKRGET, NKRSAVE
 COMMON BLOCKS : 	none
 SIDE EFFECTS : 	none
 RESTRICTIONS :	none
 PROCEDURE :
	If title keyword is set, new full window is always created.
 MODIFICATION HISTORY :
	'91/01/29   k.i.
	'92/11/12   k.i.
	'92/11/28   k.i.
	'93/02/04   k.i.
	'93/02/06   k.i.   'temp' --> outfiles(i)

 Copyright (c) 1993 P. R. Blanco, CASS, 0111, UCSD, La Jolla, CA92093-0111, USA
 and SERC, Polaris House, North Star Avenue, SN2 1ET, United Kingdom.

 Calculates an array, of dimensions specified by DIMS, of (distances)^2
 from the a central point specified by  coordinates Centers. If the center
 is not specified, it is assumed to be at the origin of the array.

 Returns:
 NDDIST2 (Float array(DIMS)) - array of distances^2.

 Input Parameters:
 DIMS (Integer , vector(NDIM)) - List of dimensions for the output array
 Centers (Float, vector(NDIM), default FLTARR(NDIM)) - Coordinates of the
 center from which the distances^2 are calculated.

 CALLED BY:
	HSI_PIXON_LOCAL_SMOOTH, LOCAL_SMOOTH
 History:
 22 Sep 93 - written by P. R. Blanco, CASS/UCSD to replace a routine of
 the same name.

 Copyright (c) 1993 P. R. Blanco                                  and
 SERC, Polaris House, North Star Ave. Swindon SN2 1ET, United Kingdom.

 Shifts an array in by Vector(1,2,3...) elements in N-dimensions
 Due to the inflexibility of IDL's SHIFT command each number of
 dimensions is handled separately. If an error occurs, the original 
 input array is returned.

 Input parameters:
 Data (array) - the array to be shifted
 Vector (1-D array, size = no. of dimensions of array Data) -  Size and 
 direction of shift. Vector may be a scalar if Data is a 1-D array.

 CALLED BY:
	wrap
 History:
 17 August 1993- new version by P. R. Blanco, CASS/UCSD.
 Handles 0 to 5 dimensions.

 NAME:
       NEAR_CONJ
 PURPOSE:
	enter Woo's ephemeris for the NEAR conjunction of 17 Feb 97
 CATEGORY:
	Yohkoh
 CALLING SEQUENCE:
       near_conj, near, /plot
 INPUTS:
       from data file (Woo's email)
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
       plot makes a plot of radial distance, with fem_grid
 OUTPUTS:
	structure with many things,
	among which near.xx is the WEST departure
                   near.yy is the NORTH departure
                   near.pp is the position angle (ecliptic coordinates):
                   "this is the angle in the plane of the sky measured with 
                   respect to ecliptic North increasing counterclockwise, 
                   so that the East limb is 90 degrees and the West limb
                   is 270 degrees" - Woo
                   
 WARNINGS
       your window if PLOT set 
 CALLS: ***
	ARROW, DAY_IN, FMT_TIMER [1], FMT_TIMER [2], anytim [1], anytim [2], anytim [3]
	anytim [4], anytim [5], fmt_tim [1], fmt_tim [2], get_rb0p [1], get_rb0p [2], rd_fem
	rd_tfile [1], rd_tfile [2]
 MODIFICATION HISTORY:
       HSH written 16-Feb-97, plated from SXT_ECLIPSE

 NAME:
     NEUPERT_PLOT
 PURPOSE:
     Make a standard plot of soft and hard X-ray light curves
 CATEGORY:
 CALLING SEQUENCE:
     neupert_plot, t1, t2
     neupert_plot,'25-aug 14:20','25-aug 15:00',/hxt,infil=infil,channel=0,/gif,back=.05
       (where infil = "hda000825.1344", made the plot in 000901.html)
 INPUTS:
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
     hxt forces HXT as the hard X-ray reference
       (default is the late lamented BATSE)
     infil=infil gets HXT counts (M1) from an HDA file
     logplot forces semilog plot
     bad, an array (e.g. 20+indgen(3)) that
       lists the HDA index entries not to be plotted
     channel, the HXT channel to be overplotted (M1 default) 
     back, an empirical background level that you have to
       enter manually if you want a pretty plot
     gif, writes a gif version (and a thumbnail)
 OUTPUTS:
 CALLS: ***
	AVE_CTS, FMT_TIMER [1], FMT_TIMER [2], KILL_INDEX [1], KILL_INDEX [2], LOADCT
	OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], OVER, RD_XDA [1], RD_XDA [2], Rd_Roadmap [2]
	UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
	UTPLOT_IO [1], UTPLOT_IO [2], UTPLOT_IO [3], WRITE_GIF, rd_gbl, rd_gxd [1]
	rd_gxd [2], rd_gxd [3], rd_roadmap [1]
 COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
 MODIFICATION HISTORY:
     31-Aug-00, written by HSH
      1-Sep-00, LF added elementary /ps, /eps switch

   Name: new_dpath

   Purpose: update data paths common area 

   Input Parameters:
      inpath - string or string array of paths to add (or remove)

   Keyword Parameters:
      remove  -  switch, if set, remove instead of add
      home    -  switch, if set, include home directory (add/remove)
      current - switch, if set, include current directory (add/remove)
      prepend - switch, if set, users new stuff comes first
      init    - switch, if set, re-initiailize before adding new path(s)
      empty   - switch, if set, empty before adding new path(s)
		 -- NOTE:initialize=system defaults, empty=empty --
      loud    - switch, if set, extra messages
      file    - two options: if switch set, read file: HOME/setup.ysdpaths
			      if string, read paths from specified file 
      save    - switch, if set, save current dpaths to HOME/setup.ysdpaths
      
   Calling Examples:
      new_dpath			   ; add current directory
      new_dpath,/current,/home		   ; same plus home directory
      new_dpath,'/usr/name/data'	   ; add named path/path arrray
      new_dpath,'/usr/name/data',/prepend ; same, but first, not last
      new_dpath,'/usr/name/data',/remove  ; remove named path/path array
      new_dpath,/init			   ; reinitialize (system defaults)
      new_dpath,/empty,/current	   ; empty and add current
      new_dpath,data_paths(),/remove	   ; same as empty
      new_dpath,/file			   ; adds whatever paths are listed
					   ; in file: HOME/setup.ysdpaths
      new_dpath,/save			   ; save current paths in file:
					   ; HOME/setup.ysdpaths

   Restrictions: proto-routine uses common block, not system variables 
		  todo - add recursive (+path) feature??

 CALLS: ***
	CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], DATA_PATHS, concat_dir [4]
	curdir [1], curdir [2], file_append [1], file_append [2], file_info2 [1]
	file_info2 [2], rd_tfile [1], rd_tfile [2], rem_elem [1], rem_elem [2], uniqo [1]
	uniqo [2]
   Common Blocks:
      data_paths - contains string array of data paths

   History:
      slf -  7-may-93
      slf -  8-may-93	add file keyword and function
      slf -  9-may-93  changed default file name to setyp.ysdpaths
			changed default-to-add-current-path restrictions
      slf - 20-may-93  improved save option (call uniqo,file_delete)
      slf - 22-may-93  add doc header terminator
      slf - 26-may-93  read of DIR_SITE_SETUP/setup.ysdpaths on init

 NAME:
   new_eff_area
 PURPOSE:
   Plot the old and new SXT effective area functions.
 MODIFICATION HISTORY:
   Written, 8-jan-93, J. R. Lemen

 NAME:
      NEW_NORMAL_POINT 
 PURPOSE:
      Generate a new OG#032 (dec) 
 CATEGORY:
      Yohkoh
 CALLING SEQUENCE:
      new_normal_point (no arguments)
      new_normal_point, ogname=['w.751','n1.76'] to force a starting point.
        This form must be used if the tohbans have not properly archived
        the current OG#032, in which case the automatic software will
        start from the wrong coordinates. The "ogname" variable is the
        two-element name of the currently running OG.
      new_normal_point, target = [500,  590] (21-Mar-00; the 
        11-Feb trial position  was [500,592])
 INPUTS:
      Current absolute coordinates, obtained from OG name
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
      target, the desired CCD coordinates of sun center (default [520,590])
      lookback, the number of SFDs in the comparison (default 20)
      nofile, disables filing the fax message
      direct, two-element array giving the SXT_CEN input directly, ie
        bypassing the automatic use of the most recent SFDs
 OUTPUTS:
 COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
      hardwired to [520,590] target (see "target" keyword)
 WARNING:
 MODIFICATION HISTORY:
      19-Nov-1999, written by Hugh Hudson 
      29-Nov-1999, more documentation (HSH)
      11-Feb-2000, added assignable target coordinates and fixed
         yet another Y2K bug (HSH, DMCK)
       3-Mar-2000, hard-wired new coordinates (HSH)
      21-Mar-2000, tweaked to [520, 590] because of Feb. 00 problems (HSH)
      10-Apr-2000, added a keyword for the look-back image number in
         the SFD reference
       7-Jun-2000, nofile and direct keywords added 

       NAME: new_win


       PURPOSE: make a new window with the next available #. Default
		 positioning is 


       METHOD: Find largest current number and add 1. The default passes
		things to wdef. With the OFFSET keyword you can position
		the new window relative to the old.


       CALLING SEQUENCE: new_num=new_win(xs,ys)[,offset=ofst]


       PARAMETERS:   	xs,ys		sizes of new window
		      	new_num		number of new window

	KEYWORDS: 	offset		offset the new window relative 
					to the last one.

 CALLS: ***
	WDEF [1], WDEF [2]
 CALLED BY:
	ar_prep, box_lc_plot, dejitter, mode_fit, rd_lonlat, sel_region
       HISTORY: Drafted March 17, 1992, A. McAllister,
		 [Code mostly copied from SXT_GRID]

   Name: newfiles

   Purpose: return file name(s) of new data 

   Calling Sequence:
	files=newfiles( [nfiles, /sfd, /sfr, /spd, /spr, /bda, /wda, /hda
			 	,oldfiles=oldfiles]
   
   Optional Keyword Parameters:

 CALLS: ***
	STR2ARR [1], STR2ARR [2], file_list [1], file_list [2], gbo_paths [1]
	gbo_paths [2], uniq_fid [1], uniq_fid [2], wmenu_sel [1], wmenu_sel [2]
 CALLED BY:
	FIRST_LIGHT [1], FIRST_LIGHT [2], LATEST_SFDS_GIF, QUICKSTRAY [2], chk_pointing
	daily_forecast [2], first_bcs, fl_goesplot [1], fl_goesplot [2], fl_goesplot [3]
	fl_sxtobsnar [1], fl_sxtobsnar [2], last_lc, lastspd, ref_term [2], sxt2mpeg
	xspr [1], xspr [2]
   History: 
	10-sep-1993 (SLF) extracted and genericized code from newsfd.pro
	 4-Mar-1994 (SLF) print warning if no files
	 9-Mar-1994 (SLF) fix bug when last keyword not used
       16-sep-1994 (SLF) expanded to GBO data
        6-apr-1995 (SLF) add ONE keyword (pass 2 wmenu_sel)
        8-aug-2000 (SLF) add sss and ssc keywords (prefix)
        10-Aug-00  (DMZ) added patch for SFD files that use week ID's
                         (e.g. sfd92_10a.12)

   Name: newfiles

   Purpose: return file name(s) of new data 

   Calling Sequence:
	files=newfiles( [nfiles, /sfd, /sfr, /spd, /spr, /bda, /wda, /hda
			 	,oldfiles=oldfiles]
   
   Optional Keyword Parameters:

 CALLS: ***
	STR2ARR [1], STR2ARR [2], file_list [1], file_list [2], gbo_paths [1]
	gbo_paths [2], uniq_fid [1], uniq_fid [2], wmenu_sel [1], wmenu_sel [2]
 CALLED BY:
	FIRST_LIGHT [1], FIRST_LIGHT [2], LATEST_SFDS_GIF, QUICKSTRAY [2], chk_pointing
	daily_forecast [2], first_bcs, fl_goesplot [1], fl_goesplot [2], fl_goesplot [3]
	fl_sxtobsnar [1], fl_sxtobsnar [2], last_lc, lastspd, ref_term [2], sxt2mpeg
	xspr [1], xspr [2]
   History: 
	10-sep-1993 (SLF) extracted and genericized code from newsfd.pro
	 4-Mar-1994 (SLF) print warning if no files
	 9-Mar-1994 (SLF) fix bug when last keyword not used
       16-sep-1994 (SLF) expanded to GBO data
        6-apr-1995 (SLF) add ONE keyword (pass 2 wmenu_sel)
        8-aug-2000 (SLF) add sss and ssc keywords (prefix)
        10-Aug-00  (DMZ) added patch for SFD files that use week ID's
                         (e.g. sfd92_10a.12)

 NAME       : newkirk (function)
 PURPOSE    :
	perform Newkirk filter operation on 10-cm coronagraph data
 CATEGORY :
	idl/lib/nkr
 CALLING SEQUENCE :
  img = newkirk(img0,norder=norder,dh=dh,nrep=nrep,verbose=ver, 	$
	flat=flat,od=od,round=round,skyint=skyint,sclbyt=sclbyt,	$
	trend=trend,com=com,sbt=sbt,sample=sample,skycoeff=skycoeff,	$
	fltcoeff=fltcoeff,i0=i0,enlarge=enlarge)
 INPUTS :
 	img0     -- image data assumed to be FIX type
 OPTIONAL INPUT PARAMETERS : 
	none
 KEYWORD PARAMETERS :
 	norder  -- order of polinomial (default=5)
 	dh      -- skin height above the limb (default=2)
	nrep    -- repeat no. of fitting procedure for lower edge of sky
			 (default=0)
 	/ver    -- display fitting curve
 	flat    -- flat image if make correction before newkirk
 	/od     -- use occulting disk edge to detect the limb
 	skyint  -- return sky intensity
	sclbyt(2) -- scale 0-255 with range of sky*[(1+sclbyt(0))~(1+sclbyt(1))]
	trend   -- remove linear trend of sky and set height from the limb
	com     -- return ox,oy,r information
	/sbt    -- make subtraction other than division
	sample  -- sampling step of sky (default=10)
	skycoeff -- return coefficients of sky fitting polinomial
	fltcoeff -- coefficients of sky fitting polinomial for flat data
 	      if sbt is set, re-correct flat inhomogeneity
	i0      -- return central intensity for calibration
	/enlarge -- operation for enlarged image ('left', 'right')
 OUTPUTS :
	none
 CALLS: ***
	CENTER, CONGRID [1], CONGRID [2], CONGRID [3], GET2EDGES, LIMBCENTER, MKROUND, NEWKIRK
	ODCENTER, ODCENTER2, POLY, POLY_FIT, SKYPOINT, TREND1
 COMMON BLOCKS : 	none
 SIDE EFFECTS :
	if /ver is set, fitting curve is displayed
 RESTRICTIONS :	none
 PROCEDURE :
 MODIFICATION HISTORY :
	K.I. 1991/12/31	
	K.I. 1992/04/22

 NAME:		NEWMOVIE
 PURPOSE:	Show a cyclic sequence of images stored in a 3D array.
               This version is controlled by buttons.
               It features a single step mode, and displays frame number
               and (optionally) frame time.
 CATEGORY:	Image display
 CALLING SEQUENCE:
	newmovie, Images [,Rate]
 INPUTS:
	Images = (n,m,nframes) byte array of image data, consisting of
		nframes images, each of n by m.  This array should be
		stored with the top row first, (order = 1) for maximum 
		efficiency.
 OPTIONAL INPUT PARAMETERS:
	Rate = initial rate, in APPROXIMATE frames per second.  If 
		omitted, the inter-frame delay is set to 0.01 second.
 KEYWORD PARAMETERS:
	Order = image ordering.  Order = 0 for images ordered bottom up,
		or = 1 for images ordered top down (default : 0).
               Order = 1 is somewhat more efficient.
       time = string vector. Contains time for each frame.
 OUTPUTS:
	No explicit outputs.
 COMMON BLOCKS:
	None.
 SIDE EFFECTS:
	Images are displayed in the lower left corner of the currently
		selected window.
       At high image rates the frame number is not updated.
 RESTRICTIONS:
	As SunView has no zoom or pan, we have to write each image to
	the display. Restricting the maximum rate, etc.  Experience has
	shown that you can count on a rate of approximately 10 frames / second
	with 192 by 192 images.  This varies according to the type of
	computer, amount of physical memory, and number of frames.

	The amount of available memory also restricts the maximum amount
	of data	that can be displayed in  a loop.
 PROCEDURE:
	Straightforward.
 MODIFICATION HISTORY:
	DMS, Nov, 1988.  Original movie procedure
       JPW, Feb, 1992.  Added display of image number and single step mode

   Name: newsfd

   Purpose: read in sfd data files from temporary reformatted data area

   Calling Sequence:
      data=newsfd(index, info)		; return most recent image
      data=newsfd(index, info, /select); menu selection of file
      data=newsfd(index, info, last=4) ; most recent 4 images are read
      data=newsfd(index, info, /sfr)	; use sfr, not sfd files
   
   Optional Keyword Parameters:
      all  - if set, all sfds in temporary area are read
      if neither keyword is set, x menu selection is used
      last - number of images to read (default=1)
      nimgs - synonym for last (default=1)
      sfds  - (output) string/string vector containg sfd files read
      sfr   - if set, use sfr files, not sfd files
      spd   - if set, use spd files

   Output Parameters:
      index - index structure array
      info  - text description 
      The function return value is the data image cube

   User Notes:
      Since this is a function, the following usage is acceptable:
         xstepper,newsfd(index,info,last=4),info
	
 CALLS: ***
	FILE2TIME, GET_INFO2, RD_SDA [1], RD_SDA [2], RD_SDA [3], Rd_Roadmap [2]
	Rd_fHeader [1], Rd_fHeader [2], Rd_fHeader [3], file_list [1], file_list [2]
	get_logenv [1], get_logenv [2], gt_filtb, gt_res, prstr [1], prstr [2]
	rd_roadmap [1], wmenu_sel [1], wmenu_sel [2]
 CALLED BY:
	background_new, lastsfd [1], lastsfd [2], sxt_summary
   Restrictions: 
      sfd files must exist on site specific DIR_SITE_NEWDATA 

   History: slf 21-Oct-92
	     slf 25-Oct-92 - rebin quarter res to 512,512
	     slf  8-may-93 - add keywords nimg, sfds, last now imgs not files
	     slf 11-may-93 - add sfr keyword (use sfr... not sfd...)
	     slf 17-may-93 - add sfm keyword
	     slf 31-aug-93 - add oldfile keyword
	     slf 23-Jun-94 - bug in select option
            slf 10-apr-95 - prob with /half but only Al.1 (typo)
            slf 22-aug-95 - env LASTSFD_DEFAULT for filter (def AlMg)
            slf 12-nov-95 - dont rebin 256 if rebin=256

 NAME:
       NEXTTIME
 PURPOSE:
       Convert the Time in the unit of HHMMSS(String*6) into Time =
     Time + 1(sec) in the unit of HHMMSS(string*6).
       example: NEXTTIME return 123141 if 123140 is inseted.
 CATEGPRY:
       PROGRAMMING
 CALLING SEQUENCE:
       NEXTTIME, time
 INPUTS:
       time = inputted time in the unit of HHMMSS(String*6)
 OUTPUTS:
       time = converted time in the unit of HHMMSS(String*6)
 COMMON BLOCKS:
       None.
 SIDE EFFECTS:
       None.
 RESTRICTIONS:
       None.
 PROCEDURE:
       Convert TIME --> TIME + 1(sec) in the unit of HHMMSS(String*6).
 PROGRAMER:
       Toshifumi Shimizu
           Department of Astronomy, University of Tokyo, JAPAN
 HISTORY:
       Ver 1.0   5-JAN-1991  made

 NAME:
       NEXTTIME
 PURPOSE:
       Convert the Time in the unit of HHMMSS(String*6) into Time =
     Time + 1(sec) in the unit of HHMMSS(string*6).
       example: NEXTTIME return 123141 if 123140 is inseted.
 CATEGPRY:
       PROGRAMMING
 CALLING SEQUENCE:
       NEXTTIME, time
 INPUTS:
       time = inputted time in the unit of HHMMSS(String*6)
 OUTPUTS:
       time = converted time in the unit of HHMMSS(String*6)
 COMMON BLOCKS:
       None.
 SIDE EFFECTS:
       None.
 RESTRICTIONS:
       None.
 PROCEDURE:
       Convert TIME --> TIME + 1(sec) in the unit of HHMMSS(String*6).
 PROGRAMER:
       Toshifumi Shimizu
           Department of Astronomy, University of Tokyo, JAPAN
 HISTORY:
       Ver 1.0   5-JAN-1991  made

 NAME:
       NEXTTIME
 PURPOSE:
       Convert the Time in the unit of HHMMSS(String*6) into Time =
     Time + 1(sec) in the unit of HHMMSS(string*6).
       example: NEXTTIME return 123141 if 123140 is inseted.
 CATEGPRY:
       PROGRAMMING
 CALLING SEQUENCE:
       NEXTTIME, time
 INPUTS:
       time = inputted time in the unit of HHMMSS(String*6)
 OUTPUTS:
       time = converted time in the unit of HHMMSS(String*6)
 COMMON BLOCKS:
       None.
 SIDE EFFECTS:
       None.
 RESTRICTIONS:
       None.
 PROCEDURE:
       Convert TIME --> TIME + 1(sec) in the unit of HHMMSS(String*6).
 PROGRAMER:
       Toshifumi Shimizu
           Department of Astronomy, University of Tokyo, JAPAN
 HISTORY:
       Ver 1.0   5-JAN-1991  made

 NAME:
       NEXTTIME
 PURPOSE:
       Convert the Time in the unit of HHMMSS(String*6) into Time =
     Time + 1(sec) in the unit of HHMMSS(string*6).
       example: NEXTTIME return 123141 if 123140 is inseted.
 CATEGPRY:
       PROGRAMMING
 CALLING SEQUENCE:
       NEXTTIME, time
 INPUTS:
       time = inputted time in the unit of HHMMSS(String*6)
 OUTPUTS:
       time = converted time in the unit of HHMMSS(String*6)
 COMMON BLOCKS:
       None.
 SIDE EFFECTS:
       None.
 RESTRICTIONS:
       None.
 PROCEDURE:
       Convert TIME --> TIME + 1(sec) in the unit of HHMMSS(String*6).
 PROGRAMER:
       Toshifumi Shimizu
           Department of Astronomy, University of Tokyo, JAPAN
 HISTORY:
       Ver 1.0   5-JAN-1991  made

 NAME       : nkrget (procedure)
 PURPOSE    :
	read header and images from norikura formatted file
 CATEGORY :
	idl/lib/nkr
 CALLING SEQUENCE :
	nkrget,file,h,img,img_no=ino,nset=nset,all=all,/quiet
 INPUTS :
	file   : filename [may be an array]
 OPTIONAL INPUT PARAMETERS : 
	none
 KEYWORD PARAMETERS :
 	h	 : header structure
 	img	 : image array
 	img_no   : specify image no. for nset>1. 
		   if img_no=-1, return only nset
 	nset     : return number of images
 	all      : extract all images in img(*,*,*)
	quiet	 : don't print any message
 OUTPUTS :
	none
 CALLS: ***
	CONGRID [1], CONGRID [2], CONGRID [3], HSWAB, NKRGET, NKRGETH, NKRHEAD, NKRREAD1
 COMMON BLOCKS : 	none
 SIDE EFFECTS : 	none
 RESTRICTIONS :	none
 PROCEDURE :
	If file is an array, the results are saved in img(*,*,*)
 MODIFICATION HISTORY :
	K.I. '91/10/31	
	K.I. '92/04/19	; pos_info
	K.I. '93/01/20	; allow '~/'
	K.I. '93/01/21	; allow float image
	K.I. '93/01/24	; multiple file OK
	K.I. '93/01/26	; fact keyword
	K.I. '93/01/31	; quiet keyword,  data_type ge 3 --> fltarr
	K.I. '93/02/07	; for h.nf > 1
	K.I. '93/04/29	; quiet keyword
	K.I. '93/10/07	; file=filea(0)

 NAME       : nkrhd_info (function)
 PURPOSE    :
	get nkrhead for SXT data from info_array(i) of test_rd  as nkrhead
 CATEGORY :
	idl/lib/nkr
 CALLING SEQUENCE :
	h = nkrhd_info(info_array)
 INPUTS :
 	info_array     -- test_rd info_array
 OPTIONAL INPUT PARAMETERS : 
	none
 KEYWORD PARAMETERS :
 OUTPUTS :
	none
 CALLS: ***
	NKRHD_INFO, NKRHEAD
 COMMON BLOCKS : 	none
 SIDE EFFECTS :	none
 RESTRICTIONS :	none
 PROCEDURE :
 MODIFICATION HISTORY :
	23 Jan.93  K.I.

 NAME       : nkrhead (function)
 PURPOSE    :
	define norikura type header structure
 CATEGORY :
	idl/lib/nkr
 CALLING SEQUENCE :
	h = nkrhead()
 INPUTS :
 CALLS: ***
	NKRHEAD
 MODIFICATION HISTORY :
	K.I. '93/01/31	

ROUTINE:
	nn_rd2_sxt_noaa
PURPOSE:
	Read a specified sfr file and create three arrays, 512x512,
	1. Al1265, MBE=16, 2. Al1265, MBE=6.
	--->  Changed to Al/Mg/Mn               
       --->  Changed back to Al1265
	--->  Changed back to Al/Mg/Mn MBE=8 and 18 
	This is a software specific to a purpose of sending images to
	NOAA and CRL routinely.
CALLING SEQUENCE:
	NN_RD2_SXT_NOAA
INPUT PARAMETER:
       INFILE: Full name of the SFR file
SIDE EFFECT:
	Four files are created on /3p/nitta/raw_data.
 CALLS:

ROUTINE:
	nn_rd3_sxt_noaa
PURPOSE:
	Read a specified sfr file and create three arrays, 512x512,
	1. Al1265, MBE=16, 2. Al1265, MBE=6.
	--->  Changed to Al/Mg/Mn               
       --->  Changed back to Al1265
	--->  Changed back to Al/Mg/Mn MBE=8 and 18 
	This is a software specific to a purpose of sending images to
	NOAA and CRL routinely.
CALLING SEQUENCE:
	NN_RD2_SXT_NOAA
INPUT PARAMETER:
       INFILE: Full name of the SFR file
SIDE EFFECT:
	Four files are created on /3p/ftp/pub.
 CALLS:

ROUTINE:
	nn_rd_sxt_noaa
PURPOSE:
	Read a specified sfr file and create three arrays, 512x512,
	1. Al1265, MBE=16, 2. Al1265, MBE=6.
	This is a software specific to a purpose of sending images to
	NOAA and CRL routinely.
CALLING SEQUENCE:
	NN_RD_SXT_NOAA
INPUT PARAMETER:
       INFILE: Full name of the SFR file
SIDE EFFECT:
	Two files are created on /3p/nitta/raw_data.
 CALLS:

 CALLS: ***
	ARR2STR [1], Arr2Str [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
	anytim2ex [1], anytim2ex [2], anytim2ints [1], anytim2ints [2], concat_dir [4]
	ex2fid [1], ex2fid [2], gt_day [1], gt_day [2], mail [1], mail [2], prstr [1]
	prstr [2], rd_tfile [1], rd_tfile [2]
HISTORY:
	Written 4-Oct-93 by M.Morrison

NAME:
	nob_name_fix
PURPOSE:
	Convert the Nobeyama filename to the Yohkoh convention.  Scale
	the image down to byte type
CALLING SEQUENCE:
	nob_name_fix, '/yd8/nobeyama_raw/s930923_14391acp.fits
INPUT:
	infil	- The input file(s).  There must be a "acp.fits"
		  or an "scp.fits"
OPTIONAL INPUT:
       outdir  - The output directory to move the files to.
                 If not passed, it calls FIND_GBO_DIR to find
                 the proper GBO weekly directory.
OPTIONAL OUTPUT:
	outfil	- The new filename for the input file(s)
 CALLS: ***
	BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2]
	CONCAT_DIR [3], FIND_GBO_DIR, RFITS [1], RFITS [2], RFITS [3], anytim2ex [1]
	anytim2ex [2], break_file [4], concat_dir [4], ex2fid [1], ex2fid [2]
	wrt_fits [1], wrt_fits [2]
HISTORY:
	Written 4-Oct-93 by M.Morrison

   Name: nobeyama_update

   Purpose: update nobeyama data base

   Keyword Parameters
      nodir - if set, do not ftp directory files (use existing or supplied)
      directory - directory file(s) - [lists of bfYYMMDD.fits names]
      mk_week_file - if set, regenerate associated nts files      

 CALLS: ***
	BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2]
	CONCAT_DIR [3], FILE_EXIST [2], break_file [4], concat_dir [4], deriv_arr [1]
	deriv_arr [2], file_exist [1], file_exist [3], ftp_copy [1], ftp_copy [2]
	get_logenv [1], get_logenv [2], nts_copy [1], nts_copy [2], prstr [1], prstr [2]
	rd_tfile [1], rd_tfile [2], tbeep [1], tbeep [2], tbeep [3], uniqo [1], uniqo [2]
	wc_where [1], wc_where [2]
   History:
      13-Apr-1994 (SLF) Written

PURPOSE:
   To compute nonparametric fits of x and y data to the function
   y = x^alpha.  Both "t-statistic" and "D-statistic" methods are
   used.  Nonparametric analyses are appropriate for data that 
   are not normally distributed (i.e., individual measurements are
   samples from parent distributions that are not Gaussian).  The
   method is described in Lee et al. (1993, ApJ, 412, 401) and Efron 
   and Petrosian (1992, ApJ, 399, 345).

   The most probable value of alpha and the 90% confidence level alphas 
   are printed to the screen.  The 90% confidence level alphas define
   the range in which the actual alpha is likely to fall (with 90% 
   confidence).

INPUT:
   xdata, ydata = x and y data arrays
   yweight = array of weighting factors for ydata
   alphamin, alphamax = minimum and maximum of the range of alpha
       values to be searched.

OUTPUT:
   alpha = array of alpha values
   tw = array of weighted "t-statistics" for the data set (x, y/x^alpha)
   d = array of "D-statistics" for the data set (x, y/x^alpha)
   probtw = array of probabilities of obtaining tw if the data are 
       uncorrelated
   probd = array of probabilities of obtaining d if the data are
       uncorrelated

 CALLS: ***
	correlation
HISTORY:
   26-apr-94, written, J. A. Klimchuk
   18-may-95, commented out print statements for D-statistic, JAK

NAME:
	norik_name_fix
PURPOSE:
	Convert the Norikura filename to the Yohkoh convention.  Also
	rotate image to put north up as default on IDL screen.
CALLING SEQUENCE:
	norik_name_fix, '/f7p/hara/930927/r1.i0', /move
INPUT:
	infil	- The input file(s).  The leading characters
		  must be "r", "g", or "y"
OPTIONAL INPUT:
       outdir  - The output directory to move the files to.
                 If not passed, it calls FIND_GBO_DIR to find
                 the proper GBO weekly directory.
OPTIONAL OUTPUT:
	outfil	- The new filename for the input file(s)
 CALLS: ***
	BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2]
	CONCAT_DIR [3], FIND_GBO_DIR, RFITS [1], RFITS [2], RFITS [3], anytim2ex [1]
	anytim2ex [2], break_file [4], concat_dir [4], ex2fid [1], ex2fid [2]
	mk_fits_head [1], mk_fits_head [2], wrt_fits [1], wrt_fits [2]
HISTORY:
	Written 1-Oct-93 by M.Morrison

 PURPOSE:  Exposure-normalize images, and byte scale them with
           a power.
 CALLING SEQUENCE:  data_norm=norm_sxt_nn(index,data,0.2,1.5)

 CALLS:

 NAME:
       NOV93_ECLIPSE	
 PURPOSE:
	enter Soma-san's calculations for the Nov 1993 eclipse
 CATEGORY:
	Yohkoh
 CALLING SEQUENCE:
	nov93_eclipse, data
 INPUTS:
	none. Uses hard-wired data file
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
       plot makes a plot of radial distance, with fem_grid
 OUTPUTS:
	structure with many things,
	among which eclipse.xx is the EAST departure
                   eclipse.yy is the NORTH departure
                   eclipse.daylight is a D/N flag
 CALLS: ***
	DAY_IN, OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], RDFILE [1], RDFILE [2], UTPLOT [1]
	UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], fem_grid [1]
	fem_grid [2], get_rb0p [1], get_rb0p [2], rd_fem
 MODIFICATION HISTORY:
	HSH written 10-nov-94 (based on NOV_ECLIPSE.PRO)

 NAME:
       NOV_ECLIPSE	
 PURPOSE:
	enter Soma-san's calculations for the Nov 1994 eclipse
 CATEGORY:
	Yohkoh
 CALLING SEQUENCE:
	nov_data, data
 INPUTS:
	none. Uses hard-wired data file
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
       plot makes a plot of radial distance, with fem_grid
 OUTPUTS:
	structure with many things,
	among which eclipse.xx is the EAST departure
                   eclipse.yy is the NORTH departure
                   eclipse.daylight is a D/N flag
 CALLS: ***
	DAY_IN, OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], RDFILE [1], RDFILE [2], UTPLOT [1]
	UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], fem_grid [1]
	fem_grid [2], get_rb0p [1], get_rb0p [2], rd_fem
 MODIFICATION HISTORY:
	HSH written 8 April 1994
	HSH added d/n and saa 10 April 1994
	HSH updated to Nov eclipse 13-May-94
       HSH added /plot 28-Sep-94
       HSH hard-wired "final" Soma predictions from 18-Oct elements

 CALLS: ***
	BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2]
	CONCAT_DIR [3], break_file [4], concat_dir [4], ex2week [1], ex2week [2]
	fid2ex [1], fid2ex [2], file_list [1], file_list [2], mail [1], mail [2]
	mk_week_file [1], mk_week_file [2], prstr [1], prstr [2], rd_tfile [1]
	rd_tfile [2], where_arr [1], where_arr [2]
 CALLED BY:
	nobeyama_update
HISTORY:
	Written 14-May-93 by M.Morrison
       14-May-93 (MDM) - Changed to use $ydb_raw and $ys
       26-May-93 (MDM) - Changed to use /ys/site/script instead of /ys/gen/script
	19-Jun-93 (MDM) - Modified to avoid file that do not begin with "bf"
			  and do not have ".fits"  as the extension
       19-Aug-93 (MDM) - Removed /self switch from Mail call
	 3-Oct-93 (MDM) - Patch since "ls *" command was not working
			  since 21-Aug for some reason
	 8-Nov-93 (MDM) - Patch since there were too many files in the
			  nobeyama_raw directory
	 3-Jan-93 (MDM) - Modified mailing list
       13-Apr-94 (SLF) - dont copy directory if files passed in.
       15-Mar-98 (SLF) - dusted off on request - simplify 'needed' logic
                         (still an  ancient routine doomed to fail...)

 CALLS: ***
	BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2]
	CONCAT_DIR [3], break_file [4], concat_dir [4], ex2week [1], ex2week [2]
	fid2ex [1], fid2ex [2], file_list [1], file_list [2], mail [1], mail [2]
	mk_week_file [1], mk_week_file [2], prstr [1], prstr [2], rd_tfile [1]
	rd_tfile [2], where_arr [1], where_arr [2]
 CALLED BY:
	nobeyama_update
HISTORY:
	Written 14-May-93 by M.Morrison
       14-May-93 (MDM) - Changed to use $ydb_raw and $ys
       26-May-93 (MDM) - Changed to use /ys/site/script instead of /ys/gen/script
	19-Jun-93 (MDM) - Modified to avoid file that do not begin with "bf"
			  and do not have ".fits"  as the extension
       19-Aug-93 (MDM) - Removed /self switch from Mail call
	 3-Oct-93 (MDM) - Patch since "ls *" command was not working
			  since 21-Aug for some reason
	 8-Nov-93 (MDM) - Patch since there were too many files in the
			  nobeyama_raw directory
	 3-Jan-93 (MDM) - Modified mailing list
       13-Apr-94 (SLF) - dont copy directory if files passed in.

 NAME:
	nuden_2nd
 PURPOSE
	To make a second-order leak correction for exposures that
	use the nuden filter.  Such exposures, with DPE= 7 or 8,
	in the second half of CY2001.
 PROCEDURE
	Average of rows 12/res through 42/res is subtracted from every
	row in the image.
 CALLING SEQUENCE
	better=nuden_2nd(instruct,indata,row=row)
 INPUT
	instruct, SXT index structure
	indata, sxt_prep'd image or data cube
 OUTPUT
	data with nuden second-order leak correction applied.
 OPTIONAL KEYWORD OUTPUT
	row, the correction vector (used by sxt_deleak)
 RESTRICTIONS
	indata must be integer, longword, or floating point.
	Correction will only be applied to images utilizing
	the neutral density filter.
 HISTORY
	 1-Apr-2004 LWA  (Not an April Fool joke!)
	1b-Apr-2004 LWA  Mod to handle nuden's with any DPE.
       16-Apr-2004 SLF Change name NuDen_2nd -> nuden_2nd.pro
	22-Apr-2004 LWA  removed NuDen's from header.
 CALLS:
 CALLED BY
	sxt_deleak [1], sxt_deleak [2]

  NAME:
	nvsstat
  PURPOSE:
	Spawn a call to the script nvsstat to get current position.
  CALLING SEQUENCE:
	frame_number = nvsstat()
	frame_number = nvsstat(cur_min,cur_sec,cur_frame,next_frame=next_frame)
	frame_number = nvsstat(cur_min,cur_sec,cur_frame,text=text,/nodisplay)
  OPTIONAL KEYWORD INPUTS:
	echo = If set, will echo the script message to the screen.
  OPTIONAL OUTPUT PARAMETERS:
	cur_min	  = Current position: minutes
	cur_sec   = Current position: sec
	cur_frame = Current position: Frame number (0-29)
	text	  = The text array from the script file call.
	next_frame= String with the time of the next frame (MMSSFF).

 CALLS: ***
	input [1], input [2]
 CALLED BY:
	go_lasdisk2 golaserdisk, go_nvs5
 RESTRICTIONS:
 	The environment variable NVSHOST should be set to 133.74.8 before
	running this routine.
 RETURNS:
	frame_number is the index of the laser disk.  A negative value
	indicates that the position could not be obtained.

 MODIFICATION HISTORY:
  Written,   1-sep-92, J. R. Lemen
  Mod, 15&16-Mar-93, GAL for frame numbers > max 2-byte integer value.
  11-oct-93, JRL, Changed to script_dir to $ys (so backup IDL software will work).
  10-oct-95, LS, changed script_dir to reflect changes made for SOHO
		  compatibilty re-arrangement; used environment variable

FUNCTION: nzmean,image,good_pixmap

The purpose of this program is to prepare an image with bad pixels
replaced by the mean of adjacent good pixels.  The good pixels
are identified by the binary array good_pixmap.  The result of the
operation is given in the new array, img.

	2/6/99, LWA
       27-aug-1999 - S.L.Freeland - subscript protection (KK)
 CALLS: