PROCEDURE:   telem_model
PURPOSE:
  Reads the Model Counts from Richard Schwartz' Model of the HESSI reponse=to
  the March 1989 Flare period (March 1-28)
  Outputs plots to the screen and IDL.PS
USAGE:
  model
INPUTS:
 Richards Mode Counts, etc.:
   ~/hessi/data/hessi_telem_sim.sav
   ~/hessi/data/counts_cutoff.sav
KEYWORDS:
CREATED BY:	David Curtis  97-11-10
 EXPLANATION



1.  If the Fast Rate > 50Kc/s/d, then increase the shutter setting

2.  If the Fast Rate << 10Kc/s/d, then decrease the shutter setting

3.  If Memory is > 25% Full, Shutter minimum setting = first in

4.  If Memory is > 50% Full, Shutter minimum setting = both in

5.  If memory is > 75% Full, Decimate Front detector events << 40KeV by 8x

6.  If Memory is > 50% Full and HESSI is in the dark, decimate Rear Background

    (I don't have a model of the background spectra as a function of time, but

     I assume I can get a 50% reduction in background event rate by

     decimation)

7.  If memory is full, stop taking data.


Memory size was set to 16Gbits, and background rate set to 2500c/s

With these rules, over the 28 days of Richards file, I am:

  10% both shutters open

  50% shutter 1 closed

  40% shutter 2 closed

  15% decimating fronts

  50% decimating rear (in the dark)

  Memory Full for about 3 hours on one occasion.


Instead of reaching Memory Full, we will probably go to strong decimation by spin of the Front Segment 
data at about 90% full to avoid the loss.


With an 8Gbyte memory (disk) The numbers are:

  10% both shutters open

  60% shutter 1 closed

  30% shutter 2 closed

  2% decimating fronts

  30% decimating rear (in the dark)

  Memory never exceeds about 85% full.


Shutter 1 switches state about 30 times, and shutter 2 about 50 times in these 28 days
 (with significant hysteresis added to minimize switching when near the memory limits on the shutter).




