[Previous]
[Next]
NAME:
page_fmt
PURPOSE:
Add form-feeds and number pages every "lines" lines in an ascii file
CALLING SEQUENCE:
page_fmt, infil, outfil
page_fmt, infil, outfil, 54
page_fmt, infil, outfil, 54, text='Page' ; Right justified
page_fmt, infil, outfil, 54, text='Page ??', /nonumber ; No page numbers
Existing <FF> (form-feed) characters are left in place and numbered
according. This allows the user to force new pages at desired locations.
\newpage commands are treated just as if it were a <FF> character.
Two blank lines are inserted after each form-feed.
If the next 2 lines following a pre-exising <FF> are blank, then no
new blank lines are inserted.
INPUTS:
infil = Input file name
outfil = Output file name
lines = Number of lines between form-feeds (default = 54)
OPTIONAL INPUT KEYWORDS:
text = Text to write (def = 'Page') -- right justified
nonumber = By default, will number the pages. If /nonumber, this will
be inhibited.
CALLS: ***
prstr [1], prstr [2], rd_tfile [1], rd_tfile [2]
MODIFICATION HISTORY:
8-feb-94, J. R. Lemen, Written.
8-sep-94, JRL, Enabled \newpage allow new pages to be inserted (so that you
don't have to know how to insert a ctl-L with the editor)
[Previous]
[Next]
NAME:
PAINT_PIXELS
PURPOSE:
Interactively correct "low 8" images. Current version assumes
64 x 64 PFI
CATEGORY:
Yohkoh analysis
CALLING SEQUENCE:
paint_pixels, image
MODIFICATION HISTORY:
HSH, Nov 1991
[Previous]
[Next]
Name:
PANALOOP
Purpose:
Play a range of frames on the Pansonic TQ-3031F laser videodisk
recorder/player or Pansonic TQ-3032F laser videodisk player
an arbitrary number of times
Calling sequence:
panaloop, f0, f1, rate_fac=rate_fac, reps=reps,
reverse=reverse, w=w, device=device
Input:
f0, f1 = start and end frames of sequence to play
Optional Input Keywords:
RATE_FAC = If positive, RATE_FAC is the multiplicative factor
applied the normal playback rate of 30 frames per sec.
If negative, RATE_FAC is the reciprocal of the multi-
plicative factor applied to the normal rate. Accepted
values are: -256,-255,-254,...,-1,1,2,3,...,10. Notice
that 0 is not accepted.
Examples: RATE_FAC = 10: 10 times normal rate, or
300 frames/sec
RATE_FAC = -10: 1/10 normal rate, or
3 frames/sec
DEFAULT = NORMAL RATE (30 FRAMES/SEC).
REPS = Number of playback repeats at each playback speed.
DEFAULT = ENDLESS (0)
REVERSE = 0 to jump from last to first frame
= 1 to reverse play from last frame to first frame
DEFAULT = 0
W = wait interval in seconds between playbacks
DEFAULT = 0
DEVICE = 0 for recorder/player (/dev/ttyd2)
= 1 for player (/dev/ttyd3)
DEFAULT = 0
Method:
Modification History:
Written, 31-aug-92 - gls
[Previous]
[Next]
NAME:
PARFIT
PURPOSE:
Force-fit a quadratic equation to three points
CALLING SEQUENCE
parfit,x,y,a,b,c
INPUTS:
X(3), Y(3)
OUTPUTS:
(a,b,c) of Y = A*X^2 + B*X + C
CALLS: ***
DETERM
HISTORY:
written by Hugh Hudson, July 2, 1992
[Previous]
[Next]
NAME:
ParsKey
Purpose:
Parse the key into an item, operator(s), and value(s).
CALLING SEQUENCE:
ParsKey, key, item, opr, val, [op2=opr2, val2=val2,
op3=opr3, val3=val3]
Input:
key string value of a search operation.
Output:
item value which corresponds to a data-vector field.
opr operator
val search value for item.
opr2&3 2nd and 3rd operators if present
val2&3 2nd and 3rd values for item
CALLS: ***
STR2ARR [1], STR2ARR [2], str_replace [1], str_replace [2]
CALLED BY:
KeySea [1], KeySea [2], KeySearch
History:
written 14-May-91 by GAL
2-Aug-91, added opr and val 1 &2.
24-oct-94, change remchar call to str_replace call
(avoid conflict with astron library routine)
[Previous]
[Next]
Name:
pass_sheets
Purpose:
Called by op_first_guess to generate the pass sheets.
Calling Sequence:
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], COVER_PAGE [1], COVER_PAGE [2], REVERSE, 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
[Previous]
[Next]
Name: path_sw
Purpose: return the path to on-line software using mapfile contents
Input Parmeters:
module - software module name or substring (.pro file, script, etc)
Optional Keyword Parameters:
nopath - if set, only return the module names (no path info)
nopro - if set, all file types in path matching module are returned
multi - if set, all occurences are returned (default is first)
(this should be enhanced to use !path to determine order)
readfile - force reading of file and regeneration of common block
swdirs - if set, return $ys directories (ignore module, if present)
Output Parameters: function returns full path specification
count - number of matches found
Calling Sequence:
fullpath=path_sw(module [,count , /nopath , multi, /nopro])
Category:
system, swmaint, util, gen
Method - on first call, reads generic file $DIR_GEN_SETUP/data/swmap.genx
into a common block - once common is established, finding module
is via. table lookup
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], RD_GENX, UNIQ [1], UNIQ [2]
UNIQ [3], concat_dir [4], str_replace [1], str_replace [2], wc_where [1]
wc_where [2]
CALLED BY:
add_pro [1], add_pro [2], fastdoc [1], fastdoc [2], oneline_event [1]
oneline_event [2], sw2tree [1], sw2tree [2], ucon_check
History: slf, 8-Aug-92
slf, 7-Jan-93 use DIR_GEN_SETUPD (=DIR_GEN_SETUP/data)
call wc_where to find pattern match
slf, 29-jan-93 use rd_genx instead of restgen (faster)
[Previous]
[Next]
NAME:
PATROL_ORDER
PURPOSE:
list the sequence of images preceding an SXT patrol image
CATEGORY:
CALLING SEQUENCE:
patrol_order, t_start, t_stop
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
CALLS: ***
FMT_TIMER [1], FMT_TIMER [2], get_info [1], get_info [2], get_info [3]
get_info [4], gt_shape_cmd, rd_obs, sort_index [1], sort_index [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
8-Aug-1997, written (HSH)
[Previous]
[Next]
NAME:
pc_tplot_grs
PURPOSE:
To analyze wbs/grs-pc in the reformat data
( under only grs_plot )
CALLS: ***
pprint [1], pprint [2], time_add1s, window_select, wmenu_sel [1], wmenu_sel [2]
CALLED BY:
grs_plot
HISTORY:
Written 20-jun-93 by T.Hiraoka
[Previous]
[Next]
NAME:
pc_tplot_rbm
PURPOSE:
To analyze wbs/rbm-pc in the reformat data
( under only grs_plot )
CALLS: ***
pprint [1], pprint [2], time_add1s, window_select, wmenu_sel [1], wmenu_sel [2]
CALLED BY:
grs_plot
HISTORY:
Written 6-jul-93 by T.Hiraoka
13-jul-93 K.Morimoto
[Previous]
[Next]
NAME:
pe_pl
PURPOSE:
Total photon energy flux of a single power law
CALLING SEQUENCE:
energy = pe_pl(k,m,e0,e1=e1,flat_to=flat_to)
INPUT:
m; spectral index
k; emission at 1keV
e0=low energy limit
OUTPUT:
energy = energy flux in photons
KEYWORDS:
e1 = the high energy limit of the calculation, default is 20MeV
flat_to=if set, the spectrum will be assumed to be flat down to
the value given in flat_to. If flat_to>e0, then flat_to=e0
HISTORY:
JMcT, 29-apr-94
[Previous]
[Next]
NAME:
PERASE
PURPOSE:
Load an empty 1280 x 1024 byte array into one of the Peritek
color planes.
CATEGORY:
Image display.
CALLING SEQUENCE:
PERASE, XSIZE, YSIZE, PLANE=PLANE, RESULT=RESULT
INPUTS:
XSIZE: X dimension of area to erase.
YSIZE: Y dimension of area to erase.
KEYWORD PARAMETERS:
PLANE: VCT plane ('red','green,','blue', or one other)
RESULT: Error message
OUTPUTS:
RESULT: Error message (optional)
PROCEDURE:
Spawns call to C routine VCTQPSEUDO which loads the array
MODIFICATION HISTORY:
10-Jul-1991 - Written by GLS
15-nov-1991 - Updated, JRL -- Added DIR_PERITEK environment
variable to reference C routines.
29-apr-1991 - Updated, GLS -- Use !version to determine which
C routine (V bus version or Q bus
version) to call.
[Previous]
[Next]
NAME:
Periph
PURPOSE:
decode the periph byte parameter and return aspect door
status, shutter mode, filter a and b positions.
CALLING SEQUENCE:
Periph, val, [aspDr=aspDr, shut=shut, filb=filb,
fila=fila]
INPUT:
val is the periph byte value.
OUTPUT/Keyword:
aspDr is the Aspect Door status: 1= open and 0=closed
shut is the shutter mode: 1=mech and 0=Frame Transfer
filb is the filter b position
fila is the filter a position
Notes:
Actual bit positions are reversed (i.e. 7 for 0) from
the documentation notes in the File Control Doc.
The comments for each bit operation will describe the
actual bit # mask and the (documented bit #).
CALLED BY:
GetStrExpLev, GetStrPeriph, LWA_TE, SXT_TE, get_active_reg, lwa_te_intact
HISTORY:
written Mar-91 by GAL
updated 12-Aug-91 to use Keywords for all formal
parameters
updated 16-Sep-91 to change input from Keyword to standard
formal parameter.
[Previous]
[Next]
NAME:
PFI_CENT
PURPOSE:
give solar coordinates for fov centers of PFI images
CALLING SEQUENCE:
nnee = pfi_cent(index) returns position(s) relative to sun
center, in arc sec N, E
nnee = pfi_cent(index,pnt_data) doesn't re-read the PNT data
nnee = pfi_cent(index,/hel) prints the heliographic coordinates
(and returns them to nnee)
INPUTS:
index entry(s) for image(s) to be located
OUTPUTS:
nnee positive coordinates (2, n)
pnt_data, if that is in the call
KEYWORDS:
/hel gives answer in heliographic coordinates (N,E)
WARNINGS:
DOES NOT HANDLE PIXEL SUMMATION MODE CORRECTLY YET! USE WITH
FULL RES ONLY PLEASE
CALLS: ***
ARCMIN2HEL, SXT_FOV_CENTER, fmt_tim [1], fmt_tim [2], gt_hxa, rd_pnt, tim2dset [1]
tim2dset [2]
HISTORY:
Hugh Hudson, Aug. 6, 1992
HSH added array entry and modified the documentation, Aug. 8, 1992
[Previous]
[Next]
NAME:
PFI_CENTER
PURPOSE:
Computes the position of the center of a PFI image in full-res SXT
pixels relative to sun center.
CATEGORY:
CALLING SEQUENCE:
position = pfi_center(index)
INPUTS:
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS
hxacode = Quality code from get_suncenter (output)
arcsec = if set, return the result in arcseconds rather than SXT pixels
heliographic = if set, return [longitude,latitude]
quiet = if set, work quietly
OUTPUTS:
position = fltarr(2,nindex) array containing the x,y coordinates of the
center of the PFI image, where nindex = number of elements in
the index array.
CALLS: ***
SXT_FOV_CENTER, conv_p2h [1], conv_p2h [2], gt_pix_size, sxt_cen [1], sxt_cen [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
Uses get_suncenter and sxt_fov_center to compute the position of the
center of the PFI image.
North and West are positive numbers.
MODIFICATION HISTORY:
T. Metcalf November, 1992
TRM June 18, 1993: CHANGED THE ORDER IN WHUCH THE HELIOGRAPHIC
COORDS ARE RETURNED FROM PFI_CENTER
TRM July 18, 1995: Changed the call to get_suncentr to sxt_cen.
[Previous]
[Next]
NAME:
pfi_composite
PURPOSE:
Prepare a composite PFI movie.
CALLING SEQUENCE:
movie = pfi_composite(index,data)
INPUTS:
index = SXT index structure
data = 3-d data cube (at least 2 images) of alternately
long and short (or short and long) exposures
If data is byte then decompress and leak_sub
If data is not byte type then assume already decompressed,
-- and background subtracted
OUTPUT:
Function result is a movie of composite frames.
plot the results as we go. The results are displayed
on a logarithmic scale.
CALLS: ***
sxt_composite [1], sxt_composite [2], sxt_composite [3]
MODIFICATION HISTORY:
LWA - August 1993
[Previous]
[Next]
NAME:
PFI_DATAGET2
PURPOSE:
Read the PFI image data from the PFI image file made by the image
generater "MKIMAGE".
This "PFI_DATAGET2" is for reading the data file being version 2
file format.
CATEGORY:
SOLAR-A/SXT Data Analysis
CALLING SEQUENCE:
image= PFI_DATAGET2(file, header)
INPUTS:
file= the file name. ex, /SXTpfi/910911/123156.dat on SXT1.
This file name must be the file that is the first ROI image file.
This variable is string type.
OUTPUTS:
image= the PFI(Matrix) image data.(ewsize * (64*nROI) bytes'
array)
header= the header of this image data.(512 * nROI bytes' array)
COMMON BLOCKS:
None.
SIDE EFFECTS:
None.
RESTRICTIONS:
This procedure can be used only in reading PFI version 2 file.
PROCEDURE:
open the file and read the header and image data. Then return
image and header.
PROGRAMMER:
Toshifumi Shimizu
Department of Astronomy, University of Tokyo, JAPAN
HISTORY:
Ver 1.0 21-SEP-1991 written
Ver 1.1 23-SEP-1991 correct wrong place
Ver 1.2 05-Oct-1992 rotate image data
[Previous]
[Next]
Name: pfi_dominant
Purpose: identify PFI dominant, return FIDS and SFR:SPR ratios
Input Parameters:
t0,t1 - start and stop time for search
Keyword Parameters:
cutoff (input) - PFI dominant ratio cutoff (default = FFI/PFI < .7 )
rd_obs (input) - switch, if set, read&return observing log / filter Flare
ss (output) - SS where FFI:PFI < cutoff
fids (output) - list of all fids within time range
ratios (output) - all FFI:PFI ratios within time range
obsrec (output) - output Observing Log records
Output Paramters:
pfistart - start times of pfi dominant modes
pfistop - stop times of pfi dominant modes (1:1 with pfistart)
Note - if /rd_obs not set, pfistart/stop are aligned with
FMS and LMS (granularity of ratio test)
Calling Sequence:
pfi_dominant, t0, t1, pfi0, pfi1 [,cutoff=NN., ss=ss, $
fids=fids, ratios=ratios, $
/rd_obs, obsrec=obsrec]
Method:
look at ratio of SFR/SPR using the xbd (tape contents) files
Optionally refine via observing log and return obsserving log records
CALLS: ***
FMT_TIMER [1], FMT_TIMER [2], Rd_TapDir, UNIQ [1], UNIQ [2], UNIQ [3]
anytim2ints [1], anytim2ints [2], fid2ex [1], fid2ex [2], fmt_tim [1], fmt_tim [2]
get_utevent [1], get_utevent [2], gt_dp_mode [1], gt_dp_mode [2], gtab_comm
str_concat [1], str_concat [2], weekfiles
History:
5-nov-1995 S.L.Freeland
7-nov-1995 S.L.Freeland - add RD_OBS and OBSLOG keywords and function
read Common table / confirm for success subset
[Previous]
[Next]
Name:
pfi_ffi
Purpose:
decode the pfi_ffi byte and return type of image (pfi
or ffi), BLS status and the count from bits 0-3.
CALLING SEQUENCE:
pfi_ffi, [val=val, itype=itype, BLS=BLS, expNo=bits03]
Input/Keyword:
val is the pfi_ffi byte-type value.
Output/Keyowrd:
itype is the image type: 1 for PFI, 2 for FFI and
3 for Patrol image
BLS is BLS status: 0 for off, and 1 for on
bits03 is the decimal value of bits 0 thru 3
Note for:
PFI-- this is the strip number within the
observing region (OR). The image size
would be n64 columns by 64 lines.
FFI-- this is the ROI (Region Of Interest)
number. (Note there is a max of 4
ROIs.)
Notes:
Actual parameter bit positions are reversed from the
positions documented in the FILE CONTROL DOC. (i.e.
bit 7 is bit 0). Comments for each bit operation list
the actual bit #s for the mask and the (documented bit
#s).
CALLED BY:
GetStrpfiffi
HISTORY:
written 26-Mar-91 by GAL
updated for ver 1.2 of the "File Contol Doc" 31-July-91
12-Aug-91, changed all formal parameters to keywords.
[Previous]
[Next]
NAME:
pfi_loc
PURPOSE:
To overlay the PFI location on an image of the sun, or to display a line
plot of the sun and show the PFI location.
SAMPLE CALLING SEQUENCE:
pfi_loc, item1
pfi_loc, '1-jun-93 1:00'
pfi_loc, index(30), /hour
pfi_loc, index(30), pfi_struct=pfi_index
pfi_loc, pfi_struct=pfi_index
pfi_loc, index, '6-nov-93 12:00',color=128
OPTIONAL INPUTS:
item1 - The index for a full frame image which was just displayed
(in which case the box location is overlayed on the image)
(OR)
A string time and date (in which case a line plot is made)
item2 - The time for the PFI images to be displayed. This option is
only relevant when "item1" is the full frame index.
The routine will correct for any commanded S/C offpoints, but
will not correct for solar rotation effects, so the delta time
between item1 and item2 should be less than a couple of days
to be very useful.
CALLS: ***
FMT_TIMER [1], FMT_TIMER [2], GET_NBYTES, anytim2ints [1], anytim2ints [2]
data_type [1], data_type [2], draw_boxcensiz [1], draw_boxcensiz [2]
fmt_tim [1], fmt_tim [2], gt_center, gt_corner, gt_res, gt_shape, plot_fov [1]
plot_fov [2], rd_obs, sxt_cen [1], sxt_cen [2], tim2dset [1], tim2dset [2]
CALLED BY:
STEPPER [3], STEPPER [4]
OPTIONAL KEYWORD INPUT:
bin - If the full frame image was rebinned, then the binning factor
must be passed in.
hours - If set to 1, then get the PFI locations +/- 1 one hour on
either side of the input time. If set to 2, then +/1 2 hours.
pfi_struct - The PFI location is determine by reading the observing
log by default. A SXT PFI roadmap or index can be passed
directly using pfi_struct.
color - color of box for pfi drawn on image, default=255
HISTORY:
Written 10-Jan-94 by M.Morrison
11-Jan-94 (MDM) - Added HOURS option
- Added second parameter option for the time of
the PFI(s) that should be displayed
21-Mar-94 (MDM) - Corrected an error which was causing PFI images
that were not FR to be displayed in the wrong
place.
- Corrected so that FFI images that are not full
frame are displayed properly.
9-Dec-94 (LWA) - Added color keyword.
[Previous]
[Next]
NAME:
PFI_TOTALIZER
PURPOSE:
quick and dirty total fluxes
CATEGORY:
Yohkoh
CALLING SEQUENCE:
pfi_totalizer,index,data,totals
INPUTS:
raw index, data
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
totals in each PFI frame, DN/pix/sec
CALLS: ***
gt_expdur [1], gt_expdur [2], sxt_decomp [1], sxt_decomp [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
HSH 6-dec-94
[Previous]
[Next]
Name:
ph2dn
Purpose:
Convert photon COUNTS into SXT DN.
Calling Sequence:
Dn = sxt_ph2dn(PH)
Dn = sxt_ph2dn(PH,Wave=Wave,gain=gain)
Inputs:
PH = Photon counts
Optional Input Keywords:
Gain = Sxt Gain constant (electrons/DN)
If not supplied, defaults to 100.
Wave = Photon wavelength (Ang). If not supplied, defaults to
10 ang.
Output:
DN = scalar or array of compressed (byte-type) DN values
or uncompressed (non-byte) DN values
CALLED BY:
UDN
Procedure:
Assumes that 3.65 eV produces 1 electron-hole pair
The equation is:
Photons = ( Wavelength (Ang) / 1.98648e-8 erg )
* 1.602192e-12 erg/eV * 3.65 eV/e * Gain (e/DN) * DN
History:
11 Apr, 1994 T.Shimizu (Univ.of Tokyo)
[Previous]
[Next]
NAME:
ph_eff_grs1
PURPOSE:
spectrum(efficiency) procedure for GRS-PHL1
CALLING SEQUENCE:
ph_eff_grs1,grs_phl1,grs_phh1,bitrate,date,doy,dp_mode,
tsec,isec,utime,utime2,utime_sec,nnn,pname
INPUTS:
(It will be written by Sato-san near future.)
CALLS: ***
ERRPLOT [1], ERRPLOT [2], dataget_grs, dtcf_ph_grs1l, header_scan, pprint [1]
pprint [2], window_select, wmenu_sel [1], wmenu_sel [2]
CALLED BY:
grs_plot
HISTORY:
22-Jun-93 (J.Sato) First written.
10-Jul-93 (K.Suga) First published.
[Previous]
[Next]
NAME:
ph_eff_grs2
PURPOSE:
spectrum(efficiency) procedure for GRS-PHL2
CALLING SEQUENCE:
ph_eff_grs2,grs_phl2,grs_phh2,bitrate,date,doy,dp_mode,
tsec,isec,utime,utime2,utime_sec,nnn,pname
INPUTS:
(It will be written by Sato-san near future.)
CALLS: ***
DTCF_PH_GRS2L, ERRPLOT [1], ERRPLOT [2], dataget_grs, header_scan, pprint [1]
pprint [2], window_select, wmenu_sel [1], wmenu_sel [2]
CALLED BY:
grs_plot
HISTORY:
22-Jun-93 (J.Sato) First written.
10-Jul-93 (K.Suga) First published.
[Previous]
[Next]
NAME:
ph_pow_grs1
PURPOSE:
Single power law spectrum procedure for GRS-PHL1 & GRS-PHH1
CALLING SEQUENCE:
ph_pow_grs1,grs_phl1,grs_phh1,bitrate,date,doy,dp_mode,
tsec,isec,utime,utime2,utime_sec,nnn,pname
INPUTS:
(It will be written by Sato-san near future.)
CALLS: ***
ERRPLOT [1], ERRPLOT [2], dataget_grs, dtcf_ph_grs1l, grs_lsm, grs_pow_fit
pprint [1], pprint [2], window_select, wmenu_sel [1], wmenu_sel [2]
HISTORY:
22-Jun-93 (J.Sato) First written.
10-Jul-93 (K.Suga) First published.
[Previous]
[Next]
NAME:
ph_pow_grs2
PURPOSE:
Single power law spectrum procedure for GRS-PHL2 & GRS-PHH2
CALLING SEQUENCE:
ph_pow_grs2,grs_phl2,grs_phh2,bitrate,date,doy,dp_mode,
tsec,isec,utime,utime2,utime_sec,nnn,pname
INPUTS:
(It will be written by Sato-san near future.)
CALLS: ***
DTCF_PH_GRS2L, ERRPLOT [1], ERRPLOT [2], dataget_grs, grs_lsm, grs_pow_fit
pprint [1], pprint [2], window_select, wmenu_sel [1], wmenu_sel [2]
HISTORY:
22-Jun-93 (J.Sato) First written.
10-Jul-93 (K.Suga) First published.
[Previous]
[Next]
NAME: photoland (landscape vers. of phototext.pro)
PURPOSE: output image and high qual. text to postscript printer.
Calling Seq:
land, img, r,g,b, [mag=mag, file=file, color=color, $
t1=t1, t2=t2, t3=t3, textoffset=textoffset, $
imgoffset=imgoffset]
INPUT:
img byte-array image
r,g,b color table vectors
OPTIONAL INPUT:
mag magnification factor for image 1, 1.5, 2, ...
file output postscript filename
color switch to specify color hardcopy
t1 1st text string to output
t2
t3
textoffset additional text offset in cm from
"bottom of page"
imgoffset additional image offset in cm from
centered position.
fillcolor name of rgb color for background fill.
CALLS: ***
CONGRID [1], CONGRID [2], CONGRID [3], PHOTOLAND, STRETCH, cfill, color2rgb
pprint [1], pprint [2]
HISTORY:
written by GAL Nov 93.
updated for color fill 1-Feb-94, gal
[Previous]
[Next]
NAME:
photoport (portrait vers. of phototext)
PURPOSE:
output image and high qual. text to postscript printer.
Calling Seq:
photoport, img, r,g,b, [mag=mag, file=file, color=color, $
t1=t1, t2=t2, t3=t3, textoffset=textoffset, $
imgoffset=imgoffset, fillcolor=fillcolor]
INPUT:
img byte-array image
r,g,b color table vectors
OPTIONAL INPUT:
mag magnification factor for image 1, 1.5, 2, ...
file output postscript filename
color switch to specify color hardcopy (default B&W)
t1 1st text string to output
t2
t3
textoffset additional text offset in cm from
"bottom of page"
imgoffset additional image offset in cm from
centered position.
fillcolor name of rgb color for background fill.
CALLS: ***
CONGRID [1], CONGRID [2], CONGRID [3], INTERPOL, PHOTOPORT, cfill, color2rgb
pprint [1], pprint [2]
HISTORY:
written by GAL Nov 93.
extensive updates 1-Feb-94, gal.
[Previous]
[Next]
NAME: PICK_CAN_PRELIM
PURPOSE: Reads strawman file and user choses passes to be deleted.
CALLING SEQUENCE: IDL> passes = pick_can_prelim(week)
INPUTS: week is the number of the week to be checked
OPTIONAL INPUTS: none
OUTPUTS: none
OPTIONAL OUTPUTS: none
CALLS: ***
Bell, CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], FILE_EXIST [2]
concat_dir [4], file_exist [1], file_exist [3], tbeep [1], tbeep [2], tbeep [3]
wmenu_sel [1], wmenu_sel [2]
CALLED BY:
CANCEL_PRELIM
RESTRICTIONS: Strawman file for the week must exist in
$DIR_GEN_ORBIT_RAW
PROCEDURE: Reads file, selects Yohkoh passes and presents them
as a menu to user who selects those to be deleted.
Those chosen are then writeen onto a JPL look-alike
form for faxing.
MODIFICATION HISTORY: Written: 2-Feb-1993 CDP
6-Feb-1993 JRL Read from $DIR_GEN_ORBIT_RAW
[Previous]
[Next]
NAME:
PICK_PIX
PURPOSE:
Get location of clicked pixel in IDL image coordinates.
CALLING SEQUENCE:
pick_pix, index
pick_pix, index, bin=2, side=9
INPUT:
data - 2-D image array already on screen
index - corresponding index
OPTIONAL KEYWORD PARAMETERS:
side - Number of pixels on a side (if not passed, the default is 5)
bin - The rebin size used to display the image
OUTPUTS:
FFI IDL image location of clicked pixel in summation mode of image.
CALLS: ***
draw_boxcorn [1], draw_boxcorn [2], get_pix_coor, gt_res
PROCEDURE:
Use get_pix_coord and index.sxt.corner_cmd to compute the position of
each image on the CCD. Add to these the values derived by clicking
on the desired feature in the image.
HISTORY:
Written by LWA 10/15/92 by copying pieces of codes by Metcalf and Morrison.
[Previous]
[Next]
NAME: PICKLAMBDA
PURPOSE:
To sum the lines at the same wavelengths
CATEGORY:
CALLING SEQUENCE:
picklambda, wvl, lflux, dummy, dumflx
INPUTS:
wvl = 1-d array of wavelengths of lines in the original Mewe file
lflux = 2-d array of line strengths with temperature (first
subscript) and wavelength (second subscript)
OPTIONAL INPUTS: none.
OUTPUTS:
dummy = 1-d array of new wavelengths after eliminating the
overlapping ones
dumflx = 2-d array of line strengths with temperature and
new wavelength (dummy)
OPTIONAL OUTPUTS: none.
CALLED BY:
linflx [1], linflx [2], xlinflx
COMMON BLOCKS: none.
SIDE EFFECTS: none.
RESTRICTIONS: none.
MODIFICATIONS: written by N.Nitta, April 1991.
J. R. Lemen, 25-jan-93 Renamed to picklambda
[Previous]
[Next]
NAME:
PINCH_PLASMA
PURPOSE:
calculates radio and X-ray fluxes of a pinched plasma
CATEGORY:
Yohkoh
CALLING SEQUENCE:
pinch_plasma, pinch, [t_in, v_in]
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
structure containing predicted radio and SXT X-ray fluxes
default initial temperature 3 MK
initial volume 6.2e27
initial density 1e9
** Structure <10c4f130>, 8 tags, length=136, refs=1:
SIGMA FLOAT 1.00000
EM FLOAT Array(3)
TEMP FLOAT Array(3)
VOLUME FLOAT Array(3)
DENSITY FLOAT Array(3)
RADIO FLOAT Array(3)
X_RAY FLOAT Array(3, 5)
CALLS: ***
sxt_flux [1], sxt_flux [2], sxt_flux [3]
COMMON BLOCKS:
SIDE EFFECTS:
WARNING:
RESTRICTIONS:
MODIFICATION HISTORY:
HSH written 5-Dec-94
[Previous]
[Next]
NAME:
PINHOLE_MAP
PURPOSE:
map out the pinhole brightness in thin Al terminators
CATEGORY:
CALLING SEQUENCE:
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
ROUTINES CALLED
FINDFILE, GT_FILTA, GT_FILTB, RD_ROADMAP, !stime
OUTPUTS:
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
HISTORY
5-Aug-97, written (HSH)
[Previous]
[Next]
NAME:
PINIT
PURPOSE:
Initialize the Peritek VCT graphics interface
(load timing, display, etc. parameters).
CATEGORY:
Image display through the Peritek VCT graphics interface.
CALLING SEQUENCE:
PINIT,RESULT=RESULT
INPUTS:
None.
KEYWORD PARAMETERS:
RESULT: Error message (optional)
OUTPUTS:
RESULT: Error message
CALLED BY:
go_lasdisk golaserdisk, go_lasdisk2 golaserdisk, video_test
PROCEDURE:
Spawns calls to two C routines which perform the necessary tasks.
MODIFICATION HISTORY:
10-Jul-1991 - Written by GLS
[Previous]
[Next]
NAME:
PIX2HEL
PURPOSE:
Converts an SXT pixel to heliographic coordinates.
CALLING SEQUENCE:
latlon = pix2hel(yy,xx)
latlon = pix2hel(yy,xx,suncenter=suncenter)
latlon = pix2hel(yy,xx,suncenter=suncenter,date=fmt_tim(index))
INPUTS:
yy = NS coordinate, full resolution pixel number (N positive)
xx = EW coordinate, full resolution pixel number (W positive)
suncenter may be obtained from hxa_suncenter. This routine
does not use the Yohkoh convention, so use [sunc(1),sunc(0)]
KEYWORDS (INPUT):
Suncenter = [center_y, center_x]. Sun center position, full pixel
units, W positiveJ; if not supplied, the default [520, 560]
will be used.
OUTPUTS:
Heliographic coordinates in decimal degrees, (N positive, W positive).
If the pixel coordinates are off the solar limb, the output is the
radial projection back to the nearest point on the limb
CONVENTIONS:
This follows the Yohkoh internal conventions:
Coordinate order is NN, WW.
Coordinates are positive in N and W directions.
CALLS: ***
conv_p2h [1], conv_p2h [2], dont_use_this [1], dont_use_this [2]
CALLED BY:
FOV_COORDS
SIDE EFFECTS:
WARNING:
Note above the treatment of points off the limb.
This has cumbersome interfaces because of the need to identify
the suncenter coordinates and date externally. Use a higher-level
routine if possible!
PROCEDURE:
Uses PB0R to calculate B0 and the solar radius.
MODIFICATION HISTORY:
HSH, October 1992 (This is just Wuelser's gbo_hel with pixel input
and the standard conventions for order of coordinates).
HSH improved the documentation, 1992 Dec. 19. Reversed the output
variables to match the convention somewhat better.
JPW Feb.93 added correction for tilt between SXT N and heliographic N.
MDM 16-Oct-93 - Removed from use - replaced code
MDM 19-Oct-93: Corrected replacement call to CONV_P2H
[Previous]
[Next]
NAME : pix_size (function)
PURPOSE :
return pix1 scale (arc sec.) for the observational data
CATEGORY :
idl/lib/nkr
CALLING SEQUENCE :
pix1 = pix_size(tel)
INPUTS :
tel : telescope name, 'ft1', 'ft3', 'c10', 'sxt', 'c25b'
OPTIONAL INPUT PARAMETERS :
none
KEYWORD PARAMETERS : none
OUTPUTS : none
CALLS: ***
PIX_SIZE
COMMON BLOCKS : none
SIDE EFFECTS : none
RESTRICTIONS : none
PROCEDURE :
MODIFICATION HISTORY :
K.I. '93/01/22
K.I. '93/05/06 append monochro
K.I. '93/09/03 append c25b
[Previous]
[Next]
NAME:
PIXEL_LC
PURPOSE:
dump out a time series for an interactively determined pixel
of a given data cube. Calls STEPPER to present an image to
interact with. Plots the results. Does no normalization or
background subtraction.
CALLING SEQUENCE:
lc = pixel_lc(index, data, pixq=pixq, same_pix=same_pix, quiet=quiet)
INPUT PARAMETERS:
data and index
OUTPUT PARAMETERS:
time series, and the pixq
KEYWORDS:
pixq= the 2d subscripts of the pixel chosen
quiet= does not ask for an ok on each pixel
SPECIAL TRICKS:
1) there is an invisible Q entry required to stop STEPPER
2) double click on center button, to exit
CALLS: ***
STEPPER [1], STEPPER [2], STEPPER [3], STEPPER [4], STEPPER [5], STEPPER [6]
STEPPER [7], UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
yesnox [1], yesnox [2]
HISTORY:
Written from old box_lc.pro, 8-22-94, jmm
[Previous]
[Next]
function pixon_map,image,func,nx,ny,npixons,top=top,resolution=resolution, $
logarithm=logarithm,imax=imax,smooth=smooth,grow=grow, $
sigmathresh=sigmathresh,noconserve=noconserve, $
guard_width=guard_width
NAME:
PIXON_MAP
PURPOSE:
Computes a pixon map for an input iage and GOF function
CATEGORY:
CALLING SEQUENCE:
pixon = pixon_map(image,func,nx,ny,npixons,top=top,resolution=resolution)
INPUTS:
image = image for which the pixon map is to be computed
func = name of function to compute GOF and gradient
nx,ny = size of image
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS
top = maximum scale
resolution = minimum scale
guard_width = width of guard ring surrounding image
/logarithm = image is the log of the true image
/smooth = smooth the pixon map by itself
/grow = grow the pixon map and then smooth the pixon map by itself
imax = when /logarithm is set, imax gives the maximum value of the log
sigmathresh = controls the noise threshold, default = 3.0
/noconserve = used to select smoothing algorithm in local_smooth.pro
OUTPUTS:
pixon = pixon map (local scale of image)
npixons = the number of pixons in the image.
CALLS: ***
DERIV, HISTERROR, LOCAL_SMOOTH, STDEV
CALLED BY:
HXTFPB, PSFPIXON
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
Assumes that there are regions with both the minimum and the maximum
scales.
PROCEDURE:
Computes the scale by oversmoothing the image and checking the gradient
to see where the smoothing is bad.
MODIFICATION HISTORY:
T. Metcalf 1994-06-27
1994-09-23 TRM Complete rewrite. The pixon map is now computed
pixels by pixel, without assuming anything about the
range of scales present in the image.
1994-10-26 TRM The pixon map is now computed using the derivative
od d(chi^2) with resolution.
[Previous]
[Next]
NAME:
PLATEAU
PURPOSE:
This function is used by CURVEFIT to fit this function and
a constant background to the data.
CATEGORY:
E2 - SURFACE FITTING.
CALLING SEQUENCE:
PLATEAU,X,A,F,PDER
INPUTS:
XY = 2-d array of values of independent variables = XY(N,2).
The second argument is 0 for values of X and 1 for values of Y.
A = PARAMETERS OF EQUATION DESCRIBED BELOW.
OUTPUTS:
F = VALUE OF FUNCTION AT EACH XY(I,K).
OPTIONAL OUTPUT PARAMETERS:
PDER = (N_ELEMENTS(X),N_ELEMENTS(Y),6) ARRAY CONTAINING THE
PARTIAL DERIVATIVES. P(I,J,K) = DERIVATIVE
AT ITH and JTH POINTS W/RESPECT TO KTH PARAMETER.
COMMON BLOCKS:
NONE.
SIDE EFFECTS:
NONE.
RESTRICTIONS:
NONE.
PROCEDURE:
X = XY(i,0) & Y = XY(i,1)
MODIFICATION HISTORY:
Written by R. Fuller, August 1992
Modified by R. Fuller, March 1993 to take 2-d array of XY
[Previous]
[Next]
NAME:
PLCT
PURPOSE:
Load the Peritek VCT graphics interface color translation
tables using the specified vectors, R,B,G.
CATEGORY:
Image display.
CALLING SEQUENCE:
VCTQ_TVLCT,R,G,B,RESULT=RESULT
INPUTS:
R,G,B: Integer arrays containing the color translation tables
in the RGB color system
KEYWORD PARAMETERS:
RESULT: Error message
OUTPUTS:
RESULT: Error message (optional)
CALLED BY:
PLOADCT, go_lasdisk golaserdisk, go_lasdisk2 golaserdisk, ptv_sxt, video_test
PROCEDURE:
Spawns call to C routine VCTQPSEUDO which loads the color tables
MODIFICATION HISTORY:
10-Jul-1991 - Written by GLS
15-nov-1991 - Updated, JRL -- Added DIR_PERITEK environment
variable to reference C routines.
30-apr-1991 - Updated, GLS -- Use !version to determine which
C routine (V bus version or Q bus
version) to call.
[Previous]
[Next]
NAME:
PLOADCT
PURPOSE:
Load Stern's predefined color tables into the Peritek
VCT graphics interface.
CATEGORY:
Image display.
CALLING SEQUENCE:
PLOADCT[,TABLE_NUMBER,RESULT=RESULT]
INPUTS:
TABLE_NUMBER: Number of table to load, from 0 to 15.
If omitted, a menu of the available tables
is printed and the user is prompted for a number.
KEYWORD PARAMETERS:
SILENT: If present and non-zero, supress the color table
message.
RESULT: Error message
OUTPUTS:
RESULT: Error message (optional)
CALLS: ***
FILEPATH, PLCT
COMMON BLOCKS:
Maintains the common block colors.
SIDE EFFECTS:
The color tables of the currently selected device are modified.
RESTRICTIONS:
Works from the file: $IDL_DIR/colors.tbl.
PROCEDURE:
The file colors.tbl is read. A maximum of 16 entries may
be defined. The titles of each table are contained in the first
16 * 32 characters/name bytes.
The colors loaded into the display are saved in the common
block COLORS, as both the current and original color vectors.
You can use the procedure modifyct to save your own color tables
in the file colors.tbl.
MODIFICATION HISTORY:
10-Jul-1991 - Written by GLS
29-apr-1992 - Updated, GLS -- Use !version to determine which
C routine (V bus version or Q bus
version) to call.
PLOADCT is a slightly modified version of LOADCT for use
exclusively with the Peritek VCT graphics display controller.
The only differences are:
1. The number of bits per color is hard-coded to eight.
2. Because the color table must be written to a temporary
file which is then read by the IDL-spawned C program,
PLOADCT calls a modified version of TVLCT (PLCT)
which writes the color table file and spawns the C
program.
[Previous]
[Next]
TITLE: plot2offset
PURPOSE:
Program to scale klik'd plot coordinates into Yohkoh offset in
degrees from present pointing.
Units: Degree changes in Yohkoh pointing.
+ means point Yohkoh E or S.
- means point Yohkoh W or N
CALLLING SEQUENCE:
plot2offset,cen [,out,blokcmd=blokcmd,yohkoh=yohkoh]
CALLED BY:
term_score2
RESTRICTIONS:
Only valid for 980x980 window and plot made with the following
idl call:
IDL> plot,[x0,x0],[y0,y0],yr=[0,1023],ystyle=1,xr=[0,1023],xstyle=1,$
title='Sun Center on CCD (92-Nov to Present)',$
xtit='CCD FR pixel',ytit='CCD FR pixel',charsize=1.4,$
psym=2,symsize=1.
where x0,y0 can be anything you like, e.g., suncenter of current
normal pointing.
INPUT:
cen, the suncenter position for present normal pointing.
CALLS: ***
POINT_HEX2 [1], POINT_HEX2 [2], WDEF [1], WDEF [2], klik
OPTIONAL KEYWORD INPUT:
/blokcmd prints out block commands for desired pointing
yohkoh, Euler angles of Yohkoh boresight.
default: yohkoh=[-0.0797222,0.0302778]
[S/C north , S/C east]
OPTIONAL OUTPUT:
out, array of Yohkoh Euler angles for offpoints.
PROGRAMS CALLED:
klik.pro, point_hex2.pro
HISTORY:
LWA 2/15/95
LWA 3/8/95, updated header.
Sometime in March '95, Weber added plot capability.
22-apr-95, JRL, If !d.window is already 980x980, don't create a new one.
Don't bomb out if no pointings are selected.
LWA 5/13/95, Changed plot3offset name back to plot2offset.
[Previous]
[Next]
TITLE: plot2offset
PURPOSE:
Program to scale klik'd plot coordinates into Yohkoh offset in
degrees from present pointing. Plots pointings.
Units: Degree changes in Yohkoh pointing.
+ means point Yohkoh E or S.
- means point Yohkoh W or N
Orbital pointing offset bias of terminators removed.
CALLLING SEQUENCE:
plot2offset,cen [,out,blokcmd=blokcmd,yohkoh=yohkoh,sft=sft,sxtf=sxtf, $
noplot=noplot,noklix=noklix,orb_bias=orb_bias,blowup=blowup]
plot2offset,cen,out,/blokcmd,/sft,/sxtf,blowup=[534,544,560,600]
plot2offset,cen,/blokcmd,/noplot
INPUT:
cen, the suncenter position for present normal pointing.
sxtf=sxtf causes plotting of sfr (or spr flare_mode) pointings.
sft=sft causes plotting or overplotting of SFT pointings.
CALLS: ***
POINT_HEX2 [1], POINT_HEX2 [2], WDEF [1], WDEF [2], klik
CALLED BY:
term_score2
OPTIONAL KEYWORD INPUT:
/blokcmd prints out block commands for desired pointing
yohkoh, Euler angles of Yohkoh boresight.
default: yohkoh=[-0.0797222,0.0302778]
[S/C north , S/C east]
noplot, /noplot assumes data already displayed, used for getting
new Yohkoh pointings from existing plot.
noklix, /noklix bypasses whole klik pointing selection process
orb_bias, offset in FR pixels of terminator offset from avg pointing.
default: orb_bias=[0.277.4.921]
blowup, array to specify expanded region, blowup=[x0,x1,y0,y1]
default: blowup=[0,1023,0,1023]
OPTIONAL OUTPUT:
out, array of Yohkoh Euler angles for offpoints.
PROGRAMS CALLED:
klik.pro, point_hex2.pro
RESTRICTIONS:
Only valid for 980x980 window and plot made with the following
idl call:
IDL> plot,[x0,x0],[y0,y0],yr=[0,1023],ystyle=1,xr=[0,1023],xstyle=1,$
title='Sun Center on CCD',$
xtit='CCD FR pixel',ytit='CCD FR pixel',charsize=1.4,$
psym=2,symsize=1.
where x0,y0 can be anything you like, e.g., suncenter of current
normal pointing.
HISTORY:
LWA 2/15/95
LWA 3/8/95, updated header.
Sometime in March '95, Weber added plot capability.
22-apr-95, JRL, If !d.window is already 980x980, don't create a new one.
Don't bomb out if no pointings are selected.
LWA 5/13/95 Changed plot3offset name back to plot2offset.
LWA 5/25/95 Separated plotting of SXTF and SFT points.
Added /noplot option and cleaned up header.
HSH 6/5/95 Added /name in call to point_hex2
LWA 11/11/95 Removed '(92-Nov to Present)' from title of plot.
Added orb_bias and blowup keywords.
Changed window call to use free_win.
Corrected error in scaling from window.
[Previous]
[Next]
NAME: PLOT_ALL
PURPOSE: plot and overplot 2D array
INPUT:
CALLING SEQUENCE: plot_all, [xvec,] data
KEYWORD PARAMETERS:
Normal - switch, true to normalize all plots
Labels - string array for line labels
Log - switch for Log Y axis
Pad - y padding as percentage of yrange
Title - main plot title (ignored if OPLOT set)
Which - integer vector of n elements for desired vectors
0:no plot, 1:plot
SIDE EFFECTS: plot to currently selected graphics device
HISTORY:
Written by S.L.Freeland, 5/21/91
[Previous]
[Next]
NAME: PLOT_ALL
PURPOSE: plot and overplot 2D array
INPUT:
CALLING SEQUENCE: plot_all, [xvec,] data
KEYWORD PARAMETERS:
Normal - switch, true to normalize all plots
Labels - string array for line labels
Log - switch for Log Y axis
Pad - y padding as percentage of yrange
Title - main plot title (ignored if OPLOT set)
Which - integer vector of n elements for desired vectors
0:no plot, 1:plot
CALLS: ***
OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], UTPLOT_ALL
SIDE EFFECTS: plot to currently selected graphics device
HISTORY:
Written by S.L.Freeland, 5/21/91
22-Nov-1993 (SLF) convert to utplot/outplot
[Previous]
[Next]
NAME:
plot_ar_pfi
PURPOSE:
To plot positions of NOAA Active Regions on a heliocentric
grid. Also plot PFI target regions over it.
CALLING SEQUENCE:
plot_ar_pfi, '15-Jun-92', pfi_index=roadmap
INPUT:
ref_item - The reference time to which all feature coordinates
will be differentially rotated before plotting on the grid.
Form can be:
(1) structure with a .time and .day field,
(2) standard 7-element external representation,
or (3) a string of the format "hh:mm dd-mmm-yy"
If no date is entered, the current date is used.
The year is not required, but if entered should be
of the form 1-Oct-91 style. The date should be entered
in string style with date first.
OPTIONAL INPUT:
st_item - The start time to search NOAA active region list for
regions. If not passed, only the time entries nearest to
REF_ITEM are plotted.
en_item - The end time to search NOAA active region list for
regions. If not passed, the NOAA data base is searched for
entries between ST_ITEM and REF_ITEM.
CALLS: ***
ADDTIME [1], ADDTIME [2], DIFF_ROT [1], DIFF_ROT [2], DRAW_GRID [1], DRAW_GRID [2]
DRAW_GRID [3], DR_PHOTO [1], DR_PHOTO [2], DR_PHOTO [3], RD_NAR, S2C [1], S2C [2]
SXT_GRID [1], SXT_GRID [2], SXT_GRID [3], WDEF [1], WDEF [2], anytim2ex [1]
anytim2ex [2], clearplot [1], clearplot [2], fmt_tim [1], fmt_tim [2]
font_size [1], font_size [2], get_pnt, get_rb0p [1], get_rb0p [2], gt_center, gt_res
gt_shape, int2secarr [1], int2secarr [2], pprint [1], pprint [2], sxt_cen [1]
sxt_cen [2], uniqo [1], uniqo [2]
OPTIONAL KEYWORD INPUT:
hc - If set, write the results to a temporary file, print
it and delete it.
OUTPUT:
HISTORY:
Written 26-Mar-93 by G. L. Slater
10-may-93 (slf) added pfi_color, pfi_lstyle keywords
1-Jun-93 (MDM) Added correction for the placement of the SXT
PFI window.
21-jun-93 (gal) quick bug fix for case where frontpts are
undefined [e.g. frontpts=where(x_coord ge 0)].
14-Jul-93 (MDM) Added draw_grid common block
16-Oct-93 (MDM) Removed call to GET_SUNCENTER and replaced with
a call to SXT_CEN
Revised 26-Feb-94 (TS) Plot PFI target regions (from roadmap) over it.
QUESTIONS, PROBLEMS, THINGS TO DO:
- Properly position window at creation, using screen size
[Previous]
[Next]
NAME:
plot_arc
PURPOSE:
To allow a user to mark an arc on an image and the intensity
along that arc will be plotted
SAMPLE CALLING SEQUENCE:
plot_arc, image
plot_arc, image, width, bin=bin, /sample
plot_arc, image, index, /enhance, /sample
plot_arc, data(*,*,1), index(1), 10, bin=4, /use_spline, /enhance
INPUT:
image - The image for which the curve wants to be plotted
OPTIONAL INPUT:
index - The index structure for the image used.
This is only necessary if /enhance is used.
width - The width in original pixels of the curve.
(ie: the number of lines to add together)
This is the width of the channel for which to get
the light curve.
CALLS: ***
DEFROI [1], DEFROI [2], DERIV, LWA_ARC, PAUSE [1], POLY, POLY_FIT, SPLINE, TVPLOT
WDEF [1], WDEF [2], deriv_arr [1], deriv_arr [2], enhancer, pause [2]
OPTIONAL KEYWORD INPUT:
enhance - Use EHNANCER to display data for increased visibity
of loops
sample - Use the /sample function of rebin to avoid smoothing.
bin - The rebinning factor.
If not passed, it will be set to whatever it takes
to make the image be 512 in the x direction.
use_spline - If set, then use the SPLINE method for smoothing
the curve, instead of the polynomial fit.
ndeg - The polynomial degree to be used for the fits. If
not passed, it will use 4th degree polynomial.
OPTIONAL KEYWORD OUTPUT:
lcur - The light curve array
METHOD:
DEFROI is used to mark the points on the image to make the
arc. The default is to make a polynomial fit to the x and
y positions separately (as a function of the distance down
the arc).
HISTORY:
Written 21-Oct-93 by M.Morrison
Enhancer and sample keywords added, LWA 14-Nov-93
[Previous]
[Next]
NAME:
plot_att_cmd
PURPOSE:
Plot the difference between att and cmd pointings.
CALLING SEQUENCE:
.run plot_att_cmd
HISTORY:
24-apr-95, J. R. Lemen, Written
V2.0 2-May-95 (MDM) - Modified to use the times of the SFM images
rather than equally spaced times
- Added options to overplot the correction
function
- Added option to plot the residuals
- Restore file /ys/gen/status/plot_att_cmd.idl
for a quick run
V2.1 2-May-95 (MDM) - Corrected sign error
V3.0 5-May-95 (MDM) - Added option to make the SHOW_PIX plot page
V3.1 6-jun-95 (SLF) - move restore file to site / dont allow restore if
not available
[Previous]
[Next]
NAME:
plot_bar_lines
PURPOSE:
Produce a line plot using histogram mode and filling
in the bins with vertical bars.
CALLING SEQUENCE:
plot_bar_lines,x,y
This routine takes all allowed plot proceedure keywords
HISTORY:
11-Jan-96, J. R. Lemen and Gan Wei Qun, written
[Previous]
[Next]
PLOT_BCS_FREE
Plots current disk usage as a function of file deletion as a function of
maximum counts in CaXIX
cdp Aug 92
CALLS:
[Previous]
[Next]
NAME:
PLOT_BSC
PURPOSE:
plot a BSC spectrum
CALLING SEQUENCE:
PLOT_BSC,BSC_INDEX,BSC_DATA
PLOT_BSC,BSC_INDEX,BSC_DATA,/NOWIDGET
INPUTS:
BSC_INDEX - BSC index structures
BSC_DATA - BSC data structure
KEYWORDS:
NOWIDGET - Don't use the widget environment
CALLS: ***
BCS_SPEC_PLOT, BSC_CHECK [1], BSC_CHECK [2], CHKARG, GT_BSC_CHAN [1]
GT_BSC_CHAN [2], GT_BSC_CRATE, OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], SEL_BSC
UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], WBSC [1]
WBSC [2], clearplot [1], clearplot [2], input [1], input [2], int2secarr [1]
int2secarr [2]
HISTORY:
Nov'92 - written by D. Zarro (ARC)
Sept'93 - added call to BCS_SPEC_PLOT for single spectrum case
30-sep-93, JRL, Check on ss= for call to BCS_SPEC_PLOT
15-oct-93, DMZ, Added ebar keyword
30-Nov-93, DMZ, added /SECOND keyword to oplot 2nd cmpt
9-Dec-93, DMZ, added VUNIT switch
13-Dec-93, DMZ, added /BLUE switch
(same effect as /second, and consistent with FIT_BSC and CAL_BSC)
24-Jan-94, DMZ, added CLEARPLOT and /PRIMARY keyword
26-Mar-94, DMZ, converted to keyword inheritance (IDL ver 3.5)
9-Jun-94, DMZ, fixed potential bug in keyword inheritance
8-Aug-94, DMZ, added call to CLEANPLOT
25-Aug-94, DMZ, fixed OUTPLOT bug
28-Sep-94, DMZ, fixed SS bug
22-Oct-94, DMZ, fixed /BLUE bug
Apr'95 , DMZ, added single spectrum channel check
16-Nov-95, DMZ, modified lightcurve time axis to use mean time
[Previous]
[Next]
NAME:
PLOT_BSC_AS2
PURPOSE:
plots BSC spectrum for use in overplot_fit.pro code.
CALLING SEQUENCE:
PLOT_BSC,BSC_INDEX,BSC_DATA,W,Z,X_LO,X_HI,/BIN,/HARD,CHAN=CHAN
INPUTS:
BSC_INDEX - BSC index structures
BSC_DATA - BSC data structure
X_LO, X_HI - lower and upper wavelength limits for plot.
KEYWORDS:
BIN - plot in bin units
OPTIONAL OUTPUTS:
W - bin or wavelength array
Z - counts/bin/sec values of plotted spectrum
CALLS: ***
BSC_CHAN_AS, BSC_CHECK [1], BSC_CHECK [2], PLOTSPEC, UTPLOT [1], UTPLOT [2]
UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], anytim2ints [1], anytim2ints [2]
fmt_tim [1], fmt_tim [2]
PROCEDURE:
Calls PLOTSPEC
HISTORY:
Modified Zarro's code! (4 Jan 93, A. Sterling)
plot_bsc_as2 version used as call from overplot_fit code. Subtracts
backgnd and changes some plotting parameters. (ACS 4 Mar 93.)
Changed bsc_chan to bsc_chan_as, as Zarro has changed his bsc_chan
version (ACS 25 Oct 93)
Written Nov'92 by D. Zarro (ARC)
[Previous]
[Next]
NAME:
PLOT_CINEMA2
PURPOSE:
To plot HXT image contours either on screens or on papers.
CALLING SEQUENCE:
plot_cinema2,images,nlevels=8
plot_cinema2,images,level=[40,60,80]
INPUT:
images - HXT image array read-in with RD_HXTIMG
CALLS: ***
SQRT2_LEV
CALLED BY:
DISP_HXTIMG
OPTIONAL KEYWORD INPUT:
nlevels - Number of contour levels displayed.
level - If set, then contour levels are set to this level (in %).
Default contour levels are:
[12.5, 17.7, 25.0, 25.4, 50.0, 70.7] %
OUTPUT:
none
OPTIONAL KEYWORD OUTPUT:
none
RESTRICTIONS:
none
SIDE EFFECTS:
Contours are displayed on your screen (default).
NOTICE:
This procedure may not be used directly; it is mainly called from
DISP_HXTIMG.
HISTORY:
version 1.0 (sometime) T.Sakao written
version 1.1 94.07.29 (Fri) @ MSSL
Options level and nlevels added.
version 1.2 98.05.03 (Sun) @ NSRO
Filename output location adjusted. This version is for the
HXT NMP release.
[Previous]
[Next]
NAME:
PLOT_CINEMA_HDR
PURPOSE:
To display header part of HXT images files.
CALLED BY:
DISP_HXTIMG
HISTORY:
version 1.0 date unknown T.Sakao written.
1.1 98.05.03 (Sun)
Y-location of the header for PS outputs modified so as
to match with the change made in PLOT_CINEMA2 ver.1.2.
This version is for the HXT NMP release.
[Previous]
[Next]
NAME:
PLOT_CLON
EXAMPLES:
PURPOSE:
Plot specified lines of Carrigton longitude
Plot a box bounded specified lines of Carrigton longitude and
heliographic latitude
Plot a Carrigton or heliographic grid
Interactively display Carrington coordinates of cursor
CATEGORY:
Image display.
CALLING SEQUENCE:
plot_clon,item,data,clon_vec=clon_vec,box_vec=box_vec, $
gridsep=gridsep,noerase=noerase,read_out=read_out, $
gridtype=gridtype,fov=fov, $
charsize=charsize,charthick=charthick,thick=thick,color=color, $
img_pos=img_pos,draw_limb=draw_limb,last_latlon=last_latlon
EXAMPLE CALLS:
plot_clon,index(0),box=[340,20,-30,30]
plot_clon,'11-jan-94',box=[24.1,24.2,15.001,15.00001],/read,gridt=2
plot_clon,index(5),clon=[130,150,180],gridt=1
REQUIRED INPUTS:
ITEM SXT index record or time in any Yohkoh format
OPTIONAL INPUTS:
POSITIONAL PARAMETERS:
DATA Image array. If passed it is displayed.
KEYWORDS PARAMETERS:
CLON_VEC Vector of Carrington longitudes to plot
BOX_VEC Vector of extrema in longitude and latitude
for box to be plotted. Form should be:
[lon_min,lon_max,lat_min,lat_max]
GRID_SEP Latitude and longitude grid separation interval
in degrees (def = 15).
NOERASE If present and nonzero, window is not erased prior
to drawing plotting (def is 1).
READ_OUT If present and non-zero, the Carrington coordinates
of the cursor will be written to a window-let.
GRIDTYPE 0 - No grid is drawn
1 - Heliocentric grid is drawn
2 - Carrington grid is drawn (def is 0).
DRAW_LIMB If present and non-zero, limb is drawn.
ULYSSES If present and non-zero, plot the extrpolated ulysse
position at the image time.
OUTPUTS:
POSITIONAL PARAMETERS:
none.
KEYWORDS PARAMETERS:
LAST_LATLON last output [during /READ_OUT mode]
[e/w, n/s] with east and south negative.
CALLS: ***
C2S, DRAW_GRID [1], DRAW_GRID [2], DRAW_GRID [3], MAKVEC, S2C [1], S2C [2], TIM2CLON
TIM2UPOS, WDEF [1], WDEF [2], anytim2ints [1], anytim2ints [2], conv_c2h
get_rb0p [1], get_rb0p [2], gt_corner, gt_res, gt_shape, restsys [1], restsys [2]
savesys [1], savesys [2], sxt_cen [1], sxt_cen [2]
COMMON BLOCKS:
None.
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
21-Jan-94 - GLS, based on SXT_GRID
18-Mar-94 - GLS, added ULYSSES keyword
TODO - Change IF to WHILE in longitude range conversions
[Previous]
[Next]
NAME:
PLOT_DIFF
PURPOSE:
Obtain the differential distribution of the "count" and
plot 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
[Previous]
[Next]
Procedure to plot dead time correction factors
Input:
roadmap - usual roadmap structure from test_rd
dt_fac - output array from DT_FACTOR
stime - start time hrs UT
etime - end time hrs UT
infil - file name for documentation
Use:
IDL> plot_dtfac,roadmap,dt_fac,stime,etime,infil
CDP March 92
CALLS:
[Previous]
[Next]
NAME:
PLOT_ERATE
PURPOSE:
Plots BCS all_cnts, lim_cnts and acc_cnts for a selected channel
and time interval
CATEGORY:
CALLING SEQUENCE:
Use YODAT to get dp_sync data
.run plot_erate
INPUTS:
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB Oct '91 Written
RDB 11-Jun-93 Made Channels 1-4
RDB 14-Jun-93 Changed to use UTPLOT
RDB 18-Jun-93 Put reset of pmulti and clearplot, etc at end
RDB 08-Jul-93 Added call to fmt_timer
[Previous]
[Next]
NAME:
plot_evn
PURPOSE:
To plot the data held in the EVN database
CALLING SEQUENCE:
plot_evn, '1-apr-92', '2-apr-92'
INPUT:
sttim - The starting date/time to plot
entim - The starting date/time to plot
CALLS: ***
RD_EVN, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ints [1], anytim2ints [2], gt_dp_mode [1], gt_dp_mode [2], int2secarr [1]
int2secarr [2], rd_fem
HISTORY:
Written 6-Aug-92 by M.Morrison
[Previous]
[Next]
NAME:
plot_fem
PURPOSE:
To plot the data held in the fem_data structure (s/c day/night,
saa, station contacts, ...)
CALLING SEQUENCE:
plot_fem, sttim
plot_fem, sttim, entim
plot_fem, '1-jun-93', '10-jun-93'
plot_fem, weeks, years
plot_fem, fem_data=fem_data
INPUT:
weeks - week(s) scalar or integer array
years - year(s) scalar or integer array
CALLS: ***
LPRINT, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ints [1], anytim2ints [2], data_type [1], data_type [2], int2secarr [1]
int2secarr [2], plot_lcur [1], plot_lcur [2], plottime [1], plottime [2], rd_fem
OPTIONAL KEYWORD INPUT:
fem_data- Pass in the FEM structure instead of reading it
fem_asca- Pass in the ASCA FEM structure instead of reading it
asca - If set, plot ASCA FEM information (default is Yohkoh)
both - If set, plot both ASCA and Yohkoh FEM data
zoom - If set, Allow interactive zooming
HISTORY:
Written 1-Mar-92 by M.Morrison
8-Apr-93 (MDM) - Added /ASCA and /BOTH options
10-Aug-93 (MDM) - Added /ZOOM option
23-Mar-95 (MDM) - Modifications to handle Wallops
2-May-95 (MDM) - Modified to trap on cases where the time
span is too short.
V3.2 15-May-95 (MDM) - Added Santiago
[Previous]
[Next]
NAME:
plot_fov
PURPOSE:
To plot the locations of the PFI fields of views (FOV).
The routine accepts either ROADMAP or INDEX.SXT
CALLING SEQUENCE:
plot_fov, index
plot_fov, roadmap, sel, ss
INPUT:
struct - Either the ROADMAP or INDEX structure
OPTIONAL OUTPUT:
selfov - If present, the user is prompted to mark a box
selecting the "window" of interest. The value
returned is an array of the same length as
the input array with the selected datasets
having a "1". The way to use this value is
something like the following:
ss = where(selfov)
ss = where(selfov and (struct.percentd eq 255))
ss - Return the list of indicies that were selected.
This should be used only if you do not want to
do other searches using the "where" function.
OPTIONAL INPUT:
box - If present, draw the size of the field of view.
The summation mode is reflected in the box size.
qstop - If present, stop at the end (for debugging)
img - If present, after plotting the line plot, display the
image and the location of the PFI images on the sun
It allows the user to check the offset being used.
nocmd - If set, use the actual pointing information from the
PNT file, not the S/C commanded value.
CALLS: ***
SETPS, font_size [1], font_size [2], pprint [1], pprint [2]
CALLED BY:
OBS_PLOT, chk_pointing, pfi_loc, sxt_uniq [1], xsxt_prep_event
HISTORY:
Written Oct-91 by M.Morrison
20-Oct-91 (MDM) - Added /box, /img, and /sun2 options
12-Nov-91 (MDM) - Added "ss" parameter option
16-Nov-91 (MDM) - Changed to use GT_FOV_CENTER so that index
can be passed directly
20-Nov-91 (MDM) - Put a temporary fix in because of GETSXTFOV
error in OR PFIs
28-Jan-93 (MDM) - Modified to use the expanded GT_FOV_CENTER
which does corrections done here before, and
corrects now for commanded pointing changes.
- Removed SUN2 option
14-Jul-93 (MDM) - Corrected IMG option
- Use GT_CENTER instead of GT_FOV_CENTER
- Allowed NOCMD option to use PNT pointing info
- Use GT_SHAPE instead of GT_SHAPE_CMD
- Replaced fixed solar radius with a call to GET_RB0P
21-Jul-93 (MDM) - Added /OBS_REGION switch to GT_SHAPE so that the
observing region shape is displayed
12-Sep-93 (SLF) - added skipfirst - if img and skipfirst, just do the
second plot (img overlayed) and dont prompt
[Previous]
[Next]
NAME:
plot_gbl
PURPOSE:
To plot the GRO/BATSE light curve data
SAMPLE CALLING SEQUENCE:
plot_gbl, '8-may-91 13:00', '8-may-91 13:15'
plot_gbl, event=1000
INPUT:
sttim - The start time to plot
entim - The end time to plot
CALLS: ***
RD_GBE, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ints [1], anytim2ints [2], rd_gbl
CALLED BY:
monthly_summary
OPTIONAL KEYWORD INPUT:
event - If set, plot the light curve for that event
HISTORY:
Written 16-Apr-93 by M.Morrison
6-May-93 (MDM) - Added YRANGE
27-Oct-94 (DMZ) - Added keyword inheritance
27-Oct-94 (SLF) - removed explicit yrange (ok via _extra)
[Previous]
[Next]
NAME:
PLOT_GRSPCH
PURPOSE:
To plot grs_pc13,14,15,16,23,24,25,16 with universal time label.
X axis range can be as large as 7 days.
CALLING SEQUENCE:
plot_grspch,index,data,pc=13,14,15,16,23,24,25 or 26
( after the run of yodat and the extract of 'data' and 'index' )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
pc - Integer indicates the energy band of pc data.
-- example.--
pc = 13 --> plot the grs_pc13.
pc = [ 13,14,15 ] --> plot the sum of grs_pc13,pc14 and pc15.
OPTIONAL INPUTS:
err - The magnitude of error bars ( in unit of sigma ) given
to the data.
-- example.--
err = 1 --> 1 sigma error bars.
err = 2 --> 2 sigma error bars.
sec - If set, count rate is 1 sec.
sum - The number of adjacent grs_pc data points to sum.
-- example.--
sum = 2 --> plot the sum of 2 adjacent points.
timerange - This gives a time range for plotting
It can be a string representation or structure.
-- example.--
timerange = [ '22:30:30 11-NOV-91','22:45 11-NOV-91' ]
timerange = [ '11-NOV-91','12-NOV-91' ]
timerange = [ index(10),index(20) ]
disp - If set, select the time range on display by mouse, and
save the time range in 'timerange' variable in 'timerange.dat'
file.
notitle - If set, no title is plotted.
PROCEDURE:
Time plotted is the middle of observation time interval.
The observation time is header time minus prestored time.
Data plotted are no deadtime-corrected.
CALLS: ***
ERRPLOT [1], ERRPLOT [2], FCHECK, Int2Ex [1], Int2Ex [2], MK_GRSPCH, UTPLOT [1]
UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], gt_day [1], gt_day [2]
gt_dp_rate [1], gt_dp_rate [2], gt_time [1], gt_time [2], wmenu_sel [1]
wmenu_sel [2]
RESTRICTIONS:
The optional input SEC is restricted to high and medium bit rates.
The optional inout ERR is restricted to high and medium bit rates.
HISTORY:
Written 22-jul-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 30-aug-93 by T.Hiraoka ; added disp option
Updated 3-dec-93 by T.Hiraoka ; used utplot
[Previous]
[Next]
NAME:
PLOT_GRSPCL
PURPOSE:
To plot grs_pc11,12,21,22 with universal time label.
X axis range can be as large as 7 days.
CALLING SEQUENCE:
plot_grspcl,index,data,pc=11,12,21 or 22
( after the run of yodat and the extract of 'data' and 'index' )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
pc - Integer indicates the energy band of pc data.
-- example.--
pc = 11 --> plot the grs_pc11.
pc = [ 11,12 ] --> plot the sum of grs_pc11 and grs_pc12.
OPTIONAL INPUTS:
err - The magnitude of error bars ( in unit of sigma ) given
to the data.
-- example.--
err = 1 --> 1 sigma error bars.
err = 2 --> 2 sigma error bars.
sec - If set, count rate is 1 sec.
sum - The number of adjacent grs_pc data points to sum.
-- example.--
sum = 2 --> plot the sum of 2 adjacent points.
timerange - This gives a time range for plotting.
It can be a string representation or structure.
-- example.--
timerange = [ '22:30:30 11-NOV-91','22:45 11-NOV-91' ]
timerange = [ '11-NOV-91','12-NOV-91' ]
timerange = [ index(10),index(20) ]
disp - If set, select the time range on display by mouse, and
save the time range in 'timerange' variable in 'timerange.dat'
file.
notitle - If set, no title is plotted.
PROCEDURE:
Time plotted is the middle of observation time interval.
The observation time is header time minus prestored time.
Data plotted are no deadtime-corrected.
CALLS: ***
ERRPLOT [1], ERRPLOT [2], FCHECK, Int2Ex [1], Int2Ex [2], MK_GRSPCL, UTPLOT [1]
UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], gt_day [1], gt_day [2]
gt_dp_rate [1], gt_dp_rate [2], gt_time [1], gt_time [2], wmenu_sel [1]
wmenu_sel [2]
RESTRICTIONS:
The optional input SEC is restricted to high and medium bit rates.
The optional inout ERR is restricted to high and medium bit rates.
HISTORY:
Written 19-jul-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 21-jul-93 by T.Hiraoka ; added timerange option
Updated 30-aug-93 by T.Hiraoka ; added disp option
Updated 3-dec-93 by T.Hiraoka ; used utplot
[Previous]
[Next]
NAME:
PLOT_GRSPHH
PURPOSE:
To plot grs_phh with universal time label.
X axis range can be as large as 7 days.
CALLING SEQUENCE:
plot_grsphh,index,data,ph=1,range=[40.0,56.0]
plot_grsphh,index,data,ph=2,channel=[8,12]
( after the run of yodat and extracting 'data' and 'index' )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
ph - Detector number of wbs grs_phh.
range - The energy range.(MeV) ( 8.00 - 104 )
-- example.--
range = [ 30.0,82.0 ] --> plot ths total of bitween each
energy.
channel - The range of channel number. ( 0 - 15 )
-- example.--
channel = 7 --> plot the channel 7 of grs_ph.
channel = [ 7,10 ] --> plot the total of channel 7 to 10.
OPTIONAL INPUTS:
err - The magnitude of error bars ( in unit of sigma ) given
to the data.
-- example.--
err = 1 --> 1 sigma error bars
err = 2 --> 2 sigma error bars
sum - The number of adjacent grs_phh data points to sum.
-- example.--
sum = 2 --> plot the sum of 2 adjacent points.
timerange - This can gives a time range for plotting.
It can be a string representation or structure.
-- example.--
timerange = [ '22:30:30 11-NOV-91','22:45 11-NOV-91' ]
timerange = [ '11-NOV-91','12-NOV-91' ]
timerange = [ index(10),index(20) ]
disp - If set, select the time range on display by mouse, and
save the time range in 'timerange' variable in
'timerange.dat' file.
notitle - If set, no title is plotted.
( title, xtitle, ytitle )
CALLS: ***
ERRPLOT [1], ERRPLOT [2], FCHECK, Int2Ex [1], Int2Ex [2], MK_GRSPHH1, MK_GRSPHH2
UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], gt_day [1]
gt_day [2], gt_dp_rate [1], gt_dp_rate [2], gt_time [1], gt_time [2], header_scan
wmenu_sel [1], wmenu_sel [2]
PROCEDURE:
Time plotted is the middle of observation time interval.
The observation time is header time minus prestored time.
Data plotted are deadtime-corrected.
PRESTRICTIONS:
The optional inout ERR is restricted to high and medium bit rates.
HISTORY:
Written 22-nov-93 by T.Hiraoka ( Rikkyo Univ.)
[Previous]
[Next]
NAME:
PLOT_GRSPHL
PURPOSE:
To plot grs_phl with universal time label.
X axis range can be as large as 7 days.
CALLING SEQUENCE:
plot_grsphl,index,data,ph=1,range=[2.22,6.23]
plot_grsphl,index,data,ph=2,channel=[8,12]
( after the run of yodat and the extract of 'data' and 'index' )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
ph - Detector number of wbs grs_phl.
range - The energy range.(MeV) ( 0.270 - 13.081 ) <- PH1
( 0.300 - 13.682 ) <- PH2
-- example.--
range = [ 4.23,6.45 ] --> plot ths total of bitween each
energy.
channel - The range of channel number. ( 0 - 127 )
-- example.--
channel = 15 --> plot the channel 15 of grs_phl.
channel = [ 15,20 ] --> plot the total of channel 15 to 20.
OPTIONAL INPUTS:
err - The magnitude of error bars ( in unit of sigma ) given
to the data.
-- example.--
err = 1 --> 1 sigma error bars
err = 2 --> 2 sigma error bars
sum - The number of adjacent grs_phl data points to sum.
-- example.--
sum = 2 --> plot the sum of 2 adjacent points.
timerange - This can gives a time range for plotting.
It can be a string representation or structure.
-- example.--
timerange = [ '22:30:30 11-NOV-91','22:45 11-NOV-91' ]
timerange = [ '11-NOV-91','12-NOV-91' ]
timerange = [ index(10),index(20) ]
disp - If set, select the time range on display by mouse, and
save the time range in 'timerange' variable in
'timerange.dat' file.
notitle - If set, no title is plotted.
( title, xtitle, ytitle )
CALLS: ***
ERRPLOT [1], ERRPLOT [2], FCHECK, Int2Ex [1], Int2Ex [2], MK_GRSPHL1, MK_GRSPHL2
UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], gt_day [1]
gt_day [2], gt_dp_rate [1], gt_dp_rate [2], gt_time [1], gt_time [2], header_scan
wmenu_sel [1], wmenu_sel [2]
PROCEDURE:
Time plotted is the middle of observation time interval.
The observation time is header time minus prestored time.
Data plotted are deadtime-corrected.
PRESTRICTIONS:
The optional inout ERR is restricted to high and medium bit rates.
HISTORY:
Written 22-nov-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 3-dec-93 by T.Hiraoka ; used utplot
[Previous]
[Next]
NAME:
plot_gsn
PURPOSE:
To plot the active region information in the SOON observing log files
SAMPLE CALLING SEQUENCE:
plot_gsn,'/ydb/gsn/sth940517.cat'
plot_gsn,'/ydb/gsn/sth940517.cat', /hc
plot_gsn,'/ydb/gsn/sth940517.cat',gif_file="soon_out.gif"
INPUT:
infil - the input file
CALLS: ***
LOADCT, MK_SQUARE2, Rd_fHeader [1], Rd_fHeader [2], Rd_fHeader [3], WDEF [1]
WDEF [2], WRITE_GIF, anytim2ints [1], anytim2ints [2], data_chk [1], data_chk [2]
draw_boxcensiz [1], draw_boxcensiz [2], fmt_tim [1], fmt_tim [2], plottime [1]
plottime [2], pprint [1], pprint [2], rd_gsn
CALLED BY:
soon_search [1], soon_search [3]
OPTIONAL KEYWORD INPUT:
hc - If set, send the output to the printer
gif_file - plot results to named gif file
HISTORY:
Written 28-Jul-94 by M.Morrison
5-Aug-94 (MDM) - Allow header to be passed in
14-jul-95 (SLF) - add gif_file keyword
[Previous]
[Next]
NAME:
plot_histo_j
PURPOSE:
To plot histograms
CALLING SEQUENCE:
plot_histo_j, f, edges, mdpts=mdpts, xrange=xrange, yrange=yrange, $
title=title, xtitle=xtitle, ytitle=ytitle, $
psym=psym, linestyle=linestyle, oplot=oplot, $
charsize = charsize, color=color
INPUT:
f = the histogram
edges= the bin edges for the histogram
KEYWORDS:
mdpts= the bin mdpts, if not set, they'll be calculated
xrange= range of x values
yrange= range of y values
title= a title for the plot
xtitle= title for the x axis
ytitle= title for the y axis
psym= psym for the plot, default is 10, histogram mode,
psym=-10 will put dashed lines on the bin edges.
linestyle= linestyle for the plot, default is 0, a solid line
oplot= if set, will overplot
charsize = the size for the characters in the titles and axis labels, 1.0 is the default
color= color for the plot
CALLED BY:
avsig_test, dist_histo
HISTORY:
20-jun-94 by JMM, dissatisfied with psym=10
13-may-1997, jmm, changed name to avoid conflict with the
not-particularly useful PLOT_HISTO in /ssw
[Previous]
[Next]
NAME:
plot_hxa
PURPOSE:
Plot the HXA scan data (loops through all data
with 2 second scan between plots)
INPUT:
hxa_scan- The hxa scan data structure
CALLS: ***
fmt_tim [1], fmt_tim [2]
HISTORY:
Written Fall '91 by M.Morrison
[Previous]
[Next]
NAME:
PLOT_HXSPC
PURPOSE:
To plot hxs_pc1 and pc2 with universal time label.
X axis range can be as large as 7 days.
CALLING SEQUENCE:
plot_hxspc,index,data,pc=1 or 2
( after the run of yodat and the extract of 'data' and 'index' )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
pc - Integer indicates the energy band of pc data.
-- example.--
pc = 1 --> plot the hxs_pc1
pc = [ 1,2 ] --> plot the sum of hxs_pc1 and hxs_pc2
OPTIONAL INPUTS:
err - The magnitude of error bars ( in unit of sigma ) given
to the data.
-- example.--
err = 1 --> 1 sigma error bars
err = 2 --> 2 sigma error bars
sec - If set, count rate is 1 sec.
sum - The number of adjacent hxs_pc data points to sum.
-- example.--
sum = 2 --> plot the sum of 2 adjacent points.
timerange - This can gives a time range for plotting.
It can be a string representation or structure.
-- example.--
timerange = [ '22:30:30 11-NOV-91','22:45 11-NOV-91' ]
timerange = [ '11-NOV-91','12-NOV-91' ]
timerange = [ index(10),index(20) ]
disp - If set, select the time range on display by mouse, and
save the time range in 'timerange' variable in
'timerange.dat' file.
notitle - If set, no title is plotted.
( title, xtitle, ytitle )
PROCEDURE:
Time plotted is the middle of observation time interval.
The observation time is header time minus prestored time.
Data plotted are deadtime-corrected.
CALLS: ***
ERRPLOT [1], ERRPLOT [2], FCHECK, Int2Ex [1], Int2Ex [2], MK_HXSPC, UTPLOT [1]
UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], gt_day [1], gt_day [2]
gt_dp_rate [1], gt_dp_rate [2], gt_time [1], gt_time [2], wmenu_sel [1]
wmenu_sel [2]
RESTRICTIONS:
The optional input SEC is restricted to high and medium bit rates.
The optional inout ERR is restricted to high and medium bit rates.
HISTORY:
Written 23-jul-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 30-aug-93 by T.Hiraoka ; added disp option
Updated 3-dec-93 by T.Hiraoka ; used utplot
[Previous]
[Next]
NAME:
PLOT_HXSPH
PURPOSE:
To plot hxs_ph with universal time label.
X axis range can be as large as 7 days.
CALLING SEQUENCE:
plot_hxsph,index,data,range=[110.5,242.4]
plot_hxsph,index,data,channel=[8,12]
( after the run of yodat and the extract of 'data' and 'index' )
INPUTS:
index - Wbs index structure from reformat data.
data - Wbs data structure from reformat data.
range - The energy range.(keV) ( 19.81 - 657.5 ) <- OS 21
( 24.90 - 826.44 ) <- OS 01
-- example.--
range = [ 234.5,345.6 ] --> plot ths total of bitween each
energy.
channel - The range of channel number. ( 0 - 31 )
-- example.--
channel = 15 --> plot the channel 15 of hxs_ph.
channel = [ 15,20 ] --> plot the total of channel 15 to 20.
OPTIONAL INPUTS:
err - The magnitude of error bars ( in unit of sigma ) given
to the data.
-- example.--
err = 1 --> 1 sigma error bars
err = 2 --> 2 sigma error bars
sum - The number of adjacent hxs_ph data points to sum.
-- example.--
sum = 2 --> plot the sum of 2 adjacent points.
timerange - This can gives a time range for plotting.
It can be a string representation or structure.
-- example.--
timerange = [ '22:30:30 11-NOV-91','22:45 11-NOV-91' ]
timerange = [ '11-NOV-91','12-NOV-91' ]
timerange = [ index(10),index(20) ]
disp - If set, select the time range on display by mouse, and
save the time range in 'timerange' variable in
'timerange.dat' file.
notitle - If set, no title is plotted.
( title, xtitle, ytitle )
CALLS: ***
ERRPLOT [1], ERRPLOT [2], FCHECK, Int2Ex [1], Int2Ex [2], MK_HXSPH, UTPLOT [1]
UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6], gt_day [1], gt_day [2]
gt_dp_rate [1], gt_dp_rate [2], gt_time [1], gt_time [2], header_scan
wmenu_sel [1], wmenu_sel [2]
PROCEDURE:
Time plotted is the middle of observation time interval.
The observation time is header time minus prestored time.
Data plotted are deadtime-corrected.
PRESTRICTIONS:
The optional inout ERR is restricted to high and medium bit rates.
HISTORY:
Written 22-nov-93 by T.Hiraoka ( Rikkyo Univ.)
Updated 3-dec-93 by T.Hiraoka ; used utplot
[Previous]
[Next]
pro plot_hxt, index0, data0, force_number=force_number, A4=A4, $
levels=level, independent=independent, ps=ps, fileps=psfile, $
noerase = noerase, title=title, sample=sample, image=image, $
rowcol=rowcol, legal=legal, nowords=nowords, nogrids=nogrids, $
noboxes=noboxes, nocontours=nocontours,limb=limb,latlon=latlon, $
bin=bin,nocloseps=nocloseps,colorbox=colorbox,colorcontour=color, $
labels=plotlabel,noscaleimg=noscaleimg,nlevels=nlevels,nolimb=nolimb
NAME:
PLOT_HXT
PURPOSE:
Plots a series of HXT images as contours.
CATEGORY:
CALLING SEQUENCE:
plot_hxt,index,data
INPUTS:
index = index with index.hxi and index.gen substructures
data = HXT image data cube
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS
rowcol = number of rows,columns to plot, [rows,cols]. rows*columns
must be GE to the number of images.
level = contour levels in fraction of max value.
default = [0.04,0.08,0.15,0.25,0.4,0.65,0.95]
nlevels = specify the number of equally spaced levels. Ignored if level is
set.
independent = Scale each image indepently. The default is to scale
the entire data cube so that relative intensities are
preserved from image to image.
/sample = use /sample in rebin
image = data cube of images to display under the contours. Must have
the same dimensions as data. However, if image is a 2-D array
with the same dimension as the images in data, the single image
is displayed under each data image.
/noscaleimg = Do no bytescale the underlying images. the default is
to display each image independently with tvscl. If
/noscaleimg is set, the images must be byte type and
will be displayed with tv rather than tvscl.
force_number = force this number of plots on the output
/limb = Plot a contour representing the solar limb. See restrictions,
below. This is now the default. Use /nolimb to turn the limb
off.
/latlon = Plot a lat/lon grid. See restrictions, below.
labels = String array giving a label for each image.
/nowords = do not annotate the plots
/nogrids = do not plot grids on images
/noboxes = do not plot boxes around images
/noerase = do not erase the display before plotting.
/nocontours = do not plot contours. This makes little sense unless you
set the image keyword.
/ps = Use PostScript. If the current graphics device is 'PS', the /ps
keyword is the default.
fileps = PostScript file, default is idl.ps
/nocloseps = Do not close the PS device after plotting. This is useful if you
want to add more annotations etc to the plot. the side effects are
that the device is left as PS after the program and, of course,
you will have to close the PS file yourself.
/legal = Use legal size paper size. Ignored except for PostScript output.
/A4 = Use A4 paper size. Ignored except for PostScript output.
title = title for plot
bin = binup factor for the images (def = 4)
OUTPUTS:
CALLS: ***
ADDTIME [1], ADDTIME [2], CONGRID [1], CONGRID [2], CONGRID [3], CONV_HXT2A
CONV_HXT2H, GET_LATLON, OCONTOUR2, PIXSQUARE, draw_box_nn [1], draw_box_nn [2]
get_rb0p [1], get_rb0p [2], gt_pix_size, gt_time [1], gt_time [2], tag_index [1]
tag_index [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
The /independent keyword only makes sense if the /nonorm keyword was
passed to hxt_multimg. If /nonorm was not used, the images are
already scaled independently.
PROCEDURE:
The data cube is byte scaled and plotted as contours. Use the
/independent keyword to scale each image independently
MODIFICATION HISTORY:
T. Metcalf 1995-06-29 Borrowed lots of code from N. Nitta's routines
1995-08-21 Added /limb and /latlon keywords
1996-11-13 Added bin keyword.
1997-01-24 Added nlevels keyword.
Modified so that the index need not be specified.
1998-03-25 Fixed bug in the lat/lon calculation.
1998-05-11 Made the limb contour smoother. Removed
restrictions on lat/lon which made the grid fail
if the FOV center was off the limb.
[Previous]
[Next]
NAME:
PLOT_INIT
PURPOSE:
CATEGORY:
CALLING SEQUENCE:
INPUTS:
POSITIONAL PARAMETERS:
KEYWORDS PARAMETERS:
OUTPUTS:
POSITIONAL PARAMETERS:
KEYWORDS PARAMETERS:
CALLED BY:
TV_SSX
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
EXAMPLE:
MODIFICATION HISTORY:
August, 1992. - Written by GLS, LMSC.
[Previous]
[Next]
NAME:
plot_int
PURPOSE
get one-dimensional distribution plot of some quantity
between two arbitrary points on a two-dimensional picture
CALLING SEQUENCE
.run plot_int.pro
INPUT
da : image data (two-dimensional array)
infss : info_array(n) corresponding to da, note that this is not array
OUTPUT
dint : one-dimensional distribution of da between two points
x1, y1 : coordinate of the first point specified by mouse
x2, y2 : coordinate of the second point specified by mouse
ix, iy : coordinate of grid points on the line between the two points
RESTRICTION
Please prepare da and infss before running this program.
Print-out command (lpr) is included in this version.
So if you do not want to print out figures, you should
delete such procedures in this program.
AUTHOR & HISTORY
written by K. Shibata, NAOJ
Ver 1.0 Apr 6, 1992
Ver 1.1 Apr 27, 1992
[Previous]
[Next]
NAME:
plot_lat_cmd
PURPOSE:
show an image of the mgram's white-light with overlay of
lat,long grid, and then calculate/ show spot positions
in units of Mm on the solar surface (heliographic coords)
CATEGORY:
data analysis
CALLING SEQUENCE:
plot_lat_cmd,mag,year,month,day,hour,bin=bin,image=image
INPUTS:
mag: magnetic field structure, *with pointing info!*
year: e.g. 1992
month: 0-12
day: 0-31
hour: decimal, eg 1.5 for 01:30
OPTIONAL INPUTS:
bin = if set, the factor to bin up for displaying the image
origin = if set, then the lat,cmd position of the input image
is at the bottom left corner - else assumed at the center
image = if named, will bytscl the named image under the grid
instead of the white-light image
KEYWORD PARAMETERS:
OUTPUTS:
CALLS: ***
POLREC3D, ROT_3D, STONY, SUN [1], SUN [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
KDL, written November 1992
[Previous]
[Next]
NAME:
plot_nar
PURPOSE:
To plot positions of NOAA Active Regions on a heliocentric
grid
CALLING SEQUENCE:
plot_nar, roadmap(i)
plot_nar, index(i)
plot_nar, '30-oct-91'
INPUT:
ref_item - The reference time to which all feature coordinates
will be differentially rotated before plotting on the grid.
Form can be:
(1) structure with a .time and .day field,
(2) standard 7-element external representation,
or (3) a string of the format "hh:mm dd-mmm-yy"
If no date is entered, the current date is used.
The year is not required, but if entered should be
of the form 1-Oct-91 style. The date should be entered
in string style with date first.
OPTIONAL INPUT:
st_item - The start time to search NOAA active region list for
regions. If not passed, only the time entries nearest to
REF_ITEM are plotted.
en_item - The end time to search NOAA active region list for
regions. If not passed, the NOAA data base is searched for
entries between ST_ITEM and REF_ITEM.
CALLS: ***
ADDTIME [1], ADDTIME [2], DIFF_ROT [1], DIFF_ROT [2], DRAW_GRID [1], DRAW_GRID [2]
DRAW_GRID [3], DR_PHOTO [1], DR_PHOTO [2], DR_PHOTO [3], RD_NAR, S2C [1], S2C [2]
SXT_GRID [1], SXT_GRID [2], SXT_GRID [3], WDEF [1], WDEF [2], anytim2ex [1]
anytim2ex [2], fmt_tim [1], fmt_tim [2], font_size [1], font_size [2]
get_rb0p [1], get_rb0p [2], gt_corner_cmd, gt_or_expnum, gt_res, gt_shape_cmd
int2secarr [1], int2secarr [2], pprint [1], pprint [2], restsys [1], restsys [2]
savesys [1], savesys [2], sxt_cen [1], sxt_cen [2], uniqo [1], uniqo [2]
CALLED BY:
fl_sxtobsnar [1], fl_sxtobsnar [2], lastsfd [1], lastsfd [2], lastspd
mk_sun_mosaic, xspr [1], xspr [2]
OPTIONAL KEYWORD INPUT:
hc - If set, write the results to a temporary file, print
it and delete it.
text - If set, list active region parameters
OUTPUT:
HISTORY:
Written 26-Mar-93 by G. L. Slater
10-may-93 (slf) added pfi_color, pfi_lstyle keywords
1-Jun-93 (MDM) Added correction for the placement of the SXT
PFI window.
21-jun-93 (gal) quick bug fix for case where frontpts are
undefined [e.g. frontpts=where(x_coord ge 0)].
14-Jul-93 (MDM) Added draw_grid common block
16-Oct-93 (MDM) Removed call to GET_SUNCENTER and replaced with
a call to SXT_CEN
3-Jul-94 (SLF) Fix bug in pfi resolution, only do 'uniq' pfis
since pfi code is slow
4-Jul-94 (SLF) use gt_corner_cmd
14-Jul-94 (SLF) Clean up, allow z-buffer, use savesys,restsys
18-Jul-94 (SLF) do uniq pfi outside this routine
29-Mar-96 (GLS) Added /nosavrest to DRAW_GRID call
Corrected wset problem
QUESTIONS, PROBLEMS, THINGS TO DO:
- Properly position window at creation, using screen size
[Previous]
[Next]
NAME: plot_new_respon
purpose: Plot and compare new and old response functions for SXT.
CALLS:
[Previous]
[Next]
NAME:
PLOT_POWERLOW
PURPOSE:
Obtain the slope of the power-law fitting curve.
INPUTS:
pos - X-axis data array
hist - Y-axis data array
OPTIONAL INPUTS:
fitrange = [low, high]
the energy range for the power-law fitting.
OUTPUTS:
CALLS: ***
PLOT_POWERLAW, POLY_FIT
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)
29 Dec, 93 Print Sigma value.
[Previous]
[Next]
atp 27-mar-1992
name: plot_ref.pro
purpose: plots and over plots a series of spectra for direct comparison
inputs: data data structure
index structure
channel - data channel you are interested in 0,1,2,3
dset_arr - array of dataset numbers - the first is plotted
the rest are overplotted.
optional inputs: title - the plot title
if not specified the channel and date from the
first index are used
decomp - do bcs decompression on the data
default is to not decompress
line - put a different line style on each plot
default is to use solid for each
scale - vector of same length as dset_arr containing
the scale factors associated with each index.
default is to not scale.
noextract - will not extract the data using the bcs_grp_plan
(useful if you have already extracted by channel)
(or are after a sub array )
This may allow use of bsd files in this program
bsd - use bsd file structures
CALLS: ***
BCS_DECOMP, BCS_GRP_PLAN, anytim2ex [1], anytim2ex [2], fmt_tim [1], fmt_tim [2]
gt_day [1], gt_day [2], gt_time [1], gt_time [2], intextime
side effects: a plot is produced on the current graphics device.
[Previous]
[Next]
NAME:
plot_sld
PURPOSE:
To plot the SXT dark current signal after the Nov-92 entrance filter
failure.
CALLING SEQUENCE:
.run plot_sld
METHOD:
Dark frame images are taken with the Be x-ray analysis filter.
After 17-Nov-92 entrance filter failure, it was noticed that the dark
current signal increase dramatically. This routine monitors that
increase (which decreased after a few months).
A full frame image is taken and broken up into 51 vertical strips
(20 columns per strip for FR, 10 for HR, and 5 for QR). The values
plotted in the first plot are strip 5 and strip 0. Strip 5 is the area
where the increase in the dark signal was seen. Strip 0 is off the
area where the sun hits the CCD and is a reference area.
The HR images are converted to the equivalent signal for a QR image.
Open /Al.1 Half Dark C 13 78.0 512x512
Open /AlMg Qrtr Dark C 13 78.0 256x256
HISTORY:
Written Winter, 1992 by M.Morrison
22-May-93 (MDM) - Added header information
V1.2 3-Sep-93 (MDM) - Added a new plot - plot each vertical strip offset
V1.3 30-Sep-93 (MDM) - Sped up the plotting steps
V2.0 30-Sep-93 (MDM) - Changed to plot in units of QR DN
V2.1 12-Oct-93 (MDM) - Modified to always make the image 800 line display
V2.2 18-Feb-94 (MDM) - Corrected labeling of the temperature plot
(it was labeled as TEC hot end and was really the
center of the metering tube).
[Previous]
[Next]
NAME:
plot_sls
PURPOSE:
To plot the SXT Scattered Light engineering data
CALLING SEQUENCE:
.run plot_sls
METHOD:
Dark frame images are taken with the narrow band analysis filters.
Since the leak after 17-Nov-92 entrance filter was intense enough
to show signal at the CCD even with the CCD closed, we use this
data to monitor the size of the aperature in the entrance filter.
The values plotted are the average signal for the whole FFI image.
The HR images are converted to the equivalent signal for a QR image.
The 70 DN "fat zero/dark current" of the QR images is still in the
values which are plotted.
There are two sets of values because of the difference in the
stopped position of the shutter blade. Two of the four positions
allow more light to leak around the blade.
NaBan/Open Half Dark C 2 1.0 512x512
NaBan/Open Qrtr Dark C 2 1.0 256x256
HISTORY:
Written Winter, 1992 by M.Morrison
22-May-93 (MDM) - Added document header
28-May-93 (MDM) - Modified header
V2.0 30-Sep-93 (MDM) - Changed to plot in units of QR DN
V2.1 25-Oct-93 (MDM) - Added option to display in image form
V2.2 15-Dec-93 (MDM) - Added code to remove the 4-NOV-93 08:21:32'
datapoint
V2.3 25-Nov-95 (MW) - Added code so image form goes past 2nd
entrance filter failure (marked by ***)
V2.4 09-Mar-98 (MW) - Even more code so image form goes past the
Jan-98 filter failure (marked by ***). Note
that the grayscale image of the data between
the 1st and 2nd failures is arbitrarily
multiplied by 2 to show up better.
[Previous]
[Next]
NAME:
plot_sls
PURPOSE:
To plot the SXT Scattered Light engineering data
CALLING SEQUENCE:
.run plot_sls
METHOD:
Dark frame images are taken with the narrow band analysis filters.
Since the leak after 17-Nov-92 entrance filter was intense enough
to show signal at the CCD even with the CCD closed, we use this
data to monitor the size of the aperature in the entrance filter.
The values plotted are the average signal for the whole FFI image.
The HR images are converted to the equivalent signal for a QR image.
The 70 DN "fat zero/dark current" of the QR images is still in the
values which are plotted.
There are two sets of values because of the difference in the
stopped position of the shutter blade. Two of the four positions
allow more light to leak around the blade.
NaBan/Open Half Dark C 2 1.0 512x512
NaBan/Open Qrtr Dark C 2 1.0 256x256
HISTORY:
Written Winter, 1992 by M.Morrison
22-May-93 (MDM) - Added document header
28-May-93 (MDM) - Modified header
V2.0 30-Sep-93 (MDM) - Changed to plot in units of QR DN
V2.1 25-Oct-93 (MDM) - Added option to display in image form
V2.2 15-Dec-93 (MDM) - Added code to remove the 4-NOV-93 08:21:32'
datapoint
[Previous]
[Next]
NAME:
plot_source_rate
PURPOSE:
Plots the countrate of BCS Cal. radioactive sources
CALLING SEQUENCE:
.run plot_source_rate
INPUTS:
OPTIONAL INPUTS:
KEYWORD INPUTS:
OUTPUTS:
OPTIONAL OUTPUTS:
RESTRICTIONS:
PROCEDURE:
reads info from file $DIR_BCS_LOGS/src_hk_summary.log
MODIFICATION HISTORY:
RDB 26-Nov-97 Written
[Previous]
[Next]
NAME:
plot_ssl
PURPOSE:
To plot the SSL results.
0. Plots of measured exposure levels
1. Plots of SXT temperatures
2. Plots of #/duration of exposures through different filters
3. Plots of percentage of data received
4. Plots of average/Max SXS1,SXS2,HXT-Low,RBM counting rates
SAMPLE CALLING SEQUENCE:
plot_ssl
plot_ssl, [0,1] ;only display the temperature plots
plot_ssl, /big
plot_ssl, /all
OPTIONAL INPUT:
code - A boolean array saying which plots to make.
code(0) = 1 says make plots of measured exposure levels
code(1) = 1 says make plots of SXT temperatures
code(2) = 1 says make plots of #/duration of exposures
through different filters
code(3) = 1 says make plots of percentage of data
received
code(4) = 1 says make plots of average SXS1,SXS2,
HXT-Low,RBM counting rates
CALLS: ***
OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], PAUSE [1], UTPLOT [1], UTPLOT [2], UTPLOT [3]
UTPLOT [4], UTPLOT [5], UTPLOT [6], gt_mbe [1], gt_mbe [2], gt_ssl_explab
gt_temp_ccd, gt_temp_hk, int2secarr [1], int2secarr [2], pause [2], plottime [1]
plottime [2], pprint [1], pprint [2], rd_ssl, yesnox [1], yesnox [2]
OPTIONAL KEYWORD INPUT:
all - If set, then make plots of everything
big - If set, then plot "Temperatures", "percent data received"
and "wbs/hxt counting rates" on one plot per page instead
of 4 plots per page.
year - If set, then put the year on the x labels
HISTORY:
Written May-92 by M.Morrison
V1.2 1-May-93 (MDM)
V2.0 25-Oct-93 (MDM) - Made into a procedure
- Added header information
13-Dec-94 (MDM) - Modified header information (described "code")
[Previous]
[Next]
Name: plot_sxg
Purpose: overplot SXT and GOES flux from SXG structure
Calling Sequence:
IDL> rd_sxg,time1, time2, sxg
IDL> [optionally filter sxg by Filter / Resolution]
IDL> plot_sxt,sxg
OR IDL> plot_sxt,sxg(ss)
CALLS:
[Previous]
[Next]
NAME:
plot_sxl
PURPOSE:
Plot SXL data (x-ray summary log of intensity of each SFD image)
SAMPLE CALLING SEQUENCE:
plot_sxl, sxl, ss, ymax, [log=log]
INPUT:
sxl - index-recs from sxl data files
ss - which indices to plot
ymax - max value to plot
log - switch to do log plot.
CALLS: ***
OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], UTPLOT [1], UTPLOT [2], UTPLOT [3]
UTPLOT [4], UTPLOT [5], UTPLOT [6], deriv_arr [1], deriv_arr [2], gt_filtb
int2secarr [1], int2secarr [2]
HISTORY:
Written by M.Morrison in 1991
18-Nov-93 (GAL) - ADDED MINIMUMAL DOCUMENTATION HEADER FOR PLOT_SXL
[Previous]
[Next]
NAME:
plot_sxl
PURPOSE:
To plot the average intensity of the SFD images over the mission
SAMPLE CALLING SEQUENCE:
.run plot_sxl
CALLS: ***
OUTPLOT [1], OUTPLOT [2], OUTPLOT [3], UTPLOT [1], UTPLOT [2], UTPLOT [3]
UTPLOT [4], UTPLOT [5], UTPLOT [6], deriv_arr [1], deriv_arr [2], gt_filtb
int2secarr [1], int2secarr [2]
HISTORY:
Written Fall, 91 by M.Morrison
V2.5 25-Oct-93 (MDM) - Changed criteria for removing bad points to be
less than 2000 (instead of less than 5000)
V2.6 18-Nov-93 (GAL) - Updated for layered structure of sxl.
18-Nov-93 (MDM) - Added /year switch to UTPLOT
[Previous]
[Next]
NAME:
plot_sxt_vig
PURPOSE:
Make the plots for SXT Cal Note 37 (Vignette function)
CALLING SEQUENCE:
plot_sxt_vig,/high
plot_sxt_vig,/high,/ps
plot_sxt_vig,/low
plot_sxt_vig,/low,/ps
plot_sxt_vig,/low,/ps,path=path ; If data files not on current directory
plot_sxt_vig,/low,/ps,/noplot ; Create plot file but don't print
INPUT KEYWORDS:
high - If set, compute the high-energy case
low - If set, compute the low-energy case
ps - If set, plot the results. Also create two files:
cal_note37_hi.ps and cal_note37_lo.ps
path - To specify the directory containing the data files.
PROCEDURE:
This is a very specialized routine to generate the final plots
of the Vignette study that was conducted by R. A. Fuller. His
subroutines are included in this file.
These routine require several data files to work. They are available,
in the on-line directories, but the user will have to copy the compressed
Tar file to his or her own scratch area. The Tar file must be uncompressed
and the files must be extracted.
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Instructions to restore the data file:
The data files are contained in a compressed Tar file called
fuller_vignette.tar.Z on the $DIR_SXT_SENSITIVE directory.
To run this program you must:
1. Copy the fuller_vignette.tar.Z file from $DIR_SXT_SENSITIVE
to a local directory. It might be best to create a new sub-directory.
Most of files will be names 'wc_*.dat' or 'wc_*.pc' and there is
a README file. On a Unix machine use the following commands:
% mkdir junk # Optional -- create a sub-directory called junk
% cd junk # Move to the sub-directory
% cp $DIR_SXT_SENSITIVE/fuller_vignette.tar.Z .
2. Make sure you are on the sub-directory which contains the compressed tar file.
% ls -l fuller_vignette*
3. At the Unix prompt type:
% uncompress fuller_vignette.tar.Z
4. To unpack the Tar file, at the Unix prompt type:
% tar xvf fuller_vignette.tar
5. Run this program in IDL. If you are not running from the same
directory as the data files, use the path=path to specify the
location of the files.
% idl
IDL> plot_sxt_vig,/lo ; For low-energy plot (see calling sequence)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
In these programs the centroid positions and the slopes have been
fixed to the best values. The routine will fit the data point to
these curves and compute the Chi square values.
The /special switch eliminates 3 data points that have partial annuli
yet are within the nominal "inner-region". This was a Fuller fix for
the fact that the inner- and outer-region cones are not concentric.
That is, the /special switch removes points that are really in the
outer region.
RESTRICTIONS/SIDE EFFECTS:
Will try to set the plot display to 'X' for diagnostic plots.
Not guaranteed on VMS.
MODIFICATION HISTORY:
2-Feb-94, J. R. Lemen, Written
[Previous]
[Next]
MODIFICATION HISTORY:
21-jul-93, R. A. Fuller, Written
29-jan-94, RAF + JRL, Fixed defaults for least scatter
run with /special to get best case
31-jan-94, JRL, Renamed to eff_area_lo; Pass in the path for the data files.
2-feb-94, JRL, Added noplot and fileout keywords
[Previous]
[Next]
MODIFICATION HISTORY
19-aug-93, R. A. Fuller, Written
29-jan-94, RAF+JRL, Modified defaults to give best results (slope and Rmax)
31-jan-94, JRL, Renamed to eff_area_hi; Pass in the path for the data files
2-feb-94, JRL, Added noplot and fileout keywords
[Previous]
[Next]
NAME:
PLOT_TEMPS2
PURPOSE:
To plot the SXT temperatures given an SXT index array
SAMPLE CALLING SEQUENCE:
.run plot_temps2
INPUT:
index - a variable with the name "index" must exist and be
an SXT index array
METHOD:
One option is to do the following:
ff = file_list('$DIR_SXT_SDC', '*')
rd_xda, ff, -1, index, /nodata
.run plot_temps2
Another option is:
get_sfm,timegrid('1-oct-91',!stime,days=1),index,/nodata
.run plot_temps2
qnot10 = 1
.run plot_temps2
HISTORY:
Written Jun-92 by M.Morrison
13-Dec-94 (MDM) - Added documentation header
V2.2 15-Feb-95 (MDM) - Added to documentation header
- Allow option to not plot 10 to a page
[Previous]
[Next]
NAME:
plot_trav
PURPOSE:
To plot the travel plans for the SXT team at ISAS. User has
ability to interactively expand the plot over a given time period.
CALLING SEQUENCE:
plot_trav ;plot all sxt people, full time range
plot_trav, /sxt ;plot all sxt people, full time range
plot_trav, /bcs ;plot all bcs people, full time range
plot_trav, '1-jan-91', '1-jan-92' ;all people for 1991
plot_trav, /lparl_ops
plot_trav, /noplot_n
plot_trav, infil='travel_loren.txt'
plot_trav, /dates
OPTIONAL INPUT:
p_st_day- The starting day to plot. If undefined, plot
from first day of data.
p_en_day- The ending day to plot. If undefined, plot to
last day of data.
CALLS: ***
CLEAR_UTPLOT [1], CLEAR_UTPLOT [2], RD_TRAVEL, TRAV_ASSIGN_GRP, UNIQ [1], UNIQ [2]
UNIQ [3], UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ex [1], anytim2ex [2], anytim2ints [1], anytim2ints [2], fmt_tim [1]
fmt_tim [2], int2secarr [1], int2secarr [2], plot_lcur [1], plot_lcur [2]
plottime [1], plottime [2], pprint [1], pprint [2]
CALLED BY:
quick_plottrav
OPTIONAL KEYWORD INPUT:
trav_arr- The structure which holds the travel information.
If undefined, read the data
lparl_ops- If set, only plot the key LPARL operations
personnel.
nozoom - If set, do not allow interactive zooming
noplot_n- If set, do not plot a summary of the number of
personnel present at ISAS.
sxt - If set, use the SXT team file
bcs - If set, use the BCS team file
infil - Read the travel information from this file
dates - If set, print the starting and ending date of the travel
on the plots
ops - If set, plot the OPS people (LPARL plus Co-I)
noclear - If set, do not call CLEAR_UTPLOT at the end of the routine
apartments- If set, print the assigned apartment number.
only_lparl - If set, only list LPARL or LTSC people
nomark - If set, don't mark the current time
invert - If set, black on white (default for PS)
charsize - character size (default=1) (scaled by existing hard-coded #)
HISTORY:
Written 29-Apr-92 by M.Morrison
22-Sep-92 (MDM) - Added /SXT and /BCS
- Added /YEAR switch to UTPLOT call
22-Feb-93 (MDM) - Added Bruner to /LPARL list
4-Mar-93 (MDM) - Added INFIL option
17-Mar-93 (MDM) - Added DATES option
1-Apr-93 (MDM) - Added /OPS option
12-Apr-93 (MDM) - Implemented the /NOPLOT_N option
- Added /NOCLEAR option
4-May-93 (MDM) - Modified DATES option to print the proper ending date
of the travel. It was off by one.
22-Sep-93 (MDM) - Minor mod (tit=' ' for # people at isas plot)
2-Nov-93 (RDB) - Small correction of name alignment on y-axis; got rid of "0" at top
8-Nov-93 (MDM) - Added Handey and Alexander to OPS and LPARL listing
7-Feb-94 (MDM) - Changed calls to UTPLOT because of changes to UTPLOT
2-Jun-94 (GAL) - Changed hardcopy calls to pprint from lprint.
8-Jun-94 (GAL) - Added apartments option and reorganized database.
3-Aug-94 (MDM) - Added /ONLY_LPARL option
3-Aug-94 (GAL) - Added Weber to OPS and LPARL listings
8-Aug-94 (MDM) - Changed plot command to not use ytickname and to draw the
labels myself (since there are over 31)
3-Nov-94 (JRL) - Fixed tick label software. Added nomark keyword. Fixed overplotting
of y-axis by the year boundary markers.
20-Feb-95 (JRL) - Added Shing to OPS and LPARL listings
28-Sep-95 (GAL) - Added Serge to OPS and LPARL listings
2-Apr-96 (GAL) - Moved Jean-Pierre to LM group.
5-Apr-96 (GAL) - Moved Tom Metcalf to LM group.
28-Jun-96 (NN) - Inserted /update in rd_travel which is used for plot
10-jul-96 (SLF) - add INVERT, keyword and function
30-Sep-96 (NN) - Eliminated Handy and Moved Alexander to Lockheed
7-May-98 (NN) - Eliminated Savy, Added Davey
11-Jan-99 (NN) - Added Fletcher (LMMS), Martens (MSU)
18-May-99 (NN) - Re-instated Handy (MSU)
29-Aug-00 (NVN) - Eliminated Weber, added Takeda as Hudson category
30-Apr-01 (NVN) - Replaced Shirts with Bartus
19-Jul-01 (NVN) - Added Sato (MSU)
[Previous]
[Next]
NAME:
plot_trav
PURPOSE:
To plot the travel plans for the SXT team at ISAS. User has
ability to interactively expand the plot over a given time period.
CALLING SEQUENCE:
plot_trav ;plot all sxt people, full time range
plot_trav, /sxt ;plot all sxt people, full time range
plot_trav, /bcs ;plot all bcs people, full time range
plot_trav, '1-jan-91', '1-jan-92' ;all people for 1991
plot_trav, /lparl_ops
plot_trav, /noplot_n
plot_trav, infil='travel_loren.txt'
plot_trav, /dates
OPTIONAL INPUT:
p_st_day- The starting day to plot. If undefined, plot
from first day of data.
p_en_day- The ending day to plot. If undefined, plot to
last day of data.
CALLS: ***
CLEAR_UTPLOT [1], CLEAR_UTPLOT [2], RD_TRAVEL, TRAV_ASSIGN_GRP, UNIQ [1], UNIQ [2]
UNIQ [3], UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ex [1], anytim2ex [2], anytim2ints [1], anytim2ints [2], fmt_tim [1]
fmt_tim [2], int2secarr [1], int2secarr [2], plot_lcur [1], plot_lcur [2]
plottime [1], plottime [2], pprint [1], pprint [2]
CALLED BY:
quick_plottrav
OPTIONAL KEYWORD INPUT:
trav_arr- The structure which holds the travel information.
If undefined, read the data
lparl_ops- If set, only plot the key LPARL operations
personnel.
nozoom - If set, do not allow interactive zooming
noplot_n- If set, do not plot a summary of the number of
personnel present at ISAS.
sxt - If set, use the SXT team file
bcs - If set, use the BCS team file
infil - Read the travel information from this file
dates - If set, print the starting and ending date of the travel
on the plots
ops - If set, plot the OPS people (LPARL plus Co-I)
noclear - If set, do not call CLEAR_UTPLOT at the end of the routine
apartments- If set, print the assigned apartment number.
only_lparl - If set, only list LPARL or LTSC people
nomark - If set, don't mark the current time
invert - If set, black on white (default for PS)
charsize - character size (default=1) (scaled by existing hard-coded #)
HISTORY:
Written 29-Apr-92 by M.Morrison
22-Sep-92 (MDM) - Added /SXT and /BCS
- Added /YEAR switch to UTPLOT call
22-Feb-93 (MDM) - Added Bruner to /LPARL list
4-Mar-93 (MDM) - Added INFIL option
17-Mar-93 (MDM) - Added DATES option
1-Apr-93 (MDM) - Added /OPS option
12-Apr-93 (MDM) - Implemented the /NOPLOT_N option
- Added /NOCLEAR option
4-May-93 (MDM) - Modified DATES option to print the proper ending date
of the travel. It was off by one.
22-Sep-93 (MDM) - Minor mod (tit=' ' for # people at isas plot)
2-Nov-93 (RDB) - Small correction of name alignment on y-axis; got rid of "0" at top
8-Nov-93 (MDM) - Added Handey and Alexander to OPS and LPARL listing
7-Feb-94 (MDM) - Changed calls to UTPLOT because of changes to UTPLOT
2-Jun-94 (GAL) - Changed hardcopy calls to pprint from lprint.
8-Jun-94 (GAL) - Added apartments option and reorganized database.
3-Aug-94 (MDM) - Added /ONLY_LPARL option
3-Aug-94 (GAL) - Added Weber to OPS and LPARL listings
8-Aug-94 (MDM) - Changed plot command to not use ytickname and to draw the
labels myself (since there are over 31)
3-Nov-94 (JRL) - Fixed tick label software. Added nomark keyword. Fixed overplotting
of y-axis by the year boundary markers.
20-Feb-95 (JRL) - Added Shing to OPS and LPARL listings
28-Sep-95 (GAL) - Added Serge to OPS and LPARL listings
2-Apr-96 (GAL) - Moved Jean-Pierre to LM group.
5-Apr-96 (GAL) - Moved Tom Metcalf to LM group.
28-Jun-96 (NN) - Inserted /update in rd_travel which is used for plot
10-jul-96 (SLF) - add INVERT, keyword and function
30-Sep-96 (NN) - Eliminated Handy and Moved Alexander to Lockheed
[Previous]
[Next]
PLOTBDA v1.5 (IDL2)
Plot bda file - Menu structure, quick look plotter, can be used to
select + extract bda datasets. Test_rd compatible.
cursor select on light curve, summation of datasets.
if you intend to change anything please read comments.
author: Andy Phillips , mssl 1991
in: filename - name of a BDA file.
out: Plots of counts vs time.
required: Rd_bsd_header.pro - read bsd header
rd_bsd_rdmap.pro - read bsd roadmap
Rd_bsd_data.pro - read bsd data sections
fmt_tim.pro - format bsd time array into a string
range.pro - do a regexp expansion.
syntax: plotbda[,bdafilename]
(else prompted for [dir]filename[*])
comments:
spectra are selected by certain range criteria.
parsing of the selection string (spectstr) is accomplished in range.pro
it is not bullet-proof. The selection string contains a comma delinited
series of numbers, and/or two wild cards.
* selects all spectra, the rest of the string is ignored.
- selects a range. eg 1-3 == 1,2,3.
bounds checking is performed, redundancy is not. it is perfectly
possible to replot the same spectrum several times eg.
1,2,3,1-2,3,3,3 will plot spectra 1,2,3,1,2,3,3,3.
Implicit assumption in range option that max(x) = x(max) , min(x) = x(0)
error handling uncertain.
borrowed mons' code for wild cards.
[Previous]
[Next]
NAME:
plotbft
PURPOSE:
Plot results of BSDFIT for BCS analysis.
CALLING SEQUENCE:
plotbft[,BSDNAM,chan='1,2,3',/plot_page,/landscape,/p2nd]
INPUTS:
BSDNAM The BSD filename in the form 'bsdyymmdd.hhmm'
OPTIONAL INPUT PARAMETERS:
CHANNEL=[2,3] Scalar or vector containing required channels
/PLOT_PAGE The presence of this parameter causes the light curve,
electron temperature, and emission measure to be plotted
on one page, and turbulent velocity (from Te-Td),
upflow emission measure and upflow velocity
time histories to be plotted on another page.
/landscape This switch should produce lanscape orientation plots
when the postscript output option is selected.
/p2nd Plot the ratio of the first and second component
emission measures.
OUTPUTS:
Only plots are produced.
CALLS: ***
ALL_VALS [1], ALL_VALS [2], BCS_GETFILE, BCS_IONVEL, BREAK_FILE [1]
BREAK_FILE [2], BREAK_FILE [3], CLEAR_UTPLOT [1], CLEAR_UTPLOT [2], ERRPLOT [1]
ERRPLOT [2], PS_LONG [1], PS_LONG [2], PS_LONG [3], PS_RESET [1], PS_RESET [2]
PS_RESET [3], RD_BSD_DATA, RD_BSD_HEADER, RD_BSD_RDMAP, READBFT, SET_UTPLOT [1]
SET_UTPLOT [2], SET_XY, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5]
UTPLOT [6], break_file [4], clearplot [1], clearplot [2], fmt_tim [1], fmt_tim [2]
int2secarr [1], int2secarr [2]
COMMON BLOCKS:
None.
SIDE EFFECTS:
1) The BSD, BTH and BFT files are openned and read.
2) Calls procedures rd_bsd_header,rd_bsd_rdmap,rd_bsd_data and READBFT
RESTRICTIONS:
BSDNAM must be a string-type scalar.
MEGAPLOT reads the BSD and FIT files. The BSD
and FIT files must exist on the same directory.
PROCEDURE:
The following plots are generated in turn for each channel:
1) Light curve - All other time plots are plotted
on the same scale.
2) Electron temperature - In M K.
3) Doppler temperature - In M K. - only when single plots on page
4) Emission measure - * 1.e50 cm-3.
5) Ion (turbulent) velocity - In km s-1 (from Te - Td).
6) Upflow emission meas - * 1.e50 cm-3
7) Uplow velocity - in km s-1
MODIFICATION HISTORY:
Nov, 1985, Written, J.R. Lemen, MSSL
Mar, 1986, RDB+JRL - Improved titles; made compatible with EOF/MSSL
May, 1986, AF - change to use new two-component FIT files
jul, 1991, atp - cut out some rubbish + transfered to idl2.
Nov, 1991 AF - made it work again
May, 1992 RDB Cleaned up routine a lot; sped it up; got rid of
redundant stuff to do with "transfer" from SMM
Jun, 1992 Atp - finished, converted to structures, much updates
much yet to do.
july 1992 vn2.1 atp - cleaned up + finished (!)
28-Feb-92 RDB Changed VEL to BCS_IONVEL
01-Mar-93 RDB Added CLEARPLOT, CLEAR_UTPLOT at end
02-Mar-92 RDB Corrected PS_LONG and /LAND interaction
Channel select now by integer
Added prompt for pause when to screen
18-May-94 RDB Fixed problem in creating tdx, and with plots
[Previous]
[Next]
NAME:
plotbft
PURPOSE:
Plot results of BSDFIT for BCS analysis.
CALLING SEQUENCE:
plotbft[,BSDNAM,chan='1,2,3',/plot_page,/landscape,/p2nd]
INPUTS:
BSDNAM = The BSD file name
OPTIONAL INPUT PARAMETERS:
/CHANNEL - String containing channels to be processed,
'*' means all, or use a range '1-3', or explicitly
'1,3,4', default is '*'
/PLOT_PAGE - The presence of this parameter causes the light curve,
electron temperature, and emission measure to be plotted
on one page, and turbulent velocity (from Te-Td),
upflow emission measure and upflow velocity
time histories to be plotted on another page.
/landscape - This switch should produce lanscape orientation plots
when the postscript output option is selected.
/p2nd - Plot the ratio of the first and second component
emission measures.
OUTPUTS:
Only plots are produced.
CALLS: ***
BCS_GETFILE, BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], ERRPLOT [1]
ERRPLOT [2], PLOTBFT_AS, RD_BSD_DATA, RD_BSD_HEADER, RD_BSD_RDMAP, READBFT
SET_UTPLOT [1], SET_UTPLOT [2], SET_XY, UTPLOT [1], UTPLOT [2], UTPLOT [3]
UTPLOT [4], UTPLOT [5], UTPLOT [6], VEL, break_file [4], int2secarr [1]
int2secarr [2], range
COMMON BLOCKS:
None.
SIDE EFFECTS:
1) The BSD, BTH and BFT files are openned and read.
2) Calls procedures rd_bsd_header,rd_bsd_rdmap,rd_bsd_data and READBFT
RESTRICTIONS:
BSDNAM must be a string-type scalar.
MEGAPLOT reads the BSD and FIT files. The BSD
and FIT files must exist on the same directory.
PROCEDURE:
The following plots are generated in turn for each channel:
1) Light curve - All other time plots are plotted
on the same scale.
2) Electron temperature - In M K.
3) Doppler temperature - In M K. - only when single plots on page
4) Emission measure - * 1.e50 cm-3.
5) Ion (turbulent) velocity - In km s-1 (from Te - Td).
6) Upflow emission meas - * 1.e50 cm-3
7) Uplow velocity - in km s-1
MODIFICATION HISTORY:
Nov, 1985, Written, J.R. Lemen, MSSL
Mar, 1986, RDB+JRL - Improved titles; made compatible with EOF/MSSL
May, 1986, AF - change to use new two-component FIT files
jul, 1991, atp - cut out some rubbish + transfered to idl2.
Nov, 1991 AF - made it work again
Jun, 1992 Atp - finished, converted to structures, much updates
much yet to do.
july 1992 vn2.1 atp - cleaned up + finished (!)
[Previous]
[Next]
NAME:
PLOTBPC
PURPOSE:
Read the BPCxxxxxx.xxxx files written by BSDCAL and the BSD data file
specified within it and plot the spectral fits as given by parameters
in the BPC file over the raw data
CALLING SEQUENCE:
PLOTBPC[,bpc_file,chan,option,start]
INPUTS:
BPC_FILE= Equal the BPC file name or I.D.
If BPC_FILE is undefined or the first character is
blank, PLOTBPC will prompt user for name.
CHAN = required channel;
OPTION = the fit option number required to display. If supplied
on the command line then only one option is allowed
START = the time of first spectrum required ( default = all)
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CALFUN, POLY, RD_BSD_HEADER
break_file [4], dom, free_in
COMMON BLOCKS:
None.
SIDE EFFECTS:
None.
MODIFICATION HISTORY:
0) Written May 1991 CDP
1) Revised to handle start time of spectra better Sep 91 CDP
2) Undo previous "bug-fix" and fix real one where rdmap variable
was incorrectly incremented between spectra Nov 91 CDP
3) Cut out bin-width scaling (not needed once BSDCAL sorted) Nov 91 CDP
4) Generalize and UNIXify handling of file names. Mar 92 CDP
5) Changed BGRND call to routine POLY. June 92, RDB
6) Use rd_bsd_head to access BSD file header July 92 CDP
[Previous]
[Next]
NAME: PLOTBSD v3.06 (IDL2)
PURPOSE: Plot bsd file - Menu structure, quick plotter
CALLING SEQUENCE: plotbsd[,bsdfilename]
(else prompted for [dir]filename[*])
INPUTS: A Valid BSD file (See Mkbsd.exe)
OUTPUTS: Plots to graphics devices.
CALLS: ***
BCS_GETFILE, BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], PLOTBSD, RD_BSD_DATA
RD_BSD_HEADER, RD_BSD_RDMAP, break_file [4], fmt_tim [1], fmt_tim [2], lcbsd
pprint [1], pprint [2], range
RESTRICTIONS: This is a procedure, unlike plotbda, so data is
not available at program termination. Use
get_bsd to read in a BSD file to IDL.
You have to explicitly print postscript files after
the program has finished. use PPRINT.
PROCEDURE: Read selected bsd file. Setup plot parameters at
main menu. Run series of plots, with behaviour
modified by plot submenu.
MODIFICATION HISTORY:
author: Andy Phillips , mssl 1991
vn 3.0: atp, 15-jun-92
[Previous]
[Next]
NAME: plotbth v3.0 (IDL2)
PURPOSE: Plot bsd file - Menu structure, quick look plotter
overplotting of theory spectra.
CALLING SEQUENCE: plotbth[,bsdfilename]
(else prompted for [dir]filename[*])
INPUTS: Valid BSD file (see MKBSD.exe)
Matching SPC file (See BSDFIT.exe)
OUTPUTS: Plots to graphics devices.
CALLS: ***
BCS_GETFILE, BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], PLOTBTH, RD_BSD_DATA
RD_BSD_HEADER, RD_BSD_RDMAP, RD_SPC_DATA, RD_SPC_HEADER, RD_SPC_RDMAP, X [1], X [2]
break_file [4], fmt_tim [1], fmt_tim [2], lcbsd, pprint [1], pprint [2], range
RESTRICTIONS: This is a procedure, unlike plotbda, so data is
not available at program termination. Use
get_bsd to read in a BSD file to IDL.
use get_bth to read in the bth (SPC) file.
You have to explicitly print postscript files after
the program has finished. use PPRINT.
PROCEDURE: Read selected BSD file. Read selected BTH file.
Setup plot parameters at main menu. Run series of
plots, with behaviour modified by plot submenu.
MODIFICATION HISTORY:
author: Andy Phillips , mssl 1991
vn 1.8: atp, 15-jun-92
vn 2.1: atp, 07-jul-92 mods for spc structures
vn 3.0: atp, 23/11/92 - rationalise with ISAS
RDB 01-Mar-93 Changed "BTH" to "bth"
RDB 26-May-93 Modified oplot calls.
[Previous]
[Next]
NAME:
PLOTCAL_ATT
PURPOSE:
Plot hxa attitude data, giro data, and (optional) sxt pixel shifts
at the same scale. Input: output from getcal_att. Units: SXT pixels
CATEGORY:
CALLING SEQUENCE:
plotcal_att,giro,tstart,tend
INPUTS:
giro = float(3,*). Pixel shift determined from giro's (calibrated)
giro(0,*) : pixel shift in x,
giro(1,*) : pixel shift in y,
giro(2,*) : time (msec).
tstart = start time. Either a string (hh:mm:ss), or a long or float
variable with time in msec.
tend = end time.
KEYWORDS (INPUT):
sxt = float(3,*). Pixel shift determined from SXT NaBan images
(e.g., with "korrel"). Same format as giro.
hxa = float(3,*). Pixel shift determined from HXA detector.
Same format as giro.
raw_hxa = float(5,*). Raw hxa data to check for fiducial marks.
Note: this keyword causes the program to plot ONLY the
raw hxa data.
OUTPUTS:
No explicit output.
KEYWORDS (OUTPUT):
CALLS: ***
TIME2MSEC
COMMON BLOCKS:
None.
SIDE EFFECTS:
Makes a plot.
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
JPW, Jan. 1992
[Previous]
[Next]
NAME:
PLOTERR_XY
PURPOSE:
Plot data points with accompanying error bars.
(See also OPLOTERR_XY.)
CATEGORY:
Plotting, two-dimensional.
CALLING SEQUENCE:
PLOTERR, X , Y , Err_X, ERR_Y [, PSYM = Psym] [, TYPE = Type]
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:
PSYM: The plotting symbol to use. The default is +7.
ITYPE: The type of plot to be produced. The possible types are:
ITYPE = 0 : X Linear - Y Linear (default)
ITYPE = 1 : X Linear - Y Log
ITYPE = 2 : X Log - Y Linear
ITYPE = 3 : X Log - Y Log
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: plotout_pha.pro
Author: Andy Phillips mssl
Purpose: Plot Pha values against chisqrd for data from readout.pro
Classification: BDA, Housekeeping, spikes correction
History: written jan 1992 atp
[Previous]
[Next]
NAME:
plots_bda
PURPOSE:
Plots the BCS spectral data. It assumes that the data
is passed in decompressed.
CALLING SEQUENCE:
plots2_bda, index, data, psym=psym
INPUT:
index - index structure
data - data array
OPTIONAL INPUT:
psym - plotting symbol. Default is 10 (histogram)
CALLS: ***
ext_bcschan [1], ext_bcschan [2], fmt_tim [1], fmt_tim [2], plottime [1]
plottime [2]
HISTORY:
written 7-Mar-92 by M.Morrison
[Previous]
[Next]
NAME:
plots_bdaf
PURPOSE:
To plot the BCS spectral data when the data mode
is "fixed extraction" (hardware mode)
CALLING SEQUENCE:
plots_bdaf, index, data
INPUT:
index - index structure
data - data array
HISTORY:
Written Fall '91 by M.Morrison
8-Mar-92 (MDM) - Renamed and made a procedure
[Previous]
[Next]
NAME:
plots_wda
PURPOSE:
Plot WBS spectral data
INPUT:
index - index structure
data - data structure
CALLS: ***
fmt_tim [1], fmt_tim [2], plottime [1], plottime [2]
HISTORY:
Written Fall '91 by Greg Slater
7-Mar-92 (MDM) - Renamed and made it a procedure
[Previous]
[Next]
NAME:
plott_bda
PURPOSE:
Plot light curve for BCS data (uses the "total_cnts" field).
The "anomalous" negative datavalues are not plotted
CALLING SEQUENCE:
plott_bda, input1, input2
plott_bda, roadmap
plott_bda, index, psym=0
plott_bda, '15-nov-91 20:00', 60
INPUT:
A single parameter call is used to pass ROADMAP, INDEX or
OBSERVING LOG in. Two parameters can be used to specify the
date/time and duration or end time.
input1 - A structure array (roadmap, index or obs log)
or
The start time to plot from
input2 - If the value is a scalar, it is the duration in
minutes to plot. The other option is to specify the
end date/time.
CALLS: ***
GET_NBYTES, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ints [1], anytim2ints [2], data_type [1], data_type [2]
gt_total_cnts [1], gt_total_cnts [2], rd_obs
OPTIONAL KEYWORD INPUT:
psym - The plotting symbol to use. Default is 3 (+)
HISTORY:
Written 7-Mar-92 by M.Morrison
14-Feb-95 (RDB) Modified to accept times in addition to an HXT
structure (index/roadmap)
(MDM mode of plott_hda of 20-Dec-93
[Previous]
[Next]
NAME:
plott_bda_dp
PURPOSE:
Plot a light curve of the DP_SYNC data. Overplots
the limited and accumulated profiles in a different color.
CALLING SEQUENCE:
plott_bda_dp, dp_sync
INPUT:
dp_sync - the dp_sync data record
CALLS: ***
UNSIGN, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
int2secarr [1], int2secarr [2]
HISTORY:
Written Fall '91 by M.Morrison
[Previous]
[Next]
NAME:
plott_hda
PURPOSE:
Plot light curve for HXT data.
CALLING SEQUENCE:
plott_hda, input1, input2
plott_hda, roadmap
plott_hda, index, psym=0
plott_hda, '15-nov-91 20:00', 60
INPUT:
A single parameter call is used to pass ROADMAP, INDEX or
OBSERVING LOG in. Two parameters can be used to specify the
date/time and duration or end time.
input1 - A structure array (roadmap, index or obs log)
or
The start time to plot from
input2 - If the value is a scalar, it is the duration in
minutes to plot. The other option is to specify the
end date/time.
CALLS: ***
GET_NBYTES, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ints [1], anytim2ints [2], data_type [1], data_type [2], 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], rd_obs
OPTIONAL KEYWORD INPUT:
psym - The plotting symbol to use. Default is 3 (+)
HISTORY:
Written 7-Mar-92 by M.Morrison
20-Dec-93 (MDM) - Modified to accept times in addition to an HXT
structure (index/roadmap)
[Previous]
[Next]
PLOTT_HI_SUM
POURPOSE:
plot the total counts(summation of all the sensors) in HXT hi
energy channel(52.7-92.8keV) versus UT time.
CATEGORY:
Yohkoh/HXT analysis
CALLING SEQUENCE:
PLOTT_HI_SUM, data, index, [fint= fint]
INPUT:
data - HXT pc data (bytarr(4, 64, 4, *))
index - HXT index structure from reformatted data
OPTIONAL INPUT:
fint - data integration in flare mode. Default = 4.
In flare mode pulse count data comes 4 times per SF.
= 4(default): Sum 4 data in SF and output value
every SF(time resolution 1SF in flare md).
= 2: time resolution 1/2 SF in flare md.
= 1: time resolution 1/4 SF in flare md.
Include all keywords available to PLOT
OUTPUT:
None.
CALLS: ***
FCHECK, HXT_HI_SUM, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5]
UTPLOT [6], dprate2sec [1], dprate2sec [2], gt_dp_rate [1], gt_dp_rate [2]
int2secarr [1], int2secarr [2]
CALLED BY:
plott_hxt
SIDE EFFECTS:
Plot a graph in the current graphic device.
PROCEDURE:
Quiet mode: No data = 0
Flare mode: Get data from pc_data (data)
These are done with HXT_HI_SUM function.
RESTRICTIONS:
This data comes only in flare mode.
MODIFICATION HISTORY:
May 27, 1992 written by Toshifumi Shimizu
(Univ. of Tokyo)
[Previous]
[Next]
PLOTT_HXT
POURPOSE:
plot the total counts(summation of all the sensors) in every HXT
4 channels versus UT time.
CATEGORY:
Yohkoh/HXT analysis
CALLING SEQUENCE:
PLOTT_HXT, data, index, [fint= fint]
INPUT:
data - HXT pc data (bytarr(4, 64, 4, *))
index - HXT index structure from reformatted data
OPTIONAL INPUT:
fint - data integration in flare mode. Default = 4.
In flare mode pulse count data comes 4 times per SF.
= 4(default): Sum 4 data in SF and output value
every SF(time resolution 1SF in flare md).
= 2: time resolution 1/2 SF in flare md.
= 1: time resolution 1/4 SF in flare md.
Include all keywords available to PLOT
OUTPUT:
None.
CALLS: ***
plott_hi_sum, plott_low_sum, plott_med1_sum, plott_med2_sum
SIDE EFFECTS:
Plot a graph in the current graphic device.
PROCEDURE:
Plot light curve of 4 HXT channels in one window or one paper,
using PLOTT_LOW_SUM etc,PLOTT_MED1_SUM,PLOTT_MED2_SUM,PLOTT_HI_SUM.
RESTRICTIONS:
MODIFICATION HISTORY:
Aug 22, 1992 written by Toshifumi Shimizu
(Univ. of Tokyo)
[Previous]
[Next]
PLOTT_LOW_SUM
POURPOSE:
plot the total counts(summation of all the sensors) in HXT low
energy channel(13.9-22.7keV) versus UT time.
CATEGORY:
Yohkoh/HXT analysis
CALLING SEQUENCE:
PLOTT_LOW_SUM, data, index, [fint= fint]
PLOTT_LOW_SUM, index
INPUT:
data - HXT pc data (bytarr(4, 64, 4, *))
index - HXT index structure from reformatted data
OPTIONAL INPUT:
fint - data integration in flare mode. Default = 4.
In flare mode pulse count data comes 4 times per SF.
= 4(default): Sum 4 data in SF and output value
every SF(time resolution 1SF in flare md).
= 2: time resolution 1/2 SF in flare md.
= 1: time resolution 1/4 SF in flare md.
Include all keywords available to PLOT
OUTPUT:
None.
CALLS: ***
FCHECK, HXT_LOW_SUM, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5]
UTPLOT [6], dprate2sec [1], dprate2sec [2], gt_dp_rate [1], gt_dp_rate [2]
int2secarr [1], int2secarr [2]
CALLED BY:
plott_hxt
SIDE EFFECTS:
Plot a graph in the current graphic device.
PROCEDURE:
(1) PLOTT_LOW_SUM, data, index
Quiet mode: Get data from HXT_DHK in index
Flare mode: Get data from pc_data (data)
(2) PLOTT_LOW_SUM, index
Get data from HXT_DHK in index every time.
These are done with HXT_LOW_SUM function.
RESTRICTIONS:
(1) PLOTT_LOW_SUM, index
Index has the only integrated data of 4 pc data in 1 SF
even in Flare mode. (Restriction of reformatted data)
(2) There is a 4 second offset buffer in the HXT electronics.
Please check it. See document of reformatted data.
MODIFICATION HISTORY:
May 27, 1992 written by Toshifumi Shimizu
(Univ. of Tokyo)
[Previous]
[Next]
PLOTT_MED1_SUM
POURPOSE:
plot the total counts(summation of all the sensors) in HXT med1
energy channel(22.7-32.7keV) versus UT time.
CATEGORY:
Yohkoh/HXT analysis
CALLING SEQUENCE:
PLOTT_MED1_SUM, data, index, [fint= fint]
INPUT:
data - HXT pc data (bytarr(4, 64, 4, *))
index - HXT index structure from reformatted data
OPTIONAL INPUT:
fint - data integration in flare mode. Default = 4.
In flare mode pulse count data comes 4 times per SF.
= 4(default): Sum 4 data in SF and output value
every SF(time resolution 1SF in flare md).
= 2: time resolution 1/2 SF in flare md.
= 1: time resolution 1/4 SF in flare md.
Include all keywords available to PLOT
OUTPUT:
None.
CALLS: ***
FCHECK, HXT_MED1_SUM, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5]
UTPLOT [6], dprate2sec [1], dprate2sec [2], gt_dp_rate [1], gt_dp_rate [2]
int2secarr [1], int2secarr [2]
CALLED BY:
plott_hxt
SIDE EFFECTS:
Plot a graph in the current graphic device.
PROCEDURE:
Quiet mode: No data = 0
Flare mode: Get data from pc_data (data)
These are done with HXT_MED1_SUM function.
RESTRICTIONS:
This data comes only in flare mode.
MODIFICATION HISTORY:
May 27, 1992 written by Toshifumi Shimizu
(Univ. of Tokyo)
[Previous]
[Next]
PLOTT_MED2_SUM
POURPOSE:
plot the total counts(summation of all the sensors) in HXT med2
energy channel(32.7-52.7keV) versus UT time.
CATEGORY:
Yohkoh/HXT analysis
CALLING SEQUENCE:
PLOTT_MED2_SUM, data, index, [fint= fint]
INPUT:
data - HXT pc data (bytarr(4, 64, 4, *))
index - HXT index structure from reformatted data
OPTIONAL INPUT:
fint - data integration in flare mode. Default = 4.
In flare mode pulse count data comes 4 times per SF.
= 4(default): Sum 4 data in SF and output value
every SF(time resolution 1SF in flare md).
= 2: time resolution 1/2 SF in flare md.
= 1: time resolution 1/4 SF in flare md.
Include all keywords available to PLOT
OUTPUT:
None.
CALLS: ***
FCHECK, HXT_MED2_SUM, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5]
UTPLOT [6], dprate2sec [1], dprate2sec [2], gt_dp_rate [1], gt_dp_rate [2]
int2secarr [1], int2secarr [2]
CALLED BY:
plott_hxt
SIDE EFFECTS:
Plot a graph in the current graphic device.
PROCEDURE:
Quiet mode: No data = 0
Flare mode: Get data from pc_data (data)
These are done with HXT_MED2_SUM function.
RESTRICTIONS:
This data comes only in flare mode.
MODIFICATION HISTORY:
May 27, 1992 written by Toshifumi Shimizu
(Univ. of Tokyo)
[Previous]
[Next]
NAME:
plott_wda
PURPOSE:
Plot light curve for WBS data.
CALLING SEQUENCE:
plott_wda, input1, input2
plott_wda, roadmap
plott_wda, index, psym=0
plott_wda, '15-nov-91 20:00', 60
INPUT:
A single parameter call is used to pass ROADMAP, INDEX or
OBSERVING LOG in. Two parameters can be used to specify the
date/time and duration or end time.
input1 - A structure array (roadmap, index or obs log)
or
The start time to plot from
input2 - If the value is a scalar, it is the duration in
minutes to plot. The other option is to specify the
end date/time.
CALLS: ***
GET_NBYTES, UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5], UTPLOT [6]
anytim2ints [1], anytim2ints [2], data_type [1], data_type [2], gt_grs1, gt_grs2
gt_hxs [1], gt_hxs [2], gt_rbmsc, gt_rbmsd, gt_sxs1 [1], gt_sxs1 [2], gt_sxs2 [1]
gt_sxs2 [2], rd_obs
OPTIONAL KEYWORD INPUT:
psym - The plotting symbol to use. Default is 3 (+)
HISTORY:
Written 7-Mar-92 by M.Morrison
20-Mar-92 (MDM) - Changed to use the "title" option of the
GT routines
20-Dec-93 (MDM) - Modified to accept times in addition to an HXT
structure (index/roadmap)
[Previous]
[Next]
plots withbroe tremps for max and solar min holes
CALLS:
[Previous]
[Next]
NAME:
ploty
PURPOSE:
To make a light curve plot of all instruments with the same time range
CALLING SEQUENCE:
ploty, bcs, hxt, sxt, wbs
ploty, bcs, w_h, sxtp, w_h
ploty, brmap, hrmap, srmap, wrmap
ploty, '8-may-92 15:00', '8-may-92 18:40'
ploty, infil=infil
ploty, infil=infil, xmin='15-nov-91 22:20', xmax='15-nov-91 22:45'
INPUT:
bcs - The BCS roadmap, index or observing log structure array
hxt - The HXT roadmap, index or observing log structure array
sxt - The SXT roadmap, index or observing log structure array
wbs - The WBS roadmap, index or observing log structure array
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CLEAR_UTPLOT [1]
CLEAR_UTPLOT [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], GET_NBYTES
LPRINT, Rd_Roadmap [2], UTPLOT [1], UTPLOT [2], UTPLOT [3], UTPLOT [4], UTPLOT [5]
UTPLOT [6], anytim2ints [1], anytim2ints [2], break_file [4], concat_dir [4]
gt_dp_mode [1], gt_dp_mode [2], gt_filta, gt_filtb, gt_hxs [1], gt_hxs [2], gt_rbmsd
gt_sum_h [1], gt_sum_h [2], gt_sum_l [1], gt_sum_l [2], gt_sum_m2 [1]
gt_sum_m2 [2], gt_sxs1 [1], gt_sxs1 [2], gt_total_cnts [1], gt_total_cnts [2]
int2secarr [1], int2secarr [2], plot_lcur [1], plot_lcur [2], plottime [1]
plottime [2], rd_obs, rd_roadmap [1], sxt_lcur
OPTIONAL KEYWORD INPUT:
zoom - Allow user to interactively zoom in on time range
option - FUTURE USE - different plotting items
flare - FUTURE USE - only plot flare portion of time range or data
sxt_filt- Integer array of 2 filters to plot the light curve for. Default
is filtb = 2,4 = thin aluminum and Be
Values 1-6 are filter B 1-6, and 11-16 are filter A 1-6
infil - input a sample file name, and it will read the
roadmaps for those files
big - If set and device is 'X' window, make it 900x900
xmin - The starting time to plot
xmax - The ending time to plot
psym - The IDL PLOT PSYM value to use
HISTORY:
Written 7-Jul-92 by M.Morrison
Ver 1.1 6-Oct-92 (MDM) - Added CLEAR_UTPLOT call and modified header info.
[Previous]
[Next]
Name: PLT_STR
Purpose: Plot Structure Fields Against Subscript
Input Parameters: Structure Type (Simple or Nested)
Optional Input: Pmulti - copied into !p.multi for page formatting
up to 5 element array [erase, colums, rows,...]
Keyword Parameter: All - If set, then all fields are plotted
Number of pages produced depends on pmulti
If not set, user is prompted (menu) for tag
Calling Sequence: PLT_STR, structure [,pmulti=pmulti, /all]
Example PLT_STR, structure, [0,4,5], /all
(plots all fields of structure, 20/page, 4col/5rows)
PLT_STR, structure
(prompts user for tags and plots until user quits)
Restrictions: Boring for non-structures
(for that matter, pretty boring for structures, too)
Note if ALL is set while running from X-Terminal and the
Number of tags > then the single page plot total as set by
pmulti, then plots are overwritten
Method: Straightforward (recursive for nested tags)
[Previous]
[Next]
NAME:
plusmark
PURPOSE
Plot a tiny cross at the specified place on a window.
CALLING SEQUENCE
plusmark,x,y[,color=color,thick=thick,big=big
INPUT
x,y - The pixel location of the desired position of the cross.
OPTIONAL KEYWORD INPUT
color - color of the cross, default is 255.
thick - line thickness of the cross, default=1.0
big - size of the cross, default = 7 pixels.
HISTORY
LWA, 22-Jun-94
LWA, added keywords thick and big. 1-Nov-94
CALLED BY
azi_calc, azimuth_avg, radial_avg, scat_avg, sector_avg
[Previous]
[Next]
NAME:
PM_CLEV
PURPOSE:
To return an array of elements which is to be used as contour levels
when calling the CONTOUR procedure. Particularly for drawing
magnetogram contours.
CALLING SEQUENCE:
clev = PM_CLEV(clevels, minus=minus)
then,
CONTOUR, ... , level=clev
INPUT:
clevels - An array which contains contour levels.
OPTIONAL KEYWORD INPUT:
plus - dummy keyword.
minus - If set, then all the elements in the CLEVELS array will
change the sign.
OUTPUT:
clev - Contour levels which are SORTed so that they can be fed to
the CONTOUR procedure with the LEVEL keyword.
CATEGORY:
Yohkoh HXT data analysis; especially for comparison with magnetograms.
HISTORY:
version 1.0 98.06.29 (Mon) T.Sakao written.
[Previous]
[Next]
NAME:
PMM
PURPOSE:
print max and min of something
CALLING SEQUENCE:
pmm, variable
CALLED BY:
HSI_EVENTLIST_TEST, HSI_EVFILE_TEST, HSI_IMAGE_TEST [2]
HSI_IMAGE_TEST__define, OBS_PLOT, READ_TRACE_FOV, VIEW_TRACE_MOSAIC, azi_calc
azimuth_teem [2], fit_gcross, fudge_scale, hsi_count_rebin, rd_hk [1], rd_hk [2]
ssw_build_trace, sxl2radiance, sxt_erg_per_dn
HISTORY:
Hugh Hudson, Dec. 18, 1992
[Previous]
[Next]
NAME:
pmnmx
PURPOSE:
1 liner to print min and max of an array
CALLING SEQUENCE:
pmnmx(array)
HISTORY:
written by A.McAllister, may-7-1992
[Previous]
[Next]
NAME:
PMODE
PURPOSE:
Change the timing parameters of the Peritek VCT graphics
interface to output NTSC, high-res-interlaced,
or high-res-non-interlaced signals
CATEGORY:
Image display through the Peritek VCT graphics interface.
CALLING SEQUENCE:
PMODE,RESULT=RESULT
INPUTS:
MODE: Video signal mode (0,1,2 = NTSC, high-res-interlaced,
high-res-non-interlaced
KEYWORD PARAMETERS:
RESULT = error message
CALLED BY:
go_lasdisk golaserdisk, go_lasdisk2 golaserdisk
PROCEDURE:
Spawns calls to two C routines which perform the necessary tasks.
MODIFICATION HISTORY:
10-Jul-1991 - Written by GLS
15-nov-1991 - Updated, JRL -- Added DIR_PERITEK environment
variable to reference C routines.
5-may-1991 - Updated, GLS -- Use !version to determine which
C routine (V bus version or Q bus
version) to call.
[Previous]
[Next]
NAME:
pn_pl
PURPOSE:
Total photon flux of a single power law
CALLING SEQUENCE:
energy = pn_pl(k,m,e0,e1=e1,flat_to=flat_to)
INPUT:
m; spectral index
k; emission at 1keV
e0=low energy limit
OUTPUT:
energy = energy flux in photons
KEYWORDS:
e1 = the high energy limit of the calculation, default is 20MeV
flat_to=if set, the spectrum will be assumed to be flat down to
the value given in flat_to. If flat_to>e0, then flat_to=e0
CALLED BY:
fsp_avsig
HISTORY:
JMcT, 29-apr-94
[Previous]
[Next]
NAME: pnn
PURPOSE: get min, max of given argument (up to ten arguments)
CATEGORY: analysis
CALLING SEQUENCE: pnn,blarch ;blarch = array, image, vector, etc...
pnn,blarch,/mm
pnn,blarch,mm=mm,/quiet ; no print
INPUTS: blarch
OPTIONAL INPUT PARAMETERS: none
KEYWORD PARAMETERS: /help, /mm
OUTPUTS: prints min, max of blarch. with mm set, return the min and max
for it too, with mm=FLTARR(2,N) for N positional parameters.
CALLS: ***
DOC_LIBRARY
COMMON BLOCKS: none
SIDE EFFECTS: none
RESTRICTIONS: none
PROCEDURE: simple min, max call
MODIFICATION HISTORY: ANM 910507. Updated 911031. Doc'd DKL 930426
[Previous]
[Next]
NAME:
pnt2atr
PURPOSE:
To take the existing PNT files and convert them into ATR files.
INPUT:
week_str- The week to process (ie: "93_24")
OPTIONAL KEYWORD OUTPUT:
run_time- The number of minutes it took to run
CALLS: ***
ATT_STRUCT [1], ATT_STRUCT [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
GEN_STRUCT [1], GEN_STRUCT [2], GET_NBYTES, OBS_STRUCT [1], OBS_STRUCT [2]
OPEN_DA_FILE, PNT_STRUCT, REF_STRUCT, Rd_NewOrb_p [1], Rd_NewOrb_p [2]
Rd_Pointer [1], Rd_Pointer [2], Rd_fHeader [1], Rd_fHeader [2], Rd_fHeader [3]
WRT1P_FH, concat_dir [4], rdwrt [1], rdwrt [2], rdwrt [3]
HISTORY:
Written 16-Jul-93 by M.Morrison
[Previous]
[Next]
NAME: pnt_exst
PURPOSE: check if there are PNT files
CALLING SEQUENCE: ans=pnt_exst(index)
PARAMETERS: times list of times
CALLS: ***
UNIQ [1], UNIQ [2], UNIQ [3], file_list [1], file_list [2], tim2orbit [1]
tim2orbit [2]
CALLED BY:
GEN_LOC
HISTORY: Cut out from SXT_CEN(MDM), 22-jun-93, AMcA.
[Previous]
[Next]
NAME:
PNT_STRUCT
PURPOSE:
Define the following PNT specific database structures
* PNT_Data_Rec
CALLING SEQUENCE:
PNT_STRUCT
CALLED BY:
mk_pnt, pnt2atr, 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, Feb 92.
[Previous]
[Next]
NAME:
POINT_HEX2
PURPOSE:
Generate an OG for a specified pointing location. The output
is in hexadecimal, appropriate for OG editing via OPOGEDT.
CATEGORY:
Yohkoh
CALLING SEQUENCE:
point_hex2, angle, og_number, bc, offset = offset
point_hex2, [-.0797+E1/3600.,.0295+E2/3600.], /name
(here the Euler angles E1 and E2 are those returned
by TERM_REVIEW,/SET - need to convert from degrees)
INPUTS:
angle [eX,eY], in degrees (note: this is a low-level
program! These target angles are the Euler angles, not the
X,Y displacements.)
og_number, the (decimal) OG number that the new commands are for
OPTIONAL (KEYWORD) INPUT PARAMETERS:
offset is the nominal pointing, default = [0.0003889,0.002111]
normal is the Euler angles of the S/C NORMAL POINT command, default
at [-.0797,.0295]
name forces the printed output of a suggested name for the OG
command, under the convention (incremental displacements in
arc min relative to the NORMAL POINT location)
OUTPUTS:
bc: an array of 8 hexadecimal numbers for entry into an OG
CALLS: ***
TO_HEX
CALLED BY:
plot2offset [1], plot2offset [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
does not create a Z-angle command (hard-wired to 00 00)
Errors as large as 2 arcsec in X (8 LSB units)
and 10 arcsec in Y (50 LSB units) have been found
in comparison with 6 previous cases.
The code has built-in limit switches to help prevent the generation
of large offpoints
WARNING:
The Euler angles are defined such that X increasing => S/C south
Y increasing => S/C east
MODIFICATION HISTORY:
HSH, written 13-Feb-1995
LWA, corrected and checked 15-Feb-95
HSH, added /name and /normal keywords, limit switches 5-Jun-95
LWA, Changed naming format to include decimals. 13-Nov-95
HSH, added print statements to match OG editor format;
added og_number input variable. 13-Nov-95
LWA, OG_NUMBER feature disabled. 13-Nov-95 (!)
SKS, added if & format statements to adjust name_x/y output. 13-Nov-95
HSH, added another example to the calling sequence. 6-aug-96
HSH, improved the documentation a bit. 22-Nov-98
HSH, added text output to hardwired directory, 19-Nov-99
HSH, added nofile to defeat text, 7-Jun-00
AkT, improved to get better approx. of dx(y)_nums from disp_x(y),
otherwise sometimes disagree with Seko-san's check, 13-May-01
[Previous]
[Next]
NAME:
POINT_HEX2
PURPOSE:
Generate an OG for a specified pointing location. The output
is in hexadecimal, appropriate for OG editing via OPOGEDT.
CATEGORY:
Yohkoh
CALLING SEQUENCE:
point_hex2, angle, og_number, bc, offset = offset
point_hex2, [-.0797+E1/3600.,.0295+E2/3600.], /name
(here the Euler angles E1 and E2 are those returned
by TERM_REVIEW,/SET - need to convert from degrees)
INPUTS:
angle [eX,eY], in degrees (note: this is a low-level
program! These target angles are the Euler angles, not the
X,Y displacements.)
og_number, the (decimal) OG number that the new commands are for
OPTIONAL (KEYWORD) INPUT PARAMETERS:
offset is the nominal pointing, default = [0.0003889,0.002111]
normal is the Euler angles of the S/C NORMAL POINT command, default
at [-.0797,.0295]
name forces the printed output of a suggested name for the OG
command, under the convention (incremental displacements in
arc min relative to the NORMAL POINT location)
OUTPUTS:
bc: an array of 8 hexadecimal numbers for entry into an OG
CALLS: ***
TO_HEX
CALLED BY:
plot2offset [1], plot2offset [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
does not create a Z-angle command (hard-wired to 00 00)
Errors as large as 2 arcsec in X (8 LSB units)
and 10 arcsec in Y (50 LSB units) have been found
in comparison with 6 previous cases.
The code has built-in limit switches to help prevent the generation
of large offpoints
WARNING:
The Euler angles are defined such that X increasing => S/C south
Y increasing => S/C east
MODIFICATION HISTORY:
HSH, written 13-Feb-1995
LWA, corrected and checked 15-Feb-95
HSH, added /name and /normal keywords, limit switches 5-Jun-95
LWA, Changed naming format to include decimals. 13-Nov-95
HSH, added print statements to match OG editor format;
added og_number input variable. 13-Nov-95
LWA, OG_NUMBER feature disabled. 13-Nov-95 (!)
SKS, added if & format statements to adjust name_x/y output. 13-Nov-95
HSH, added another example to the calling sequence. 6-aug-96
HSH, improved the documentation a bit. 22-Nov-98
HSH, added text output to hardwired directory, 19-Nov-99
HSH, added nofile to defeat text, 7-Jun-00
[Previous]
[Next]
NAME:
POLY_MOVIE
PURPOSE:
performs the poly_2d function on a movie (3-D array)
CATEGORY:
CALLING SEQUENCE:
out = poly_movie(in,pp)
INPUTS:
in = byte(*,*,*). Input movie.
pp = float(2,2,*). Coordinate pairs for coregistration. Created
by setpts procedure. Optional.
KEYWORDS (INPUT):
pp = float(2,2,*). Same as positional parameter pp.
The pp keyword is an alternate way of supplying pp,
mostly for consistency of syntax with other programs.
The pp keyword is ignored if the positional parameter is present.
key = transformation type. Used in conjunction with pp.
See "transtype" keyword in CALTRANS for details.
Default is a general linear transform.
index = vector of frame numbers to be used for output movie.
xy = float(2,*). x/y offsets for each frame of input movie in units
of input pixels. Also indexed by "index", just as input movie.
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.
/nointerp = Nearest neighbor instead of bilinear interpolation.
The program then corrects for the 0.5 pixel offset to keep the
non-interpolated images co-aligned with the interpolated ones.
missing = scalar. Intensity to be used for missing data points.
OUTPUTS:
out = output movie.
KEYWORDS (OUTPUT):
CALLS: ***
CALTRANS
CALLED BY:
mk_coal_movie
COMMON BLOCKS:
None.
SIDE EFFECTS:
RESTRICTIONS:
Movies that show several spectral channels (option of rmovie)
cannot be processed directly. They have to be split up into
simple movies beforehand.
PROCEDURE:
MODIFICATION HISTORY:
JPW, Jan. 1992
JPW, 21-jun-94 added missing keyword
JPW, 21-jun-94 pp also as keyword now.
JPW, 22-nov-94 added "key" keyword
[Previous]
[Next]
NAME:
POSN_PLOT
PURPOSE:
CATEGORY:
CALLING SEQUENCE:
.run posn_plot
INPUTS:
Read files "posn_log.dat" on current directory
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
VMS specific at the moment
PROCEDURE:
MODIFICATION HISTORY:
RDB 22-Mar-94 Written
[Previous]
[Next]
NAME:
PP_MWLT2SXT
PURPOSE:
Transforms a pp vector pointing to MWLT data into a pp vector
pointing to SXT data.
CATEGORY:
CALLING SEQUENCE:
pp_gbo_sxt = pp_mwlt2sxt(pp_gbo_mwlt,mwlt_data,sxt_index)
INPUTS:
pp_gbo_sxt = point array for transforming ground based data onto
the grid of Mees White Light Telescope images.
mwlt_data = MWLT image, or alternatively, 5-element vector containing
[xcen,ycen,rad,oblat,ob_angle], e.g., from find_limb.
sxt_index = index structure of the reference SXT image.
KEYWORDS (INPUT):
time_mwlt = time of MWLT image, either in any Yohkoh format,
or as a time string 'hh:mm:ss'. In this latter case,
the mwlt_date keyword must also be present.
The time information is used to correct for solar
rotation between the time of the MWLT image and the
Yohkoh image.
Default: time of Yohkoh image.
date_mwlt = date string 'dd-mm-yy'. Use only if mwlt_time given,
but NOT in a Yohkoh format.
roll_mwlt = roll angle of MWLT image. CCD N (up) measured eastward
from geocentric N (in radian).
If not supplied, then the MWLT roll database is read.
/xreverse = flag for mirrored image (true = reversed).
Note: raw MWLT images are reversed (N=up, W=left),
JPEG compressed MWLT images are not (N=up, W=right).
Examples: non-reversed: rollm=0 -> N up, W right
rollm=pi/2 -> E up, N right
reversed: rollm=0 -> N up, W left
rollm=pi/2 -> E up, N left
k_in = input "key" parameter, setting constraints for transformation.
See the description of the "transtype" keyword in CALTRANS.
/old = flag for data from old MWLT telescope. If not supplied, then
the info from the roll database is used.
/verbose = flag to print some limb fitting parameters.
OUTPUTS:
pp_gbo_sxt = point array for transforming ground based data onto
the grid of SXT images.
KEYWORDS (OUTPUT):
k_out = output "key" parameter. k_in is modified for use with SXT.
CALLS: ***
CONV_A2H [1], CONV_A2H [2], CONV_H2A [1], CONV_H2A [2], FITIM2INTS, MWLT_CENTER
anytim2ints [1], anytim2ints [2], get_mwlt_roll, get_rb0p [1], get_rb0p [2]
get_roll, gt_corner, gt_pix_size, gt_res, int2secarr [1], int2secarr [2]
sxt_cen [1], sxt_cen [2]
COMMON BLOCKS:
None.
PROCEDURE:
RESTRICTIONS:
Image must be dark outside the solar limb (no annotations!)
MODIFICATION HISTORY:
JPW, 25-mar-94 for mirrored (E-W reversed) MWLT images
JPW, 17-nov-94 program now assumes non-mirrored image
added keyword /xreverse
JPW, 25-nov-94 added k_in, k_out
JPW, 2-dec-94 fixed the default pp vector
[Previous]
[Next]
NAME:
PP_NABAN2SXT
PURPOSE:
Adds the pointing offset between the SXT and the Aspect Telescope
to a transformation point array pp.
CATEGORY:
CALLING SEQUENCE:
pp_gbo_sxt = pp_naban2sxt(pp_gbo_naban)
INPUTS:
pp_gbo_naban = point array for transforming ground based data onto
the grid of aspect telescope images.
KEYWORDS (INPUT):
resolution = 1, 2, or 4. Binning factor of the SXT/NaBan images.
1=full, 2=half, 4=quarter resolution. Default=1.
sxt2gbo = 1. Set flag if input pp array represents a transformation
from the aspect telescope grid to the GBO (ground based)
grid (instead of GBO grid to aspect grid).
The pointing offset is then applied so that the output
pp array properly represents the transformation from the
SXT grid to the GBO grid.
OUTPUTS:
pp_gbo_naban = point array for transforming ground based data onto
the grid of SXT soft X-ray images.
COMMON BLOCKS:
None.
PROCEDURE:
MODIFICATION HISTORY:
JPW, May 1992.
JPW, Aug 92. Added keywords "resolution" and "sxt2gbo".
[Previous]
[Next]
NAME:
pp_reverse
PURPOSE:
Calculate the "pp" vector and the "transtype" for the inverse
transform
CATEGORY:
CALLING SEQUENCE:
pp_reverse,pp_in,pp_out
pp_reverse,pp_in,pp_out,ttyp_in,ttyp_out
INPUTS:
pp_in = fltarr(2,2,*). Point array defining transform.
ttyp_in = string, float element, or structure. Sets constraints.
See "transtype" keyword in CALTRANS for details.
KEYWORDS (INPUT):
OUTPUTS:
pp_out = fltarr(2,2,*). Point array defining inverse transform.
ttyp_out = string or structure. Constraints for inverse transf.
KEYWORDS (OUTPUT):
CALLED BY:
OPLOT_B
COMMON BLOCKS:
None.
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
JPW, 22-nov-94
[Previous]
[Next]
NAME:
PP_SXT2ANY
PURPOSE:
Transforms a pp vector pointing to SXT data into a pp vector
pointing to any image with given pointing and scaling info.
CATEGORY:
CALLING SEQUENCE:
pp_in_any = pp_sxt2any(pp_in_sxt,sxtindex=index,pnt=pnt)
INPUTS:
pp_in_sxt = point array for transforming input data onto
the grid of SXT soft X-ray images. If undefined or
an element variable, then input=sxt is assumed.
KEYWORDS (INPUT):
sxtindex = index structure of the reference SXT image. REQUIRED!
sizpix = (x,y). target image pixel size in arcseconds. default (1,1).
rotang = pos. angle of solar N measured eastward from image N
in degrees. Default=0
/geo : rotang relative to geographic N instead of solar N.
/geo without rotang keyword sets rotang = P_angle.
refpix = reference pixel (x,y), first pixel of array = 0,0 = default!
pnt = solar coordinates for reference pixel (ew,ns),
in arcseconds from suncenter, W and N positive.
/helio : reference pixel coordinates in heliographic system,
instead of angular offsets from suncenter
OUTPUTS:
pp_in_any = point array for transforming input data onto
the target grid.
CALLS: ***
CONV_H2A [1], CONV_H2A [2], get_rb0p [1], get_rb0p [2], get_roll, gt_corner
gt_pix_size, gt_res, sxt_cen [1], sxt_cen [2]
COMMON BLOCKS:
None.
PROCEDURE:
RESTRICTIONS:
Change to work for half res. data not well tested
MODIFICATION HISTORY:
JPW, Dec. 1993.
JPW, Nov. 1996, derived from pp_sxt2vla, not much testing
[Previous]
[Next]
NAME:
PP_SXT2VLA
PURPOSE:
Transforms a pp vector pointing to SXT data into a pp vector
pointing to VLA data.
CATEGORY:
CALLING SEQUENCE:
pp_gbo_sxt = pp_sxt2vla(pp_gbo_sxt,sxt_index=index, $
sun_radec=s_radec,vla_radec=v_radec,vla_pixref=v_pix)
INPUTS:
pp_gbo_sxt = point array for transforming ground based data onto
the grid of SXT soft X-ray images.
KEYWORDS (INPUT):
sxt_index = index structure of the reference SXT image.
sun_radec = celest. coord. of suncenter (['hh:mm:ss.sss','dd:mm:ss'])
vla_radec = celestial coord. of VLA ref. pixel
vla_pixref = VLA reference pixel (x,y), first pixel of array = 0,0 !
vla_spix = VLA pixel size in arcseconds
vla_dpix = center of de-rotation back to solar north.
OUTPUTS:
pp_gbo_vla = point array for transforming ground based data onto
the grid of VLA images.
CALLS: ***
anytim2ints [1], anytim2ints [2], get_rb0p [1], get_rb0p [2], get_roll, gt_corner
gt_pix_size, gt_res, sxt_cen [1], sxt_cen [2]
COMMON BLOCKS:
None.
PROCEDURE:
RESTRICTIONS:
Solar and VLA coordinates must both be in the same system,
e.g., apparent geocentric coordinates.
Change to work for half res. data not tested yet!
Does not correct for diurnal parallax. Is usually taken
into account by the VLA (output values are geocentric).
MODIFICATION HISTORY:
JPW, Dec. 1993.
[Previous]
[Next]
NAME:
PPHA
PURPOSE:
Extract BCS PHA data from the DP_SYNC data stream
plots PHA spectra of selected time interval
CALLING SEQUENCE:
.run PPHA
INPUTS:
OUTPUTS:
Plot of PHA spectra
OPTIONAL OUTPUTS:
hardcopy is optional
RESTRICTIONS:
Only works on and X-terminal
This programme should be run after PP_CARTEX or PLOT_ERATE
PROCEDURE:
The user is allowed to select a region from the existing
light curve. The PHA data corresponding to this interval
is extracted from the DP_SYNC data and is thrown up on the
screen for the user to further select the required time interval
The PHA spectra of the sum of this period is plotted and can
optionally be sent to the hardcopy.
MODIFICATION HISTORY:
RDB 1991/92
[Previous]
[Next]
NAME: Ppha.pro
PURPOSE: plot a contour map of pha data, sum and display pha data
optionally fit and append results to file in readout format.
CALLING SEQUENCE: .run test_rd (-888,y) |.run plot_erate |.run ppha
INPUTS: structures from test_rd
OUTPUTS: plots of pha,
OPTIONAL OUTPUTS: appends to output.dat results of fit.
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
written -- rdb 1991
modified for fitting -- atp 7/3/1992
Classification: BDA, Housekeeping, PHA, Spikes (fixed patterning) correction
[Previous]
[Next]
NAME: PR_CAN_PRELIM
PURPOSE: To produce the Preliminary pass cancellation form
CALLING SEQUENCE: IDL> pr_can_prelim, pass [,sheet,tot_sheet,file=file]
INPUTS: pass The string array containing details
of the passes taken from the strawman
file for the relevant week
OPTIONAL INPUTS: file if specified the Postscript graphics
are dumped to that file, otherwise
default is IDL.PS
sheet the number of the current sheet
tot_sheet the total numbers of sheets
flare1 If set, try to print to Flare1 printer
OUTPUTS: none
OPTIONAL OUTPUTS: none
CALLS: ***
PRINT_FLARE1, pprint [1], pprint [2]
CALLED BY:
CANCEL_PRELIM
RESTRICTIONS: none
PROCEDURE: The procedure formats a sheet into the form of the
JPL change request form and adds the pass information
where appropriate.
MODIFICATION HISTORY: Written 2-Feb-93 CDP
5-Feb-93 JRL Call print_flare1 to print plot
[Previous]
[Next]
pro diff_angle,yazim,yelev,aazim,aelev,ang
f=!pi/180.
yx=cos(f*yelev(0))*cos(f*yazim(0))
yy=cos(f*yelev(0))*sin(f*yazim(0))
yz=sin(f*yelev(0))
ax=cos(f*aelev(0))*cos(f*aazim(0))
ay=cos(f*aelev(0))*sin(f*aazim(0))
az=sin(f*aelev(0))
ang=acos(yx*ax+yy*ay+yz*az)/f ; unit:(degree)
end
CALLS:
[Previous]
[Next]
NAME: pr_dark_time
CALLING SEQUENCE: y=pr_leak_time(index)
INPUT: index - SXT index with the his structure indicating processed data.
can be a vector.
OUTPUT: time(s) (string(s)) in the Yohkoh style, e.g. '21-AUG-92 10:30:00'
CALLS: ***
GET_INFO_DARK, fmt_tim [1], fmt_tim [2]
HISTORY: Written by Nariaki Nitta, 22-NOV-97
[Previous]
[Next]
NAME:
pr_dates_warm
PURPOSE:
To list the dates that SXT CCD was warm
SAMPLE CALLING SEQUENCE:
pr_dates_warm, ssl
pr_dates_warm, ssl, outfil='dum.dum'
OPTIONAL INPUT/OUTPUT:
ssl - The SSL structure. If undefined, it calls RD_SSL and
reads it. It is an optional parameter since it takes so
long to read it, if you want to call the routine twice,
it won't have to re-read the SSL files.
CALLS: ***
gt_day [1], gt_day [2], gt_temp_ccd, rd_ssl
OPTIONAL KEYWORD INPUT:
outfil - An output file name
HISTORY
Written 16-Dec-93 by M.Morrison
22-Dec-93 (MDM) - Added guide option
[Previous]
[Next]
NAME:
pr_evn
PURPOSE:
To list the times when there is Yohkoh Data.
CALLING SEQUENCE:
pr_evn ;contacts for today
pr_evn, '1-jan', /flare
pr_evn, '1-jan', /flare, mindur=2
pr_evn, '30-oct-91', '2-nov-91', /hc
INPUT:
st_date - The JST date to search for. If no date is entered,
the current date is used. The year is not
required, but if entered should be of the
form 1-Oct-91 style. The date should be
entered in string style with date first.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
CALLS: ***
CHECK_TIME, DPRINT, Ex2Int [1], Ex2Int [2], Int2Ex [1], Int2Ex [2], RD_EVN, RD_GBE
RD_GEV, anytim2ex [1], anytim2ex [2], anytim2ints [1], anytim2ints [2]
fmt_tim [1], fmt_tim [2], gt_dp_mode [1], gt_dp_mode [2], gt_time [1], gt_time [2]
int2secarr [1], int2secarr [2], tag_index [1], tag_index [2], tim2orbit [1]
tim2orbit [2]
CALLED BY:
FIRST_LIGHT [1], FIRST_LIGHT [2]
OPTIONAL KEYWORD INPUT:
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
flare - If set, only display the flare modes
mindur - Only list events which have a duration longer than
this value (in minutes)
counts - If set, print the count rates instead of # of datasets
allgoes - If set, list all goes events, even if there was not a
Yohkoh event
gro - If set, list the BATSE GRO events
OUTPUT:
The times of the S/C start of day and start of night.
Also lists the station contacts during the orbit (not the times)
and whether there is an SAA passage. List the FILEID also.
HISTORY:
Written 1-May-92 by M.Morrison
6-Jul-92 (MDM) - Added "allgoes" option
28-Aug-92 (MDM) - Added INFIL option
29-Oct-92 (MDM) - Adjusted the header to be aligned
5-Jan-93 (MDM) - Corrected error of when using /ALLGOES option to
print the proper DP mode when not in flare mode
21-Jan-93 (MDM) - Added the /GRO option
6-Apr-93 (MDM/SLF) - Modified call to RD_EVN to not use the "flag"
keyword - use "status" instead
14-May-93 (MDM) - Added INDIR option (removed INFIL option)
27-May-93 (MDM) - Modified to accept any time format as input and
to use the hours of the input as requested
18-Jun-93 (MDM) - Modified to use RD_GBE instead of RD_BATSE
[Previous]
[Next]
NAME:
pr_fem
PURPOSE:
To list the times of the S/C day and night
CALLING SEQUENCE:
pr_fem ;contacts for today
pr_fem, '1-jan
pr_fem, '30-oct-91', '2-nov-91', /hc
INPUT:
st_date - The JST date to search for. If no date is entered,
the current date is used. The year is not
required, but if entered should be of the
form 1-Oct-91 style. The date should be
entered in string style with date first.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
CALLS: ***
CHECK_TIME, DPRINT, Ex2Int [1], Ex2Int [2], Int2Ex [1], Int2Ex [2], anytim2ints [1]
anytim2ints [2], fmt_tim [1], fmt_tim [2], gt_time [1], gt_time [2], rd_fem
timstr2ex [1], timstr2ex [2]
OPTIONAL KEYWORD INPUT:
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
list_saa- If set, list the times that the S/C is in SAA.
asca - If set, use the ASCA spacecraft FEM files
OUTPUT:
The times of the S/C start of day and start of night.
Also lists the station contacts during the orbit (not the times)
and whether there is an SAA passage. List the FILEID also.
HISTORY:
Written 23-Apr-92 by M.Morrison
15-Jan-93 (MDM) - Added LIST_SAA option
7-Apr-93 (MDM) - Added ASCA option
[Previous]
[Next]
NAME:
pr_gbe
PURPOSE:
To list the GRO/BATSE events
CALLING SEQUENCE:
pr_gbe, '1-jan
pr_gbe, '30-oct-91', '2-nov-91', /hc
pr_gbe, event=3944
INPUT:
st_date - The JST date to search for. If no date is entered,
the current date is used. The year is not
required, but if entered should be of the
form 1-Oct-91 style. The date should be
entered in string style with date first.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
CALLS: ***
DPRINT, Ex2Int [1], Ex2Int [2], Int2Ex [1], Int2Ex [2], RD_GBE, anytim2ints [1]
anytim2ints [2], fmt_tim [1], fmt_tim [2], gt_time [1], gt_time [2], timstr2ex [1]
timstr2ex [2]
OPTIONAL KEYWORD INPUT:
above - If passed, then only list events that have a counting
rate above this value
event - If passed, only list event with this BATSE event number.
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
OUTPUT:
The times and classes of the GOES events
HISTORY:
Written 13-Apr-93 by M.Morrison
RDB 6-Aug-93 Added /event keyword
[Previous]
[Next]
NAME:
pr_gev
PURPOSE:
To list the GOES events
CALLING SEQUENCE:
pr_gev, '1-jan
pr_gev, '30-oct-91', '2-nov-91', /hc
INPUT:
st_date - The JST date to search for. If no date is entered,
the current date is used. The year is not
required, but if entered should be of the
form 1-Oct-91 style. The date should be
entered in string style with date first.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
CALLS: ***
DPRINT, Ex2Int [1], Ex2Int [2], Int2Ex [1], Int2Ex [2], RD_GEV, anytim2ints [1]
anytim2ints [2], fmt_tim [1], fmt_tim [2], gt_time [1], gt_time [2], timstr2ex [1]
timstr2ex [2]
CALLED BY:
BCS_CREATE_CAT, chk_flares [1], chk_flares [2]
OPTIONAL KEYWORD INPUT:
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
OUTPUT:
The times and classes of the GOES events
HISTORY:
Written 6-Jun-92 by M.Morrison
12-Apr-93 (MDM) - Used STATUS keyword with RD_GEV to see if any
events are available for the time period
16-May-93 (MDM) - Modified format in order to fit on 80 column terminal
14-May-93 (HSH) - Added the GEV_DATA output and put a /quiet keyword
28-Jun-2002 - S.L.Freeland - tweak justify for records with NOAA but no
location (this routine could use a rewrite using
gt_tagval and strjustify, probably becomes a five
liner vectorized (no for loop), but I'm busy right now...
27-oct-2004 - S.L.Freeland - allow five digit NOAA AR (finally...)
[Previous]
[Next]
NAME:
PR_GRPPLN
PURPOSE:
Splits BCS grouper plan into its components and prints it
CATEGORY:
CALLING SEQUENCE:
pr_grppln,plans,jp
INPUTS:
plans grouper plans structure
jp plan to print
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLED BY:
pr_plan_area [1], pr_plan_area [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB 5-May-94 Written
RDB 17-May-94 Changed .grpid to .planid for reformatter compatability
[Previous]
[Next]
NAME:
pr_gsn
PURPOSE:
To display the information in the SOON observing log files
SAMPLE CALLING SEQUENCE:
pr_gsn,'/ydb/gsn/sth940517.cat', /summary
pr_gsn,'/ydb/gsn/sth940517.cat', array, /summary
INPUT:
infil - the input file
OPTIONAL OUTPUT:
array - Instead of printing the listing information, return it
as an array
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], CONCAT_DIR [1], CONCAT_DIR [2]
CONCAT_DIR [3], CONV_A2H [1], CONV_A2H [2], DPRINT, anytim2ints [1]
anytim2ints [2], break_file [4], concat_dir [4], data_type [1], data_type [2]
fmt_tim [1], fmt_tim [2], rd_gsn
CALLED BY:
soon_search [1], soon_search [3]
OPTIONAL KEYWORD INPUT:
outfil - If defined, save results to this file
summary - If set, display a summary line for each tape file
hc - If set, send the output to the printer
tapefile - If set, then do the per image listing only on the
tapefile passed in.
Optionally, tapefile can be a string array of the sort:
' sth940517 7 '
which is from the PR_GSN array output
HISTORY:
Written Dec-93 by M.Morrison
28-Jul-94 (MDM) - Modified summary output to include the image types
3-Aug-94 (MDM) - Added TAPEFILE option
5-Aug-94 (MDM) - Added GSN_HEAD option
- Added COMMENTS option
8-Aug-94 (MDM) - Modified to display deviation from reference location
in arcminutes instead of arcsec
23-Sep-94 (MDM) - Added ARRAY output option
- Expanded the TAPEFILE option to be able to pass a
string which has the tape name and file number in
a string
- Changed RD_GSN call to not read obs if doing
summary listing
2-Aug-95 (MDM) - Corrected how image listing works for multiple
file inputs
[Previous]
[Next]
NAME:
pr_his_index
PURPOSE:
To print out the information from the SXT history index
SAMPLE CALLING SEQUENCE:
pr_his_index, index
pr_his_index, index, /longform
INPUT:
index - The index
CALLS: ***
ARR2STR [1], Arr2Str [2], GET_INFO2, HIS_EXIST, MORE [1], MORE [2], STR2ARR [1]
STR2ARR [2], anytim [1], anytim [2], anytim [3], anytim [4], anytim [5]
anytim2ints [1], anytim2ints [2], fmt_tim [1], fmt_tim [2], gt_expdur [1]
gt_expdur [2], gt_res, gt_tagval [1], gt_tagval [2], rd_obs, ssw_deltat
strarrcompress, strjustify
CALLED BY:
convert YohkohSXT history record more legible form, sxt_orhis
OPTIONAL KEYWORD INPUT:
longform - If set, then give full GET_INFO2 information on
each of the dark and leak images used
Keyword OUTPUT -
The following keywords are OUTPUT booleans which indicate whether
implied processing was (1) or was not (0) done for each input index:
dark, leak, synleak, second, despiked, destreaked, vignette, descatter
registered, ssc_norm res_norm
ncomposite (output) - number of images if composite
only_flags - (input switch) - if set, dont bother which text or tty
output - ie, only output state flags are of interest
(will run faster for large number of input index records)
Restriction : This version adds the following line without an exception,
for the Yohkoh Legacy data;
'Intensity units: DN/Image-Pixel/Exposure-Duration'.
HISTORY:
Written 1994 by M.Morrison
20-Mar-95 (MDM) - Modified to report information on Vignette
~1-Mar-99 S.L.Freeland - add despike/destreak, some formatting
22-sep-99 S.L.Freeland - add DC ORBIT CORRECT check
[high bit = 2^15 in index.HIS.q_dark_sub
23-sep-99 S.L.Freeland - Colum Blem history check
[low bit = 1 in index.HIS.q_extra(1,nn) ]
10-Jan-2000 S.L.Freeland - 2nd order leak correction check
16-jan-2004 S.L.Freeland - Scattering his (sxt_decon...)
14-Jul-2004 S.L.Freeland - exposure normalization flg/algorithm
29-Sep-2004 S.L.Freeland - externalize info via keyword flags
and add /ONLY_FLAGS keyword and funct
(for use by sxt_index2filename for
example)s
5-oct-2004 S.L.Freeland - enable leak type (synthetic vs terminator)
13-dec-2004 S.L.Freeland - check/verbalize LOW8 restor
20-oct-2005 Aki Takeda - omit some blanks in output strings, to avoid
truncated info found in FITS header.
21-oct-2005 Aki Takeda - carry milisec info for SSC history.
15-dec-2005 Aki Takeda - omit 'Exposure normalization' information.
- add unit in the 'Exposure duration' line
- add a new line, 'Intensity units: .....'
22-jan-2006 Aki Takeda - modify 'Intensity units' proper for movies.
[Previous]
[Next]
NAME:
pr_his_index
PURPOSE:
To print out the information from the SXT history index
SAMPLE CALLING SEQUENCE:
pr_his_index, index
pr_his_index, index, /longform
INPUT:
index - The index
CALLS: ***
ARR2STR [1], Arr2Str [2], GET_INFO2, HIS_EXIST, MORE [1], MORE [2], STR2ARR [1]
STR2ARR [2], anytim [1], anytim [2], anytim [3], anytim [4], anytim [5]
anytim2ints [1], anytim2ints [2], fmt_tim [1], fmt_tim [2], gt_expdur [1]
gt_expdur [2], gt_res, gt_tagval [1], gt_tagval [2], rd_obs, ssw_deltat
strarrcompress, strjustify
CALLED BY:
convert YohkohSXT history record more legible form, sxt_orhis
OPTIONAL KEYWORD INPUT:
longform - If set, then give full GET_INFO2 information on
each of the dark and leak images used
HISTORY:
Written 1994 by M.Morrison
20-Mar-95 (MDM) - Modified to report information on Vignette
~1-Mar-99 S.L.Freeland - add despike/destreak, some formatting
22-sep-99 S.L.Freeland - add DC ORBIT CORRECT check
[high bit = 2^15 in index.HIS.q_dark_sub
23-sep-99 S.L.Freeland - Colum Blem history check
[low bit = 1 in index.HIS.q_extra(1,nn) ]
10-Jan-2000 S.L.Freeland - 2nd order leak correction check
16-jan-2004 S.L.Freeland - Scattering his (sxt_decon...)
[Previous]
[Next]
NAME:
pr_his_index
PURPOSE:
To print out the information from the SXT history index
SAMPLE CALLING SEQUENCE:
pr_his_index, index
pr_his_index, index, /longform
INPUT:
index - The index
CALLS: ***
ARR2STR [1], Arr2Str [2], GET_INFO2, HIS_EXIST, MORE [1], MORE [2], STR2ARR [1]
STR2ARR [2], anytim [1], anytim [2], anytim [3], anytim [4], anytim [5]
anytim2ints [1], anytim2ints [2], fmt_tim [1], fmt_tim [2], gt_expdur [1]
gt_expdur [2], gt_res, gt_tagval [1], gt_tagval [2], rd_obs, ssw_deltat
strarrcompress, strjustify
CALLED BY:
convert YohkohSXT history record more legible form, sxt_orhis
OPTIONAL KEYWORD INPUT:
longform - If set, then give full GET_INFO2 information on
each of the dark and leak images used
HISTORY:
Written 1994 by M.Morrison
20-Mar-95 (MDM) - Modified to report information on Vignette
[Previous]
[Next]
NAME:
pr_image
PURPOSE:
To print the average DN level of a box to be marked by the
user. It can also display the raw values.
CALLING SEQUENCE:
pr_image, data
pr_image, data, 10
pr_image, data, 20, /nomark, /raw, bin=4
pr_image, data, index=index
INPUT:
data - 2-D data image array
OPTIONAL INPUT:
side - Number of pixels on a side (if not passed, the
default is 9
nomark - If present, the location on the image which
is selected is NOT marked (default is to mark it)
raw - If present, print out the raw data values also
bin - The rebin size used to display the image
index - The index structure that goes with the data image.
CALLS: ***
draw_boxcorn [1], draw_boxcorn [2], get_pix_coor
HISTORY:
Written 20-Nov-91 by M.Morrison
7-May-92 (MDM) - Added capability to print out the CCD address
(need to pass in the index for that to work)
15-Oct-92 (MDM) - Modified call to DRAW_BOXCORN to add 1 to x1 and y1
to encircle the selected pixels more properly
- Corrected the sub-image selected when "side" was
an even number (it used to give 17x17 when side
was set to 16)
[Previous]
[Next]
NAME:
pr_nar
PURPOSE:
To list the NOAA Active Region list
CALLING SEQUENCE:
pr_nar, '1-jan
pr_nar, '30-oct-91', '2-nov-91', /hc
INPUT:
st_date - The JST date to search for. If no date is entered,
the current date is used. The year is not
required, but if entered should be of the
form 1-Oct-91 style. The date should be
entered in string style with date first.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
CALLS: ***
DPRINT, RD_NAR, anytim2ex [1], anytim2ex [2], anytim2ints [1], anytim2ints [2]
fmt_tim [1], fmt_tim [2], gt_day [1], gt_day [2], timstr2ex [1], timstr2ex [2]
OPTIONAL KEYWORD INPUT:
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
noprint - If set, do not print the results to the screen
nearest - Passed on to RD_NAR
noaa - List only data for a specified NOAA region number
OPTIONAL KEYWORD OUTPUT:
outstring- The results of the NAR listing in a string array
OUTPUT:
The times and classes of the GOES events
HISTORY:
Written 28-Jan-93 by M.Morrison
8-Apr-93 (MDM) - Added OUTSTRING and NOPRINT options
- Added NEAREST option
17-Aug-00 (HSH) - Added NOAA keyword
[Previous]
[Next]
NAME:
pr_nel
PURPOSE:
To list the Nobeyama Radio Observatory events
CALLING SEQUENCE:
pr_nel, '1-jan
pr_nel, '30-oct-91', '2-nov-91', /hc
INPUT:
st_date - The UT date to search for.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
CALLS: ***
DPRINT, Ex2Int [1], Ex2Int [2], Int2Ex [1], Int2Ex [2], RD_NEL, anytim2ints [1]
anytim2ints [2], fmt_tim [1], fmt_tim [2], timstr2ex [1], timstr2ex [2]
OPTIONAL KEYWORD INPUT:
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
OUTPUT:
The times and classes of the GOES events
HISTORY:
Written 19-Apr-93 by M.Morrison
[Previous]
[Next]
NAME:
PR_NODE
PURPOSE:
Returns information about the specified node.
CALLING SEQUENCE:
pr_node,nodename
INPUTS:
nodename name of node about which information is needed
OUTPUTS:
prints to screen
OPTIONAL OUTPUTS:
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], concat_dir [4]
RESTRICTIONS:
PROCEDURE:
Does a search or grep (depending on o/s) on the file
$DIR_BCS_DOC/node_list.dat
MODIFICATION HISTORY:
RDB 02-Jul-92 Written
RDB 21-Jul-92 UNIX check and fix
[Previous]
[Next]
NAME:
pr_orbit_sum
PURPOSE:
Generate a listing of the predicted Yohkoh daylight times.
Can optionally exclude SAA times or provide a full listing showing
day/night, SAA and station contacts.
CALLING SEQUENCE:
pr_orbit_sum,orbit,outfile,[/full,orbit_params=orbit_params,/nosaa,/echo]
INPUTS:
orbit = Strcture of predicted orbital events.
outfil = File name of output file.
OPTIONAL INPUT PARAMETERS:
echo If set to 1, echo the listing to the terminal with /full.
full If set to 1, will only print day/night, SAA, station contacts.
nosaa If set to 1 and not /full, exclude times with SAA.
No effect if /full switch is set.
OUTPUTS:
None.
CALLS: ***
fmt_tim [1], fmt_tim [2]
CALLED BY:
mk_orbit_sum
MODIFICATION HISTORY:
5-oct-91, written, MDM/JRL
26-oct-91, updated, JRL -- Changed timarr to standard format
(h,m,s,ms,d,m,y)
6-mar-92, Updated, JRL -- Changed the format of the output
[Previous]
[Next]
NAME:
PR_PATROLS
PURPOSE:
print details of the patrol images for a given time
CATEGORY:
CALLING SEQUENCE:
pr_patrols, t1
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
CALLS: ***
anytim [1], anytim [2], anytim [3], anytim [4], anytim [5], gtab_patrols
tim2dset [1], tim2dset [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
19-sep-01, written by HSH
[Previous]
[Next]
NAME:
pr_plan_area
PURPOSE:
Unplacks the bytes in the BCS plans parameter area
CATEGORY:
CALLING SEQUENCE:
pr_plan_area,plan_area [,default_area]
INPUTS:
plan area plane parameter area (512 bytes long, starting 3800h)
OPTIONAL INPUT PARAMETERS:
default area area that hold default sequeneces 0,1 and flops 0,1
KEYWORD PARAMETERS:
full breaks mode lists within sequence down
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLS: ***
ALL_VALS [1], ALL_VALS [2], PR_GRPPLN, RD_BCS_GRPDTB, iword
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB 11-May-94 Written, based on pre-launch routine
rdb 21-Nov-97 More obvious flag if default are not supplied
[Previous]
[Next]
NAME:
pr_plan_area
PURPOSE:
Unplacks the bytes in the BCS plans parameter area
CATEGORY:
CALLING SEQUENCE:
pr_plan_area,plan_area [,default_area]
INPUTS:
plan area plane parameter area (512 bytes long, starting 3800h)
OPTIONAL INPUT PARAMETERS:
default area area that hold default sequeneces 0,1 and flops 0,1
KEYWORD PARAMETERS:
full breaks mode lists within sequence down
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
CALLS: ***
ALL_VALS [1], ALL_VALS [2], PR_GRPPLN, RD_BCS_GRPDTB, iword
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB 11-May-94 Written, based on pre-launch routine
[Previous]
[Next]
NAME:
pr_pnt_hist
PURPOSE:
To list the Yohkoh S/C commanded pointing history
SAMPLE CALLING SEQUENCE:
pr_pnt_hist
pr_pnt_hist, outfil='pr_pnt_hist.txt'
pr_pnt_hist, /outfil
CALLS: ***
fmt_tim [1], fmt_tim [2], int2secarr [1], int2secarr [2], prstr [1], prstr [2]
rd_pnt_hist, sel_timrange [1], sel_timrange [2]
OPTIONAL KEYWORD INPUT:
outfil - The name of the output file to be written to. If it
is an integer, it will use the file 'pr_pnt_hist.txt'
HISTORY:
Written Jul-93 by M.Morrison
5-Aug-93 (MDM) - Modified to print the header information when running
interactively (not just when writing to a file)
21-nov-94 (SLF) - add start and stop time input parameters
buffer data
[Previous]
[Next]
NAME:
pr_qs_hxi
PURPOSE:
To print the contents of the quasi-static HXI file.
CALLING SEQUENCE:
pr_qs_hxi,qs_structure
pr_qs_hxi,qs_structure,/hc
pr_qs_hxi,qs_structure,/hc,outfile=outfile
INPUT:
qs_structure The quasi-static structure returned by the
routine rd_qs
CALLS: ***
DPRINT
OPTIONAL KEYWORD INPUT:
outfile If a filename is passed, the results will also
be written to that file.
hc If set, write the results to a temporary file,
print it and delete it.
OUTPUT:
Results are sent tty. The format is similar to that produced by
Facom mainframe calibration file.
MODIFICATION HISTORY:
written, 12-aug-92, J.R. Lemen, LPARL
[Previous]
[Next]
Name: pr_sfc
Purpose: list sfc file information
CALLS:
[Previous]
[Next]
NAME:
pr_sxt_term
PURPOSE:
To summarize times of data from the SXT terminator database (sft)
CALLING SEQUENCE:
pr_sxt_term ; List terminators for today
pr_sxt_term, '1-jan
pr_sxt_term, '30-oct-91', '2-nov-91', /hc, /by_point
pr_sxt_term, '30-oct-91', '2-nov-91', outfile=outfil, /cmd
INPUT:
st_date - The UT date to search for. If no date is entered,
the current date is used. The year is not
required, but if entered should be of the
form 1-Oct-91 style. The date should be
entered in string style with date first.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
OPTIONAL INPUT KEYWORDS:
all - Normally summarize images which have:
dpe_lo <= dpe <= dpe_hi and (=5.3s)
10s <= tim2night <= 25s (prior to D/N)
/all returns all images in the sft database.
dpe_hi - Default = 25 (5.3 sec): Sets max dpe value to include
dpe_ho - Default = 25 (5.3 sec): Sets min dpe value to include
tim2night_lo- Default = 10 sec. Sets min time prior to D/N to include.
tim2night_hi- Default = 25 sec. Sets max time prior to D/N to include.
by_point- If set, will summarize results by unique pointings and
by gt_filtb values (includes ND mesh and non-ND mesh cases).
cmd - If set, will call gt_center with /cmd to get command pointings
(less accurate, but considerably faster)
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
OPTIONAL OUTPUT KEYWORDS:
rmap - Roadmaps of the terminator images
index - Vector of indexes of the terminator images
dset - Vector of image numbers
infil - Vector of sft file names
coord - (2,N) array of results from gt_center(index,/angle)
CALLS: ***
DPRINT, RD_SDA [1], RD_SDA [2], RD_SDA [3], Rd_Roadmap [2], UNIQ [1], UNIQ [2]
UNIQ [3], fmt_tim [1], fmt_tim [2], gt_center, gt_dp_mode [1], gt_dp_mode [2]
gt_dp_rate [1], gt_dp_rate [2], gt_dpe, gt_filta, gt_filtb, gt_res, gt_shape_cmd
prstr [1], prstr [2], rd_roadmap [1], sft_files, tim2orbit [1], tim2orbit [2]
CALLED BY:
CHKTERM, GET_POINTS, term_score2
HISTORY:
14-Oct-94, J. R. Lemen, Written.
[Previous]
[Next]
NAME:
pr_sxtobs
PURPOSE:
Read the SXT observing log and list the images taken
CALLING SEQUENCE:
pr_sxtobs, '1-jan', /ffi
pr_sxtobs, '30-oct-91', '2-nov-91', /hc, /ffi
pr_sxtobs, '1-jan-92', '1-jan-92 14:00', outfil='jan92.sxtobs'
INPUT:
st_date - The UT date to search for. If no date is entered,
the current date is used. The year is not
required, but if entered should be of the
form 1-Oct-91 style. The date should be
entered in string style with date first.
OPTIONAL INPUT:
en_date - The ending date to list data for. If no date is passed,
it only searches for "st_date"
CALLS: ***
GET_INFO2, GET_NBYTES, anytim2ints [1], anytim2ints [2], fmt_tim [1], fmt_tim [2]
int2secarr [1], int2secarr [2], prstr [1], prstr [2], rd_obs
OPTIONAL KEYWORD INPUT:
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
ffi - If set, then only search FFI. Default is to list both PFI/FFI
pfi - If set, then only search PFI. Default is to list both PFI/FFI
notit - If set, then do not put the title on the output
OPTIONAL KEYWORD OUTPUT:
list - The "get_info" array information (string array)
obs - The SXT observing log data
HISTORY:
Written 28-Oct-93 by M.Morrison
[Previous]
[Next]
NAME:
PR_TABLES
PURPOSE:
prints table contents
CATEGORY:
CALLING SEQUENCE:
a = pr_tables(serno,/hard)
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
/hard prints to file *.tmp
OUTPUTS:
CALLS: ***
gtab_comm, gtab_entry, gtab_ffi, gtab_pfi, gtab_roi
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
HSH, written May 14, 1993
[Previous]
[Next]
NAME:
pr_teem
PURPOSE
To interactively select regions, and obtain temp and em for.
CALLING SEQUENCE
pr_teem,index,data,comp=comp
INPUT
Index : Index ¿
Data : Data Cube contining two images through separate SXT
filters.
OUTPUT
Interactive temperatures and emission measures.
OPTIONAL OUTPUT
none
OPTIONAL KEYWORD INPUT
comp : display images on a log scale.
RESTRICTION
some-- delta version 0.1415927
HISTORY
Written CRF 12 Nov 95 : caf@mssl.ucl.ac.uk
CALLS:
[Previous]
[Next]
NAME:
pr_trav_summ
PURPOSE:
To print out a summary of the YOHKOH personnel travel
SAMPLE CALLING SEQUENCE:
pr_trav_summ, '1-jan-92', '31-dec-92'
pr_trav_summ, '1-jan-92', '31-dec-92', outfil='trav.1992'
pr_trav_summ, '1-jan-92', '31-dec-92', /hc
INPUT:
st_tim - The starting time
en_tim - The ending time
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], DPRINT, Ex2Int [1], Ex2Int [2]
RD_TRAVEL, anytim2ex [1], anytim2ex [2], concat_dir [4], gt_day [1], gt_day [2]
int2secarr [1], int2secarr [2]
OPTIONAL KEYWORD INPUT:
outfil - The output file to put the results into
hc - If set, write results to a file and send to the printer
HISTORY:
Written 1992 by M.Morrison
23-Mar-93 (MDM) - Added /HC option
30-May-93 (MDM) - Added INFIL option
24-Aug-93 (MDM) - Modified to write to the home directory if
outfil is not specified
[Previous]
[Next]
Name: pr_visible
Purpose: print KSC orbit data, optionally make hardcopy
Calling Statement:
pr_visible ; Print the current week
pr_visible,startt,stopt ; Start and Stop times (JST)
pr_visible,week='96_04 ; Print week 96_04
pr_visible,file='outfile.dmp' ; Send output to a file
Input Parameters:
startt, stopt - start and stop times (JST!!)
week - week# (default is current week)
Keyword Parameters
hc - if set, send to printer
week - week number for full week listing other than current week
noprint - dont echo info to terminal (usually with outarry keyword)
file - Name of output file.
outarray - (output) string array of requested information
Calling Examples:
pr_visible [,/hc] ; print current week
pr_visible, day ; specified day (+24 hours)
pr_visible, t0,t1 ; specified time range
pr_visible, t0,t1,outarr=outarr,/noprint ; return info, no print
CALLED BY:
FIRST_LIGHT [1], FIRST_LIGHT [2]
History:
3-may-1995 (SLF)
10-jan-1996 (SLF) allow look at last years data...
22-Jan-1996 (JRL) Improved header and added file= for output file (via prstr)
11-dec-1996 (SLF) Added 1 hour to end of week determination (JST/GMT)
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], STR2ARR [1], STR2ARR [2]
anytim2weekid [1], anytim2weekid [2], concat_dir [4], cont2time [1]
cont2time [2], data_chk [1], data_chk [2], fmt_tim [1], fmt_tim [2], gt_day [1]
gt_day [2], gt_time [1], gt_time [2], jst, pr_week2tim [1], pr_week2tim [2]
prstr [1], prstr [2], rd_tfile [1], rd_tfile [2], sel_timrange [1]
sel_timrange [2], strsplit, tbeep [1], tbeep [2], tbeep [3], timegrid, where_arr [1]
where_arr [2]
Restrictions: does not span weekly files
[Previous]
[Next]
Name: prcols
Purpose: print up to 10 array parameters in aligned columns w/optional header
Input Parameters:
up to 10 arrays - number of elements must match
Keyword Paramters:
right/left/center - specify alignment/justification (see strjustify)
box - if set, draw box the whole thing (see strjustify)
header - optional header string to align with data - comma or blank delimited string will
(number of fields should match array sizes)
noprint - switch - if set, dont print to terminal (output via OUTARR keyword)
coldelim=coldelim - if string, delimiter to insert between columns
if number, number of blanks to insert between cols (def=2 blanks)
outarr - the aligned text array
Method: call strjustify, execute
CALLS: ***
MORE [1], MORE [2], PRCOL, STR2ARR [1], STR2ARR [2], data_chk [1], data_chk [2]
strjustify
History:
11-sep-1995 (SLF) - provide terminal front end to strjustify
[Previous]
[Next]
Name: pref_deldir
Purpose: manage directories for pref_super
Input Parameters:
prestage - stage number ABOUT to happen
Restrictions - hard coded for standard ref directories
CALLS:
[Previous]
[Next]
Name: pref_info
Purpose: return info to simplify pref_super management
CALLS:
[Previous]
[Next]
Name: pref_super
Purpose: build and submit a 'go_ref_super' job with parallel
processing for stages 1 and 2
History:
1-Oct-1994 (SLF)
CALLS: ***
ARR2STR [1], Arr2Str [2], CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3]
anytim2ex [1], anytim2ex [2], concat_dir [4], data_chk [1], data_chk [2]
ex2fid [1], ex2fid [2], file_append [1], file_append [2], get_logenv [1]
get_logenv [2], go_batch [1], go_batch [2], pr_status [1], pr_status [2]
set_logenv [1], set_logenv [2], str_replace [1], str_replace [2], tbeep [1]
tbeep [2], tbeep [3]
Restrictions:
prototype - not all ref_super keywords available
assume standard disks used for reformatting
[Previous]
[Next]
NAME:
PreficCk
Purpose:
return the a list of requested reformatted file prefixes.
CALLS: ***
PREFIXCK
HISTORY:
written 14-Nov-91, gal
updated, 17-Feb-92, added noobs
updated, 22-May-92, gal to sort returned prefix into
order on tape... if we add any files this will
simple sorting will need to be updated.
updated, 30-May-92, gal bug fix
[Previous]
[Next]
NAME:
PREFLARE_SEARCH
PURPOSE:
find displacements between flare location and preflare data
CATEGORY:
Yohkoh
CALLING SEQUENCE:
preflare_search, t1, t2, data_out
INPUTS:
start time, stop time
OPTIONAL (KEYWORD) INPUT PARAMETERS:
ROUTINES CALLED
GT_CENTER, RD_OBS
OUTPUTS:
data_out, structure with info
CALLS: ***
fmt_tim [1], fmt_tim [2], gt_center, gt_dp_mode [1], gt_dp_mode [2], gt_time [1]
gt_time [2], rd_obs
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
(HSH) written 26-Jun-96
[Previous]
[Next]
Name: prep_dat
Purpose: prep. weekly log file data (eg sfd) for CD rom writes.
Parameters:
dpath path to data
cdmax total space available on CD for data
archdir returned archive dir. for ydb files.
mvfiles list of files to move to /cdr_base
logdiscna String with disc name, etc. for logging
dest destination for CD
select not used yet.
CALLED BY:
mk_cd [1], mk_cd [2]
History: 26-May-95, GAL, written
26-Mar-96, GAL, mod for modvolume_inf.
12-Jan-97, PGS, changed value used for 'MB' to: 1048576.
[Previous]
[Next]
Name: prep_dat
Purpose: prep. weekly log file data (eg sfd) for CD rom writes.
Parameters:
dpath path to data
cdmax total space available on CD for data
archdir returned archive dir. for ydb files.
mvfiles list of files to move to /cdr_base
logdiscna String with disc name, etc. for logging
dest destination for CD
select not used yet.
CALLED BY:
mk_cd [1], mk_cd [2]
History: 26-May-95, GAL, written
26-Mar-96, GAL, mod for modvolume_inf.
[Previous]
[Next]
NAME: prep_gendat
PURPOSE: prepare general data for writing to CD.
Parameters:
dpath path(s) to data
flist list of files to write to CD
cdmax max. data size in Mb
archdir user defined archive directory
mvfiles fully qual. list of file names to mv /cdr_base
mvcmd move "mvfiles" command when dir-trees are used
logdiscna strings for logging and disc names
dest destination for CD
select switch to allow interactive selection of files
cdpath maps one to one to dpath list giving
user specified root directories.
CALLS: ***
CONCAT_DIR [1], CONCAT_DIR [2], CONCAT_DIR [3], concat_dir [4], file_size [1]
file_size [2], get_dirtree, input [1], input [2], modvolume_inf [1]
modvolume_inf [2], wmenu_sel [1], wmenu_sel [2]
CALLED BY:
mk_cd [1], mk_cd [2]
History: written by GAL, 26-May-95
19-Sep-95, GAL, bug fix for upcase in volume.inf file
13-Feb-96, GAL, corrected Mb conversion.
26-Mar-96, GAL, used modvolume_inf.
[Previous]
[Next]
NAME:
PREP_HXTDAT
PURPOSE:
Prepare a data structure from which images are synthesized by the
MOVIE program. The data structure is created according to a card file.
CALLING SEQUENCE:
hxt = prep_hxtdat(index,data,ssss,infil,bgd,'flare.card')
INPUT:
index - HXT index
data - HXT data
ssss - HXT ssss (automatically given by YODAT)
infil - HXT infil (same as above)
bgd - Background data (structure) prepared by MK_HXTDATA
program.
cardfile - A card file according to which output data is
prepared. Content of this file should be as
follows:
hda920811.2136
1090 0 1097 3
1098 0 1103 3
1104 0 1106 3
1107 0 1111 3
hda... - Reformat file name
1090, 1097, ... - MF number contained in ssss
0, 3, ... - Specifies portion of a MF
(0,1,2,3 available)
CALLS: ***
MK_HXT_SBT, MK_HXT_SUM, RD_CARD, WRT_HXTDAT, input [1], input [2], yesnox [1]
yesnox [2]
OPTIONAL KEYWORD INPUT:
none
OUTPUT:
hxt - An HXT data structure from which images are synthesized.
OPTIONAL KEYWORD OUTPUT:
none
RESTRICTIONS:
May be powerful (= detailed analyses available) than HXTPRO at this
moment, however need more preparations before running this program.
Contact HXT people for more details on the program.
HISTORY:
version 1.0 94.07.27 (Tue) T.Sakao written at MSSL.
[Previous]
[Next]
Name: prep_week
PURPOSE: prep Yohkoh archive tape data for a CD archive.
dpath if data is already on disk this is the path to it.
cdmax max. size of archive
archdir Returned archive dir. name for weekly files.
mvfiles Returned list of files to mv to /cdr_base and archive
when tape is used this list comes from xbd files.
logdiscna Ret: String with disc name, for logging.
dest Input: destination for CD, used to build logging string
rdtap_script String-array with commands for dumping tape
weekid Input: tape Archive week to dump.
stage_disk staging area, default = /sxt8data5/stage_b1
type_comp type of compression: none or 1 = Unix compress
2 = gzip
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], PREFIXCK, RdTap [1], RdTap [2]
RdTap [3], Rd_TapDir, break_file [4], ck, comp_sfr_arch, def_tapd [1], def_tapd [2]
def_tapd [3], get_afile_size [1], get_afile_size [2], input [1], input [2]
modvolume_inf [1], modvolume_inf [2], weekfiles, yo_arch_size
CALLED BY:
mk_cd [1], mk_cd [2]
History: written, Sep 95, GAL
26-Mar-96, GAL, mod for modvolume_inf.
09-Jan-98, PGS, mod modvolume_inf, type_comp,
changed flow of program, expanded call to
compression algorithm, added type_comp, stage_disk
changed defn. of mb, updated mvfiles, sfr_comp_arch call.
[Previous]
[Next]
Name: prep_week
PURPOSE: prep Yohkoh archive tape data for a CD archive.
dpath if data is already on disk this is the path to it.
cdmax max. size of archive
archdir Returned archive dir. name for weekly files.
mvfiles Returned list of files to mv to /cdr_base and archive
when tape is used this list comes from xbd files.
logdiscna Ret: String with disc name, for logging.
dest Input: destination for CD, used to build logging string
rdtap_script String-array with commands for dumping tape
weekid Input: tape Archive week to dump.
CALLS: ***
BREAK_FILE [1], BREAK_FILE [2], BREAK_FILE [3], PREFIXCK, RdTap [1], RdTap [2]
RdTap [3], Rd_TapDir, break_file [4], ck, comp_sfr_arch, def_tapd [1], def_tapd [2]
def_tapd [3], get_afile_size [1], get_afile_size [2], input [1], input [2]
modvolume_inf [1], modvolume_inf [2], weekfiles, yo_arch_size
CALLED BY:
mk_cd [1], mk_cd [2]
History: written, Sep 95, GAL
26-Mar-96, GAL, mod for modvolume_inf.
[Previous]
[Next]
NAME:
PRFITS
PURPOSE:
Prints a string array (e.g., a FITS header) on the screen, or
(optionally) on a printer
CALLING SEQUENCE:
prfits, stringarray
INPUTS:
stringarray = string array to be printed
OPTIONAL (KEYWORD) INPUT PARAMETERS:
printer = integer or string. If present and not of type string,
output is directed to the default printer. If string,
output is sent to named printer.
OUTPUTS:
No explicit outputs.
COMMON BLOCKS:
None.
SIDE EFFECTS:
Screen or printer output.
MODIFICATION HISTORY:
JPW, Nov. 1991
[Previous]
[Next]
NAME:
print_ogsel
PURPOSE:
FTP the output of OPOGEDT from FACOM and print it
CALLING SEQUENCE:
print_ogsel
print_ogsel,/browse,/print
CALLS: ***
anytim2ex [1], anytim2ex [2], ftp_from_facom, rd_tfile [1], rd_tfile [2]
HISTORY:
97/08/08 T. Sakurai
97/09/02 H. Hudson
[Previous]
[Next]
NAME:
PRINT_OP
PURPOSE:
FTP the OP file created by OPOGEDT from the FACOM, print it, and archive it
CALLING SEQUENCE:
print_op
print_op, /browse
OPTIONAL KEYWORD INPUT
/browse prints to screen and not to paper.
RESTRICTIONS
This program can only be run from flare20 because of ftp restrictions.
CALLS: ***
anytim2ex [1], anytim2ex [2], ftp_from_facom, rd_tfile [1], rd_tfile [2]
HISTORY:
97/08/05 T. Sakurai
97/08/06 H. Hudson, converted to procedure with keywords. Extracted
OP name for archiving.
[Previous]
[Next]
NAME:
PRINT_TEXT
PURPOSE:
Print Text Section from SXT Generic Files
Recursively Parses for line break characters
CALLING SEQUENCE:
PRINT_TEXT, text [,device]
data structure
INPUTS:
TEXT - string, generally extracted from SXT generic file
DEVICE - [Optional], Logical unit for output
COMMON BLOCKS;
NONE
CALLED BY:
INQ_GEN [1], INQ_GEN [2]
MODIFICATION HISTORY:
Version 0 - SLF, 3/15/91
[Previous]
[Next]
Name: printing
Purpose: toggle state of ys_noprint environmental/logical
CALLS: ***
set_logenv [1], set_logenv [2]
History:
7-Nov-1993 (SLF)
Calling Sequence
printing, /off
printing [,/on] (default)
[Previous]
[Next]
NAME:
pro chkop,file,ps=ps
KEYWORDS: Yohkoh Operation
PURPOSE:
show OGs in OP on AOSLOSGD.
CALLING SEQUENCE:
chkop,'op61021d.txt'
chkop,'op61021d.txt',ps='filename.ps'
INPUTS:
file: text filename of OP automatically generated by OP_FIRST_GUESS
CALLS: ***
CHKOP, INTERPRETFEM, LOADOPFEM, PLOTOGONALGD, PSSETUP, PSSHUT
HISTORY:
Ver 1.00 96-10-29 Prelim. Ver. Created by TY
[Previous]
[Next]
NAME:
pro disp2dimtable,pos,xr,yr,charsize=chars,ccol=ccol, $
xtitle=xtitle,ytitle=ytitle,ticklen=ticklen
KEYWORDS: Tool, Image Processing
PURPOSE:
display 2 dim. color table.
CALLING SEQUENCE:
ty2dimcolor
wdef,0
disp2dimtable,[30,30,100,100],[0,3],[2,5],ccol=14*14-1
INPUTS:
pos: position [left,bottom,right,top]
xr: x range
yr: y range
charsize: charactor size
ccol: color of frame and text (default 14*14-1)
ticklen : tick length (default -0.1)
CALLS: ***
DISP2DIMTABLE, TYTV
HISTORY:
rev. 1.0 19-Sep-96 T. Yoshida Created
[Previous]
[Next]
NAME:
pro dispcolortable,pos,vrg,crg,chars=chars,ccol=ccol,title=title, $
ticklen=ticklen
PURPOSE:
Display color table.
CALLING SEQUENCE:
dispcolortable,[10,10,20,100],[3.2,4.5],[0,!p.color],chars=0,5, $
ccol=!p.color*0.5,title='erg/s/cm!e-2!n'
INPUTS:
pos : position : 4-dimensional array
vrg : value range : 2 elements array
crg : color rabge : 2 elements array
KEYWORDS:
chars : size of characters
ccol : color of outline
title : 'erg/s/cm!e-2!n'
ticklen : length of tick (default -0.1)
CALLS: ***
DISPCOLORTABLE
HISTORY:
rev. 1.0 : Creation
rev. 2.0 30-Mar-95 : General procedure.
rev. 2.1 05-Apr-95 : Add the parameter 'ticklen'
[Previous]
[Next]
NAME:
pro interpretfem,fem,femtime,femflg,femstr=femstr,verbose=verbose
KEYWORDS: Yohkoh Operation
PURPOSE:
Interpret fem and return sorted times and event flags
CALLING SEQUENCE:
loadopfem,'op61021d.txt',time,og,fem,/verbose
interpretfem,fem,femtime,femflg,/verbose
INPUTS:
fem: fem data
OUTPUTS:
femtime: sorted time array of events
femflg: event flag (byte)
0:DAY 1:NIGHT 2:START SAA 3:END SAA
128:START KSC 129:START DSN-M 130:START DSN-C
131:START DSN-G 132:START DSN-W 133:START DSN-S
192:END KSC 193:END DSN-M 194:END DSN-C
195:END DSN-G 196:END DSN-W 197:END DSN-S
OPTIONAL OUTPUT:
femstr: COMMENT of the event flags (STRING ARRAY 256)
CALLS: ***
INTERPRETFEM, anytim2ints [1], anytim2ints [2], gt_day [1], gt_day [2]
gt_time [1], gt_time [2], sort_index [1], sort_index [2]
HISTORY:
Ver 1.0 96-10-23 Created by TY
[Previous]
[Next]
NAME:
pro loadopfem,file,time,OG,fem,verbose=verbose
KEYWORDS: Yohkoh Operation
PURPOSE:
load OP from the file generated by OP_FIRST_GUESS,
and fem by using rd_fem.
CALLING SEQUENCE:
loadopfem,'op61021d.txt',time,og,fem
print,gt_day(time,/str)
print,gt_time(time,/str)
print,og
INPUTS:
file: file name
OUTPUTS:
time: DAY and TIME (array of STRUCT = -> ANYTIM2INTS)
og: OG ID (byte array)
fem:
CALLS: ***
LOADOPFEM, LOADOPTEXT, anytim2ints [1], anytim2ints [2], gt_day [1], gt_day [2]
gt_time [1], gt_time [2], rd_fem
HISTORY:
Ver 1.0 96-10-23 Created by TY
[Previous]
[Next]
NAME:
pro loadoptext,file,time,OG,verbose=verbose
KEYWORDS: Yohkoh Operation
PURPOSE:
load OP from the file generated by OP_FIRST_GUESS.
CALLING SEQUENCE:
loadoptext,'op61021d.txt',time,og
print,gt_day(time,/str)
print,gt_time(time,/str)
print,og
INPUTS:
file: file name
OUTPUTS:
time: DAY and TIME (array of STRUCT = -> ANYTIM2INTS)
og: OG ID (byte array)
CALLS: ***
LOADOPTEXT, anytim2ints [1], anytim2ints [2], gt_day [1], gt_day [2], gt_time [1]
gt_time [2], rd_tfile [1], rd_tfile [2]
HISTORY:
Ver 1.0 96-10-22 Created by TY
Ver 1.01 96-10-25 Modify some bugs
Ver 1.1 97-02-21
[Previous]
[Next]
NAME:
pro optimizedbdr,file,showparam=showparam
KEYWORDS: Yohkoh Operation
PURPOSE:
load OP from the file generated by OP_FIRST_GUESS,
and show optimized bitrate sequences.
CALLING SEQUENCE:
optimizedbdr,'op61021d.txt',/showparam
INPUTS:
file: op text file name
OPTIONAL KEYWORDS:
showparam: show parameters
CALLS: ***
CNT_DSNAOS, CNT_DSNLOS, CNT_HIGHANDMED, CNT_KSCAOSD, CNT_KSCAOSN, CNT_KSCLOS
CNT_NIGHT, CNT_TOTBDR, INTERPRETFEM, LOADOPFEM, OPTBDRGTEVT, OPTBDRGTOGFEM
OPTBDRMAIN, OPTBDRSHOWSEQ, OPTIMIZEDBDR, OROR, TIMEINTERVAL, gt_time [1]
gt_time [2], sort_index [1], sort_index [2]
HISTORY:
Ver 1.0 96-10-25 Created by TY
[Previous]
[Next]
NAME:
pro plotogonalgd,ogtime,og,femtime,femflg,starttime=starttime,endtime=endtime,ps=ps
KEYWORDS: Yohkoh Operation
PURPOSE:
show og on X window
CALLING SEQUENCE:
loadopfem,'op61021d.txt',time,og,fem,/verbose
interpretfem,fem,femtime,femflg,/verbose
plotogonalgd,time,og,femtime,femflg
INPUTS:
ogtime: time array
og: og ID
femtime: time array
femflg: event flag
OPTIONAL INPUTS:
starttime: start time of the plot
endtime: end time of the plot
CALLS: ***
PLOTOGONALGD, SHOWALGDPOS, SHOWAOSLOSGD, TIMEINTERVAL, anytim2ints [1]
anytim2ints [2], rd_tfile [1], rd_tfile [2]
HISTORY:
Ver 0.50 96-10-28 Prelim. Ver. Created by TY
[Previous]
[Next]
NAME:
pro pssetup,land=land,range=range,offset=offset,color=color,nofont=nofont $
,xoffset=xoffset,yoffset=yoffset,xsize=xsize,ysize=ysize $
,filename=filename,epsi=epsi,quiet=quiet
KEYWORDS: Tool
PURPOSE:
Open PostScript Device
(See also 'psshut.pro')
CALLING SEQUENCE:
pssetup
pssetup,offset=[0.8,0.8],range=[19,27]
pssetup,/color,/land
INPUTS KEYWORDS:
xoffset : X Offset
yoffset : Y Offset
offset : [xoffset,yoffset] (cm)
(default [1.905,12.7](port) [1.905,26.035](land))
xsize : X canvas size
ysize : Y canvas size
range : [xsize,ysize] (cm)
(default [17.78,12.7](port) [24.13,17.78](land))
filename : default 'idl.ps' or 'idl.eps'
OPTIONAL KEYWORDS:
land : set landscape (default portrait)
color : set color postscript. (default monochro)
nofont : Use idl original fonts. (default postscript fonts)
epsi : set epsi mode.
(Call ps2epsi (one of the 'ghostscript' commands),
in the program 'psshut.pro' when shut the output device.)
quiet : quiet mode
COMMON
common psutil,fname,eps,epsfname
HISTORY
95-09-20 ver 2.0 T.Yoshida
CALLS:
[Previous]
[Next]
NAME:
pro psshut,quiet=quiet
KEYWORD: Tool
PURPOSE:
Close Post Script Device
(See also 'pssetup.pro')
CALLING SEQUENCE:
psshut
OPTIONAL KEYWORDS:
quiet : quiet mode
COMMON
common psutil,fname,eps,epsfname
HISTORY
95-09-20 ver 2.0 T.Yoshida
CALLS:
[Previous]
[Next]
NAME:
pro showalgdpos,sec,psec,pclm,timestr=timestr
KEYWORDS: SUBROUTINE for showaoslosgd.pro
CALLS:
[Previous]
[Next]
NAME:
pro showaoslosgd,femtime,femflg, $
starttime=starttime,endtime=endtime,no10minreset=no10minreset
KEYWORDS: Yohkoh Operation
PURPOSE:
show aoslosgd on X window
CALLING SEQUENCE:
loadopfem,'op61021d.txt',time,og,fem,/verbose
interpretfem,fem,femtime,femflg,/verbose
showaoslosgd,femtime,femflg
INPUTS:
femtime: time array
femflg: event flag
OPTIONAL INPUTS:
starttime: START TIME of the plot
CALLS: ***
SHOWALGDBOX, SHOWALGDEVENT, SHOWALGDFRAME, SHOWALGDINIT, SHOWALGDPOS
SHOWALGDTBLINIT, SHOWAOSLOSGD, TIMEINTERVAL, anytim2ints [1], anytim2ints [2]
gt_day [1], gt_day [2], gt_time [1], gt_time [2]
HISTORY:
Ver 1.0 96-10-28 Created by TY
[Previous]
[Next]
NAME:
pro ty2dimcolor,g=gau
KEYWORDS: Tool, Image Processing
PURPOSE:
Set 2 dimensional color table.
CALLING SEQUENCE:
ty2dimcolor
ty2dimcolor,g=2.0
ty2dimcolor,g=0.5
INPUTKEYWORDS:
g: default 1.0
CALLS: ***
RETRUN, TY2DIMCOLOR
HISTORY:
rev. 1.0 19-Sep-96 T. Yoshida Created
[Previous]
[Next]
NAME:
pro ty2dimtv,x,y,xrange=xr,yrange=yr,range=rr, $
xpos,ypos,xsize=xsize,ysize=ysize
KEYWORDS: Tools, Image Processing
PURPOSE:
display overlap image
CALLING SEQUENCE:
ty2dimcolor
ty2dimtv,img1,img2,xr=[0,1000],yr=[100,2000],100,100,xsize=500,ysize=300
INPUTS:
x,y: images to overlap
xrange: Range of the first image
yrange: Range of the second image
range: Range of the two images
xpos: The first parameter of the IDL procedure 'TV'
ypos: The second parameter of the IDL procedure 'TV'
INPUT KEYWORDS:
xsize,ysize: size of the image to be displayed
CALLS: ***
NORMIMGC, TY2DIMTV, TYTV
HISTORY:
rev. 1.0 19-Sep-96 T. Yoshida Created
[Previous]
[Next]
NAME:
pro tyoshidacolor,ro=ro,go=go,bo=bo, $
tgamma=tgamma,toffset=toffset,trange=trange, $
rgamma=rgamma,rrange=rrange
KEYWORDS: Tool, Image Processing
PURPOSE:
This is a subroutine for the program 'xtyoshidact.pro'
Red Temperature$B$H!"(BRainbow$B$r$^$<$?(Bcolor table$B$N:n@.!#(B
CALLING SEQUENCE:
tyoshidacolor
tyoshidacolor,tgamma=0.7,toffset=35,trange=[10,80]
INPUTS KEYWORDS:
tgamma :$BFz%F!<%V%k$N(BGamma correction factor
toffset :$BFz%F!<%V%k$NCf?4(Bcolor$BHV9f(B
[trange(0)<toffset<trange(1)]
trange :$BFz%F!<%V%k$NHO0O(B [bottom,top]
color$B?t$O(B100$B$J$N$G!"(B[0:99]$B$NHO0O$,(Breal
rgamma :$B@V%F!<%V%k$N(BGamma correction factor
rrange :$B@V%F!<%V%k$NHO0O(B(%) [bottom,top]
OUTPUTS KEYWORDS:
ro,go,bo : $B%+%i!<%F!<%V%k(B
CAUTIONS:
0$BHV$KGr!"(B[1:100]$B$KFz%F!<%V%k(B
[101:!d.table_size-2]$B$K@V%F!<%V%k!"(B!d.table_size-1$B$K9u(B
$BFz%F!<%V%k$K(B100$B?'!"(B0$B$KGr!":GBg$K9u$N(B
$B9g7W(B102$B?'0J>e$G$kC<Kv$G$J$$$H%(%i!<$r5/$3$9!#(B
CALLS: ***
REVERSE, TYOSHIDACOLOR
HISTORY:
95-11-22 $B5HED(B
95-11-28 $B5HED(B $B:G=i$N?'<h$j%P%0$N=$@5(B
[Previous]
[Next]
NAME:
pro xtyoshidact,param,nox=nox
PURPOSE:
Set the color table in which Red-Temperature table and
Rainbow table are included.
$BFz%F!<%V%k$H(BRed Temperature$B%F!<%V%k$N:.:_$7$?%F!<%V%k$N?'9g$$(B
$BJQ99!#(B
CALLING SEQUENCE:
xtyoshidact
xtyoshidact,[ 0, 100, 0, 0, 99, 50, 0]
xtyoshidact,[ 0, 100, 0, 0, 99, 50, 0],/nox
INPUTS:
param : Parameter array
which will be showed (printed out)
at the end of this program.
$B%Q%i%a!<%?(Barray
$B=*N;$7$?$H$-$K=*N;;~$N%Q%i%a!<%?$,I=<($5$l$k!#(B
$B8e$K:F8=$7$?$$>l9g$KMQ$$$k(B
OPTIONAL KEYWORDS:
nox : Run this program without X-widget.
This parameter should be set under the device 'PS'.
$B;XDj$5$l$?(Bparam$B$b$7$/$OD>A0$N(Bparam$B$G%+%i!<%F!<%V%k$r(B
$B@_Dj$7(BX-Widget$B$r3+$+$:$K=*N;!#(B
$B0lEY?'9g$$$r@_Dj$7$?8e!"(BPS$B$J$I$K=PNO$9$k$H$-$K;H$&!#(B
CALLS: ***
NORMIMGC, TYOSHIDACOLOR, XMANAGER, XTYCTCHGPARAM, XTYCTCMTOUT, XTYCTGETPARAM
XTYCTRESETTBL, XTYCTSETCOLOR, XTYCTSETPARAM, XTYCT_EVENT, xtyoshidact
HISTORY:
95-11-21 $B5HED(B
95-11-28 $B5HED(B $B:G=i$N?'<h$j%P%0$N=$@5(B
[Previous]
[Next]
NAME:
pro_list
PURPOSE:
To return a list of the IDL procedures (.PRO files)
that are in the IDL path (!path)
SAMPLE CALLING SEQUENCE:
pros = pro_list()
pros = pro_list(/hc, /land)
CALLS: ***
DPRINT, REVERSE, SPLIT_FILES, STR2ARR [1], STR2ARR [2], STR_LASTPOS [1], UNIQ [1]
UNIQ [2], UNIQ [3], file_list [1], file_list [2], str_lastpos [2]
OPTIONAL KEYWORD INPUT:
outfile - If a filename is passed, the results will also
be written to that file
hc - If set, write the results to a temporary file, print
it and delete it.
OUTPUT:
Return a list of procedure/function/program names
OPTIONAL KEYWORD OUTPUT:
dirs - the directory associated with the IDL procedure
HISTORY:
Written 23-Apr-92 by M.Morrison
29-Apr-93 (MDM) - Added sample calling sequence to the header
[Previous]
[Next]
NAME:
PROCESS_BCS_CALS
PURPOSE:
Uses source log to process all BCS cals in log
CATEGORY:
CALLING SEQUENCE:
.run process_bcs_cals
INPUTS:
Source of input file (local or on $DIR_BCS_CALS) requested.
User must specify if stims are to be processed (def=source cal)
OPTIONAL INPUT PARAMETERS:
Can create log and plot files
set PLOT = 1 for plot (set_plot to PS for PostScript files)
set LOG = 1 for logging of the results to default files
Debug mode assumed by default (just checks existance of files!)
set DEBUG = 0 for normal running
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
Creates log and plot files if requested
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB 19-May-94 Written from seperate routines
RDB 18-Aug-94 Some improvment to documentation
RDB 09-Feb-95 Allowed use of general input file
RDB 29-Nov-97 Clarified some questions; specific dirs at ISAS.
[Previous]
[Next]
NAME:
PROCESS_BCS_CALS
PURPOSE:
Uses source log to process all BCS cals in log
CATEGORY:
CALLING SEQUENCE:
.run process_bcs_cals
INPUTS:
Source of input file (local or on $DIR_BCS_CALS) requested.
User must specify if stims are to be processed (def=source cal)
OPTIONAL INPUT PARAMETERS:
Can create log and plot files
set PLOT = 1 for plot (set_plot to PS for PostScript files)
set LOG = 1 for logging of the results to default files
Debug mode assumed by default (just checks existance of files!)
set DEBUG = 0 for normal running
KEYWORD PARAMETERS:
OUTPUTS:
OPTIONAL OUTPUT PARAMETERS:
Creates log and plot files if requested
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
PROCEDURE:
MODIFICATION HISTORY:
RDB 19-May-94 Written from seperate routines
RDB 18-Aug-94 Some improvment to documentation
RDB 09-Feb-95 Allowed use of general input file
[Previous]
[Next]
NAME:
PROM_ON_DISK
PURPOSE:
overlay prominence positions onto a disk image
CATEGORY:
Yohkoh
CALLING SEQUENCE:
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
CALLS: ***
LIMB2DISK, fmt_tim [1], fmt_tim [2], get_rb0p [1], get_rb0p [2], gt_pix_size
int2secarr [1], int2secarr [2], sxt_cen [1], sxt_cen [2]
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
24-nov-97, ad hoc version (HSH)
[Previous]
[Next]
NAME:
PROTRACTOR
PURPOSE:
CATEGORY:
CALLING SEQUENCE:
protractor, index, data
INPUTS:
OPTIONAL (KEYWORD) INPUT PARAMETERS:
OUTPUTS:
CALLS: ***
HXA_SUNCENTER, LIMB_TRACE, fmt_tim [1], fmt_tim [2], gt_pix_size
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
MODIFICATION HISTORY:
HSH written 12-Sep-94
[Previous]
[Next]
NAME:
prt_afile_size
Purpose:
List files and sizes in Kbytes by file-types.
CALLING SEQUENCE:
prt_afile_size, infn, fn_sz
infn is a list of file names
fn_sz is a list of file sizes in bytes
CALLS: ***
PREFIXCK, UNIQ [1], UNIQ [2], UNIQ [3]
CALLED BY:
go_rdtap [1], go_rdtap [2]
History:
written, 17-Sept-92, gal
updated, 1-Dec-92, by gal, to output 0 byte requests
and success count.
[Previous]
[Next]
Name: prt_flst
Purpose: Print a list of files on /cdr_base to be pre-mastered
or re-mastered.
Parameters:
dpath path to data
cdmax total space available on CD for data
logdiscna String with disc name, etc. for logging
dest destination for CD
select not used yet.
CALLED BY:
mk_cd [1], mk_cd [2]
History: 03-Aug-95, GAL, written
12-Jan-98, PGS, changed value used for MB to: 1048576.
----------------------------------------------------------------
[Previous]
[Next]
Name: prt_flst
Purpose: Print a list of files on /cdr_base to be pre-mastered
or re-mastered.
Parameters:
dpath path to data
cdmax total space available on CD for data
logdiscna String with disc name, etc. for logging
dest destination for CD
select not used yet.
CALLED BY:
mk_cd [1], mk_cd [2]
History: 03-Aug-95, GAL, written
----------------------------------------------------------------
[Previous]
[Next]
NAME:
PS_LONG
PURPOSE:
Setup to fill portait page on PostScript printer
CALLING SEQUENCE:
PS_LONG
INPUTS:
none
OUTPUTS:
none
CALLED BY:
ANAL_STIMS, BCS_24HR_PLOT [1], BCS_24HR_PLOT [3], BCS_MULTI [1], BCS_MULTI [2]
LCBDA, PLOT_WINDOWS, lcbsd, lcobs, plot_bcs_free, plotbft [1]
RESTRICTIONS:
PROCEDURE:
Issues command: DEVICE,YSIZ=25.,YOFFSET=3
If Environment variable US_FUNNY_PAPER exists then the
command is: DEVICE,YSIZ=23.,YOFFSET=3
MODIFICATION HISTORY:
RDB Oct-92
RDB 22-Jun-92 Different size for US_FUNNY_PAPER
RDB 5-Jan-94 Return without action if YOMAN exists
[Previous]
[Next]
NAME:
PS_RESET
PURPOSE:
Reset device plot size to PostScript normal
This is the antidote to PS_LONG; also clears /landscape
CALLING SEQUENCE:
PS_RESET
INPUTS:
none
OUTPUTS:
none
CALLED BY:
ANAL_STIMS, BCS_24HR_PLOT [1], BCS_24HR_PLOT [3], BCS_MULTI [1], BCS_MULTI [2]
LCBDA, lcbsd, plot_bcs_free, plotbft [1]
RESTRICTIONS:
PROCEDURE:
Issues the command:
DEVICE,/INCHES,YSIZE=5.,YOFFSET=5.,/portrait
MODIFICATION HISTORY:
RDB Oct 91
RDB 2-Mar-93 Added /protrait switch
RDB 5-Jan-94 Return without action if YOMAN exists
[Previous]
[Next]
function psfpixon,datain,psfin,partan=partan,continued=continued,nosnr=nosnr, $
nochi2=nochi2,snr=snr,pseudoscale=pseudoscale
NAME:
PSFPIXON
PURPOSE:
Deconvolves a point spread function from an image using the fractal
pixon basis of Pina and Puetter.
CATEGORY:
CALLING SEQUENCE:
new_image = psfpixon(image,psf)
INPUTS:
image = image to be cleaned
psf = point spread function (centered in array of arbitrary size).
PSF must have the spatial scale of the image *after* it is rebinned
by pseudoscale, i.e. if pseudoscale is not 1, image and
psf must have different spatial scales.
OPTIONAL INPUT PARAMETERS:
KEYWORD PARAMETERS
/snr = Use SNR weighting in the iteration
/nochi2 = Use the GOF from Pina and Puetter rather than chi^2
/partan = Use parallel tangent iteration rather than conjugate gradient
/continued = Use continued parallel tangent iteration rather than
conjugate gradient
pseudoscale = Magnification of pseudo image (required). Integer GE 1.
OUTPUTS:
new_image = cleaned image
CALLS: ***
CHI_SQR, CONGRID [1], CONGRID [2], CONGRID [3], EXPERIMENTAL, LOCAL_SMOOTH
PIXON_MAP, PSF_GOF_FUNC, conj_grad2, minF_partan
CALLED BY:
SXTPIXON
COMMON BLOCKS:
SIDE EFFECTS:
RESTRICTIONS:
The edges of the image are not deconvolved properly so the images
should be centered in the FOV.
PROCEDURE:
Uses FPB Pixon deconvolution of Pina and Puetter.
MODIFICATION HISTORY:
T. Metcalf June 23, 1994
[Previous]
[Next]
NAME:
PTV
PURPOSE:
Load an array into the Peritek VCT graphics interface for
display as a pseudocolor or true color image.
CATEGORY:
Image display.
CALLING SEQUENCE:
PTV,A,FN=FN,RESULT=RESULT
INPUTS:
A = A vector, two or three dimensional array to be loaded as an image.
KEYWORD PARAMETERS:
FN = Name of file to save image in C-readable format
PLANE = VCT plane ('red','green,','blue', or one other)
RESULT = Error message
OUTPUTS: Result (optional)
CALLED BY:
ptv_sxt, wr_peritek
PROCEDURE:
Spawns call to C routine VCTQPSEUDO or VCT*TRUE which displays
the image
MODIFICATION HISTORY:
10-Jul-1991 - Written by GLS
15-nov-1991 - Updated, JRL -- Added DIR_PERITEK environment
variable to reference C routines.
29-apr-1991 - Updated, GLS -- Use !version to determine which
C routine (V bus version or Q bus
version) to call.
17-aug-1992 - Updated, gal -- uses dimensions on input image to
determine which C routine to use.
if img is 2-d pseudo color is used
and 3-d uses true color 24 bit color.
[Previous]
[Next]
Name:
ptv_sxt
Purpose:
Construct and load two image arrays and an annotation array into the
Peritek graphics display interface.
Calling Sequence:
ptv_sxt,index0,data0,index1,data1[,nimg=nimg, $
comment0=comment0,comment1=comment1,color_tab0=color_tab0, $
color_tab1=color_tab1,num=num
Input Parameters:
index0 = SDA file dataset index for large image.
data0 = Array containing large image.
index1 = SDA file dataset index for small image.
data1 = Array containing small image.
Optional Input Parameters:
None.
Optional Input Keywords:
nimg = Number of images to display (default = 1).
comment0 = First string which is put into the text box of the display.
comment1 = Second string which is put into the text box of the display.
color_tab0 = three dimension byte array (3,256) of r,g,b colors for
data0.
color_tab1 = three dimension byte array (3,256) of r,g,b colors for
data1.
charsize = Character size of displayed text.
num = Number displayed in text box.
delta_sec = Number seconds the 2nd image is w.r.t. to the primary image
If present, this will be display as hours on the image.
CALLS: ***
CONGRID [1], CONGRID [2], CONGRID [3], PLCT, PTV, gt_filta, mk_title [1], mk_title [2]
mk_title [3], sxt_decomp [1], sxt_decomp [2], wr_peritek
Side effects:
mk_title will create window,1 for displaying the titles.
These are read back using tvrd.
Modification History:
Written, 21-may-92, G. L. Slater
Modified, 26-may-92, JRL: Check the result keyword of PTV for errors
Modified, 13-Aug-92, GAL: added logic for using 24 bit color.
Modified, 17-Aug-92, GAL: added new parameters for titles
Modified, 9-Jan-92, GAL, mk_title change
[Previous]
[Next]
NAME:
PTVSCL
PURPOSE:
Byte scale an array and load it into the Peritek VCT graphics
interface for display as a pseudocolor image.
CATEGORY:
Image display.
CALLING SEQUENCE:
PTVSCL,A,FN=FN,RESULT=RESULT
INPUTS:
A: A vector or two dimensional array to be loaded as an image.
KEYWORD PARAMETERS:
FN: Name of file to save image in C-readable format
PLANE: VCT plane ('red','green,','blue', or one other)
RESULT: Error message
OUTPUTS:
RESULT: Error message (optional)
CALLED BY:
video_test
PROCEDURE:
Spawns call to C routine VCTQPSEUDO which displays the image
MODIFICATION HISTORY:
10-Jul-1991 - Written by GLS
15-nov-1991 - Updated, JRL -- Added DIR_PERITEK environment
variable to reference C routines.
29-apr-1991 - Updated, GLS -- Use !version to determine which
C routine (V bus version or Q bus
version) to call.