ibclrf $SMEI/ucsd/sat/idl/toolbox/ibclrf.pro
[Previous] [Next]
 NAME:
	ibclrf
 PURPOSE:
	Takes a longword (32-bit) integer, and clears (sets to zero) a bit
 CATEGORY:
	Bits and bytes
 CALLING SEQUENCE:
	Result = ibclrf(bitset, bitpos)
 INPUT:
	bitset	longword    A longword integer, one of it
		    bits will be clear or zero
	bitpos	longword    The position of the bit that
		    will be cleared; 0 <= bitpos <= 31
 OUTPUT:
	Result	longword    Same as bitset, (bitwise) but
		    at bit position, bitpos, it will
		    definitely be a zero
 INCLUDE:
	@compile_opt.pro    ; On error, return to caller
 CALLED BY:
	ice_pack
 SEE ALSO:
	ibsetf, ibtestf
 MODIFICATION HISTORY:
	8/23/00 	Kevin Nguyen


ibsetf $SMEI/ucsd/sat/idl/toolbox/ibsetf.pro
[Previous] [Next]
 NAME:
	ibsetf
 PURPOSE:
	Sets a bit in a longword (32-bit) integer
 CATEGORY:
	Bits and bytes
 CALLING SEQUENCE:
	Result = ibsetf(bitset, bitpos)
 INPUT:
	bitset longword integer 32 bits longword integer.
		    Its bit at bit position bitpos
		    will be set at 1
	bitpos longword integer The bit position of bit to be
		    set as 1 (0 <= bitpos <= 31)
 OUTPUT:
	Result longword integer Same at bitset bitwise except
		    at bit position bitpos
 INCLUDE:
	@compile_opt.pro    ; On error, return to caller
 CALLED BY:
	ice_pack, ice_unpack, smei_frm_flag
 SEE ALSO:
	ibclrf, ibtestf
 MODIFICATION HISTORY:
	8/23/00     Kevin Nguyen


ibtestf $SMEI/ucsd/sat/idl/toolbox/ibtestf.pro
[Previous] [Next]
 NAME:
	ibtestf
 PURPOSE:
	Tests if the bit at bitpos of a longword (32-bit) integer is set
 CATEGORY:
	Bits and bytes
 CALLING SEQUENCE:
	Result = ibtestf(bitset, bitpos)
 INPUT:
	bitset	longword integer  will be used to test if
		      bit position, bitpos is zero or one
	bitpos	longword integer  the bit postion, bit pos, being tested
 OUTPUT:
	ibtestf     scalar; type:  byte
		1B if bit is set;  0B bit is clear
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLED BY:
	ice_pack, ice_unpack, smei_frm_flag
 SEE ALSO:
	ibclrf, ibsetf
 MODIFICATION HISTORY:
	8/23/00     Kevin Nguyen