VERSION:  1.
LAST MODIFICATION:  09-11-10

 PROJECT:
	HESSI

 NAME: 
       Telemetry To FITS

 PURPOSE: 
       Converts Telemetry Binary Source Packet files into
       FITS binary tables for packet data and lookup table.

 CATEGORY:
       Utilities (hessi/util)

 CALLING SEQUENCE:
       HSI_Telem2Fits, telemFile [, fitsFile]

 CALLS: ***
	BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], DOC_MENU
	HSI_LOOKUPTABLEINSERTENTRY, HSI_LOOKUPTABLEWRITE, HSI_PACKETFILECLOSE
	HSI_PACKETFILEOPEN, HSI_PACKETWRITE, HSI_STRUCT_DEF [1], HSI_STRUCT_DEF [2]
	HSI_TELEM2FITS, break_file [4]
 INPUTS: 
       telemFile:  the name of the file containing the 
		    telemetry source packets

 OPTIONAL INPUT:
       fitsFile: the name of the FITS file. If not set, the name of
                 the fits file will be the same as the name of the
                 telemetry file, but with extension '.fits' 

 KEYWORDS:
    FORCEREMOVE: If set, the fits file to create is removed without
                 prior asking for confirmatin. It has no effect with
                 the keyword APPEND.
    VERBOSE: Comments the progress of the procedure

 SIDE EFFECTS:
	Creates a FITS file with the same basename as the input file.
       The files contain two binary tables, one for the packet data
       and one for the lookup table (further information in http://xxx )

 RESTRICTIONS:
	none

 PROCEDURE:
	none

 MODIFICATION HISTORY:
       Release 1 dev, csillag@ssl.berkeley.edu, Dec 98
	Release 0, A.Csillaghy, csillag@ssl.berkeley.edu, Oct 98

first correct for att_state = 4, to do this you locate the
att_state=1 nearest to the offending 4 state, and reset the count
rate and the att_state, jmm, 19-feb-2003
rates=(nchannel, ntimes)
Added atten_state =2 to calculation
Pro temp_atten421, rates, atten_state
 CALLED BY
	hsi_obs_summ_flag__define, hsi_obs_summ_rate__define

 This is different than temp_rear_decimation_correct it has the
 offset set to zero on the front end, 1 on the back..., and npts=1
 This is now used for temp_rear_decimation, for post FSW 2.4, since
 there are now rear decimation flags, jmm, 2003-03-26
 Added check and correction for isolated 0's which should never
 happen, but there is one early on 24-mar-2003
 Added check for data gaps, to correct 21-apr-2002 problem
 changed to avoid interpolating over datagaps, 29-nov-2003, jmm
 Added lt1 keyword, 7-sep-2004, jmm
 CALLS:
 CALLED BY
	hsi_obs_summ_flag__define, hsi_obs_summ_rate__define
	test_front_decimation_correct, test_rear_decimation_correct

 PROJECT:
       HESSI
 NAME:
       TEMP_BCK_TIME_INTV()
 PURPOSE: 
       Looks 1 orbit ahead, 1 orbit back, 15 ahead, 15 back, to get
       appropriate time intervals
 CALLING SEQUENCE: 
       tr = temp_bck_time_intv(time_range)
 INPUTS:
       time_range = the input time range
 OUTPUTS:
       tr=a 2X4 array of output times
 KEYWORDS:
       delta= absolute distance in km between the ephemeris at the
       midpoint of the input time range and the ephemeris at the
       midpoint of the background time range
       n_orbits = just take the background at this many orbits away
 CALLS: ***
	HSI_TIME_N_ORBITS, YOU_REALLY_WANT_A_TBEEP, anytim [1], anytim [2], anytim [3]
	anytim [4], anytim [5], yesnox [1], yesnox [2]
 HISTORY:
       Version 1, September 19, 2002, 
           jmm, jimm@ssl.berkeley.edu

  Returns an estimated upper limit for the count rate in the obs_summ
  energy channels, for a given particle rate, p. This may be of use
  in flare-listing. Only returns results for channel 2 (12-25 keV)
  and above. 20-nov-2003, jmm
  Added minp keyword, 29-nov-2003,jmm, LE minp the rate is set to
  zero, the default is 2.0
 CALLS:
 CALLED BY
	hsi_flare_list_fill [1], hsi_flare_list_fill [2]

	Pro twod_subs,ij,ipix,jpix,i,j
 Returns (i,j) the two-d subscripts of position ij in an (ipix,jpix) matrix

NAME:
 temp_random_fn
PURPOSE:
 returns random deviates with an arbitrary input probability
 distribution
CALLING SEQUENCE:
 yy = temp_random_fn(seed, n, y, fy, dy, n_pts = n_pts)
INPUT:
 seed= a seed variable, either passed into RANDOMU or generated by it
 n= the number of random numbers desired
 y= the independent variable, for which you want the random #'s
 fy = a function of y, n_elements(y), used as the probability distribution,
      but it does not need to be normalized, that is done inside the program
 dy = the grid size for y's either a scalar or a vector of n_elements(y)
KEYWORDS:
 n_pts = the number of points used for the interpolation of the inverse
         of the probability integral, the default is 1001
HISTORY:
 jmm, 14-oct-1998, from the method described in Numerical Recipes, ch 7.1
 jmm, 2-sep-1999, restricted to output only in the input range

NAME:
 temp_random_fn_a
PURPOSE:
 returns random deviates with an arbitrary input probability
 distribution
CALLING SEQUENCE:
 yy = temp_random_fn_a(seed, n, y, fy, dy, n_pts = n_pts)
INPUT:
 seed= a seed variable, either passed into RANDOMU or generated by it
 n= the number of random numbers desired
 y= the independent variable, for which you want the random #'s
 fy = a function of y, n_elements(y), used as the probability distribution,
      but it does not need to be normalized, that is done inside the program
 dy = the grid size for y's either a scalar or a vector of n_elements(y)
KEYWORDS:
 n_pts = the number of points used for the interpolation of the inverse
         of the probability integral, the default is 1001
 CALLS: ***
	INTERPOL, hsi_sim_random_2d, interp2d [1], interp2d [2], interp2d [3]
 CALLED BY:
	hsi_sim_random_2d
HISTORY:
 jmm, 14-oct-1998, from the method described in Numerical Recipes, ch 7.1
 jmm, 2-sep-1999, restricted to output only in the input range
 ras, 15-feb-2005, fixed problem with windows where normalized cumulative
                    distrbution ty sometimes does not reach 1.0
 jmm, 13-apr-2005, fixed bug which would cause non-uniform Y grids to
                   be handles incorrectly

NAME:
 temp_redo_catalog
PURPOSE:
 Redoes the Hessi qlook catalog for a file, deletes and rewrites the file
CALLING SEQUENCE:
 temp_redo_catalog, filex_fits, ignore_orbit=ignore_orbit, simulated_data = simultaed_data
INPUT:
 filex_fits = the file name, full path please
OUTPUT:
 None, the file is rewritten
KEYWORDS:
 ignore_Orbit = if set, ignore the orbital information
 simulated_data = set this if the data is simulated,
 obs_time_interval = the obs_time_interval for the file
 CALLS: ***
	DATATYPE [1], DATATYPE [2], DATATYPE [3], DELVARX [1], DELVARX [2], DELVARX [3]
	DELVARX [4], delvarx [5]
HISTORY:
 19-June-2000, jmm, jimm@ssl.berkeley.edu

Pro temp_split_flare, rate, st_ss, en_ss, split_range = split_range, $
                      _extra = _extra
 Splits a solar flare up into individual peaks and valleys, using
 local_max and local_min, with range=split_range. This is designed to
 be run in hsi_excess_intervals, or in hsi_flare_list_fill
 CALLS:
 CALLED BY
	Hsi_excess_intervals [1], Hsi_excess_intervals [2], local_max, local_min

 Returns the start and end indices of intervals where a condition
 applies
 CALLED BY
	HSI_DO_QLOOK_PLOTS, HSI_QLOOK_ASPECT, HSI_QLOOK_SUMMARY_PAGE__DEFINE
	HSI_VCFILE_PACKET_RATE, Hsi_excess_intervals [1], Hsi_excess_intervals [2]
	hsi_flare_list_fill [1], hsi_flare_list_fill [2], hsi_test_packet_rate
	old_rear_decimation_correct, temp_atten_correct

 Attempts to correct countrate for front_decimation
 paases _extra through to temp_atten_correct, 7-sep-2004
 CALLS:
 CALLED BY
	hsi_obs_summ_flag__define, hsi_obs_summ_rate__define

 Attempts to correct countrate for rear_decimation
 Fixed to deal with 1d rate, 25-feb-2004, jmm
 Passes _extra to temp_atten_correct, for /lt1 case, 7-sep-2004, jmm
Function test_rear_decimation_correct, rate, wt, ch, tim_arr, $
                                       ebands = ebands, $
                                       allchan = allchan, $
                                       _extra = _extra
 CALLS:
 CALLED BY
	hsi_obs_summ_flag__define, hsi_obs_summ_rate__define

 PROJECT:
       HESSI

 NAME:
       TIME_AXIS__DEFINE

 PURPOSE:
		This defines a generic time axis that can be specified in several ways.

 CATEGORY:
		Utilities

 CALLING SEQUENCE:
       o = time_axis()

 METHODS:
 		o->set, parameter = value: sets a parameter from the list below.
       var = o->get( /parameter ): gets the value of a parameter from the list below.

 PARAMETERS:
		time_range: the time range, in anytim format. A 2-element array.
		time_resolution the resolution of the sub_intervals between time_range

 EXAMPLES:
		see time_axis_test below

 CALLS: ***
	EXIST, TIME_AXIS::GET, TIME_AXIS::SET, TIME_AXIS_TEST, anytim [1], anytim [2]
	anytim [3], anytim [4], anytim [5]
 SEE ALSO:
	binning__define
 HISTORY:
       apr-2004 --- acs further development for the new hsi_image object
       dec-2003 --- acs, separate display parameters from axis def. parameters
       Version 1, August 21, 2003,
           A Csillaghy, csillag@ssl.berkeley.edu