[Previous]
[Next]
NAME:
obs_coord
PURPOSE:
To print out the OR table coordinate for a region clicked on
a FFI image
CALLING SEQUENCE:
obs_coor, data [, sum=sum ,side=side, nomark=nomark, $
suncenter=suncenter, date=date. bin=bin, mktv=mktv]
obs_coor, data {simplest case for HR FFI}
obs_coor, data, sum=4, side=128, date='2-mar-92', bin=2, /mktv
INPUT:
data - 2-D data image array
sum - summation mode of data (FR=1, HR=2, QR=4)
side - size of square in FR pixels to draw on image (if not
passed the default is 64)
KEYWORDS (INPUT):
suncenter = [center_x, center_y]. Sun center position
(in Yohkoh coordinates).
If not supplied, the default [123,145] will be used.
date = string 'dd-mmm-yy', e.g., '23-feb-92'.
If not supplied, today's date will be used.
bin = rebin factor of tvscl image.
mktv - If present, data is tvscled with rebin factor bin.
nomark - If present, the location on the image which
is selected is NOT marked (default is to mark it)
CALLS: ***
GBO_HEL, draw_boxcorn [1], draw_boxcorn [2]
HISTORY:
Written 20-Nov-91 by M.Morrison
Derived 2-Apr-92 by L.Acton
Modified to accept rebinned images of any resolution. LWA 6-Apr-92
Added keyword 'suncenter' to call for gbo_hel. LWA 17-Apr-92
Corrected error in use of keyword SUM. LWA 21-Apr-92
Changed default suncenter to [123,145]. LWA 9-Dec-92
[Previous]
[Next]
NAME:
obs_coord2
PURPOSE:
To print out the OR table coordinate for a region clicked on
a HR FFI image
CALLING SEQUENCE:
obs_coor, data [, sum=sum ,side=side, nomark=nomark, $
suncenter=suncenter, date=date. bin=bin, mktv=mktv]
obs_coor, data {simplest case for HR FFI}
obs_coor, data, sum=4, side=128, date='2-mar-92', bin=2, /mktv
INPUT:
data - 2-D data image array
sum - summation mode of data (FR=1, HR=2, QR=4)
side - size of square in FR pixels to draw on image (if not
passed the default is 64)
KEYWORDS (INPUT):
suncenter = [center_x, center_y]. Sun center position
(in Yohkoh coordinates).
If not supplied, the default [140,119] will be used.
date = string 'dd-mmm-yy', e.g., '23-feb-92'.
If not supplied, today's date will be used.
bin = rebin factor of tvscl image.
mktv - If present, data is tvscled with rebin factor bin.
nomark - If present, the location on the image which
is selected is NOT marked (default is to mark it)
CALLS: ***
G_H, draw_boxcorn [1], draw_boxcorn [2]
HISTORY:
Written 20-Nov-91 by M.Morrison
Derived 2-Apr-92 by L.Acton
Modified to accept rebinned images of any resolution. LWA 6-Apr-92
Derived from obs_coord to give lwa control of suncenter again!
18-Apr-92
Updated suncenter default and changed order of NS,EW input.
[Previous]
[Next]
NAME:
OBS_DEFAULTS
PURPOSE:
Loads default values into common and optionally widgets for WOBS_PLOT
CATEGORY:
CALLING SEQUENCE:
obs_defaults,bcs=bcs,hxt=hxt,sxt=sxt,wbs=wbs,goes=goes
obs_defaults,bcs=bcs,hxt=hxt,sxt=sxt,wbs=wbs,goes=goes,/widget
INPUTS:
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
widget if set, loads widgets...
bcs see obs_plot...
hxt
sxt
wbs
goes
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLED BY:
OBS_EVENT, OBS_PLOT, OBS_WIDGET
COMMON BLOCKS:
obs_tbl, obs_cm2 and obs_btn
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB Feb-95 Written
RDB 10-May-95 Added SFM option
RDB 13-Nov-97 Some simplification of GOES stuff
[Previous]
[Next]
NAME:
OBS_EVENT
PURPOSE:
Event processor called by widget created by OBS_WIDGET
CATEGORY:
CALLING SEQUENCE:
INPUTS:
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLS: ***
ADDTIME [1], ADDTIME [2], OBS_DEFAULTS, OBS_PLOT, fmt_tim [1], fmt_tim [2]
pprint [1], pprint [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
Responds to interrupts and sets values. Plot button calls OBS_PLOT
MODIFICATION HISTORY:
RDB Feb-95 Written
RDB 24-Apr-95 Added hardcopy code
RDB 10-May-95 Added SFM option
RDB 13-Nov-97 New GOES spacecraft added - droplist used
RBM added to WBS; error messages passed back to widget
now has start and stop time.
[Previous]
[Next]
NAME:
OBS_PLOT
PURPOSE:
Plots an orbit of data given file ID.
CATEGORY:
CALLING SEQUENCE:
obs_plot,tstart
obs_plot,tstart,tstop
obs_plot,tstart,dur=40,bcs=[4,3]
obs_plot,tstart,bcs=3,/hxt,/sxt
obs_plot,tstart,/hxt,sxt=[2,4]
obs_plot,'16-dec-91 04:40',dur=30,/hxt,sxt=[2,4],/goes,/line
INPUTS:
tstart Start time as a date/time string (e.g. '1-dec-92 19:45')
or in File ID format (e.g. '921201.1945')
OPTIONAL INPUT PARAMETERS:
tstop Stop time as a date/time string (e.g. '1-dec-92 20:00')
or in File ID format (e.g. '921201.2000')
- if present, tstop overrides duration
KEYWORD PARAMETERS:
bcs BCS channel number in range 1-4:
1 = Fe XXVI, 2 = Fe XXV, 3 = Ca XIX, 4 = S XV
Array of channels allowed. (/bcs == Fe XXVI)
hxt HXT channel number in range 1-4:
1 = Sum L, 2 = Sum M1, 3 = Sum M2, 4 = Sum H
Only single channel allowed. (/hxt == Sum L)
wbs WBS channel number in range 1-7:
1 = SXS1, 2 = SXS2, 3 = HXS, 4 = GRS1, 5 = GRS2
6 = RBMSC, 7 = RBMSD
Only single channel allowed. (/wbs == SXS1)
sxt SXT channel number in the range 1-7
1 = All Filters, 2 = Thin Al, 3 = Al/Mg/Mn
4 = Be, 5 = Thick Al, 6 = Mg, 7 = White Light
Array of filters can be selected. (/sxt == All Filters)
goes GOES channel number (as listed below)
60 = GOES6 Lo, 61 = GOES6 Hi, 70 = GOES7 Lo, 71 = GOES7 Hi
80 = GOES8 Lo, 81 = GOES8 Hi, 90 = GOES9 Lo, 91 = GOES9 Hi
Only single channel allowed.
(/goes == GOES7 Lo or GOES9 Lo, selected by date)
batse plot combined BATSE (CGRO) lightcurve
fov if set, plot 2nd window with centre of SXT FOV's
sfm if set, plos 3rd window with SXT image from SFM file (if present)
duration duration of time interval plotted (def=75 mins)
exact plot time axis exactly as defined
fem if set, overplot yohkoh ephemeris data
normalize normalizes all lightcurved to 1.0
line if set, joints dots on X-windows plot
keep inhibits clear_utplot so that overplotting possible
over allows overplotting on by subsequent call(s)
nosetup if set, does not create the widget common blocks (WIDGET only)
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLS: ***
ADDTIME [1], ADDTIME [2], ALL_VALS [1], ALL_VALS [2], BCS_TVEC, CLEAR_UTPLOT [1]
CLEAR_UTPLOT [2], GOOD_MAX [1], GOOD_MAX [2], HXT_TVEC, LINCOL, LOADCT, OBS_DEFAULTS
OBS_SETUP, OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], PMM [3], SXT_TVEC, UTPLOT [1]
UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], WBS_TVEC, anytim2ex [1]
anytim2ex [2], anytim2ints [1], anytim2ints [2], anytim2weeks [1]
anytim2weeks [2], bcs_channels, ex2fid [1], ex2fid [2], fem_grid [1], fem_grid [2]
fid2ex [1], fid2ex [2], fmt_tim [1], fmt_tim [2], get_sfm, gt_filta, gt_filtb, gt_grs1
gt_grs2, gt_hxs [1], gt_hxs [2], gt_rbmsc, gt_rbmsd, gt_sum_h [1], gt_sum_h [2]
gt_sum_l [1], gt_sum_l [2], gt_sum_m1 [1], gt_sum_m1 [2], gt_sum_m2 [1]
gt_sum_m2 [2], gt_sxs1 [1], gt_sxs1 [2], gt_sxs2 [1], gt_sxs2 [2], plot_fov [1]
plot_fov [2], pmm [1], pmm [2], rd_gbl, rd_gxd [1], rd_gxd [2], rd_gxd [3], rd_obs
CALLED BY:
OBS_EVENT
COMMON BLOCKS:
obs_plot_save Carries some parameters for /over
SIDE EFFECTS:
On X-terminal, loads a colour table for line plotting into bottom
16 levels. Revert to normal colour table with "loadct".
RESTRICTIONS:
Can plot a maximum of two instruments at a time (except SXT)
Hardcopy will only copy the timelime window (0)
PROCEDURE:
Plots with UTPLOT. If several channels are requested for an
instrument, then they are overplotted. First instrument uses
left-hand axis and second uses right-hand axis.
If no keywords set, defaults to the following:
BCS CaXIX, HXT L, SXT All Filters, /fem, /fov
MODIFICATION HISTORY:
RDB 20-Jan-94 Written
RBD 21-Jan-94 Allowed more than one BCS channel, added duration
RDB 29-Jan-94 Added normalize keyword, implemented HXT SUM_L
and multiple instrument plots
RDB 2-Feb-94 Allowed normal date string for start time
RDB 18-Feb-94 Added other HXT and all WBS channels
RDB 16-Jun-94 Corrected error if single BCS channel is array
Added /exact, goes and /batse keywords
Color lines if on X-terminal
Added /over and markers for SXT images
RDB 27-Jun-94 Added definition of usersym for plots
Added /fov and /line keywords
modified adding of osp/osf structs.
RDB 6-Jul-94 Allow array of channels for sxt
Correct test for date/fid at start
RDB 6-Aug-94 Corrected file used, batse array,
removed trap n_exp for batse and goes
Added tstop parameter, corrected missing channel test
RDB 23-Aug-94 Corrected SXT time vector for exposure
Added /keep keyword - i.e. clear_utplot normally issued
RDB 29-Aug-94 Added /fem keyword
RDB 17-Nov-94 Added /widget and /default keywords
RDB 24-Feb-95 Removed /widget; created WOBS_PLOT as a widget front
end to obs_plot - allows repeated modifications
Uses set of defaults if not keywords provided
- previously defaulted to BCS Ca XIX...
RDB 24-Apr-95 Supress FOV plot if !d.name is 'PS'
RDB 10-May-95 Added SFM option and harcopy related stuff
RDB 26-Jan-96 Added RBMSC and RBMSD to WBS channels
RDB 16-Oct-96 Added GOES7/stuff - modifed WOBS_plot as well
RDB 13-Nov-97 Finished off GOES 7/9 select by date.
Fixed a few bugs
[Previous]
[Next]
NAME:
OBS_SETUP
PURPOSE:
Creates the arrays in common needed by OBS_WIDGET
CATEGORY:
CALLING SEQUENCE:
obs_setup
INPUTS:
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLED BY:
OBS_PLOT, WOBS_PLOT
COMMON BLOCKS:
obs_tbl and obs_cm2
SIDE EFFECTS:
Most things filled with zero, duration set to 40 to allow modification
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB Feb-95 Written
RDB 11-Nov-97 Changed GOES setup (goes8/9 added)
[Previous]
[Next]
NAME:
OBS_STRUCT
PURPOSE:
Define the following observing log database structures
* Obs_NewOrb_Hd_Rec
* Obs_NewOrbit_Rec
* Obs_FileID_Rec
* Obs_OrbitSol_Rec
* Obs_WBSHXT_Rec
* Obs_BCS_Rec
* Obs_BCS_Status_Rec
* Obs_HXT_Status_Rec
* Obs_WBS_Status_Rec
* Obs_SXT_Rec
CALLING SEQUENCE:
OBS_STRUCT
CALLED BY:
GET_SDA, Rd_NewOrb_p [1], Rd_NewOrb_p [2], mk_att [1], mk_att [2], mk_evn [1]
mk_evn [2], mk_gbl, mk_gsn_obs, mk_gx, mk_pnt, mk_week1_file, mk_week_file [1]
mk_week_file [2], pnt2atr, rd_gsn, rd_obs, rd_old_obs, rd_week_file [1]
rd_week_file [2], rd_week_file [3], rd_week_file [4], weekid [2], ydb_exist [2]
HISTORY:
written by Mons Morrison, Fall 90.
[Previous]
[Next]
NAME:
obs_summary
PURPOSE:
Given a short time span, read the observing log and print (a) the
number of datasets for each instrument, (b) the peak flux in HXT-LOW
and SXS2, (c) the GOES flare event level, and (d) the SXT PFI image
summary (#images for each filter/resolution, DP mode/rate, heliocentric
location, NOAA active region #)
SAMPLE CALLING SEQUENCE:
obs_summary, sttim, entim
obs_summary, '1-jun-93', '1-jun-93 2:00', outfil='limb_search.txt'
obs_summary, evn(i), anytim2ints(evn(i), off=evn(i).duration), /hc
INPUT:
sttim - The starting time to search
entim - The ending time to search
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], DPRINT, FILE_EXIST [2], GET_NBYTES
RD_GEV, SEQ_SUMMARY, anytim2ints [1], anytim2ints [2], concat_dir [4]
data_type [1], data_type [2], file_exist [1], file_exist [3], fmt_tim [1]
fmt_tim [2], get_most_comm [1], get_most_comm [2], gt_center, gt_day [1]
gt_day [2], gt_sum_l [1], gt_sum_l [2], gt_sxs2 [1], gt_sxs2 [2], gt_time [1]
gt_time [2], rd_obs, tim2match [1], tim2match [2], tim2match [3]
CALLED BY:
sxtpnt_sum
OPTIONAL KEYWORD INPUT:
outfil - If set, save the output to the file name (default is to
display to the screen). If you call with /OUTFIL then the
output file name is "obs_summary.txt"
append - If set, append the results to the output file if it exists
hc - If set, send the results to the printer
HISTORY:
Written 28-Jan-93 by M.Morrison
19-Aug-93 (MDM) - Modified to use /EXTRA_OUT switch on SEQ_SUMMARY
1-Oct-93 (MDM) - Added /APPEND option
- Modified the summary line to hold the PFI location
in heliocentric coordinates and the GOES event level
4-Oct-93 (MDM) - Modified so that the default is to come to the screen
- Added /HC option
- Added documentation header
[Previous]
[Next]
NAME:
OBS_WIDGET
PURPOSE:
Sets up widget for WOBS_PLOT
CATEGORY:
CALLING SEQUENCE:
obs_widget
INPUTS:
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLS: ***
OBS_DEFAULTS, XMANAGER, XMENU [1], XMENU [2]
CALLED BY:
WOBS_PLOT
COMMON BLOCKS:
obs_tbl, obs_cms, obs_btn, obs_times
SIDE EFFECTS:
RESTRICTIONS:
should only be called as part of wobs_plot
PROCEDURE:
MODIFICATION HISTORY:
RDB Feb-95 Written
RDB 24-Apr-95 Added Hardcopy button and code
Does not plot the FOV window in this mode...
RDB 10-May-95 Added hardcopy and SFM options, modifed widget layout
RDB 13-Nov-97 New Widget Layout; now has start AND stop time;
extra WBS channels; NEW GOES spacecraft, etc.
[Previous]
[Next]
NAME:
ocenter
PURPOSE:
obtain the center of the Sun from optical full frame image
DESCRIPTION:
makes a 'silhouette image' consisting of 1's and 0's, then
does a centroid to determine x and y. The area is used to
determine the image radius. The input image can be compressed
or restored, and in any summation mode.
CATEGORY:
image processing
CALLING SEQUENCE:
ocenter,image,x_center,y_center, radius, cut = cut
INPUTS:
image
KEYWORDS;
cut = fraction of image maximum to define the limb. Default 0.2.
OUTPUTS:
x_center = the Sun center in the direction of EW
y_center = the Sun center in the direction of NS
radius = solar radius
CALLED BY:
mk_sot
NOTE: all outputs in full-resolution pixels
CALLS: ***
CENTROIDW [1], CENTROIDW [2], sxt_decomp [1], sxt_decomp [2]
COMMON BLOCKS:
none
HISTORY & AUTHOR:
Ver. 1.0 Jun 16, 1992
written by Hirohisa Hara, Univ. of Tokyo
Ver. 2.0 Jun 21, 1992, Hugh Hudson
- read image size automatically
- make decompression an option (allow use of low-8 images)
- threshold keyword added
- radius determination
- output in full-resolution pixels, adjusted for summation mode
[Previous]
[Next]
NAME:
ocontour
PURPOSE:
To overlay a contour on top of an image
CALLING SEQUENCE:
ocontour, img
ocontour, data(*,*,5), bin=6
ocontour, img, 100, 100
INPUT:
img - the image to be contoured on top of the
image aleady displayed with "TV"
OPTIONAL INPUT:
x0 - If the origin of the image is not (0,0) then it is
necessary to pass the left corner
y0 - see "x0"
bin - If the image was displayed using REBIN or CONGRID
then the routine needs to know that so that it scales
properly. This is the rebinning factor that was used
to display the first image.
OPTIONAL INPUT:
levels - IDL contour command
max_value- IDL contour command
color - IDL contour command
c_thick - IDL contour command
c_linestyle - IDL contour command
nlevels - IDL contour command
c_labels- IDL contour command
c_color - IDL contour command
c_charsize - IDL contour command
c_annotation - IDL contour command
c_orientation - IDL contour command
c_spacing - IDL contour command
CALLS: ***
OCONTOUR2
CALLED BY:
FAXABLE_SFD, boundary, eitoversxt, go_teem_nn, go_teem_t, lcur_plotr [1]
lcur_plotr [2], makimg0, rest_low8_cube [1], rest_low8_cube [2], rest_low8_ffi
sector_avg, sel_img, xy_lwa
HISTORY:
Written 28-Oct-91 by M.Morrison
[Previous]
[Next]
NAME:
ocontour
PURPOSE:
To overlay a contour on top of an image
CALLING SEQUENCE:
ocontour, img
ocontour, data(*,*,5), bin=6
ocontour, img, 100, 100
INPUT:
img - the image to be contoured on top of the
image aleady displayed with "TV"
OPTIONAL INPUT:
x0 - If the origin of the image is not (0,0) then it is
necessary to pass the left corner
y0 - see "x0"
bin - If the image was displayed using REBIN or CONGRID
then the routine needs to know that so that it scales
properly. This is the rebinning factor that was used
to display the first image.
OPTIONAL INPUT:
levels - IDL contour command
max_value- IDL contour command
color - IDL contour command
c_thick - IDL contour command
c_linestyle - IDL contour command
nlevels - IDL contour command
c_labels- IDL contour command
c_color - IDL contour command
c_charsize - IDL contour command
c_annotation - IDL contour command
c_orientation - IDL contour command
c_spacing - IDL contour command
CALLS: ***
OCONTOUR2
CALLED BY:
FAXABLE_SFD, boundary, eitoversxt, go_teem_nn, go_teem_t, lcur_plotr [1]
lcur_plotr [2], makimg0, rest_low8_cube [1], rest_low8_cube [2], rest_low8_ffi
sector_avg, sel_img, xy_lwa
HISTORY:
Written 28-Oct-91 by M.Morrison
[Previous]
[Next]
NAME:
ocontour
PURPOSE:
To overlay a contour on top of an image
EXPLANATION:
Note that the SOHO-CDS version of OCONTOUR.PRO was displaced by
this version. To get the same functionality of the SOHO-CDS
version, use the CONTOUR command with /OVERPLOT switch.
CALLING SEQUENCE:
ocontour, img
ocontour, data(*,*,5), bin=6
ocontour, img, 100, 100
INPUT:
img - the image to be contoured on top of the
image aleady displayed with "TV"
OPTIONAL INPUT:
x0 - If the origin of the image is not (0,0) then it is
necessary to pass the left corner
y0 - see "x0"
bin - If the image was displayed using REBIN or CONGRID
then the routine needs to know that so that it scales
properly. This is the rebinning factor that was used
to display the first image.
OPTIONAL INPUT:
levels - IDL contour command
max_value- IDL contour command
color - IDL contour command
c_thick - IDL contour command
c_linestyle - IDL contour command
nlevels - IDL contour command
c_labels- IDL contour command
c_color - IDL contour command
c_charsize - IDL contour command
c_annotation - IDL contour command
c_orientation - IDL contour command
c_spacing - IDL contour command
nodata - IDL contour command
tv2 - If set, then the image display is being done with TV2 and
if the output device is PS, then set the position a bit
differently.
CALLS: ***
WCONTOUR, tbeep [1], tbeep [2], tbeep [3]
CALLED BY:
FAXABLE_SFD, boundary, eitoversxt, go_teem_nn, go_teem_t, lcur_plotr [1]
lcur_plotr [2], makimg0, rest_low8_cube [1], rest_low8_cube [2], rest_low8_ffi
sector_avg, sel_img, xy_lwa
HISTORY:
Written 28-Oct-91 by M.Morrison
12-May-94 (MDM) - Corrected typo on "c_labels" parameter
12-May-94 (MDM) - Added /TV2 option
6-Sep-94 (MDM) - Included a check that the window is the
proper size (large enough to fit the image)
16-Feb-95 (MDM) - Added /NODATA option
9-Mar-95 (MDM) - Fixed a typo on the window size check
(it was using "x0" instead of "y0")
25-May-95 (LWA) - Added /c_orientation and /c_spacing options.
21-Jun-95 (MDM) - Updated header information (SOHO-CDS note)
18-Jul-95 (MDM) - Added /NORMAL to the contour call for /TV2 option
[Previous]
[Next]
NAME : odcenter (procedure)
PURPOSE :
obtain center of occulting disk & its radius
CATEGORY :
idl/lib/nkr
CALLING SEQUENCE :
odcenter,img,ox,oy,r,rxy=rxy
INPUTS :
img -- 10-cm image data
OPTIONAL INPUT PARAMETERS :
none
KEYWORD PARAMETERS :
rxy -- ration between x & y-radius
OUTPUTS :
ox,oy,r -- center and radius
CALLS: ***
GET2EDGES, ODCENTER
COMMON BLOCKS : none
SIDE EFFECTS :
if /ver is set, fitting curve is displayed
RESTRICTIONS : none
PROCEDURE :
MODIFICATION HISTORY :
K.I. 23 Jan.1993
[Previous]
[Next]
NAME:
OFF_AXIS
PURPOSE:
Compute the off-axis sensitivity reponse function for the SXT.
The function is normalized to 1.0 on-axis and drops off as the
edge of the SXT field of view is reached.
CALLING SEQUENCE:
Result_array = OFF_AXIS(x0,nx,y0,ny[,sum_mode=sum_mode,version=version])
In order to correct an SXT image for off-axis sensitivity, use the
following expression:
IMG = IMG / OFF_AXIS(x0,nx,y0,ny)
In this case, it is assumed that IMG is an nx by ny array whose
starting coordinates are at x0 and y0 (in CCD pixel coordinates).
INPUTS:
X0,Y0 = The starting coordinates in CCD pixel numbers.
X refers to CCD column numbers and Y to row numbers.
(The first non-BLS column is considered X0 = 0.)
NX,NY = The number of columns and rows.
OUTPUTS:
Functional result is the off-axis sensitivity response of SXT.
The result is a 2-d floating array normalized to 1.0 at the center
of the SXT field of view.
OPTIONAL KEYWORD PARAMETERS:
VERSION = (Output) String containing the version number of the
response function.
SUM_MODE =(Input) The assumed summation mode. 0=1x1, 1=2x2, 2=4x4
If not present, SUM_MODE = 0 is assumed.
COMMON BLOCKS;
None.
FILE I/O:
None.
METHOD:
Uses a 2-d polynomial representation to compute the SXT response
function for the desired position in the field of view.
CALLS: ***
POLY
MODIFICATION HISTORY:
Version 1.0 - JRL, 4-sep-91
VERSION = 'V1.0 4-sep-91'
[Previous]
[Next]
NAME:
ofr2fre
PURPOSE:
To convert the OFR (Old Full Resolution) pixel results to the
new FRE (Full Resolution Equivalent) pixel coordinates
SAMPLE CALLING SEQUENCE:
ofr2fre, index, sot
ofr2fre, index, sxc
INPUT:
prefix - The prefix being processed ('SOT' or 'SXC')
index - The index structure
INPUT/OUTPUT:
data - The SOT or SXC data structure
CALLS: ***
gt_res
CALLED BY:
rd_sot, rd_sxc
HISTORY:
Written 14-Jul-93 by M.Morrison
[Previous]
[Next]
NAME:
old_raster
PURPOSE:
To make a raster pattern of a time series of SXT images
to display time evolution of something or other.
The images start at upper left and go to lower right,
and the formatted time is written at the bottom of each
image.
CALLING SEQUENCE:
input data, index, and rebin factor (e.g. 2 to double)
this routine opens IDL window 1 with the right size and
writes to it.
COMMENT:
this is the old version of XY_RASTER, which has some additional
features. This one uses TVSCL and allows an arbitrarily large
number of images to fit into its window.
CALLS: ***
fmt_tim [1], fmt_tim [2]
HISTORY:
HSH, written May, 1992
[Previous]
[Next]
NAME: OMOVIE
PURPOSE: Show a cyclic sequence of images stored in a 3D array.
Omovie loads a special color table for overlay movies.
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:
omovie, 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 procedure movie
JPW, Feb, 1992. Added display of image number, single step mode,
and overlay color table
[Previous]
[Next]
NAME:
one_ftpt
PURPOSE:
Make up a two footpoint source, to test imaging programs on
CALLING SEQUENCE:
image = one_ftpt(btot, npix=npix, mid_pt=mid_pt, gaussian=gaussian, cyl=cyl, $
moff = moff, r1 = r1, point_source=point_source)
INPUT:
btot= total counts in image
OUTPUT:
A one footpoint image
KEYWORDS:
npix = Image will be npixXnpix, default = 64
mid_pt = midpoint of the footpoint, default is
[(npix-1)/2, (npix-1)/2]
gaussian= gaussian profile for the footpoint, the default
cyl= if set, the fp is cylindrical.
moff = if set, use cylindrical moffat fns, f=c/(1+r^2/r0^2)^moff(0)
r1 = radius for each footpoint in pixels, the default is 1
point_source= the footpoint is a point source
CALLS: ***
cyl_2d, cyl_moff, gauss_2d
CALLED BY:
find_ftpts
HISTORY:
11-aug-94, jmm
[Previous]
[Next]
Name: oneline
Purpose: widget front end to fastdoc
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], DELVARX [1], DELVARX [2]
DELVARX [3], DELVARX [4], XMANAGER, XMENU [1], XMENU [2], concat_dir [4]
delvarx [5], rd_tfile [1], rd_tfile [2]
CALLED BY:
fastdoc [1], fastdoc [2]
History:
10-Jan-1994 (SLF)
[Previous]
[Next]
Name: oneline_event
Purpose: event driver for oneline
CALLS: ***
DELVARX [1], DELVARX [2], DELVARX [3], DELVARX [4], EVENT_NAME, GET_WUVALUE
GET_WVALUE, XMANAGER, delvarx [5], path_sw [1], path_sw [2], rd_tfile [1]
rd_tfile [2], wc_where [1], wc_where [2]
History:
10-Jan-1994 (SLF)
[Previous]
[Next]
NAME:
OP_ADD_SPECIAL
PURPOSE:
Add the special commands for op_first_guess
CALLING SEQUENCE:
event2 = op_add_special( event2, og_table, op_rules, spec_times, spec_og )
CALLS: ***
ADDTIME [1], ADDTIME [2], BSORT [1], BSORT [2], BSORT [3], BSORT [4], DATATYPE [1]
DATATYPE [2], DATATYPE [3], Ex2Int [1], Ex2Int [2], anytim2ex [1], anytim2ex [2]
save_event [1], save_event [2]
CALLED BY:
make_32 [1], op_edit [1], op_edit [2], op_time_sort
HISTORY:
6-oct-94, J. R. Lemen, Written
[Previous]
[Next]
NAME:
op_bdr_opt
PURPOSE:
optimize the filling of the BDR
CALLING SEQUENCE:
op_bdr_opt, event, og_table, qbdr_opt=qbdr_opt
INPUTS:
event - op Event structure
og_table - OG table structure
qbdr_opt - 0: Don't optimize -- just summarize
1: Only move QT/H command (don't insert new commands)
2: Insert QT/H commands as necessary
CALLS: ***
ADDTIME [1], ADDTIME [2], MORE [1], MORE [2], OP_BDROPTIMIZE, OP_BDR_CNT
OP_START_END, gt_time [1], gt_time [2], save_event [1], save_event [2]
RESTRICTIONS:
The .dn flag (day/night) in the event structure must be
set by calling get_dn_rate first
HISTORY:
26-feb-97, J. R. Lemen, Written
7-mar-97, J. R. Lemen, Fixed slight error in the way OGs were counted
[Previous]
[Next]
NAME:
op_edit
PURPOSE:
Select the orbital events for the Yohkoh op command list.
CALLS: ***
ADDTIME [1], ADDTIME [2], BSORT [1], BSORT [2], BSORT [3], BSORT [4], Ex2Int [1]
Ex2Int [2], OP_ADD_SPECIAL, contact_sum, fmt_pass, op_get_event [1]
op_get_event [2], op_same_time [1], op_same_time [2], op_terminator [1]
op_terminator [2], save_event [1], save_event [2], valid_pass, valid_pass2
valid_pass3, valid_pass4 [1], valid_pass4 [2]
CALLED BY:
make_32 [1], op_time_sort
History:
15-nov-91, Written, J. Lemen
10-dec-91, Updated, J. Lemen -- Retain KSC los at the end; add pass_sum
11-jan-92, Modified,J. Lemen -- Bug fixed to remove Night og's
6-jul-92, Modified, JRL -- Fixed bug where 1st og was sometimes
day. Now it forces NOP.
13-oct-93, JRL, Added terminator image option
30-nov-93, JRL, Moved code to eliminate pre- and post-KSC events
to valid_pass4. Added call to op_same_time
14-jun-94, JRL, Added qrate_change keyword
29-jun-94, JRL, Pass og_table to op_same_time
30-jul-94, JRL, Pass ksc_str_remove to valid_pass4\
6-oct-94, JRL, Add spec_times, spec_og for special OG
28-mar-95, JRL, Call op_get_event twice to get Wallops passes
29-aug-95, JRL, Added a time sort before putting in the opstop command
25-feb-97, JRL, Added qbdr_opt switch
8-mar-97, JRL, Changed terminator keyword to term_max; added day_og keyword
[Previous]
[Next]
NAME:
op_edit
PURPOSE:
Select the orbital events for the Yohkoh op command list.
CALLS: ***
ADDTIME [1], ADDTIME [2], BSORT [1], BSORT [2], BSORT [3], BSORT [4], Ex2Int [1]
Ex2Int [2], OP_ADD_SPECIAL, contact_sum, fmt_pass, op_get_event [1]
op_get_event [2], op_same_time [1], op_same_time [2], op_terminator [1]
op_terminator [2], save_event [1], save_event [2], valid_pass, valid_pass2
valid_pass3, valid_pass4 [1], valid_pass4 [2]
CALLED BY:
make_32 [1], op_time_sort
History:
15-nov-91, Written, J. Lemen
10-dec-91, Updated, J. Lemen -- Retain KSC los at the end; add pass_sum
11-jan-92, Modified,J. Lemen -- Bug fixed to remove Night og's
6-jul-92, Modified, JRL -- Fixed bug where 1st og was sometimes
day. Now it forces NOP.
13-oct-93, JRL, Added terminator image option
30-nov-93, JRL, Moved code to eliminate pre- and post-KSC events
to valid_pass4. Added call to op_same_time
14-jun-94, JRL, Added qrate_change keyword
29-jun-94, JRL, Pass og_table to op_same_time
30-jul-94, JRL, Pass ksc_str_remove to valid_pass4\
6-oct-94, JRL, Add spec_times, spec_og for special OG
28-mar-95, JRL, Call op_get_event twice to get Wallops passes
29-aug-95, JRL, Added a time sort before putting in the opstop command
[Previous]
[Next]
NAME:
op_facom_ftp
PURPOSE:
Transfer the op_first_guess file to the FACOM
Assumes s88239 account.
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], anytim2ex [1], anytim2ex [2]
break_file [4], facom_ftp
MODIFICATION HISTORY:
5-may-93, J. Lemen, Written.
21-Feb-96, JRL, Updated for new mainframe GS8400/10
[Previous]
[Next]
ROUTINE: op_first_guess
PURPOSE: Generate a first guess for the OP editor tables.
Reads the output file of AOSLOS (as run on FACOM
mainframe).
Reads a lookup of table of OG commands.
Reads a files containing "rules" for computing
og times.
For diagnostic mode, set noprint=1 before running to prevent printing.
MODIFICATION HISTORY:
15-oct-91, JRL, V0.0 Written
21-oct-91, JRL, V0.2 Added Hex output option
28-oct-91, JRL, V1.0 Added an information message
14-nov-91, JRL, V1.1 DSN commanding added.
22-nov-91, JRL, V1.2 Default output path to /3p/guest/opog
10-dec-91, JRL, V1.3 Summarize KSC and DSN contacts
11-dec-91, JRL, V1.4 Improved wording of one question
14-jan-92, JRL, V1.5 Add DSN los night/day op rules
23-jan-92, JRL, V1.6 Add dsn_switch op rule.
6-jul-92, JRL, V1.7 Changed the re-start sequence.
5-aug-92, JRL, V1.8 Changed print-out of station summary
11-aug-92, JRL, V1.9 Added call to pass_sheets; Changed
default output directory for isass1,isass4
11-sep-92, JRL, V2.0 Fixed pass_sheets, Read RASM file;
6-feb-93, JRL, V2.1 Changed call to print_flare1
27-apr-93, JRL, V2.2 Change default directory for opog output
Added call to cover_page in pass_sheets
5-may-93, JRL, V2.3 Changed DSN Rep logic at night. Changed output names.
Added ftp code.
6-may-93, JRL, V2.4 Fixed the output name.
14-oct-93, JRL, V2.5 Align times to 32sec boundaries; Add SXT terminator images
30-nov-93, JRL, V2.6 Change the way the start time is handled.
If AOS Ant is different from actual on pass sheet, put the
actual antenna in ( )'s to try to prevent misunderstanding.
Added call to op_same_time (from op_edit)
Changed KSC pre-antenna algorithm to catch all cases
where antenna change occurs before 5 degrees.
3-dec-93, JRL, V2.7 Fixed a bug in pass_sheet involved with
position of Day OG (error made in V2.6 mod)
6-dec-93, JRL, V2.8 Fixed a bug in AOS selection (V2.6).
Force output times '000000' to be '240000'
16-jun-94, JRL, V2.9 Added code to switch to Quiet Med Rate during SAA
Added new rule for op_terminator (must be 10s before D/N)
19-jun-94, JRL, V3.0 Don't ask about FTP to mainframe if .hex file is not created.
25-jun-94, JRL, V3.1 Display number of OG slots used before asking about .txt file
30-jun-94, JRL, V3.2 Changed op_same_time to display the correct conflic entry numbers
2-jul-94, JRL, V3.3 Fixed a bug with pre-contact antenna change
30-jul-94, JRL, V3.4 Fixed a bug with QT Med Rate during SAA
Allow user to delete any pass (for ASCA conflicts)
7-oct-94, JRL, V3.5 Added op_special options
16-nov-94, JRL, V3.6 Make 'No' and 'Yes' responses to questions work as 'N' and 'Y'
Fixed a bug op_saa_med
9-dec-94, JRL, V3.7 Changed timing of CMPN OBS OFF to 6 * 32s after DAY OG
28-mar-95, JRL, V3.8 Add Wallops
22-oct-95, JRL, V3.9 Changed printing procedure
5-feb-96, JRL, V4.0 Changed the AOS algorithm -- new op_first_guess.rule also
26-feb-97, JRL, V4.1 Added BDR optimization code. Updated terminator selection
(term_rule4 and rule 6). Changed op_saa_med so that changes
are not made when all slots in the op are used up.
8-mar-97, JRL, V4.2 Terminators inserted by changing briefly from QT/M to QT/H.
Simplified some questions.
11-mar-97, JRL, V4.3 Improved method to set 2nd day KSC passes to QT/M
Don't allow terminators between last 2 KSC contacts
[Previous]
[Next]
routine: op_first_guess
purpose: Generate a first guess for the OP editor tables.
Reads the output file of AOSLOS (as run on FACOM
mainframe).
Reads a lookup of table of OG commands.
Reads a files containing "rules" for computing
og times.
For diagnostic mode, set noprint=1 before running to prevent printing.
Modification History:
15-oct-91, JRL, V0.0 Written
21-oct-91, JRL, V0.2 Added Hex output option
28-oct-91, JRL, V1.0 Added an information message
14-nov-91, JRL, V1.1 DSN commanding added.
22-nov-91, JRL, V1.2 Default output path to /3p/guest/opog
10-dec-91, JRL, V1.3 Summarize KSC and DSN contacts
11-dec-91, JRL, V1.4 Improved wording of one question
14-jan-92, JRL, V1.5 Add DSN los night/day op rules
23-jan-92, JRL, V1.6 Add dsn_switch op rule.
6-jul-92, JRL, V1.7 Changed the re-start sequence.
5-aug-92, JRL, V1.8 Changed print-out of station summary
11-aug-92, JRL, V1.9 Added call to pass_sheets; Changed
default output directory for isass1,isass4
11-sep-92, JRL, V2.0 Fixed pass_sheets, Read RASM file;
6-feb-93, JRL, V2.1 Changed call to print_flare1
27-apr-93, JRL, V2.2 Change default directory for opog output
Added call to cover_page in pass_sheets
5-may-93, JRL, V2.3 Changed DSN Rep logic at night. Changed
output names. Added ftp code.
6-may-93, JRL, V2.4 Fixed the output name.
14-oct-93, JRL, V2.5 Align times to 32sec boundaries;
Add SXT terminator images
30-nov-93, JRL, V2.6 Change the way the start time is handled.
If AOS Ant is different from actual on pass sheet, put the
actual antenna in ( )'s to try to prevent misunderstanding.
Added call to op_same_time (from op_edit)
Changed KSC pre-antenna algorithm to catch all cases
where antenna change occurs before 5 degrees.
3-dec-93, JRL, V2.7 Fixed a bug in pass_sheet involved with
position of Day OG (error made in V2.6 mod)
6-dec-93, JRL, V2.8 Fixed a bug in AOS selection (V2.6).
Force output times '000000' to be '240000'
16-jun-94, JRL, V2.9 Added code to switch to Quiet Med Rate during SAA
Added new rule for op_terminator (must be 10s before D/N)
19-jun-94, JRL, V3.0 Don't ask about FTP to mainframe if .hex file is not created.
25-jun-94, JRL, V3.1 Display number of OG slots used before asking about .txt file
30-jun-94, JRL, V3.2 Changed op_same_time to display the correct conflic entry numbers
2-jul-94, JRL, V3.3 Fixed a bug with pre-contact antenna change
30-jul-94, JRL, V3.4 Fixed a bug with QT Med Rate during SAA
Allow user to delete any pass (for ASCA conflicts)
7-oct-94, JRL, V3.5 Added op_special options
16-nov-94, JRL, V3.6 Make 'No' and 'Yes' responses to questions work as 'N' and 'Y'
Fixed a bug op_saa_med
9-dec-94, JRL, V3.7 Changed timing of CMPN OBS OFF to 6 * 32s after DAY OG
28-mar-95, JRL, V3.8 Add Wallops
22-oct-95, JRL, V3.9 Changed printing procedure
5-feb-96, JRL, V4.0 Changed the AOS algorithm -- new op_first_guess.rule also
[Previous]
[Next]
NAME:
op_get_event
PURPOSE:
Get orbital events for Yohkoh
INPUTS:
orbit = Strcture of predicted orbital events.
og_table = Structure of og's
op_rules = Structure of rules for generating the op load.
OPTIONAL INPUT PARAMETERS:
DSN = If set, will compute length of DSN contacts based on
actual AOS and LOS. Otherwise, the BDR rep, pause and
los will be based on parameters in the op_rules structure.
FIELD = To select which of contacts to use (0 or 1)
OUTPUTS:
outarr = Structure array of og commands.
CALLS: ***
ADDTIME [1], ADDTIME [2], int2secarr [1], int2secarr [2], save_event [1]
save_event [2]
CALLED BY:
make_32 [1], op_edit [1], op_edit [2], op_time_sort, save_event [1]
save_event [2]
MODIFICATION HISTORY:
7-oct-91, written, JRL
28-oct-91, Updated, JRL -- Enabled the night OG
12-nov-91, Updated, JRL -- Add DSN rep and pause og's
14-nov-91, Updated, JRL -- Changed AOS / LOS alogorithm to detect change
by using the " "/"+" or "+"/" " transitions.
12-dec-91, Modified, JRL -- Fixed bug with antenna change during DSN
15-jan-92, Modified, JRL -- Added dsn_los_d/dsn_los_n distinction
Changed NOP alogorithm
23-jan-92, Modified, JRL -- Fixed bug: Assign DN_flag incase DSN
station encountered before D/N flag
1-jul-92, Modified, JRL -- Fixed bugs: Fixed BDR pause before DSN REP.
Fixed undefined entry_ant for one case.
6-jul-92, Modified, JRL -- Select different KSC AOS antenna if the change
is at less than 5 degs.
13-jul-92, Modified, JRL -- Select KSC AOS as NT if the planned AOS comes
before the DAY OG.
23-sep-92, Modified, JRL -- Check for D/N between AOS and AOS+1.5 min
4-may-93, JRL, Added DN flag to the DSN BDR Rep event
30-nov-93, JRL, Fixed KSC AOS if antenna changes bewteen '+' and 'U'
6-dec-93, JRL, Fixed bug in previous modification
29-mar-94, JRL, Fixed crashed caused by orbit(0).station = '+'
2-jul-94, JRL, Still trying to make antenna changes between '+' and '<anything>' work
17-dec-94, JRL, Decoupled 10 deg criteria and Min KSC antenna value
28-mar-95, JRL, Make routine work with new style orbit structure.
Renamed to op_get_event
8-jun-95, JRL, Make sure timarr_day is initialized to some value.
3-feb-96, JRL, Changed the AOS antenna algorithm.
8-Mar-97, JRL, Changed when the Pass ID is defined (at the '' to '+' transition)
[Previous]
[Next]
NAME:
op_get_event
PURPOSE:
Get orbital events for Yohkoh
INPUTS:
orbit = Strcture of predicted orbital events.
og_table = Structure of og's
op_rules = Structure of rules for generating the op load.
OPTIONAL INPUT PARAMETERS:
DSN = If set, will compute length of DSN contacts based on
actual AOS and LOS. Otherwise, the BDR rep, pause and
los will be based on parameters in the op_rules structure.
FIELD = To select which of contacts to use (0 or 1)
OUTPUTS:
outarr = Structure array of og commands.
CALLS: ***
ADDTIME [1], ADDTIME [2], int2secarr [1], int2secarr [2], save_event [1]
save_event [2]
CALLED BY:
make_32 [1], op_edit [1], op_edit [2], op_time_sort, save_event [1]
save_event [2]
MODIFICATION HISTORY:
7-oct-91, written, JRL
28-oct-91, Updated, JRL -- Enabled the night OG
12-nov-91, Updated, JRL -- Add DSN rep and pause og's
14-nov-91, Updated, JRL -- Changed AOS / LOS alogorithm to detect change
by using the " "/"+" or "+"/" " transitions.
12-dec-91, Modified, JRL -- Fixed bug with antenna change during DSN
15-jan-92, Modified, JRL -- Added dsn_los_d/dsn_los_n distinction
Changed NOP alogorithm
23-jan-92, Modified, JRL -- Fixed bug: Assign DN_flag incase DSN
station encountered before D/N flag
1-jul-92, Modified, JRL -- Fixed bugs: Fixed BDR pause before DSN REP.
Fixed undefined entry_ant for one case.
6-jul-92, Modified, JRL -- Select different KSC AOS antenna if the change
is at less than 5 degs.
13-jul-92, Modified, JRL -- Select KSC AOS as NT if the planned AOS comes
before the DAY OG.
23-sep-92, Modified, JRL -- Check for D/N between AOS and AOS+1.5 min
4-may-93, JRL, Added DN flag to the DSN BDR Rep event
30-nov-93, JRL, Fixed KSC AOS if antenna changes bewteen '+' and 'U'
6-dec-93, JRL, Fixed bug in previous modification
29-mar-94, JRL, Fixed crashed caused by orbit(0).station = '+'
2-jul-94, JRL, Still trying to make antenna changes between '+' and '<anything>' work
17-dec-94, JRL, Decoupled 10 deg criteria and Min KSC antenna value
28-mar-95, JRL, Make routine work with new style orbit structure.
Renamed to op_get_event
8-jun-95, JRL, Make sure timarr_day is initialized to some value.
3-feb-96, JRL, Changed the AOS antenna algorithm.
[Previous]
[Next]
NAME:
op_get_special
PURPOSE:
Get the times for Special commands and return
list of special command sequence in a structure.
CALLING SEQUENCE:
op_get_special, op_special_file, spec_times, spec_og
INPUTS:
op_special_file - File name (if file does not exist, -1 is returned)
OUTPUTS:
spec_times = Start time(s) in external format (7,N)
spec_og = Structure containing the sequence of OGs and relative times
CALLS: ***
DELVARX [1], DELVARX [2], DELVARX [3], DELVARX [4], FILE_EXIST [2], anytim2ex [1]
anytim2ex [2], delvarx [5], file_exist [1], file_exist [3], fmt_tim [1]
fmt_tim [2], input [1], input [2], op_rd_special, yesnox [1], yesnox [2]
HISTORY:
7-oct-94, J. R. Lemen, Written
[Previous]
[Next]
NAME:
OP_MF2WS
PURPOSE:
reformat the ASCII version of an print file into single column
CALLING SEQUENCE:
ws_file = op_mf2ws(mf_file)
OPTIONAL KEYWORD INPUT
RESTRICTIONS
HISTORY:
97/10/20 H. Hudson
[Previous]
[Next]
Name:
op_pass_sheets
Purpose:
Called by op_first_guess to generate the pass sheets.
Calling Sequence:
op_pass_sheets,contacts,version=version,comment=comment, $
orbit_params=orbit_params [,/qtest,file=file]
Optional Input Keywords:
file = The IDL plot file. If not specified, defaults to idl.ps
This routine will create the plot file, but not print it.
CALLS: ***
ADDTIME [1], ADDTIME [2], OP_COVER_PAGE, REVERSE, clearplot [1], clearplot [2]
gt_time [1], gt_time [2], tbeep [1], tbeep [2], tbeep [3]
Modification History:
12-aug-92, written, JRL
23-sep-92, Modified, JRL. Fixed some bugs; Added TMX H REP for Night entry
27-apr-93, JRL, Call Dave Pike's cover_page routine. Removed MOD CTL AT
30-nov-93, JRL, Always force the entry comment (KSC AOS) to be printed
if the AOS antenna doesn't match the actual antenna, put
the actual antenna in ( ) (Kosugi san's suggestion).
3-dec-93, JRL, Fixed the DAY OG position (bec of previous change)
6-dec-93, JRL, Fixed the "non-printing" of SXT CHK DIS when AOS is -NT.
30-jul-94, JRL, Pass ksc_pass_num and weekend keywords to cover_page
22-oct-95, JRL, Renamed, op_pass_sheets. Do a clearplot first
8-mar-97, JRL, Added day_og_evn keyword
8-oct-97, JRL, Fixed a bug when the Day OG event is printed
and small format update.
[Previous]
[Next]
Name:
op_pass_sheets
Purpose:
Called by op_first_guess to generate the pass sheets.
Calling Sequence:
op_pass_sheets,contacts,version=version,comment=comment, $
orbit_params=orbit_params [,/qtest,file=file]
Optional Input Keywords:
file = The IDL plot file. If not specified, defaults to idl.ps
This routine will create the plot file, but not print it.
CALLS: ***
ADDTIME [1], ADDTIME [2], OP_COVER_PAGE, REVERSE, clearplot [1], clearplot [2]
gt_time [1], gt_time [2], tbeep [1], tbeep [2], tbeep [3]
Modification History:
12-aug-92, written, JRL
23-sep-92, Modified, JRL. Fixed some bugs; Added TMX H REP for Night entry
27-apr-93, JRL, Call Dave Pike's cover_page routine. Removed MOD CTL AT
30-nov-93, JRL, Always force the entry comment (KSC AOS) to be printed
if the AOS antenna doesn't match the actual antenna, put
the actual antenna in ( ) (Kosugi san's suggestion).
3-dec-93, JRL, Fixed the DAY OG position (bec of previous change)
6-dec-93, JRL, Fixed the "non-printing" of SXT CHK DIS when AOS is -NT.
30-jul-94, JRL, Pass ksc_pass_num and weekend keywords to cover_page
22-oct-95, JRL, Renamed, op_pass_sheets. Do a clearplot first
[Previous]
[Next]
NAME:
op_rd_special
PURPOSE:
Read special block of commands for op_first_guess
If filnam does not exist, return -1, else return a structure
CALLING SEQUENCE:
special = op_rd_special(filnam)
CALLS: ***
FILE_EXIST [2], STR2ARR [1], STR2ARR [2], file_exist [1], file_exist [3]
rd_tfile [1], rd_tfile [2]
CALLED BY:
op_get_special
HISTORY:
6-oct-94, J. R. Lemen, Written
[Previous]
[Next]
Name:
op_saa_med
Purpose:
Called by valid_pass4 (called by op_edit from op_first_guess)
Switch SAA times into MED RATE
Will switch to MED RATE at RDCT ON + op_rules.med_rate.a
Will switch to HIGH RATE at RDCT OFF - op_rules.med_rate.a
Don't make the switch if:
1. SAA is entirely in Night
2. Already in MED rate at RD HV RDCT ON
3. Length of Daylight time in MED rate is less
than op_rules.saa_min_time
4. Any DSN operations occur within the SAA
CALLS: ***
ADDTIME [1], ADDTIME [2], fmt_tim [1], fmt_tim [2], get_dn_rate [1]
get_dn_rate [2], save_event [1], save_event [2]
CALLED BY:
get_dn_rate [1], get_dn_rate [2], make_32 [2], valid_pass4 [1], valid_pass4 [2]
RESTRICTIONS:
Assumes that event is time-ordered on input.
Modification History:
14-jun-94, J. R. Lemen, Written
24-jun-94, JRL, Changed the diagnostic statement
16-nov-94, JRL, Fixed logic for rule 3
28-mar-95, JRL, The pass_num for the rate change is now set to -1
8-mar-97, JRL, Try to estimate total OGs and quit adding if we reach (128-1)
[Previous]
[Next]
Name:
op_saa_med
Purpose:
Called by valid_pass4 (called by op_edit from op_first_guess)
Switch SAA times into MED RATE
Will switch to MED RATE at RDCT ON + op_rules.med_rate.a
Will switch to HIGH RATE at RDCT OFF - op_rules.med_rate.a
Don't make the switch if:
1. SAA is entirely in Night
2. Already in MED rate at RD HV RDCT ON
3. Length of Daylight time in MED rate is less
than op_rules.saa_min_time
4. Any DSN operations occur with the SAA
CALLS: ***
ADDTIME [1], ADDTIME [2], fmt_tim [1], fmt_tim [2], get_dn_rate [1]
get_dn_rate [2], save_event [1], save_event [2]
CALLED BY:
get_dn_rate [1], get_dn_rate [2], make_32 [2], valid_pass4 [1], valid_pass4 [2]
RESTRICTIONS:
Assumes that event is time-ordered on input.
No time ordering performed on the output array
Modification History:
14-jun-94, J. R. Lemen, Written
24-jun-94, JRL, Changed the diagnostic statement
16-nov-94, JRL, Fixed logic for rule 3
28-mar-95, JRL, The pass_num for the rate change is now set to -1
[Previous]
[Next]
NAME:
op_same_time
PURPOSE:
Called by op_edit (by op_first_guess) to resolve conflicting OG times.
CALLS: ***
ADDTIME [1], ADDTIME [2], BSORT [1], BSORT [2], BSORT [3], BSORT [4], Ex2Int [1]
Ex2Int [2], tbeep [1], tbeep [2], tbeep [3]
CALLED BY:
make_32 [1], op_edit [1], op_edit [2], op_time_sort
PROCEDURE:
There are six rules:
1. If the conflict involes OP STOP and anything,
move the OP STOP later by 32 sec
2. If the conflict involves the QUIET RATE HIGH command & an LOS command,
change the LOS to QT/H and remove the QUIETT RATE HIGH OG.
3. If the conflict involves the QUIET RATE HIGH, move earlier 32 sec
(unless rules 1 or 2 apply).
4. If the conflict involves the OP DAY, move later 32 sec (unless
rules 1 or 2 apply).
5. If the conflict involves the QUIET RATE MED and an RB HV RDCT command,
move the QUIET RATE MED command 32 sec later.
6. If the conflict involves the DSN REPRO command and an antenna change
command, back the antenna change command up by 32 sec and make it a
"ANT CHNG TO X" command.
7. If the conflict involves any other type of conflict, just warn the user.
HISTORY:
30-nov-93, JRL, Written
30-jun-94, JRL, When warning the user, give the same event number as
shows up in the printout.
27-jul-94, JRL, Make the warning more prominent. Added OP DAY case.
26-feb-97, JRL, Added rule for QUIET RATE HIGH. This was made necessary because
BDR optimization makes conflicts more likely.
8-mar-97, JRL, Added additional rule for QUIET RATE HIGH and for QUIET RATE
MED, which was made necessary by changes to terminator selection.
[Previous]
[Next]
NAME:
op_same_time
PURPOSE:
Called by op_edit (by op_first_guess) to resolve conflicting OG times.
CALLS: ***
ADDTIME [1], ADDTIME [2], BSORT [1], BSORT [2], BSORT [3], BSORT [4], Ex2Int [1]
Ex2Int [2], tbeep [1], tbeep [2], tbeep [3]
CALLED BY:
make_32 [1], op_edit [1], op_edit [2], op_time_sort
PROCEDURE:
There are three rules:
1. If the conflict involes OP STOP and anything,
move the OP STOP later by 32 sec
2. If the conflict involves the OP DAY, move later than 32 sec (unless
rule 1 applies).
3. If the conflict involves the DSN REPRO command and an antenna change
command, back the antenna change command up by 32 s and make it a
"ANT CHNG TO X" command.
4. If the conflict involves any other type of conflict, just warn the user.
HISTORY:
30-nov-93, JRL, Written
30-jun-94, JRL, When warning the user, give the same event number as
shows up in the printout.
27-jul-94, JRL, Make the warning more prominent. Added Case 3 above.
[Previous]
[Next]
NAME:
op_term_score
PURPOSE:
Determine how well op_first_guess is scheduling terminators
CALLING SEQUENCE:
op_term_score
op_term_score, time1, time2
op_term_score, buff=buff ; Return string output array in buff
op_term_score, outfile=outfile ; Write results to a file for printing
Reads the file on $ydb/orbit/opog/term_times.txt . This file is
updated by a cron job which reads recent op*txt files.
OPTIONAL INPUTS:
time1, time2 - Specify the range times
outfile - If filename exists, write text to file
OPTIONAL OUTPUTS:
buff - String array with results
files_sfc - Names of the sfc files
index_accidental - Indicies of the accidentals
If time1 and time are not specified, then the first and last time
of the term_times.txt file is assumed.
CALLS: ***
ADDTIME [1], ADDTIME [2], BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3]
DELVARX [1], DELVARX [2], DELVARX [3], DELVARX [4], FIND_FILE, MORE [1], MORE [2]
RD_SFC, RD_TERM_TIMES, Rd_Roadmap [2], STDEV, STR2ARR [1], STR2ARR [2], anytim [1]
anytim [2], anytim [3], anytim [4], anytim [5], break_file [4], delvarx [5]
fmt_tim [1], fmt_tim [2], gt_filta, gt_filtb, gt_res, int2secarr [1], int2secarr [2]
prstr [1], prstr [2], rd_roadmap [1], rd_tfile [1], rd_tfile [2], sxt_cen [1]
sxt_cen [2], tim2dset [1], tim2dset [2], tim2orbit [1], tim2orbit [2]
HISTORY:
27-Feb-97, V1.0, J. R. Lemen, Written
12-Mar-97, V1.1, JRL, Improved the header information
20-Jun-99, V1.2, H. Hudson, added NOVERB and XY keywords
30-Jan-00, V1.3, DMCK and HSH patched up the listing
[Previous]
[Next]
NAME:
op_term_sum
PURPOSE:
Called by op_first_guess routines to find the times of
the programmed terminator images. Returns an information string
CALLING SEQUENCE:
Text_string = op_term_sum(event2,num)
INPUTS:
event2 = The op_edit event data structure
OPTIONAL OUTPUTS:
num = The number of terminator images.
CALLS: ***
fmt_tim [1], fmt_tim [2]
CALLED BY:
op_terminator [1], op_terminator [2]
MODIFICATION HISTORY:
14-oct-93, J. R. Lemen (LPARL), Written
[Previous]
[Next]
NAME:
op_terminator
PURPOSE:
Try to put the command sequences for the terminator
images into the op command.
Called by op_edit.
CALLING SEQUENCE:
event2 = op_terminator(event,og_table,num)
INPUTS:
event = The op_first_guess (from get_event) data structure
og_table = The op_first_guess og_table structure
OPTIONAL INPUT KEYWORDS:
term_max = Specify an upper limit to the number of terminators.
If -1, all possible terminators are returned.
RETURNED:
The event structure is modified with all possible terminators
inserted.
OUTPUTS:
num = the number of terminators inserted.
Terminator Rules:
1. Do not insert if in SAA
2. Do not insert terminators on the second KSC day
3. If terminator will occur at less than 21s (Term_rule2) prior to D/N
terminator, do not select. [Added 15-jun-94/mod 26-feb-97]
4. Do not insert if a DSN or KSC contact occurs during or before
the next DAY OG + Term_rule1*32 sec [Changed 9-dec-94]
5. Do not insert if QT/M OG occurs during or before the next
the next DAY OG + Term_rule1*32 sec [Added 8-Mar-97]
6. If 2 D/N's occur between AOS's, take the first one.
7. If the OGs exceed 128, quit adding terminators
If the number of terminators exceeds term_max, quit adding terminators
Command Sequence:
Time Command
D/N - 3*32s CNT MAN (3=Term_rule0 below)
D/N - 2*32s Campaign Obs On
D/N - 1*32s CNT AT
DAY + 6*32s Campaign Obs Off (4=Term_rule1 below)
CALLS: ***
ADDTIME [1], ADDTIME [2], Ex2Int [1], Ex2Int [2], fmt_tim [1], fmt_tim [2]
op_term_sum, prstr [1], prstr [2], save_event [1], save_event [2], tim2orbit [1]
tim2orbit [2]
CALLED BY:
make_32 [1], op_edit [1], op_edit [2], op_time_sort
MODIFICATION HISTORY:
14-oct-93, J. R. Lemen, Written
15-jun-94, JRL, Added logic for New rule 3 (less than 10s prior to D/N)
This based on a good suggestion by Barry Labonte.
30-jun-94, JRL, Added a debug statment for local diagnostics
9-dec-94, JRL, Changed value of Term_rule1 to 6 units.
26-feb-97, JRL, Changed value of Term_Rule2 to 21s;
Don't add terminator if no slots left in OP (Rule 7).
8-mar-97, JRL, Changed logic to change MED to HIGH to obtain terminators
Added the term_max keyword.
11-mar-97, JRL, Don't put terminator between last two KSC passes
[Previous]
[Next]
NAME:
op_terminator
PURPOSE:
Try to put the command sequences for the terminator
images into the op command.
CALLING SEQUENCE:
event2 = op_terminator(event,num)
INPUTS:
event = the op_first_guess /get_event data structure
RETURNED:
The event structure is modified with all possible terminators
inserted.
OUTPUTS:
num = the number of terminators inserted.
Terminator Rules:
1. Do not insert if in SAA
2. Do not insert if in Medium rate
3. Do not insert if a DSN or KSC contact occurs during or before
the next DAY OG + 6*32 sec [Changed 9-dec-94]
4. If terminator will occur at less than 10s prior to D/N
terminator, do not select. [Added 15-jun-94]
5. If 2 D/N's occur between AOS's, take the first one.
Command Sequence:
Time Command
D/N - 3*32s CNT MAN (3=Term_rule0 below)
D/N - 2*32s Campaign Obs On
D/N - 1*32s CNT AT
DAY + 6*32s Campaign Obs Off (4=Term_rule1 below)
CALLS: ***
ADDTIME [1], ADDTIME [2], Ex2Int [1], Ex2Int [2], fmt_tim [1], fmt_tim [2]
op_term_sum, prstr [1], prstr [2], save_event [1], save_event [2], tim2orbit [1]
tim2orbit [2]
CALLED BY:
make_32 [1], op_edit [1], op_edit [2], op_time_sort
MODIFICATION HISTORY:
14-oct-93, J. R. Lemen, Written
15-jun-94, JRL, Added logic for New rule 4 (less than 10s prior to D/N)
This based on a good suggestion by Barry Labonte.
30-jun-94, JRL, Added a debug statment for local diagnostics
9-dec-94, JRL, Changed value of Term_rule1 to 6 units.
[Previous]
[Next]
NAME:
op_time_sort
PURPOSE:
Time order the event structure
Also, call make_32 to align to 32s boundaries
INPUTS:
event = The event structure from op_get_event
RETURNS:
A new event structure with modified times.
OPTIONAL INPUT KEYWORDS:
no32 = If set, don't call make_32
CALLS: ***
ADDTIME [1], ADDTIME [2], BSORT [1], BSORT [2], BSORT [3], BSORT [4], Ex2Int [1]
Ex2Int [2], OP_ADD_SPECIAL, contact_sum, fmt_pass, op_edit [1], op_edit [2]
op_get_event [1], op_get_event [2], op_same_time [1], op_same_time [2]
op_terminator [1], op_terminator [2], save_event [1], save_event [2], valid_pass
valid_pass2, valid_pass3, valid_pass4 [1], valid_pass4 [2]
MODIFICATION HISTORY:
25-feb-97, JRL, Changed calls to addtime to be with the /sec switch set
[Previous]
[Next]
Name:
op_times
Purpose:
Get start and stop times to run op_first_guess
INPUTS:
filein_aoslos = Aoslos file name or lun
OUTPUTS:
start_time,end_time= Start and end times
start_rec,end_rec = Start and end record numbers
ksc_first = First KSC contact time
OPTIONAL INPUT KEYWORDS:
stime = To specify the the current system time
(For diagnostic purposes)
CALLS: ***
ADDTIME [1], ADDTIME [2], DELVARX [1], DELVARX [2], DELVARX [3], DELVARX [4]
GET_AOSREC, anytim2ints [1], anytim2ints [2], contacts [1], contacts [2]
delvarx [5], fmt_tim [1], fmt_tim [2], input [1], input [2], rd_fem, rd_orbit_event
tbeep [1], tbeep [2], tbeep [3], timstr2ex [1], timstr2ex [2], yesnox [1]
yesnox [2]
Modification History:
Written, 16-oct-91, J. Lemen
Updated, 11-dec-91, J. Lemen, A little error checking on input added
Updated, 14-sep-92, JRL, Format of default start time improved.
Updated, 27-apr-93, JRL, Force start time to be 1 hour prior to actual
KSC contact. calls rd_fem to get info.
27-nov-93, JRL, Rewritten. Now back off 10 min from user entered time and then search
ahead for the next KSC contact. Then back up 1 hour.
16-jun-94, JRL, Display selected passes using contacts
12-dec-99, JRL, Updated to fix y2k bugs
[Previous]
[Next]
NAME:
OPEN_SSH_WIN
SELECT_SHIFT
CLOSE_SSH_WIN
PURPOSE:
SELECT SHIFT quantity BY MOUSE. (called from ALIGN_MAN)
CATEGORY:
CALLING SEQUENCE:
W=OPEN_SSH_WIN()
DV=SELECT_SHIFT(W)
CLOSE_SSH_WIN(W)
INPUTS:
OUTPUTS:
CALLS: ***
CLOSE_SSH_WIN, DEL_DIGIT, INCR_X, INCR_Y, IS_LARGE_DN, IS_LARGE_LEFT, IS_LARGE_RIGHT
IS_LARGE_UP, IS_OK, IS_SMALL_DN, IS_SMALL_LEFT, IS_SMALL_RIGHT, IS_SMALL_UP
SELECT_SHIFT
CALLED BY:
ALIGN_MAN
COMMON BLOCKS:
SIDE DEFFECTS:
RESTRICTIONS:
INLUDE THIS FILE AT THE TOP OF THE SOURCE FILE,
ADDING THE FOLLOWING
@SELECT_SHIFT
MODIFICATION HISTORY:
Goro Kawai 25 Jan., 1993
EXAMPLES:
[Previous]
[Next]
NAME: oplin
PURPOSE: draw a line on an image or cube
CALLING SEQUENCE: ss=oplin(image,x0,y0,x1,y1,color)
PARAMETERS: image is the image (or cube) to be drawn on
x0,y0 are the beginning point coordinates
x1,y1 are the end point coordinates
color value to be given to the line pixels
ss return the indices of the points on the line
CALLS: LNPNT,CUBE_DEPTH
HISTORY: drafted nov. 1994, A. McAllister
[Previous]
[Next]
NAME:
OPLOT_B
PURPOSE:
Overlay a (vector) magnetogram on top of an image. The magnetogram
is displayed as contours and line segments. Optionally, a linear
transformation can be performed on the magnetogram to have it match
an underlying image.
CATEGORY:
CALLING SEQUENCE:
oplot_b,bl,bt,ba
oplot_b,bz,bx,by,/bzxy
INPUTS:
bl = longitudinal B-field. 2-D array.
bt = transverse B-field strength. Optional.
ba = azimuth angle of transverse B-field. Angle counterclockwise
from (geocentric) N. (-> N-E-S-W)
With /bzxy option:
bz = vertical B-field. 2-D array.
bx = transverse B-field x-component. Optional.
by = transverse B-field y-component. Optional.
KEYWORDS (INPUT):
/bzxy = input vector field supplied as z/x/y components. Use 2nd
form of call. This option not well tested! Check proper
direction of transverse field!
pp = float(2,2,*). Point array for transformation. Refer to the
procedure "setpts" for details.
key = transformation type. Used in conjunction with pp.
See "transtype" keyword in CALTRANS for details.
Default is a general linear transform.
xy = float(2). x/y offset for magnetogram in units of input pixels.
("dejitter" offset.)
zoom = integer. Zoom factor of the transformed output.
size = integer. Size of (zoomed) output image.
offset = integer. Output image offset in pixel units of output image.
xoff = offset of plot in display window
yoff = offset of plot in display window
levels = vector of contour levels to be displayed.
/nodash = if set, negative contours will solid instead of dashed.
/c_spline = cubic spline interpolation of contours.
c_line = line style of contours. Overrides /nodash. Separate styles
for pos. and neg. contours can be given as a 2-elements vector.
c_color = color of contours. Separate colors for pos. and neg.
contours can be given as a 2-elements vector.
CALLS: ***
CALTRANS, pp_reverse
CALLED BY:
cont_movie
COMMON BLOCKS:
None.
SIDE EFFECTS:
RESTRICTIONS:
The transverse field vectors will be distorted if the original
magnetic field image has a pixel aspect ratio other than 1.
PROCEDURE:
MODIFICATION HISTORY:
JPW, modification of an older program (plotbnew), Feb. 1992
JdLB, simplified plotting of negative contour levels, Oct. 1992
JPW, added c_spline, c_line, and c_color keywords.
JPW, added xy keyword, Jan. 1993
JPW, added "key" keyword, and call to pp_reverse, 22-nov-94
[Previous]
[Next]
NAME:
OPLOT_DIFF
PURPOSE:
Obtain the differential distribution of the "count" and
oplot it on the current device.
INPUTS:
count - the data (1-d)
bin - the bin size (10^bin)
interval - the total observing time
OPTIONAL INPUTS:
fitrange = [low, high]
the energy range for the power-law fitting.
nofit= no fitting
OUTPUTS:
CALLS: ***
FCHECK, OPLOTERR, PLOT_POWERLAW
SIDE EFFECTS:
1, Input the fitting range by the mouse.
2, Plot on the current device.
REGISTRATIONS:
HISTORY:
20 July, 93 written by T.Shimizu (Univ. of Tokyo, JAPAN)
07 Aug, 93 Add the uncertainty assummed Possion statistics.
28 Dec, 93 Revise /errbar
20 Feb, 94 Add /nofit
20 Feb, 94 "oplot_diff" was created from "plot_diff"
[Previous]
[Next]
NAME:
oplot_nts
PURPOSE:
To overplot the Nobeyama correlation coefficient time series (NTS)
SAMPLE CALLING SEQUENCE:
oplot_nts
oplot_nts, .6
oplot_nts, .6, .5
OPTIONAL INPUT:
yfst - The starting location on the plot in fractions (0 to 1)
If not passed, it will use 0.75
yfh - The height of the NTS plot in fractions (0 to 1)
If not passed, it will use 0.23
OPTIONAL KEYWORD OUTPUT:
nts - The NTS data structure
CALLS: ***
OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], RD_NTS, anytim2ints [1], anytim2ints [2]
deriv_arr [1], deriv_arr [2], fmt_tim [1], fmt_tim [2], int2secarr [1]
int2secarr [2]
HISTORY:
Written 17-May-93 by M.Morrison
[Previous]
[Next]
NAME:
OPLOTERR_XY
PURPOSE:
Overplot data points with accompanying error bars.
(See also PLOTERR_XY.)
CATEGORY:
Plotting, two-dimensional.
CALLING SEQUENCE:
PLOTERR, X , Y , Err_X, ERR_Y [, linestyle = linestyle]
INPUTS:
X: The array of abcissae.
Y: The array of Y values.
Err_X: The array of error-bar values in X.
Err_Y: The array of error-bar values in Y.
Format of Err_X & Err_Y
[Err_X_min, Err_X_max]: Specified Low & high values.
Err_X (same dimmension as X) : Specified the width of errors.
Err_Y is same as above.
OPTIONAL KEYWORD PARAMETERS:
linestyle: The line style to use. The default is +0.
COMMON BLOCKS:
None.
SIDE EFFECTS:
None.
RESTRICTIONS:
Arrays cannot be of type string. There must be enough points to
plot.
PROCEDURE:
A plot of X with error bars drawn from X - ERR_X to X - ERR_X
versus Y with those drawn from Y - ERR to Y + ERR is written
to the output device.
MODIFICATION HISTORY:
William Thompson Applied Research Corporation
July, 1986 8201 Corporate Drive
Landover, MD 20785
DMS, April, 1989 Modified for Unix.
24 Dec, 1993 Made with the base of PLOTERR.
by T.Shimizu
[Previous]
[Next]
NAME:
ORB_EXTRAP
PURPOSE:
anticipate orbit times for future events
CATEGORY:
Yohkoh
CALLING SEQUENCE:
orb_extrap, target_tim
INPUTS:
target time (Yohkoh string format)
OPTIONAL (KEYWORD) INPUT PARAMETERS:
tt_arr is the output; leave this in place to re-do the plot quickly
OUTPUTS:
tt_arr(3,20): 0, the epoch of the elements
1, the first n/d terminator time
2, the next d/n terminator time
CALLS: ***
ORBIT_DN, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
WDEF [1], WDEF [2], fmt_tim [1], fmt_tim [2], gt_time [1], gt_time [2], rd_tfile [1]
rd_tfile [2]
COMMON BLOCKS:
SIDE EFFECTS:
grabs IDL 0 and makes a plot
calls ORBIT_DN, which makes a temporary ASCII file
in the current directory and then deletes all files
named '*ascii' <== be careful
RESTRICTIONS:
it won't work unless the target time is in orbit night, sorry
MODIFICATION HISTORY:
HSH, written 17-jun-97
[Previous]
[Next]
NAME:
ORB_STRUCT
PURPOSE:
Structure containing information about the orbital position of
the spacecraft, typically every 300 secs.
Defines orb_record
CALLING SEQUENCE:
ORB_STRUCT
CALLED BY:
MK_ORB_WEEK, rd_week_file [1], rd_week_file [2], rd_week_file [3]
rd_week_file [4], weekid [2], ydb_exist [2]
HISTORY:
Written 28-Aug-96 by R.D.Bentley
[Previous]
[Next]
NAME:
ORBIT_DN
PURPOSE:
find D/N transitions
CATEGORY:
yohkoh
CALLING SEQUENCE:
INPUTS:
orbit file in the form $DIR_GEN_ORBIT_SOL/orbit*
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
CALLS: ***
AOS2EX, EX2AOS, anytim [1], anytim [2], anytim [3], anytim [4], anytim [5]
fmt_tim [1], fmt_tim [2], rd_tfile [1], rd_tfile [2], run_dsnfil
CALLED BY:
ORB_EXTRAP
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
HSH, written 16-jun-97
[Previous]
[Next]
NAME:
ORBIT_INPUT_TIME
PURPOSE:
compute a string for MK_DSNFIL entry
CATEGORY:
CALLING SEQUENCE:
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
CALLS: ***
anytim [1], anytim [2], anytim [3], anytim [4], anytim [5]
CALLED BY:
GET_YOHKOH_POS
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
HSH, written 29-Jan-97
[Previous]
[Next]
Name: orcent
Purpose: return fov center in macro pixels for observing region
Input Parameters:
x,y - observing region size (NX,NY) - (units of 64 pixels)
res - resolution of observing region desired
Keyword Paramters:
half - if set, set res=2
quart - if set, set res=4
Calling Example:
(Assume FFI displayed in current window)
CALLS: ***
BOX_CURSOR, REVERSE, STR2ARR [1], STR2ARR [2]
History:
13-May-93 (SLF)
8-Nov-93 (SLF) add erase,thick,linestyle
[Previous]
[Next]
NAME:
OUT_GRSPHL1
PURPOSE:
To output grs_phl1 data structure.
CALLING SEQUENCE:
ph = out_grsphl1( index, data, stime = stime )
ph = out_grsphl1( index, data, /disp )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
stime - Signal time. for example, stime = ['22:37 11-NOV-91',
'22:38:32.114 11-NOV-91']
disp - If set, you can select the data on display by mouse.
OUTPUTS:
ph - WBS_GRSPHL1 structure.
OPTIONAL INPUTS:
int - Integration time of signal. for example, int = 4 then
integration time is 4 seconds. If present, divide the signal
time by int.
pname - String containing the name of printer.
If not present, the printer name is 'lp'.
CALLED BY:
GRS_SPEFF, TEST_OUTG
PROCEDURE:
Both Sttime and Edtime are the observation times. These times are
header time minus prestored time.
Data are deadtime-corrected.
CALLS: ***
DATAGET3_GRS, Ex2Int [1], Ex2Int [2], HXT_DECOMP, Int2Ex [1], Int2Ex [2]
anytim2ex [1], anytim2ex [2], dtcf_ph_grs1l, gt_day [1], gt_day [2]
gt_dp_rate [1], gt_dp_rate [2], gt_time [1], gt_time [2], header_scan
RESTRICTIONS:
Data are restricted to high and medium bit rates.
HISTORY:
Written 7-dec-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 13-dec-93 by T.Hiraoka ; added energy in output structure.
[Previous]
[Next]
NAME:
OUT_GRSPHL2
PURPOSE:
To output grs_phl2 data structure.
CALLING SEQUENCE:
ph = out_grsphl2( index, data, stime = stime )
ph = out_grsphl2( index, data, /disp )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
stime - Signal time. for example, stime = ['22:37 11-NOV-91',
'22:38:32.114 11-NOV-91']
disp - If set, you can select the data on display by mouse.
OUTPUTS:
ph - WBS_GRSPHL2 structure.
OPTIONAL INPUTS:
int - Integration time of signal. for example, int = 4 then
integration time is 4 seconds. If present, divide the signal
time by int.
pname - String containing the name of printer.
If not present, the printer name is 'lp'.
CALLED BY:
GRS_SPEFF, TEST_OUTG
PROCEDURE:
Both Sttime and Edtime are the observation times. These times are
header time minus prestored time.
Data are deadtime-corrected.
CALLS: ***
DATAGET3_GRS, DTCF_PH_GRS2L, Ex2Int [1], Ex2Int [2], HXT_DECOMP, Int2Ex [1]
Int2Ex [2], anytim2ex [1], anytim2ex [2], gt_day [1], gt_day [2], gt_dp_rate [1]
gt_dp_rate [2], gt_time [1], gt_time [2], header_scan
RESTRICTIONS:
Data are restricted to high and medium bit rates.
HISTORY:
Written 7-dec-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 13-dec-93 by T.Hiraoka ; added energy in output structure.
[Previous]
[Next]
NAME:
OUT_HXSPH
PURPOSE:
To output hxs_ph data structure.
CALLING SEQUENCE:
ph = out_hxsph( index, data, stime = stime )
ph = out_hxsph( index, data, /disp )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
stime - Signal time. for example, stime = ['22:37 11-NOV-91',
'22:38:32.114 11-NOV-91']
disp - If set, you can select the data on display by mouse.
OUTPUTS:
ph - WBS_HXSPH structure.
OPTIONAL INPUTS:
int - Integration time of signal. for example, int = 2 then
integration time is 2 seconds. If present, divide the signal
time by int.
pname - String containing the name of printer.
If not present, the printer name is 'lp'.
CALLED BY:
TEST_OUTH
PROCEDURE:
Both Sttime and Edtime are the observation times. These times are
header time minus prestored time.
Data are deadtime-corrected.
CALLS: ***
DATAGET3_HXS, Ex2Int [1], Ex2Int [2], HXT_DECOMP, Int2Ex [1], Int2Ex [2]
anytim2ex [1], anytim2ex [2], dtcf_ph_hxs, gt_day [1], gt_day [2], gt_dp_rate [1]
gt_dp_rate [2], gt_time [1], gt_time [2], header_scan
RESTRICTIONS:
Data are restricted to high and medium bit rates.
HISTORY:
Written 7-dec-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 13-dec-93 by T.Hiraoka ; added energy in output structure.
[Previous]
[Next]
Overlays a pair of images with alternating colortables.
PRO over_color, im1, im2
INPUT PARAMETERS:
im1, im2 = images to be overlaid.
METHOD:
Chooses image value for each pixel as image1 > image2.
Forces image1 to even values, image2 to odd values.
Uses actual number of available colors.
Try OVERLAY for a different approach.
RESTRICTIONS:
Images must be same size.
Scaling the image values must be done first.
Coalignment and sizing must be done first.
HISTORY:
Written May 17, 1994 Barry LaBonte
[Previous]
[Next]
NAME:
OVER_MOVIE
PURPOSE:
Overlays 2 movie sequences. Use "omovie,out" to run it.
movie_a is displayed in red, movie_b in blue.
Important: the result replaces movie_a.
CATEGORY:
CALLING SEQUENCE:
over_movie,movie_a,movie_b
INPUTS:
movie_a = first movie. Is destroyed and replaced by output.
movie_b = second movie. Must be of same size as ima. Not modified.
KEYWORDS (INPUT):
amax = intensity of ima will be scaled evenly between 0 and amax.
If not supplied, image scaled according to brightest pixel.
bmax = intensity of imb will be scaled evenly between 0 and bmax.
OUTPUTS:
movie_a. Replaces first input movie.
CALLED BY:
mk_coal_movie
COMMON BLOCKS:
None.
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
The idea is to use the lower 4 bit of the image display for one
image, and the upper 4 bit for the other image. Except that this
version only uses 15 * 15 = 225 colors instead of 16 * 16 = 256.
MODIFICATION HISTORY:
JPW, Jan. 92
[Previous]
[Next]
NAME:
OVER_THE_LIMB
PURPOSE:
Displays NOAA data and estimates over-the-limb flare locations
CATEGORY:
Yohkoh helpful stuff
CALLING SEQUENCE:
over_the_limb,t1,t2,ar,flare,otl_out
INPUTS:
t1, t2 search times (string format)
ar = NOAA active region number
flare = time of the 'unknown' flare (string form)
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
makes a plot of the residuals and throws in some comments
otl_out structure containing guesses about the 'unknown' flare
- number of flares
- extrapolated longitude corresponding to entered time
- scatter of points in linear fit to flare longitudes
- height of tangent ray (+- one sigma)
- apparent rotation rate of flare locations (deg/day)
- mean latitude of flares in the selected AR
CALLS: ***
POLY_FIT, RD_GEV, STDEV, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5]
UTPLOT [6], anytim2ints [1], anytim2ints [2], fmt_tim [1], fmt_tim [2], gt_day [1]
gt_day [2], gt_time [1], gt_time [2]
MODIFICATION HISTORY:
HSH written 30-July-1993
[Previous]
[Next]
NAME:
OVERLAY
PURPOSE:
Overlays 2 images and displays them on the screen. One image is
displayed in red, the other one in blue.
CATEGORY:
CALLING SEQUENCE:
overlay,ima,imb,xoffset,yoffset
INPUTS:
ima = first image. Must be a 2-D array.
imb = second image. Must be of same size as ima.
xoffset = int, optional. Image is displayed with offset in x.
yoffset = int, optional. Image is displayed with offset in y.
KEYWORDS (INPUT):
amax = intensity of ima will be scaled evenly between 0 and amax.
If not supplied, image scaled according to brightest pixel.
bmax = intensity of imb will be scaled evenly between 0 and bmax.
coffset = lowest intensity value for color map.
cgain = intensity gain for color map. Make shure that
coffset + cgain * 14 < 256. Optimum values for coffset and
cgain depend upon the display hardware.
OUTPUTS:
No explicit outputs.
COMMON BLOCKS:
None.
SIDE EFFECTS:
An image is displayed, and the color map is modified.
RESTRICTIONS:
Requires a color display with at least 225 colors.
PROCEDURE:
The idea is to use the lower 4 bit of the image display for one
image, and the upper 4 bit for the other image. Except that this
version only uses 15 * 15 = 225 colors instead of 16 * 16 = 256.
MODIFICATION HISTORY:
JPW, Jan. 92
[Previous]
[Next]
NAME:
OVERPLOT_FIT.PRO
PURPOSE:
OVERPLOTS SYNTHETIC FE25 SPECTRA ONTO OBSERVED BCS SPECTRA. Assumes
that synthtic spectrum is in an ascii file with wavelength in first
column and spectrum values in second column. You will be propted to
input the filename.
CALLING SEQUENCE:
.RUN OVERPLOT_FIT
INPUTS:
KEYWORDS:
OPTIONAL OUTPUTS:
PROCEDURE:
Must have BSC_INDEX and BSC_DATA structures available. Produce these
using GO_BSC, or RD_BSC, etc. Must specify lower and upper plot
limits. Good default values for Fe25 are X_LO=1.85 & X_HI=1.89. Then
do .run overplot_fit. Use cursor to select from the plotted lightcurve
the spectrum to be fit. You will then be asked to click on the
location of the peak of the reasonance line. The code scales the fit
spectrum to the most intense point in the observed spectrum. In some
cases (such as when the resonace line is not the strongest feature)
this scaling may not be what is desired. You can rescale the synthetic
spectrum mannually if desired by multiplying the spectral array, y,
by the factor and overplotting, eg, IDL> oplot,x,y*factor
To make hardcopies, set_plot,'ps' and rerun overplot_fit.
HISTORY:
8 Mar 93 A. Sterling (D. Zarro)
Light curve plot feature of this code created
by modifying Zarro's PLOT_BSC code.