ice_pack $SMEI/ucsd/sat/idl/util/ice_pack.pro
[Previous] [Next]
 NAME:
	ice_pack
 PURPOSE:
	Compresses an integer array using a modified Rice algorithm
 CATEGORY:
	Bits and bytes
 CALLING SEQUENCE:
	FUNCTION ice_pack, Orig, Pack, kmax=kmax, kshift=kshift, lenbit=lenbit, sign=sign, perm=perm
 INPUTS:
	Orig	array; type: long integer
		array to be compressed
 OPTIONAL INPUTS:
	kmax	scalar; type: integer; default: 8
		cutoff for applying compression (see PROCEDURE) (kmax > 0)
	kshift	scalar; type: integer; default: 0
		minimum bits copied for each difference
	lenbit	scalar; type: integer; default: 16
		# bits encoded when full value instead of a difference
		is used (see PROCEDURE)
	sign	scalar; type: integer; default: 1
		-1: all values in Orig < zero
		+1: all values in Orig >= zero
		 0: both positive and negative values are present
	perm	array[0:kmax-kshift]; type: integer
		permutation table (see PROCEDURE)
 OUTPUTS:
	ibit	scalar; type: long integer
		# bits in Pack used to hold the compressed data;
		on failure ibit = 0 is returned; this happens only if
		the compressed array Pack is bigger than the input
		array Orig.
	Pack	array; type: long integer
		compressed data;
		The # elements in Pack is 1+(ibit-1)/nbit.
		Bits 0..ibit-1 are filled with compressed data.
		For the last element only bits 0 .. ib-1 are
		used [ib = 1+ ((ibit-1) mod nbit)].
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	InitVar, ibclrf, ibsetf, ibtestf
 SEE ALSO:
	ice_unpack
 PROCEDURE:
	The technique is based on the algorithm described by Michael W. Richmond and Nancy E.
	Ellman, in March 1996 paper which we have in preprint form.
	COMPRESSION ALGORITHM
	---------------------
	The Rice algorithm stores an integer in the following bit pattern:
	    0.......0 1 x.......x
	    |_______|	|_______|
	      ipad	itop + 1
	a string of ipad 0-bits, followed by a terminating 1 bit, followed by a string of itop + 1
	bits encoding the integer value.

	The actual number encoded is either the full value from the input array, or the difference
	with the previous element in the array (!!! for the first element in the array, the 'previous
	element' is assumed to be zero).

	The two cases are distinguished  using the kmax value: if the absolute difference is less than
	2^(kmax - 1) then the difference is encoded; if not then the full value is encoded

	CASE A: Non-zero absolute difference less than 2^(kmax - 1)
	-----------------------------------------------------------
	2^(kmax -1) has only bit kmax - 1 set. I.e. for absolute differences less than 2^(kmax - 1) this
	bit, and all higher bits, are NOT set, i.e. only the first kmax - 1 bits 0..kmax - 2 might be set.

	Let itop be the highest 1-bit of the absolute difference. For a non-zero difference
	0 <= itop <= kmax-2. Since bit itop is by definition set, the absolute difference
	is fully described by the first itop bits 0...itop-1. An additional bit is needed
	to encode the sign of the difference. So itop+1 bits are needed to store the difference:
	first ipad=itop+1 0-bits, followed by a terminator 1-bit, indicate how many bits are
	used. The terminator bit is then followed by the first itop bits of the absolute difference
	followed by a sign bit. The total number of bits used is 2*itop+3.

	    0..........0 1 x..........x
	    |__________|   |__________|
	    ipad=itop+1       itop+1

	CASE B: Zero difference
	-----------------------
	No bit at all is needed to store a zero, except for a terminator bit. Not that this
	is CASE A with ipad = 0, itop = -1.

	CASE C:  Absolute difference greater/equal 2^(kmax-1)
	-----------------------------------------------------
	In this case the full value instead of the difference is encoded.
	For case A the maximum number of leading 0-bits is kmax-1. A string of kmax 0-bits
	is used to indicate that a full value is encoded. After the terminator bit follow
	the first lenbit bits. The total number of bits used is kmax+1+lenbit

	    0..........0 1 x..........x
	    |__________|   |__________|
	    ipad=kmax	      lenbit

	MINIMUM NUMBER OF BITS
	----------------------

	The above describes the the compression algorithm for the case kshift=0.
	Setting kshift to non-zero value modifies the way differences are stored: for each
	difference the kshift lowest bits are always stored. The differences are now encoded
	as follows:

	Zero difference (ipad = 0)
	    1 0
	The 1-bit is the terminating bit (i.e. no leading 0-bits), followed by a 0-bit to
	indicate the presence of a zero difference.

	Difference needing 1..kshift bits (1 <= ipad <= kshift)
	    1 1 s x.....x
		  |_____|
		   kshift
	The leading 1-bit is the terminating bit (i.e. no leading 0-bits). It is followed by
	a 1-bit to be able to distinguish it from a zero difference. Then a bit follows
	storing the sign of the difference (s-bit). Then follow the lowest kshift bits.

	Differences needing kshift+1...kmax-1 bits (kmax-1 <= ipad <= kmax-1 above)
	    0.......0 1 x........x
		|_______|   |________|
	       ipad-kshift     ipad

	Full value (ipad=kmax)
	    0.......0 1 x........x
		|_______|   |________|
	       kmax-kshift    lenbit


	PERMUTATION TABLE
	-----------------

	The number of leading 0-bits varies from izero=0 to izero=kmax-kshift. These differences
	serve merely as identifiers for the type of information stored after the terminating
	1-bit. Rather than using them directly as described above, it is more efficient to
	use a permutation table based on the histogram of the number of bits needed for to store
	the differences. E.g. if for some sky image 3-bit difference occur most, followed by 2-bit,
	then 1-bit, then 0-bit, then 4-bit, 5-bit, etc. then the permutation table would be
	    perm = 3,2,1,0,4,5...kmax-kshift
 MODIFICATION HISTORY:
	AUG-2000, Paul Hick (UCSD; pphick@ucsd.edu)
	    Rewrite of Andy Buffington's ricearc.for.
	    Main modifications:
	    - large intermediate byte array for storing individual bits is removed
	    - # leading zeroes used to store full values was reduced from kmax+1 to kmax
	    - use of kshift allows for slightly better compression
	8/24/2000, Kevin Nguyen
	    - Rewrite the code from fortran to idl


ice_read $SMEI/ucsd/sat/idl/util/ice_read.pro
[Previous] [Next]
 NAME:
	ice_read
 PURPOSE:
	Reads a compressed file written by ice_write into an integer array
 CATEGORY:
	Bits and bytes
 CALLING SEQUENCE:
	FUNCTION ice_read, cFile, Orig, info=info
 INPUTS:
	cFile	    scalar; type: string
		    file name of compressed file
 OUTPUTS:
	Result	    scalar; type: integer
		    0: failure (an error message is displayed)
		    1: success
	Orig	    array; type: long integer
 OPTIONAL OUTPUT PARAMETERS:
	info=info   scalar; type: string
		    trailer string extracted from end of file
 INCLUDE:
	@compile_opt.pro    ; On error, return to caller
 CALLS: ***
	IsType, ice_unpack
 SEE ALSO:
	ice_pack, ice_write
 MODIFICATION HISTORY:
	AUG-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	AUG-2000, Kevin Nguyen translated from Fortran to IDL


ice_unpack $SMEI/ucsd/sat/idl/util/ice_unpack.pro
[Previous] [Next]
 NAME:
	ice_unpack
 PURPOSE:
	Deompresses an integer array compressed using ice_pack
 CATEGORY:
	Bits and bytes
 CALLING SEQUENCE:
	FUNCTION ice_unpack, ibit, Pack, Orig, kmax=kmax, kshift=kshift, lenbit=lenbit, sign=sign, perm=perm
 INPUTS:
	ibit	    scalar; type: integer
		    # bits in Pack containing compressed data.

 OPTIONAL INPUTS:
	kmax	    scalar; type: integer; default: 8
		    cutoff for applying compression
	kshift	    scalar; type: integer; default: 0
		    minimum # bits copied for each difference
	lenbit	    scalar: type: integer; default: 16
		    # bits encoded  when the full value instead
		    of a difference is used
	sign	    scalar; type: integer
		    -1: all values in Orig are <  zero
		    +1: all values in Orig are >= zero
		     0: both positive and negative values are present
	perm	    array[0:kmax-kshift]; type: integer; default: lindgen(kmax-kshift+1)
		    permutation table
	Pack	    array; type:long integer
		    compressed data
 OUTPUTS:
	Result	    scalar; type: long integer
		    # bits in Pack used to hold the compressed data
		    as determined during the decompression. This
		    should match the input value ibit.
	Orig	    array; long integer
		    decompressed array
 INCLUDE:
	@compile_opt.pro    ; On error, return to caller
 CALLS: ***
	InitVar, ORIG, ibsetf, ibtestf
 CALLED BY:
	ice_read
 SEE ALSO:
	ice_pack
 PROCEDURE:
 >	The input values for ibit, kmax, lenbit, sign, perm must be the same as for
	the ice_pack call that created the compressed array Pack.
 >	Inverts the compression of ice_pack.
 MODIFICATION HISTORY:
	AUG-2000, Paul Hick (UCSD; pphick@ucsd.edu)


IDL_postop_linux $SMEI/pro/idl_postop_linux.pro
[Previous] [Next]
 NAME:
	IDL_postop_linux
 PURPOSE:
	Final startup chores
 CATEGORY:
	Startup
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	Reset_Colors, hide_env, set_page
 CALLED BY:
	IDL_startup, IDL_vm_startup
 PROCEDURE:
	JUL-2005, Paul Hick (UCSD/CASS)
	    Added check for env. var IDL_PRINTER


IDL_postop_vms $SMEI/pro/idl_postop_vms.pro
[Previous] [Next]
 NAME:
	IDL_postop_vms
 PURPOSE:
	Initialize miscellaneous IDL V2 facilities
 CATEGORY:
	utility
 CALLING SEQUENCE:
	invoked when IDL session is started
 CALLS: ***
	Reset_Colors, set_page
 CALLED BY:
	IDL_startup, IDL_vm_startup
 PROCEDURE:


IDL_postop_win $SMEI/pro/idl_postop_win.pro
[Previous] [Next]
 NAME:
	IDL_postop_win
 PURPOSE:
	Executed during startup after the IDL path has been established.
 CATEGORY:
	startup
 CALLING SEQUENCE:
	IDL_postop_win
 CALLS: ***
	CheckDir, FILEPATH, Reset_Colors, set_page
 CALLED BY:
	IDL_startup, IDL_vm_startup
 PROCEDURE:
	Sets defaults for foreground and background color.
	Initializes printer (with set_page).
	Moves to working directory


IDL_startup $SMEI/pro/idl_startup.pro
[Previous] [Next]
 NAME:
	IDL_startup
 PURPOSE:
	Initialize IDL session
 CATEGORY:
	Startup
 CALLING SEQUENCE:
	Invoked when IDL session is started
 INPUTS:
	(none)
 OUTPUTS:
	(none)
 CALLS:
	IDL_postop_linux, IDL_postop_vms, IDL_postop_win, defsysv_smei
 PROCEDURE:
 >	Environment variables:
	SMEI:
	    The environment variable 'SMEI' MUST be defined, and point to the top
	    of the SMEI software tree.
	sys, exe, com:
	    Will be defined by this procedure if they don't exist yet. Both are
	    subdirectories in 'SMEI' with the same name as the environment variable.
	TUB, DAT:
	    Temporary and data directory. Both are best defined externally.
	    No attempt is made to define them here

	VMS:	The startup file should be defined as a logical at LOGIN
		'$DEFINE IDL_STARTUP DSK_NAME::[DIR_NAME]IDL_STARTUP.PRO'

	LINUX:	Set startup file in the IDL development environment to this file.
		Also set the working directory to the $SMEI directory
		at the top of the software tree (where this file is located).

	WIN	Same as Linux
 MODIFICATION HISTORY:
	FEB-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IDL_startup_gen $SMEI/user/phick/pro/idl_startup_gen.pro
[Previous] [Next]
 NAME:
	IDL_startup_gen
 PURPOSE:
	Initialize IDL session
 CATEGORY:
	Startup
 CALLING SEQUENCE:
	Invoked when IDL session is started
 PROCEDURE:
	Only puts $SMEI/gen/idl in the path. Used to test integrity of /gen/idl tree.
	File needs to be put in env var IDL_STARTUP.
 MODIFICATION HISTORY:
	JAN-2002, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IDL_startup_pure $SMEI/user/phick/pro/idl_startup_pure.pro
[Previous] [Next]
 NAME:
	IDL_startup_pure
 PURPOSE:
	Initialize IDL session
 CATEGORY:
	Startup
 CALLING SEQUENCE:
	Invoked when IDL session is started
 PROCEDURE:
	Only puts $SMEI/gen/idl and $SMEI/sat/idl in the path.
	Used to test integrity of /sat/idl tree.
	File needs to be put in env var IDL_STARTUP.
 MODIFICATION HISTORY:
	JAN-2002, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IDL_startup_sat $SMEI/user/phick/pro/idl_startup_sat.pro
[Previous] [Next]
 NAME:
	IDL_startup_sat
 PURPOSE:
	Initialize IDL session
 CATEGORY:
	Startup
 CALLING SEQUENCE:
	Invoked when IDL session is started
 PROCEDURE:
	Only puts $SMEI/gen/idl and $SMEI/sat/idl in the path.
	Used to test integrity of /sat/idl tree.
	File needs to be put in env var IDL_STARTUP.
 MODIFICATION HISTORY:
	JAN-2002, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IDL_vm_startup $SMEI/pro/idl_vm_startup.pro
[Previous] [Next]
 NAME:
	IDL_vm_startup
 PURPOSE:
	Initialize IDL session
 CATEGORY:
	Startup
 CALLING SEQUENCE:
	Invoked when IDL session is started
 INPUTS:
	(none)
 OUTPUTS:
	(none)
 CALLS: ***
	CheckDir, FILEPATH, IDL_postop_linux, IDL_postop_vms, IDL_postop_win, defsysv_smei
 PROCEDURE:
 >	Environment variables:
	SMEI:
	    The environment variable 'SMEI' MUST be defined, and point to the top
	    of the SMEI software tree.
	SYS, EXE, COM:
	    Will be defined by this procedure if they don't exist yet. Both are
	    subdirectories in 'SMEI' with the same name as the environment variable.
	TUB, DAT:
	    Temporary and data directory. Both are best defined externally.
	    No attempt is made to define them here

	VMS:	The startup file should be defined as a logical at LOGIN
		'$DEFINE IDL_STARTUP DSK_NAME::[DIR_NAME]IDL_STARTUP.PRO'

	LINUX:	Set startup file in the IDL development environment to this file.
		Also set the working directory to the $SMEI directory
		at the top of the software tree (where this file is located).

	WIN	Same as Linux
 MODIFICATION HISTORY:
	FEB-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


img_read $SMEI/ucsd/gen/idl/util/img_read.pro
[Previous] [Next]
 NAME:
	img_read
 PURPOSE:
	Read image from file
 CATEGORY:
	I/O
 CALLING SEQUENCE:
	FUNCTION img_read, img_file, img_data, Red, Green, Blue,    $
	    pseudocolor     = PseudoColor   ,	$
	    truecolorimage  = TrueColorImage,	$
	    rgbplane	    = RGBPlane	    ,	$
	    truedimension   = truedimension ,	$
	    exten_no	    = exten_no	    ,	$
	    silent	    = silent	    ,	$
	    _extra	    = _extra	    ,	$
	    img_info	    = img_info	    ,	$
	    errormessage    = errormessage
 INPUTS:
	img_file    scalar; type: string
			image file to be read
			If no file is specified then the pickfile_dialog is used to
			select a file name
 OPTIONAL INPUT PARAMETERS:
	/truedimension
		    scalar; type: integer: default: 1
			dimension used for color interleave (1,2 or 3)
			(works only for JPEG?)
	/pseudocolor	(used only for .BMP and .TIFF files)
			If this keyword is set and the file stores a truecolor image then it
			is converted to pseudocolor using the IDL color_quan function.
			The color table is returned in the Red, Green and Blue arrays
 OUTPUTS:
	status	    0: some error occurred
		    1: file properly read
	img	    array[n,m] or array[3,n,m]; type: depends
			if a truecolor image [3,n,m] is returned then the keyword TrueColor
			returns 1 and RGBPlane returns the indices for the color planes
			!!! if no file name is specified, and no file is selected with
			from the file dialog then the scalar (img=-1) is returned.
 OPTIONAL OUTPUT PARAMETERS:
	red=Red, green=Green, blue=Blue
		    array[*]; type: byte (?)
			red, green and blue color arrays, defining the color table for the image
			Should be loaded with tvlct, Red, Green, Blue
			The arrays are always returned for GIF files (256 colors) and pseudocolor
			(8-bit) .TIFF and .BMP files. For truecolor .TIFF or .BMP files they are
			only returned if the /pseudocolor keyword is set.
	truecolor=TrueColor
		    scalar; type: integer
			1 if true color image is returned; otherwise 0
	rgbplane=RGBPlane
		    array[3]; type: integer; always [0,1,2] (see below)
			if a true color image is returned then RGBPlane stores the indices for
			the red, green and blue color plane, respectively.
			(!!! currently the planes are rearranged if necessary, so that
			RGBPlane =[0,1,2])
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	GetFileSpec, InitVar, IsType, MagnifyArray, READFITS, READ_BMP, READ_GIF, SubArray
	SuperArray, TimeGet, TimeUnit, bin_read, do_file, flt_read, grd_read, gunzip_file
	hide_env, smei_filename, smei_frm_cvhdr, smei_frm_read, smei_property
 CALLED BY:
	RemoteView_Display3D, even_light, put_logo, qImage_Pick, qView_FileFilter
	qView_GetData, qView_PickFiles, view
 PROCEDURE:
	The file name type is used to determine how to read the file:
	.gif: gif files are read using IDL procedure read_gif
	.tif: tiff files are read using IDL procedure read_tiff

	All remaining file types deal with grey-scale images.

	.fits and .fts: read using readfits function (from SOLARSOFT library ???)
	.txt: read using the flt_read function

	There are several types of simple binary files created using bin_read, which
	can also be read again by bin_read (the dimensions are stored in a header):
	'.pph': status = bin_read ( img_file, img )
		return short integer, long integer and floating point arrays.

	Two file types refer to CCD images:
	.img: these should be images from Andy's old Photometrics CCD camera. They are read using

	    status = bin_read(img_file,img,/aint,off=160,nx=384)
	    img = rotate(img,3)

	    i.e. a 160 byte header is skipped; then rows of 384 short integers are read
	    untill the end of file is reached, followed by clockwise rotation over 90 degree.

	'.nic': these are the images from the SMEI CCD camera. They are read using

	    status = bin_read (img_file, img)

	    bin_read returns a long integer array of 1280x600 for these files.
 MODIFICATION HISTORY:
	FEB-2000, Paul Hick (UCSD/CASS)
	    Added documentation
	OCT-2002, Paul Hick (UCSD/CASS)
	    Added keyword truedimension
	MAR-2003, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Added detection of .gz files. These are now decompressed. The decompressed
	    file is then processed in the usual, and is then deleted.


InfiniteValue $SMEI/ucsd/sat/idl/toolbox/infinitevalue.pro
[Previous] [Next]
 NAME:
	InfiniteValue
 PURPOSE:
	Provide the infinity value for any data type
 CATEGORY:
	sat/idl/toolbox
 CALLING SEQUENCE:
	FUNCTION InfiniteValue, X
 INPUTS:
	X	    variable of any type
 OUTPUTS:
	Result	    scalar; type: as implied by 'size' vector
			NaN value
 INCLUDE:
	@compile_opt.pro    ; On error, return to caller
 CALLS: ***
	IsType
 CALLED BY:
	ThomsonLOSFar, ThomsonLOSRomb
 PROCEDURE:
	For integers the number zero is returned.
	For floating point the appropriate value from !values is returned.
 MODIFICATION HISTORY:
	MAR-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


InitVar $SMEI/ucsd/gen/idl/toolbox/initvar.pro
[Previous] [Next]
 NAME:
	InitVar
 PURPOSE:
	Initialize a variable
 CATEGORY:
	gen/toolbox
 CALLING SEQUENCE:
	InitVar, X, Xinit [, /keyword_var], set=set
 INPUTS:
	X		any variable
	Xinit		any variable; used to initialize X
 OPTIONAL INPUT PARAMETERS:
	/keyword_var	if set, X is assumed to be a keyword the initialization
			then is done with the command: X = keyword_set(X)
	set=set 	instead of initializing X, the variable 'set' is initialized.
			The content of X is transferred to 'set' if X exists.
			(this keyword is ignored if /keyword_var is set).
	count=count scalar; type: integer; default: 0
			# elements in X that will trigger the initialization
			(ignored if /keyword_var is set).
 OUTPUTS:
	X		any variable; the initialized variable
 INCLUDE:
	@compile_opt.pro	    ; On error, return to caller
 CALLED BY:
	AlignHeliosphere, AngleRange, AngleUnits, ArrayLocation, CheckDir, ColorPolarBox
	ColorSkybox, CombineRotations, CvPrecess, CvRotation, CvSky, Distance2Sun
	EarthSky3DLoc, EarthTransit3DLoc, Find2DGlitch, FindAllFiles, FindAllSubDirs
	FindPeaks, FindStars, FishEye, GeographicInfo, GetColors, GetFileSpec, GroupPixels
	GuessRecs, InSitu, InsituTimeSeries, IntegrateLOS, InterpolateHeliosphere
	IsBadTime, IsDisplay, IsPrinter, IsTime, IsType, KeplerOrbit, LocalExtrema
	MagnifyArray, NewcombSun, PA_Pole, PlotCurve, PlotEarthSkymap, PlotPlanarCut
	PlotPolarSkymap, PlotSolarDisk, PlotSynopticMap, PlotUserstring, PutFileSpec
	RGBO_Project, ReadSynopticMap, RemoteView, RemoteView_BodyLoc
	RemoteView_Display2D, RemoteView_Display3D, RemoteView_FOV, RemoteView_FOV_loc
	RemoteView_FOV_xyz, RemoteView_Init, RemoteView_Init_Display
	RemoteView_Init_Matrix, RemoteView_Init_View, RemoteView_ZEclipticPlane
	RemoteView_rgbo, Reset_Colors, SetFileSpec, SetRange, SubArray, TMO_skymotion
	TagArray, ThomsonElectron, ThomsonElectronFar, ThomsonLOSFar, ThomsonLOSRomb
	ThomsonLOSStep, ThomsonMidpoint, ThomsonSolarFlux, TimeArray, TimeGet
	TimeInterpol, TimeLInterpol, TimeLimits, TimeMonth, TimeOp, TimeOrigin, TimePieces
	TimePosn_test, TimeSet, TimeShift, TimeSplit, TimeString, TimeSystem, TimeUnit
	TimeXAxis, TimeYDate, TimeYDoy, ToSolarRadii, anicenum, bargraph, big_body
	big_elongation, big_eph, big_orbit, bin_read, bin_write, boost, coord3to2, cvsmei
	cvsmei_init, do_file, eclipsed_area, edit_smei, editsmei, even_light
	even_light_corrections, even_light_figures, even_light_info
	even_light_pedestal, even_light_photometry, even_light_plot
	even_light_registration, fancy, findfile_fix, flat_centerofmass, flt_clean
	flt_read, flt_string, forecast, forecast_cfg, forecast_env, forecast_ftp
	forecast_html, forecast_ice, forecast_info, forecast_movie, forecast_movie_cp
	get_page, getnagoyasources, getootyasources, getsmeiasources, grd_read, gridfill
	gridgen, gridgen1d, gunzip_file, gzip_file, hide_env, ice_pack, ice_unpack, img_read
	is_running, jpl_body, jpl_close, jpl_eph, jpl_init, jpl_mag, jpl_parallax, jpl_phase
	jpl_sizeofsun, jpl_state, jpl_test, lsqLinearFit, lsqNormalFit, lstree, maygeometry
	mk_celias, mk_flick, mpc_comets, mpc_eph, mpc_minor_planets, mpc_orbit_eph
	nagoya_glevel, nrZbrac, nso_fe_plot, nso_fe_read, nso_fe_start, os_separator
	plot3darc, plot3dcube, plot3dline, plot3dtext, put_logo, qBar, qEphem, qEphem_State
	qImage_Send, qImage_cw, qImage_cw_Box, qImage_cw_BoxCosine, qImage_cw_BoxImage
	qImage_cw_BoxZoom, qImage_cw_Ellipse, qImage_cw_Ephem, qImage_cw_MinMax
	qImage_cw_Property, qImage_cw_Show, qImage_cw_Tool, qImage_cw_Transform
	qImage_cw_Update, qImage_cw_Wedge, qImage_cw_WedgeSection, qImage_cw_ZWedge
	qImage_cw_smei_frm, qLine, qLine_YAction, qLoadCT, qRemoteView
	qRemoteView_Calculate, qRemoteView_Histogram, qRemoteView_Time, qSave2File
	qSave2File_Pick, qShow, qTool, qTool_State [1], qTool_State [2], qTool_Zoom [1]
	qTool_Zoom [2], qTool_rebin [1], qTool_rebin [2], qView, qView_ApplyGain
	qView_GetImage, qView_Image, qView_ImageInfo, qView_PlotTrack, qView_Sensitive
	qView_SubtractBase, qView_TMO_tracksky, qView_UpdateActive, qView_XYZ_Update
	qnagoya, qset_page, qslider_cw, qsmei_hdr, qsmei_sky, qvu, scalarproduct, set_page
	sgp4_eph, sgp4_tlm, show_wso, skyd_cat, skyd_equ, skyd_version, smei_base_testcase
	smei_buf, smei_buf_get, smei_buf_getframe, smei_buf_mget, smei_buf_prep
	smei_buf_read, smei_cam2angle, smei_camera, smei_camera_gain, smei_ccd2sky
	smei_coriolis, smei_filename, smei_filepath, smei_findcrazy, smei_findpnt
	smei_frm_base, smei_frm_cp, smei_frm_cvhdr, smei_frm_drive, smei_frm_eclipse
	smei_frm_findpoint, smei_frm_flatfield, smei_frm_hbar, smei_frm_info
	smei_frm_read, smei_frm_rebin, smei_frm_smoothdark, smei_frm_summary
	smei_frm_track, smei_frm_update, smei_frm_where, smei_fts_read, smei_getfile
	smei_hdr_get, smei_hdr_make, smei_hdr_plot, smei_hdr_update, smei_mksidereal
	smei_orbit_get, smei_orbits_stat, smei_property, smei_radial2theta
	smei_setup_roi, smei_sgp4_orbits, smei_shutterwrong, smei_sky, smei_sky2cam
	smei_sky2ccd, smei_sky_cleanedge_fov, smei_sky_cleanedge_map, smei_sky_field
	smei_sky_hdr2range, smei_sky_read, smei_sky_testcase, smei_sky_track
	smei_star_corepsf, smei_star_filename, smei_star_fit, smei_star_fitone
	smei_star_formatpnt, smei_star_info, smei_star_list, smei_star_lsqfit
	smei_star_readpnt, smei_star_remove, smei_star_show, smei_star_standard
	smei_star_update, smei_star_writepnt, smei_zld_model, smei_zld_remove
	smeidb_mounted, sphere_distance, sphere_great_arc, sphere_smooth, stardistance
	statpos, stopwatch, strposn, timeposn, twin, txt_read, unexpected_event, unhide_env
	unitvectors, usno_body, usno_eph, usno_init, usno_test, vectorproduct, view, vox_read
	vox_write, vu_check, vu_correlate, vu_cvgrid, vu_earthskymap, vu_extract
	vu_filename, vu_fnc, vu_get_page, vu_getdata, vu_gettime, vu_header, vu_image
	vu_insitu, vu_insitu_raw, vu_insitucurve, vu_is_sequence, vu_lineofsight
	vu_localskymap, vu_mean, vu_movie, vu_nagoyasourcemap, vu_new_time, vu_planarcut
	vu_point_source, vu_prefix, vu_read, vu_remoteview, vu_select, vu_set, vu_solardisk
	vu_synopticmap, vu_timeseries, vu_type_insitu, vu_type_skymap, vu_update_hours
	vu_update_marker, vu_vox_drawsphere, vu_vox_read, vu_vox_write, vu_whatis
	vu_write, wedge_bounding_box, wedge_content, where_common, who_am_i, wso_read
	wso_write, www_help, www_help_called_by, www_help_change_tabs, www_help_files
	www_help_get_header, www_help_mailto, www_help_make, www_help_rsi, www_help_smei
	www_help_style, www_help_tree, xhcs, xyoff
 PROCEDURE:
	If neither X nor Xinit exist, program is terminated
 MODIFICATION HISTORY:
	JUL-2002, Paul Hick (UCSD/CASS)
	JAN-2005, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Added keyword set=set.


Inside_Wedge $SMEI/ucsd/sat/idl/toolbox/math/inside_wedge.pro
[Previous] [Next]
 NAME:
	Inside_Wedge
 PURPOSE:
	Tests whether image locations are inside a wedge
 CATEGORY:
 CALLING SEQUENCE:
	result = Inside_Wedge(p_box,p)
 INPUTS:
	p_box	    array[2,2]; type: float
			two opposite corners of the wedge in polar coordinates
			in the form [ [angle1,radius1],[angle2,radius2] ]
	p	    array[2,n,m,...]; type: float
			image locations in polar coordinates
			(same form as p_box)
 OPTIONAL INPUT PARAMETERS:
	exclude_p_box=exclude_p_box
		    array[2,2]; type: float
			two opposite corners of another wedge in polar coordinates
			in the form [ [angle1,radius1],[angle2,radius2] ]
			This second wedge must lie entirely inside the first one.
			If specified the coordinate array p is tested for inside
			p_box and outside exclude_p_box.
 OUTPUTS:
	result	    array[n,m,...]; type: byte
			0 for locations p outside the wedge
			1 for locations inside the wedge
 CALLS: ***
	IsType, SubArray
 CALLED BY:
	TMO_tracksky, cvsmei, even_light, qGlitch_Run, qImage_cw_BoxCosine
	qImage_cw_SmeiMask, qImage_cw_ZEllipse, qView_PlotSeries, smei_ccd2sky
	smei_sky2ccd, stardistance, wedge_content
 RESTRICTIONS:
	If p contains !values.f_nan the value 1 is returned
 PROCEDURE:
 >	Formerly called qImage_cw_WTest
 >	The array p can have more than two dimensions, e.g. [2,n,m].
 MODIFICATION HISTORY:
	FEB-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


InSitu $SMEI/ucsd/sat/idl/util/insitu.pro
[Previous] [Next]
 NAME:
	InSitu
 PURPOSE:
	Plots in-situ data from various spacecraft
 CATEGORY:
 CALLING SEQUENCE:
	PRO InSitu, T,	$
	    delt    = delt,	$
	    weight  = weight,	$

	    source  = source,	$
	    step    = step,	$

	    xtime   = xtime,	$
	    xlng    = xlng,	$
	    xlat    = xlat,	$

	    ylng    = ylng,	$
	    ylat    = ylat,	$
	    ydis    = ydis,	$
	    yden    = yden,	$
	    yvel    = yvel,	$
	    ybr     = ybr ,	$
	    ybt     = ybt ,	$
	    ybn     = ybn ,	$
	    ybb     = ybb ,	$
	    ykp     = ykp ,	$
	    ydst    = ydst,	$
	    ypdyn   = ypdyn,	$
	    yvlng   = yvlng,	$
	    yvlat   = yvlat,	$

	    traceback=traceback,$
	    xrange  = xrange,	$
	    charsize=charsize,	$

	    _extra  = _extra
 INPUTS:
	T	    array[2], type: standard time structure
			    time range for time series
 OPTIONAL INPUT PARAMETERS:
	/xtime, /xlng,	/xlat
			select one to set the quantity plotted on the x-axis: time, heliographic longitude
			or heliographic latitude. Only one can be specified; if none is specified
			then /xtime is assumed
	source=source
		scalar; type: integer
			identifies in situ instrument (see function Instrument)
	/ylng,	/ylat,	/ydis,	/yden,	/yvel, /ybr, /ybt, /ybn, /ybb, /ykp, /ydst, /ypdyn, /yvlng, yvlat
			detemines quantity plotted on y-axis. Multiple keywords can be specified.
			If none is specified then velocity and density is plotted.
	traceback=traceback
			scalar; type: any; default:0 (no traceback)
			heliocentric traceback distance in AU; adjusts time and heliographic longitude to
			give a 'time and location of origin' at the traceback distance
 OUTPUTS:
	Output to screen
 OPTIONAL OUTPUT PARAMETERS:
	(none)
 INCLUDE:
	@compile_opt.pro		; On error, return to caller
 CALLS: ***
	Carrington, InitVar, InsituTimeSeries, Instrument, IsType, PlotCurve, TimeSet
	TimeUnit, TimeXAxis
 PROCEDURE:
 MODIFICATION HISTORY:
	SEP-1999, Paul Hick (UCSD/CASS)
	DEC-2001, Paul Hick (UCSD/CASS)
	    added call to Carrington to allow input of time as Carrington variable
	SEP-2006, Paul Hick (UCSD/CASS)
	    Added /ypdyn keyword to process Mars Global Surveyor dynamic pressure data
	SEP-2006, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Added /yvlng and /yvlat keywords to plot velocity angles from the
	    SOMNI and EOMNI databases.


InsituTimeSeries $SMEI/ucsd/sat/idl/util/insitutimeseries.pro
[Previous] [Next]
 NAME:
	InsituTimeSeries
 PURPOSE:
	Extract in-situ data in specified interval, and average if necessary
 CALLING SEQUENCE:
	FUNCTION InsituTimeseries,$
	    data_request	, $
	    t_request		, $
	    t_given		, $
	    data_given		, $
	    delt    = delt	, $
	    weight  = weight	, $
	    traceback=traceback , $
	    source  = source_	, $
	    silent  = silent
 INPUTS:
	data_request	array, type: string
			    determines which data are extracted.
			    Identifiers are:
			    'lng','lat','dis','den','vel','br','bt','bn','bb','pdyn','vlng','vlat'
			    'lng'   = heliographic longitude (deg)
			    'lat'   = heliographic latitude (deg)
			    'dis'   = heliocentric distance (AU)
			    'den'   = density (cm^-3)
			    'vel'   = velocity (km/s)
			    'br'    = radial magnetic field (nT)
			    'bt'    = tangential magnetic field (nT)
			    'bn'    = normal magnetic field (nT)
			    'bb'    = field strength (nT)
			    'pdyn'  = dynamic pressure (nPa)
			    'vlng'  = azimuthal flow angle
			    'vlat'  = elevation flow angle
	t_request	array[2],array[N] type: standard time structure
			    If array[2]: data are extracted between t_request[0] and t_request[1]
				at times extracted from s/c insitu files
			    If array[N]: data are extracted at times t_request
			    Should be in chronological order
 OPTIONAL INPUT PARAMETERS:
	source=source	scalar; type: integer or string
			    identifies insitu instrument either by name or integer
			    as provided by function Instrument.
	delt=delt	scalar; type: time difference structure, or float; default: 0
			    if present, and non-zero, time window used for averaging.
			    in-situ data inside [t-delt/2,t+delt/2] are averaged.
			    If a scalar (float or integer) is specified, then it is
			    assumed to be the time window in days.
	/weight 	if set, standard deviations are used to weight the averages (if available
			    for selected instrument and data)
			    if not set, an unweighted average is calculated
	traceback=traceback
			scalar; type: any; default:0
			    trace back distance (AU)
			if set and nonzero then the in situ data are traced back to a distance of
			'traceback' AU, i.e. the time and heliographic longitude reflect the time of
			origin at the traceback distance.
			This option explicitly assumes that the distances in the in situ files are
			in AU, the velocities in km/s, and the heliographic longitudes in degrees.
			If this option is used date_request MUST contain 'lng','dis','vel'
 OUTPUTS:
	Result		scalar; type: integer
			0: if something went wrong
			1: time series succesfully set up
			(first check this status indicator before using the output arrays!!!)
	t_given 	array[*]; type: time structure
			    times for output time series
			    If delt is set to a positive value and t_request contains more than two
			    elements then t_given=t_request and the time series contains data
			    averaged over delt days.
			    If delt is zero, or t_request contains only 2 elements then
			    t_given contains all times for which spacecraft data were found
			    (averaged over delt days if delt > zero).
			    (if /traceback is used then the time range is shifted relative to
			    the input range)
	data_given	array[*,n_elements(data_request)], type: float
			    data_given[*,i] are the data for the quantity indicated by data_request[i]
			    missing data are marked by the value !values.f_nan
 OPTIONAL OUTPUT PARAMETERS:
	source=source	scalar; type: integer
			    if not set on input, the return value is set to whatever default
			    is set by the function Instrument (currently /somni).
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	BadValue, CvSky, FILEPATH, InitVar, Instrument, IsTime, IsType, TimeFixYear, TimeGet
	TimeLimits, TimeOp, TimeSet, TimeUnit, UNIQ, big_body, big_eph, flt_read, jpl_body
 CALLED BY:
	InSitu, vu_insitu
 RESTRICTIONS:
	Needs system variables: !sun.spiral and !sun.au
	The input time array should be in chronological order
 PROCEDURE:
 >	Input files for density and velocity are hourly-averages. The filename is:
	    fil = filepath(root=getenv('SSW_SMEI_DAT'),sub='insitu','sw'+spacecraft+'_'+strcompress(year,/rem)+'.hravg')
	    e.g. file = 'imp8_2000.hravg'
 >	The files are assumed to be chronologically ordered.
 >	The file structure is described internally as column numbers for all quantities.
 > Heliographic longitudes are returned in monotonic decreasing order by adding a multiple of 360
	degrees when the timeseries crosses from one Carrington rotation to the next.
 MODIFICATION HISTORY:
	AUG-1998, Paul Hick (UCSD/CASS)
	JUL-1999, Paul Hick (UCSD/CASS); documentation added
	SEP-2001, Paul Hick (UCSD/CASS)
	    Added GSE & GSM to RTN conversion; currently this is only used
	    for the real-time ACE data.
	JUL-2002, Paul Hick (UCSD/CASS)
	    Fixed a bug with the delt keyword. When delt was input as a
	    difference time structure it was converted to hours, instead of days.
	DEC-2003, Paul Hick (UCSD/CASS)
	    Changed /omni to /somni (accesses the files with B in RTN)
	    Added /eomni (accesses the files with B in GSM)
	    The insitu data files can now also be located in a subdirectory
	    of $SSW_SMEI_DATA/insitu, and the filenames can be more general
	    than the default template described in PROCEDURE.
	SEP-2006, Paul Hick (UCSD/CASS)
	    Added Dana Crider Mars Orbiter data
	OCT-2006, Paul Hick (UCSD/CASS)
	    All ephemeris calculations now done by big_eph.
	    Added 'vlng','vlat' to pull the off-radial flow angles of the solar
	    wind speed out of the OMNI databases. Note that ONMI2 (EOMNI) and
	    modified OMNI (SOMNI) use different coordinates systems and different
	    definitions for the flow angles.
	FEB-2007, Paul Hick (UCSD/CASS)
	    Replaced "from" keyword by "source" keyword. "from" is still accepted
	    but is now considered obsolete.
	AUG-2007, Paul Hick (UCSD/CASS)
	    Input keyword can now also be a string (instead of just an integer)
	    In addition made the order of the various spacecraft independent
	    of the order in which function Instruments sets them up.
	    Removed keyword 'from'.
	SEP-2007, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Added instrument vhm to plot Ulysses magnetic field data.
	    Added instrument windb to plot Wind magnetic field data
	    Rename instrument wind to windsw (plots n,v data)


Instrument $SMEI/ucsd/sat/idl/util/instrument.pro
[Previous] [Next]
 NAME:
	Instrument
 PURPOSE:
	Provides a means for consistently dealing with all the different instrument for which
	in situ observations are available
 CATEGORY:
	Auxilliary
 CALLING SEQUENCE:
	FUNCTION Instrument, which, $
	    imp8    = imp8	, $
	    celias  = celias	, $
	    windsw  = windsw	, $
	    windb   = windb	, $
	    acesw   = acesw	, $
	    swace   = swace	, $
	    aceb    = aceb	, $
	    somni   = somni	, $
	    eomni   = eomni	, $
	    helios1 = helios1	, $
	    helios2 = helios2	, $
	    stereoa = stereoa	, $
	    stereob = stereob	, $
	    swoops  = swoops	, $
	    vhm     = vhm	, $
	    mgs     = mgs	, $

	    name    = name	, $
	    label   = label	, $

	    nearearth = nearearth,$
	    L1	    = L1	, $
	    list    = list
	Result = Instrument( which [, /name, /label, /nearearth )
	Result = Instrument( /imp8 [, /name, /label, /nearearth )
 INPUTS:
	which	    scalar; string or integer
			integer ID of primary name of instrument
 OPTIONAL INPUT PARAMETERS:
	The following instrument keywords are available
	(instead of these keywords the argument 'which' can be set to 'imp8','celias', etc.)
	/imp8
	/celias
	/windsw
	/windb
	/acesw
	/aceb
	/somni
	/helios1
	/helios2
	/stereoa	Stereo Ahead
	/stereob	Stereo Behind
	/swoops 	Ulysses plasma data
	/vhm	    Ulysses magnetic field data
	/mgs		Mars global surveyor

	/name		retrieves primary name
	/label		retrieves alternative name (usually same as primary name)
	/nearearth	identifies instrument as near-earth or deep-space
	/L1		identifies instrument as at L1.
 OUTPUTS:
	Result		scalar; integer, byte or string
			/name set: string containing the primary name
			/label set: string containing the alternative name
			/nearearth set: 1B if instrument is a nearearth instrument;
			    0B if it is a deep-space instrument
			no keyword set: integer identifier for the instrument
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLED BY:
	InSitu, InsituTimeSeries, run_forecast_insitu, vu_insitu, vu_movie
 PROCEDURE:
	Each of the in situ instruments is associated with an integer identifier.
	Currently listed by name are
	    ['imp8','celias','windsw','somni','helios1','helios2','stereoa','stereob','swoops','acesw','aceb']
 MODIFICATION HISTORY:
	FEB-2000, Paul Hick (UCSD/CASS)
	FEB-2000, Paul Hick (UCSD/CASS)
	    Changed /omni to /somni; added /eomni
	FEB-2000, Paul Hick (UCSD/CASS)
	    Added Mars Global Surveyor
	AUG-2007, John Clover (UCSD/CASS; jclover@ucsd.edu)
	    Added STEREO A/B
	AUG-2007, Paul Hick (UCSD/CASS)
	    Added keyword list to pull out list of all instruments.
	SEP-2007, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Added instrument 'vhm'


IntegrateLOS $SMEI/ucsd/sat/idl/util/integratelos.pro
[Previous] [Next]
 NAME:
	IntegrateLOS
 PURPOSE:
	Calculate a sky map by integrating through a 3D heliospheric matrix.
	Two forms of integrated expressions are implemented:
	(1) A weighted integration (used for g-level and IPS velocity):
	    F = Integral(F*W ds)/Integral(W ds)
	(2) Straight integration (used for Thomson scattering):
	    F = Integral(F ds)
	See PROCEDURE for more details.
 CATEGORY:
	sat/idl/util
 CALLING SEQUENCE:
	FUNCTION IntegrateLOS,	R, FF, WW,	$

	    ; Keywords passed to InterpolateHeliosphere

	    R3D     = R3D	, $
	    dR3D    = dR3D	, $
	    xcrange = xcrange	, $
	    cv2carrington=cv2carrington,$
	    xlrange = xlrange	, $
	    opengrid= opengrid	, $
	    xcgrid  = xcgrid	, $
	    xlgrid  = xlgrid	, $
	    rrgrid  = rrgrid	, $
	    fillbad = fillbad	, $

	    degrees = degrees	, $
	    minRR   = minRR	, $

	    f3dfnc  = f3dfnc	, $
	    f3darg  = f3darg	, $
	    w3dfnc  = w3dfnc	, $
	    w3darg  = w3darg	, $
	    rr_earth= rr_earth	, $
	    ut_earth= ut_earth	, $
	    pa_earth= pa_earth	, $
	    elo_earth=elo_earth , $
	    elo_sun = elo_sun	, $
	    silent  = silent
 INPUTS:
	R	array[3,*,n]; type: float
		    Locations in spherical coordinates for all points along all
		    lines of sight. '*' can be absent or represent more than one dimension.
		    E.g. for a typical skymap there will be two dimensions l,m
		    representing longitude and latitude. 'n' represents the
		    'depth' dimension along the lines of sight.

		    R[0,*,n] contains longitude angles in [0,360]
		    R[1,*,n] contains latitudes in [-90,90]
		    R[2,*,n] contains heliocentric distances in AU

	The last dimension is the 'depth' dimension along the line of sight over
	which the integration is performed.

	'Densities' and 'Weights' can be specified in two ways.

	If R3D and xcgrid are NOT specified then FF and WW directly specify densities
	and weights along all lines of sight:

	FF	    array[*,n]; type: float
			volume data at locations 'R'.
	WW	    array[*,n]; type: float
			weights at locations 'R'.

	If R3D or xcgrid are specified than calls to InterpolateHeliosphere
	are made with arrays FF and WW to determine the function values and
	weight along each line of sight.

	In this case the following keywords are passed to InterpolateHeliosphere
	    R
	    FF or WW
	    R3D     = R3D
	    dR3D    = dR3D
	    xcrange = xcrange
	    cv2carrington=cv2carrington
	    xlrange = xlrange
	    opengrid= opengrid
	    xcgrid  = xcgrid
	    xlgrid  = xlgrid
	    rrgrid  = rrgrid
	    fillbad = fillbad

	See InterpolateHeliosphere for more information.
	NOTE: R[0,*] MUST BE longitude angles, even though InterpolateHeliosphere
	    in principle allows both R[0,*] and xcrange to be Carrington variables.
	    R[0,*] must be a longitude angle because R is also passed as
	    argument to f3dfnc and w3dfnc, where it is used in combination
	    with the longitude angles in rr_earth

 OPTIONAL INPUT PARAMETERS:
	/degrees    if set then all angles are in degrees (default: radians)

	The integrands are set up using a number of additional keywords:

	ut_earth=ut_earth
	rr_earth=rr_earth
		position Earth in same coordinates as R.
	pa_earth=pa_earth
		array[*,n]; type: float
		    line of sight position angle measured counterclockwise
		    from ecliptic north
	elo_earth=elo_earth
		array[*,n]; type: float
		    angle Sun-observer-segment on line of sight (i.e. the elongation)
		    for all line of sight segments (note that array[*,0] = array[*,1] etc.
	elo_sun=elo_sun
		array[*,n]; type: float
		    angle observer-Sun-segment on line of sight
 OUTPUTS:
	Result	array[*]; type: float
		    sky map with line-of-sight integrated quantities
 INCLUDE:
	@compile_opt.pro		; On error, return to caller
 CALLS: ***
	BadValue, InitVar, InterpolateHeliosphere, IsType
 CALLED BY:
	RemoteView_Display2D, vu_earthskymap, vu_lineofsight
 PROCEDURE:
 > Two types of expressions are evaluated

	(1) <F> = Integral[ w3dfnc(s,W(s))*f3dfnc(s,F(s)) ds]/ Integral[ w3dfnc(s,W(s)) ds]
	and
	(2) <F> = Integral[ w3dfnc(s,W(s))*f3dfnc(s,F(s)) ds]

	The line of sight integrations is done for all lines of sight specified.
	The integration dimension always is the last dimension in the input arrays.
	The integration is performed simply by summing over the integration dimension,
	implicitly assuming a constant stepsize ds:

	(1) <F> = Sum[ w3dfnc(s,W(s))*f3dfnc(s,F(s)) ]/ Sum[ w3dfnc(s,W(s)) ]
	and
	(2) <F> = Sum[ w3dfnc(s,W(s))*f3dfnc(s,F(s)) ]

	Note that for expression (1) the stepsize ds will cancel. For expression (2)
	it is the users responsibility to absorb the stepsize ds into one of the
	input arrays F or W.

 >	W(s) and F(s) represent two physical quantities known on a heliographic grid
	(for example solar wind velocity and density, and the s-dependence reflects some
	geometrical aspect of the integration.

	In all applications I can think of the dependence of f3dfnc on F(s) will be linear:
	f3dfnc(s,t) = f3dfnc(s)*t. In that case the integration looks like:
	    <F> = Integral[ w3dfnc(s,W(s))*f3dfnc(s)*F(s) ds ]/ Int[ w3dfnc(s,W(s)) ds ]
	This form applies to the calculation of IPS velocities, where the weight factor
	depends on the density: w3dfnc is 'IPS_velocity_w', W(s) = n(s), and
	the perpendicular velocity is integrated: f3dfnc(s) = sin(angle(s)) (is
	'IPS_velocity') and F(s) = V(s).

	In its simplest form w3dfnc(s,t) = w3dfnc(s) and f3dfnc(s,t) = t. In that case:
	    <F> = Integral[ w3dfnc(s)*F(s) ds ]/ Int[ w3dfnc(s) ds ]
	i.e. <F> is a weighted average of F(s) along the line of sight. This form applies
	to the calculation of IPS g-level with a weight factor that depends only on the
	line of sight geometry: w3dfnc is 'IPS_hlevel_w', and the quantity F(s) is
	the 'gamma function' (introduced in the tomography manifesto).

	A special case occurs when there is no denominator:
	    <F> = Integral[ w3dfnc(s)*F(s) ds ]
	This applies to the Thomson scattering brightness where w3dfnc is the scattered
	intensity from a single electron ('ThomsonBrightness'), and F(s) is the
	solar wind density n(s). The stepsize ds is absorbed into ThomsonBrightness.

 >	Bad values in the input array F3D are detected with the IDL 'finite' function.
	This is done InterpolateHeliosphere. The treatment of these bad values depends on
	the setting of 'fillbad'.

 MODIFICATION HISTORY:
	AUG-1999, Paul Hick (UCSD/CASS)
	AUG-2004, Paul Hick (UCSD/CASS)
	    Bug fix: w3dfnc was not be used if WW did not exist. This is wrong
	    when w3dfnc only uses R and not WW (as for the g-level determination).
	    Bug fix: if neither w3dfnc nor WW is defined (as for Thomson
	    scattering the average rather than the total along the line of
	    sight was calculated.
	    Also rearranged some of the arrangement. Now the presence of R3D
	    triggers a call to InterpolateHeliosphere.
	JUN-2006, Paul Hick (UCSD/CASS)
	    Fixed bug in the summation expressions. It was explicitly set to
	    the third dimension. Now it is set to the last dimension present.
	    in F and W. Added grid keywords xcgrid, xlgrid and rrgrid to call to
	    InterpolateHeliosphere.
	FEB-2007, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Added test for existence of xcgrid to trigger call to
	    InterpolateHeliosphere.


InterpolateHeliosphere $SMEI/ucsd/sat/idl/util/interpolateheliosphere.pro
[Previous] [Next]
 NAME:
	InterpolateHeliosphere
 PURPOSE:
	Get interpolated function values on scalar heliospheric matrix at specified
	locations in the heliosphere.
 CATEGORY:
	sat/idl/util
 CALLING SEQUENCE:
	FUNCTION InterpolateHeliosphere, R, F3D, R3D, dR3D,	$
	    xcrange	= xcrange   , $
	    xlrange	= xlrange   , $
	    is_carrington=is_carrington, $
	    cv2carrington=cv2carrington, $
	    opengrid	= opengrid  , $
	    periodic	= periodic  , $
	    xcgrid	= xcgrid    , $
	    xlgrid	= xlgrid    , $
	    rrgrid	= rrgrid    , $
	    fillbad	= fillbad   , $
	    degrees	= degrees

	(1) Matrix specified on implicit regular grid in spherical coordinates:
	    (usually heliographic coordinates)

	F = InterpolateHeliosphere(R, F3D, R3D, dR3D,	$
	    [xcrange=xcrange, xlrange=xlrange, /cv2carrington, /opengrid, /periodic, /fillbad, /degrees])

	(2) Matrix on explicitly specified grid in spherical coordinates:
	    (usually heliographic coordinates)

	F = InterpolateHeliosphere(R, F3D, xcgrid=xcgrid, xlgrid=xlgrid, rrgrid=rrgrid,     $
	    [/fillbad, /degrees])
 INPUTS:
	All distances (R[2,*], R3D, dR3D) must have the same units (usually AU).
	The coordinate system used for R must be consistent with the grid on which
	F3D is defined. For the tomography heliographic coordinates are used, but
	other coordinate systems (e.g. ecliptic) are allowed.

	(1) Matrix specified on regular grid in spherical coordinates:

	R	    array[3,*]; type: float
			heliospheric locations where interpolated values are needed
			R[0,*] longitudes (angle in radians or degrees, or Carrington variables)
			R[1,*] latitudes (angle in radians or degrees)
			R[2,*] heliocentric distance
	F3D	    array[L,M,N,K]; type: float
			3D matrix of function values of scalar heliospheric quantity
			1 dim: longitude angles or Carrington variables
			    The range is set by keyword xcrange. xcrange covers 360 degrees.
			    For longitude angles: l=0,L-1 covers [xcrange[0],xcrange[1]]
				where xcrange[1]-xcrange[0] is 360 degrees.
			    For Carrington variables: l=0,L-1 covers [xcrange[1],xcrange[0]]
				where xcrange[1]-xcrange[0] is 1
			2 dim: latitude angles
			    The range is set by keyword xlrange.
			    m=0,M-1 covers [xlrange[0],xlrange[1] degrees

			For both angles the grid is closed (with the outer grid points on
			the edges of the ranges. If /opengrid is set the angular
			grids are open (with the outer grid points half a grid spacing
			away from the edges of the ranges).

			3 dim: heliocentric distance: n=0,N-1 covers R3D+i*dR3D
			4 dim: absent for scalar field; 3 for the 3 components for a vector field
	R3D	    scalar; type: float
			heliocentric distance of inner boundary of F3D
	dR3D	    scalar; type: float
			heliocentric distance resolution of F3D

	See keywords xcrange, xlrange, opengrid, periodic, is_carrington and
	cv_carrington for more control over the implicit grid for F3D.

	(2) Matrix on explicitly specified grid in heliographic coordinates:

	R	    array[3,*]; type: float
			heliospheric locations where interpolated values are needed
			R[0,*] longitude angles (radians or degrees) or Carrington variables
			R[1,*] latitude angles (radians or degrees)
			R[2,*] heliocentric distance
	F3D	    array[L,M,N,K]; type: float
			3D matrix of function values of scalar or vector heliospheric quantity

	xcgrid=xcgrid
		    array[L]; type: any
			grid of longitude angles or Carrington variables, matching
			R[0,*] (covering 1st dimension of F3D)
	xlgrid=xlgrid
		    array[M]; type: any
			grid of latitude angles matching R[1,*]
			(covering 2nd dimension of F3D)
	rrgrid=rrgrid
		    array[N]; type: any
			grid of heliocentric distances matching R[2,*]
			(covering 3rd dimension of F3D)

 OPTIONAL INPUT PARAMETERS:
	/degrees    if set then the longitude and latitude are assumed to be in
			degrees (default: radians)
	/fillbad    if set then GridFill is called first to fill in all 'bad' values

	(1) Matrix specified on implicit regular grid in spherical coordinates:

	/is_carrington
		    if SET then xcrange and R[0,*] are assumed to be
		    Carrington variables
	/cv2carrington
		    if SET then xcrange is assumed to be specified as
			as Carrington variable and R[0,*] is assumed to
			be heliographic longitude.
			The longitudes are converted to Carrington variables using
			    Carrington( mean(xcrange), R[0,*], /get_variable )
			These values will be within 0.5 of mean(xcrange) if the longitudes
			are in the range [0,360]
			(the input R is NOT MODIFIED!!)

	If neither /is_carrington nor /cv2carrington is set then xcrange and R[0,*]
	are assumed to be longitude angles

	xcrange=xcrange
		    scalar, array[2]; type: float; default: [0,360]
			(or [0,1] if /is_carrington or /cv2carrington is SET).
			If xcrange is a scalar then xcrange+[0,360] is assumed.
			(or xcrange+[0,1] if /is_carrington or /cv2carrington is SET).
			Range of longitude angles or Carrington variables covered by
			array F3D.
			!!! the default [0,1] is useful only if /cv2carrington is
			    set and F3D corresponds to an integer rotation (longitude
			    range [0,360])
	xlrange=xlrange
		    scalar; array[2]; type: float; default: [-90,90]
			If xlrange is a scalar then xlrange*[-1,1] is assumed.
			Range of heliographic latitudes covered by array F3D
	/opengrid   if set then the angular grids (longitude and latitude) are assumed
			to be open. By default the grids are assumed closed.
	/periodic   maps all R[0,*] values back into the range defined by xcrange
			Only useful if F3D represents exactly one rotation.
 OUTPUTS:
	F	    array[*]; type: float
			interpolated function values
 INCLUDE:
	@compile_opt.pro	    ; On error, return to caller
 CALLS: ***
	BadValue, Carrington, InitVar, IsType, MEAN, REVERSE, ToRadians, boost, gridfill
 CALLED BY:
	IntegrateLOS, qvu_draw, vu_atlocation, vu_timeseries
 PROCEDURE:
 >	Bad values are identified with the IDL 'finite' function
 >	The location R are translated into index values. The IDL function
	interpolate is used to obtain an interpolated value. 'Missing' values
	are set to BadValue(F3D).

	(1) Matrix specified on implicit regular grid in spherical coordinates:

 >	If /is_carrington is SET then R[0,*] and xcrange MUST contain Carrington variables
 >	If /cv2carrington is SET then R[0,*] MUST contain heliographic longitudes and
	    xcrange MUST contain Carrington variables
 >	If neither are SET then R[0,*] and xcrange MUST contain longitude angles in
	    a spherical coordinate system.

 MODIFICATION HISTORY:
	OCT-1999, Paul Hick (UCSD/CASS)
	SEP-2002, Paul Hick (UCSD/CASS)
	    Added check for bad entries in positions R. If any one of the three
	    coordinates is bad the matching entry in F is set to BadValue(F3D)
	NOV-2003, Paul Hick (UCSD/CASS)
	    Fixed problem with replicate statement (didn't work in IDL 5.3)
	JUN-2006, Paul Hick (UCSD/CASS)
	    Changed order of arguments (F3D is now 2nd instead of 4th argument)
	    Added keywords xcgrid, xlgrid and rrgrid to handle matrices defined
	    on irregular grids.
	OCT-2006, Paul Hick (UCSD/CASS)
	    Modified to allow interpolation on vector fields.
	FEB-2007, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Generalized to work with spherical coordinate systems other the
	    just heliographic. Introduces keywords /is_carrington and /cv2carrington
	    to support the special case of heliographic coordinates
	    (used by the heliospheric tomography programs).


IPS_hlevel_w $SMEI/ucsd/sat/idl/toolbox/ips/ips_hlevel_w.pro
[Previous] [Next]
 NAME:
	IPS_hlevel_w
 PURPOSE:
	Calculates weights for integrating IPS g-level
 CATEGORY:
	sat/idl/toolbox
 CALLING SEQUENCE:
	FUNCTION IPS_hlevel_w, rr_earth, R, F, w3darg,	$
	    ut_earth	= ut_earth  , $
	    pa_earth	= pa_earth  , $
	    elo_earth	= elo_earth , $
	    elo_sun	= elo_sun   , $
	    degrees	= degrees
 INPUTS:
	rr_earth    array[3]		for 'plain' sky map
		    array[3,N]		for 'transit' sky map
					location of Earth for all lines of sight
					in same spherical coordinates as R.
	R	    array[3,N,L,M]	(N,L may be 1 or absent)
					locations of all segments for all lines of sight
					in same spherical coordinates as rr_earth
					(usually heliographic)
					This can be a grid of NxL lines of sight with M
					segments along each line of sight
					R[0,*,*,*] : longitude
					R[1,*,*,*] : latitude
					R[2,*,*,*] : heliocentric distance
	F				not used
	w3darg	    array[4]		w3darg[0] = observing freq (MHz)
					w3darg[1] = source size (arcsec)
					w3darg[2] = extra source size scale factor
					w3darg[3] = beta_r
 OPTIONAL INPUT PARAMETERS:
	/degrees			if set all input angles should be in degrees (default: radians)
	ut_earth			not used
	pa_earth			not used
	elo_earth			not used
	elo_sun 			not used
 OUTPUTS:
	R	    array[N,L,M]	weight factors
 INCLUDE:
	@compile_opt.pro	    ; On error, return to caller
 EXTERNAL BY:
	vu_earthskymap, vu_lineofsight
 CALLS: ***
	CV_COORD, IPS_WeightFnc, SubArray, SuperArray
 CALLED BY:
	IPS_velocity_w
 RESTRICTIONS:
	For each line of sight segment the weight is the product of a factor
	which is a function heliocentric distance and a factor which is a
	function of the geocentric distance. The second factor is calculated
	for one line of sight only on the assumption that it's the same for
	all lines of sight.
 PROCEDURE:
 MODIFICATION HISTORY:
	SEP-1999, Paul Hick (UCSD/CASS)
	FEB-2007, Paul Hick (UCSD/CASS; pphick@ucsd.edu
	    Modified to work with list of los specified in any configuration
	    (not just as a 2D skymap)


IPS_params $SMEI/ucsd/sat/idl/toolbox/ips/ips_params.pro
[Previous] [Next]
 NAME:
	IPS_params
 PURPOSE:
	Sets parameters for calculating IPS velocities and g-levels
 CATEGORY:
	Physics
 CALLING SEQUENCE:
	R = IPS_params( /nagoya )
	R = IPS_params( freqmhz=FreqMHz, arcsec=ArcSec, sourcescale=SourceScale)
 OPTIONAL INPUT PARAMETERS:
	/nagoya 	if set, returns [327 , 0.1, 1.0]
	/cambridge	if set, returns [81.5, 0.3, 1.0]
	/sandiego	if set, returns [73.8, 0.3, 1.0]
 OUTPUTS:
	R	    array[3]; type: float; default: [327, 0.1, 1.0] (Nagoya values)
			R[0]: observing frequency in MHz
			R[1]: source size (arcsec)
			R[2]: source scale (multiplicative factor to be applied to source size)
 CALLED BY:
	IPS_WeightFnc, vu_earthskymap, vu_lineofsight
 MODIFICATION HISTORY:
	JUN-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IPS_velocity $SMEI/ucsd/sat/idl/toolbox/ips/ips_velocity.pro
[Previous] [Next]
 NAME:
	IPS_velocity
 PURPOSE:
	Calculates component of solar wind velocity perpendicular to line of sight
 CATEGORY:
	sat/idl/toolbox
 CALLING SEQUENCE:
	FUNCTION IPS_velocity, rr_earth, R, F, w3darg,	$
	    ut_earth	= ut_earth  , $
	    pa_earth	= pa_earth  , $
	    elo_earth	= elo_earth , $
	    elo_sun	= elo_sun   , $
	    degrees	= degrees
 INPUTS:
	rr_earth			not used
	R				not used
	F	    array[N,L,M]	solar wind velocity matrix
	w3darg				not used
 OPTIONAL INPUT PARAMETERS:
	/degrees			if set all input angles are in degrees
	ut_earth			not used
	pa_earth			not used
	elo_earth			angle Sun-Earth-LOS segment
	elo_sun 			angle Sun-LOS segment-Earth
 OUTPUTS:
	R	    array[N,L,M]	weight factors
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 EXTERNAL BY:
	vu_earthskymap, vu_lineofsight
 CALLS: ***
	ToRadians
 MODIFICATION HISTORY:
	SEP-1999, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IPS_velocity_w $SMEI/ucsd/sat/idl/toolbox/ips/ips_velocity_w.pro
[Previous] [Next]
 NAME:
	IPS_velocity_w
 PURPOSE:
	Calculates weights for integrating IPS g-level
 CATEGORY:
	Physics: IPS
 CALLING SEQUENCE:
	FUNCTION IPS_velocity_w, rr_earth, R, F, w3darg,    $
	    ut_earth	= ut_earth  , $
	    pa_earth	= pa_earth  , $
	    elo_earth	= elo_earth , $
	    elo_sun	= elo_sun   , $
	    degrees	= degrees
 INPUTS:
	rr_earth    array[3]		for 'plain' sky map
		    array[3,N]		for 'transit' sky map
					location of Earth for all lines of sight
					in same spherical coordinates as R.
	R	    array[3,N,L,M]	(N,L may be 1 or absent)
					locations of all segments for all lines of sight
					in same spherical coordinates as rr_earth
					(usually heliographic)
					This can be a grid of NxL lines of sight with M
					segments along each line of sight
					R[0,*,*,*] : longitude
					R[1,*,*,*] : latitude
					R[2,*,*,*] : heliocentric distance
	F	    array[N,L,M]	(N,L may be 1 or absent)
					g-values at all los segments
	w3darg	    array[4]		w3darg[0] = observing freq (MHz)
					w3darg[1] = source size (arcsec)
					w3darg[2] = extra source size scale factor
					w3darg[3] = betar
 OPTIONAL INPUT PARAMETERS:
	/degrees			if set all input angles should be in degrees (default: radians)
	ut_earth			not used
	pa_earth			not used
	elo_earth			not used
	elo_sun 			not used
 OUTPUTS:
	R	    array[N,L,M]	weight factors
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 EXTERNAL BY:
	vu_earthskymap, vu_lineofsight
 CALLS: ***
	IPS_hlevel_w
 RESTRICTIONS:
	For each line of sight segment the weight is the product of a
	factor which is a function heliocentric distance and a factor
	which is a function of the geocentric distance. The second
	factor is calculated for one line of sight only on the assumption
	that it's the same for all lines of sight.
	These weights are calculated by ips_hlevel_w, and are here
	multiplied by the F array.
 MODIFICATION HISTORY:
	SEP-1999, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IPS_WeightFnc $SMEI/ucsd/sat/idl/toolbox/ips/ips_weightfnc.pro
[Previous] [Next]
 NAME:
	IPS_WeightFnc
 PURPOSE:
	Calculates weight function along the line of sight for IPS
 CATEGORY:
	Physics
 CALLING SEQUENCE:
	W = IPS_WeightFnc(ZAu, params)
 INPUTS:
	ZAu	    array[N]; type float
			distance from Earth along line of sight (AU)
	params	    array[3]; type: float
			params[0]: observing freq in MHz
			params[1]: source size (arcsec)
			params[2]: extra scale factor for source size
 OUTPUTS:
	W	    array[N]; type: float
			weight factors
 CALLS: ***
	IPS_params, SPECTRALFNC, nrSimpson
 CALLED BY:
	IPS_hlevel_w
 PROCEDURE:
	CALCULATE INTEGRAL FOR WEIGHTING FUNCTION OF Z.
	SIMPSON METHOD DOUBLE INTEGRAL
	SpectralFnc = q*fresnel(q)*size(q)*spectr(q)
	fresnel= fresnel filter term
	size   = source size term
	spectr = spectrum fluctuation term
 MODIFICATION HISTORY:
	MAR-1999, Paul Hick (pphick@ucsd.edu)
	    Adapted from STELab Fortran code


is_running $SMEI/ucsd/gen/idl/environment/is_running.pro
[Previous] [Next]
 NAME:
	is_running
 PURPOSE:
	Check whether specified process is running
 CATEGORY:
	gen/idl/environment
 CALLING SEQUENCE:
	FUNCTION is_running, cmd, silent=silent, count=count
 INPUTS:
	cmd
 OPTIONAL INPUT PARAMETERS:
	/silent     suprresses informational messages
 OUTPUTS:
	Result	    0: 'cmd' NOT running
		    1: 'cmd' IS running
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	InitVar
 CALLED BY:
	run_map, smei_hdr_update
 PROCEDURE:
	Preliminary version. Works only on Unix/Linux
	Check 'cmd' against output of
	    ps -eo args
 MODIFICATION HISTORY:
	OCT-2006, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IsBadTime $SMEI/ucsd/gen/idl/toolbox/time/isbadtime.pro
[Previous] [Next]
 NAME:
	IsBadTime
 PURPOSE:
	Defines and checks for bad times
 CATEGORY:
	Tell time
 CALLING SEQUENCE:
	IsBad = IsBadTime(T)
 INPUTS:
	T	array[*]; type: time structure
		    time origin
 OUTPUTS:
	IsBad	array[*]
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	InitVar, TimeArray, TimeOp, TimeUnit
 CALLED BY:
	vu_new_time, vu_select, vu_update_hours
 PROCEDURE:

 MODIFICATION HISTORY:
	OCT-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IsDisplay $SMEI/ucsd/gen/idl/environment/isdisplay.pro
[Previous] [Next]
 NAME:
	IsDisplay
 PURPOSE:
	Check whether !d.name is the terminal display or the printer
 CATEGORY:
	Environment
 CALLING SEQUENCE:
	FUNCTION IsDisplay, dev
 INPUTS:
	(none)
 OUTPUTS:
	Result	    scalar; type: byte
			0: it's a printer device
			1: it's the display
 INCLUDE:
	@compile_opt.pro    ; On error, return to caller
 CALLS: ***
	InitVar
 CALLED BY:
	even_light_corrections, even_light_photometry, even_light_registration
	get_page, qBar, qLine_Curve, qnagoya_skymap, twin, view, vu_image
 PROCEDURE:
	Returns true if !d.name is !TheTerminal. !TheTerminal is set during IDL startup
 MODIFICATION HISTORY:
	FEB-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IsPrinter $SMEI/ucsd/gen/idl/environment/isprinter.pro
[Previous] [Next]
 NAME:
	IsPrinter
 PURPOSE:
	Check whether !d.name is the printer
 CATEGORY:
	Environment
 CALLING SEQUENCE:
	FUNCTION IsPrinter, dev
 INPUTS:
	(none)
 OUTPUTS:
	Result	    scalar; type: byte
			1: it's a printer device
			1: it's not the printer
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	InitVar
 CALLED BY:
	get_page
 PROCEDURE:
	Returns true if !d.name is !ThePrinter or !d.name = 'PS'.
	!TheTerminal is set during IDL startup
 MODIFICATION HISTORY:
	FEB-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IsSymLink $SMEI/ucsd/gen/idl/toolbox/files/issymlink.pro
[Previous] [Next]
 NAME:
	IsSymLink
 PURPOSE:
	Checks whether a file or directory is a symbolic link
 CATEGORY:
	Environment
 CALLING SEQUENCE:
	B = IsSymLink(file)
 INPUTS:
	file	    scalar; type: string
			name of file or directory to be tested
 OUTPUTS:
	B	    scalar; type: integer
			0: not a symbolic link
			1: it's symbolic link
		    On non-Unix systems, always 0 is returned
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLED BY:
	FindAllSubDirs
 RESTRICTIONS:
 >	Has only been tested on Linux (bash shell) under IDL 5.3 and 5.5
 >	The islink script is needed for IDL version <= 5.4 and must be in the Linux path
 >	Does not work on a path that contains a symbolic link, e.g.
	if /usr/local/rsi is a symbolic link then /usr/local/rsi/lib will
	return zero.
 PROCEDURE:
 >	Only useful on Unix/Linux systems.
	If !version.os_family is not 'unix' then always 0 is returned.

 >	If !version.os_family is 'unix' then:
	- for IDL 5.3 and earlier the script islink is spawned. The
	    output from this script (0 or 1) is captured in the 'result'
	    argument of spawn.
	- for IDL 5.4 and later, a 'test -L' is spawned, and the
	    exit_status keyword is used directly (this makes the test
	    independent of the islink script).
 MODIFICATION HISTORY:
	JULY-2001, Paul Hick (UCSD/CASS; pphick@ucsd.edu)


IsTime $SMEI/ucsd/gen/idl/toolbox/time/istime.pro
[Previous] [Next]
 NAME:
	IsTime
 PURPOSE:
	Check whether argument is a time structure
 CATEGORY:
	smei/gen/idl/toolbox/time
 CALLING SEQUENCE:
	FUNCTION IsTime, T, orbit=orbit
 INPUTS:
	T	array; type: anything
 OPTIONAL INPUTS:
	/orbit
 OUTPUTS:
	Result	scalar; type: byte
		    0 if argument is not a time structure
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	InitVar, IsType
 CALLED BY:
	Carrington, CvSky_RTN, GeographicInfo, InsituTimeSeries, PlotCurve, TimeDay
	TimeEaster, TimeXAxis, TimeYDate, TimeYDoy, cvsmei, eclipsed_area, forecast_cfg
	forecast_ice, jpl_eph, jpl_parallax, jpl_sizeofsun, lsqLinearFit, mpc_eph
	mpc_orbit_eph, nagoya_glevel, qImage_cw_Where, qLine, qLine_Curve, qLine_FitPlot
	qLine_XBase, qRemoteView_ChangeTimes, qRemoteView_Time, qView_ApplyGain
	qView_PlotSeries, qView_PlotTrack, qView_UpdateTime, qsmei_hdr, qsmei_sky_pick
	smei_buf_get, smei_buf_getframe, smei_buf_read, smei_cam_quaternion
	smei_coriolis, smei_filename, smei_filepath, smei_frm_cp, smei_frm_hbar
	smei_frm_mask, smei_frm_summary, smei_frm_where, smei_getfile, smei_hdr_make
	smei_star_fit, smeidb_mounted, usno_eph, vu_gettime, vu_insitu, vu_insitu_raw
	vu_insitucurve, vu_localskymap, vu_movie, vu_select, vu_set_time_entry
 PROCEDURE:
	Keys on structure name
 MODIFICATION HISTORY:
	JAN-2004, Paul Hick (UCSD/CASS)
	FEB-2008, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
	    Added keyword /orbit


IsType $SMEI/ucsd/gen/idl/toolbox/istype.pro
[Previous] [Next]
 NAME:
	IsType
 PURPOSE:
	Check for type of variable, or get type code
 CATEGORY:
	Toolbox: generic
 CALLING SEQUENCE:
	R = IsType(X, /floating_x)	    tests whether X is of type 'floating'
	R = IsType(/floating_x) 	    returns type code for type 'floating'
	R = IsType(X)			    returns type code for X
 INPUTS:
	X	any variable
 OPTIONAL INPUT PARAMETERS:
	Only one of these should be set:
	/Byte_X
	/Integer_X
	/Short_Integer_X	(same as Integer_X)
	/Longword_X
	/Floating_X
	/Double_X
	/Complex_floating_X
	/String_X
	/Structure_X
	/Complex_double_X
	/Pointer_X
	/Object_reference_X
	/Unsigned_Integer_X
	/Unsigned_Short_Integer_X   (same as Unsigned_Integer_X
	/Unsigned_Longword_X
	/Integer_64_bit_X
	/Unsigned_Integer_64_bit_X
 OUTPUTS:
	IsTime		scalar; type: long integer
			if argument X is set, in combination with one of the type keywords:
			    0 if argument is not of specified type; 1 if it is
			otherwise
			    type code of X, or type code for specified type keyword
 OPTIONAL OUTPUT PARAMETERS:
	type =type	scalar; type: integer
			    code of specified type
	bytes=bytes	scalar; type: integer
			    # bytes in scalar of specified type
	name =name	scalar; type: string
			    name of specified type
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	InitVar
 CALLED BY:
	AngleRange, AngleUnits, BadValue, Carrington, ColorSkybox, CombineRotations
	CvPointOnLos, CvRotation, CvSky, EarthTransit3DLoc, FindAllFiles, FindPeaks
	FindStars, FishEye, FitKeplerOrbit, GeographicInfo, GetColors, GetFileSpec
	GroupPixels, HammerAitoff, InSitu, InfiniteValue, Inside_Wedge, InsituTimeSeries
	IntegrateLOS, InterpolateHeliosphere, IsTime, MagnifyArray, MercatorProj
	PlotCurve, PlotEarthSkymap, PlotPlanarCut, PlotPolarSkymap, PlotSolarDisk
	PlotSynopticMap, PlotUserstring, RGBO_DepthCue, RGBO_Project, RemoteView
	RemoteView_CubeFrame, RemoteView_CurrentSheet, RemoteView_Display2D
	RemoteView_Display3D, RemoteView_FOV, RemoteView_FOV_xyz, RemoteView_Init
	RemoteView_Init_Matrix, RemoteView_Init_View, RemoteView_Matte
	RemoteView_rgbo, SetFileSpec, SubArray, SuperArray, TagArray, ThomsonElectron
	ThomsonLOSDensity, ThomsonLOSFar, ThomsonLOSRomb, ThomsonLOSStep
	ThomsonRadialFilter, ThomsonSetupLOS, TimeArray, TimeGet, TimeInterpol
	TimeLInterpol, TimeLimits, TimeMonth, TimeOp, TimeOrigin, TimePieces, TimeSet
	TimeShift, TimeSplit, TimeString, TimeSystem, TimeUnit, TimeXAxis, TimeYDate, TimeYDoy
	WhatIs, WhatIs0, anicenum, arrow3d, bargraph, big_body, big_elongation, big_eph
	big_eph_clean, big_eph_short, bin_read, bin_write, boost, cvsmei, cvsmei_init, do_file
	eclipsed_area, even_light_info, even_light_photometry, flat_centerofmass
	flt_read, forecast, forecast_html, forecast_movie, get_page, getnagoyasources
	getootyasources, getsmeiasources, grd_read, gridgen, gridgen1d, hide_env, ice_read
	img_read, jpl_body, jpl_eph, jpl_mag, lsqLinearFit, lsqQuadraticFit, maygeometry
	mk_celias, mk_flick, mpc_body, mpc_comets, mpc_eph, mpc_minor_planets, mpc_orbit_eph
	nagoya_glevel, nso_fe_read, plot3darc, plot3dcube, plot3dline, put_logo, qBar, qEphem
	qEphem_State, qImage, qImage_FileInfo, qImage_Pick, qImage_Send
	qImage_SendDestination, qImage_TriggerSend, qImage_cw, qImage_cw_BoxCosine
	qImage_cw_BoxImage, qImage_cw_DrawCross, qImage_cw_DrawEphem
	qImage_cw_Ellipse, qImage_cw_Ephem, qImage_cw_EphemKill [1], qImage_cw_MinMax
	qImage_cw_Property, qImage_cw_Set_Value, qImage_cw_Transform, qImage_cw_Update
	qImage_cw_Wedge, qImage_cw_ZEllipse, qImage_cw_ZWedge, qImage_cw_ctable
	qImage_cw_smei_frm, qLine, qLine_XBase, qLoadCT, qRemoteView_Kill
	qRemoteView_List, qRemoteView_Time, qSave2File, qSave2File_Save, qShow, qTool
	qTool_State [1], qTool_State [2], qView, qView_ApplyGain,