keiko-6: ssh -X saku@tomo
saku@tomo's password:
Last login: Thu May 22 14:03:58 2008 from keiko
Warning - Computer Use Policy in Effect
Unauthorized access and use of this computer is a criminal violation
of the Computer Fraud and Abuse Act (18 U.S.C. Sec. 1030). The
responsibilities of all users of this system are described in
Smithsonian Directive 931 (http://www.cfa.harvard.edu/do/sd931.pdf).
Use of this computer and any other CfA IT assets implies an acceptance
of the terms and conditions stated in SD 931. Any violation of these
policies should be reported to syshelp@head.cfa.harvard.edu
Badly placed ()'s.
tomo-1: cd /data/saku/Keck.dir/xte1118.dir/dan
tomo-2: molly
molly -- spectrum analysis. If you are new to molly type 'help molly' for
information about the program. Type 'help' or '?' for a command list.
This version was linked at 14:03:11, 14/08/2006
Startup file found.
Opened1: /data/tomo/Software/TRMSOFT/F77/molly/startup
Maximum spectrum length set = 4201
Maximum number of spectra = 2856
Stored user command number 1
Command name: splot
Command def: plot $01 $01 \
Stored user command number 2
Command name: gops
Command def: device /PS
Stored user command number 3
Command name: goxw
Command def: device /XS
Files will be written even if a file of the same name exists.
Closed1: /data/tomo/Software/TRMSOFT/F77/molly/startup
No local startup file found
molly> mxpix 26836
All old spectra will be lost.
Do you want to proceed ? [n]: y
Maximum spectrum length set = 26836
Maximum number of spectra = 447
molly> @molly.in
molly> edit 1 74
Commands:
clear -- Clear headers
file -- Read in headers from special file
any -- Read in a column from a file
term -- Entry from terminal.
del -- Delete header item.
ren -- Rename header item.
quit -- Quit
edit> file
Enter filename [header]: molly.info
Variable 1 name: Object , type: C
Variable 2 name: Record , type: I
Variable 3 name: Day , type: I
Variable 4 name: Month , type: I
Variable 5 name: Year , type: I
Variable 6 name: UTC , type: D
Variable 7 name: Dwell , type: R
Variable 8 name: RA , type: D
Variable 9 name: DEC , type: D
Variable 10 name: Equinox , type: D
edit> Quit
(to compute HJDs, heliocentric velocities, etc:)
molly> hfix
First slot to fix [ 1]:
Last slot to fix [ 1]: 74
Telescope ('list' for list, 'skip' to ignore) [INT]:
molly> hfix
First slot to fix [ 1]:
Last slot to fix [ 1]: 74
Telescope ('list' for list, 'skip' to ignore) [INT]: list
Please identify the telescope
Palomar 200in
Wilson
Campanas
Lemmon
WHT
INT
JKT
UKIRT
Kitt Peak
AAT
CTIO
McDonald
MMT
VLT
SAAO 1.9m
or 'elsewhere'
telescope> elsewhere
What is the name of the telescope ?
molly> ax a m
molly> trail
First slot for trail [ 1]:
Last slot for trail [ 1]: 74
Level for lowest colour index [ 0.000000 ]:
Level for highest colour index [0.2000000 ]:
Enter X axis type [W]:
Start wavelength [ 3884.960 ]:
End wavelength [ 10399.55 ]:
Number of pixels [ 26836]:
molly> trail
First slot for trail [ 1]:
Last slot for trail [ 74]:
Level for lowest colour index [ 0.000000 ]:
Level for highest colour index [0.2000000 ]: 0.08
Enter X axis type [W]:
Start wavelength [ 3884.960 ]: 5800
End wavelength [ 10399.55 ]: 6300
Number of pixels [ 26836]:
molly> splot
First slot to plot [ 1]:
Last slot to plot [ 1]:
1, obj: XTEJ1118 , run: 1, exp: 300.0 , time: 14/02/04, 07:52:04
mask> plot
Plot which slot? [ 1]:
1, obj: XTEJ1118 , run: 1, exp: 300.0 , time: 14/02/04, 07:52:04
Plot in P(ixels), W(avelength) or V(elocity)? [w]:
Enter plot limits (Left,Right, Bottom,Top) :
3559. 1.0741E+04 -6.668 10.11
0, 0 for old pair, A(bort), to proceed
5800 6300 0 0.1
Enter plot limits (Left,Right, Bottom,Top) :
5800. 6300. 0.000 0.1000
0, 0 for old pair, A(bort), to proceed
mask> mask
Set mask via C(ursor) or T(erminal)? [C]: t
Mask in P(ixels), W(avelength) or V(elocity)? [w]:
Enter angstrom range (0 0 to quit): 5850 5950
Enter angstrom range (0 0 to quit): 6100 6200
Enter angstrom range (0 0 to quit): 0 0
mask> show
Currently masked regions
Masked Angstroms: 5850. to 5950.
Masked Angstroms: 6100. to 6200.
mask> quit
Cycle 0 RMS = 1.35939705 (sigma) Rejects = 1139
Cycle 1 RMS = 1.33863544 (sigma) Rejects = 1222
Cycle 2 RMS = 1.33455098 (sigma) Rejects = 1241
Cycle 3 RMS = 1.33350205 (sigma) Rejects = 1246
Cycle 4 RMS = 1.3335017 (sigma) Rejects = 1246
101, obj: XTEJ1118 , run: 1, exp: 300.0 , time: 14/02/04, 07:52:04
molly>
load flux.mol (has list of all files)
need orbital phases
For Bowen blend 4686:
splot
sfit (to get rid of continuum)
select region around line
put in velocity space
write out for all files
put in an ephemeris to place date
modmap> V
rest wavelength: "4685.7"
gamma velocity:
number of poixels in velocity grid: 10
Maximum velocity in tomogram: 1000
FWHM of local line profile: "3.8"
modmap> d
sfit: spline fit
pfit: polynomial fit
Gaussfit: to fit lines
st spectrum to fit [ 74]: 1
First slot for fits [ 75]: 2
Number of splines [ 3]:
Lower reject threshold [-3.000000 ]:
Upper reject threshold [ 3.000000 ]:
Number of reject cycles [ 10]:
Weight factor (0=no change, 1=unit) [ 0.000000 ]:
Set mask? [YES]: no
No regions masked
Cycle 0 RMS = 1.36208713 (sigma) Rejects = 234
Cycle 1 RMS = 1.3415755 (sigma) Rejects = 319
Cycle 2 RMS = 1.33680844 (sigma) Rejects = 342
Cycle 3 RMS = 1.33680773 (sigma) Rejects = 342
2, obj: XTEJ1118 , run: 1, exp: 300.0 , time: 14/02/04, 07:52:04
molly> splot
First slot to plot [ 1]:
Last slot to plot [ 74]: 1
1, obj: XTEJ1118 , run: 1, exp: 300.0 , time: 14/02/04, 07:52:04
molly> exit
Are you sure ? [n]
molly> trail
First slot for trail [ 1]:
Last slot for trail [ 74]:
Level for lowest colour index [ 0.000000 ]:
Level for highest colour index [0.8000000E-01]:
Enter X axis type [W]:
Start wavelength [ 5800.000 ]:
End wavelength [ 6300.000 ]:
Number of pixels [ 26836]:
molly> mgfit
First slot to fit [ 1]:
Last slot to fit [ 74]:
First slot to store fits (0 to ignore) [ 701]:
Input file defining fit [gaussfit.seven]:
No polynomial will be fitted.
Read 2 gaussian definitions.
FWHM resolution has not been specified
Output file for results? to ignore [ ]:
Lower sigma reject threshold [-4.000000 ]:
Upper sigma reject threshold [ 4.000000 ]:
Number of sub-intervals/exposure [ 1]:
Fit A(ll) or I(ndividually)? [I]:
Reset to initial values for each spectrum? [yes]:
For XTEJ1118+480 spectra with larger errors(?) are 29,31,34,39,47,50,52,68
molly> plot 0 0
Enter list of spectra to plot
Slots selected currently are:
25
pick> p
Enter series of slot ranges
N,-M for M copies of spectrum N; 0,0 to stop
25,-1
N,-M for M copies of spectrum N; 0,0 to stop
725,-1
N,-M for M copies of spectrum N; 0,0 to stop
0,0
Entry finished
pick> q
2 slots chosen.
25, obj: XTEJ1118 , run: 25, exp: 300.0 , time: 14/02/04, 10:27:44
725, obj: XTEJ1118 , run: 25, exp: 300.0 , time: 14/02/04, 10:27:44
molly>
molly> wasc
Output file name [fit40.dat]: spectrum40.dat
Slot to dump [ 140]: 40
X units [ANGSTROMS]:
Y units [MJY]:
Dumped XTEJ1118 to disc. 932 points.
molly>