supermongo script.
http://www.dur.ac.uk/physics.astrolab/sm.html
-----------------------------------------------------------------------------------------
crontab example:
33 16 9 3 * /data/lgou/LMCX-1/rxte/DataAnalysis/exampl_sm/script 0 gou.ps examp /data/lgou/LMCX-1/rxte/DataAnalysis/exampl_sm >& /dev/null
-------------------------------------------------------------------------------------------
csh script file which call for supermongo:
#! /bin/csh
#
# specify time on screen, select models, the script filename, and the working directory.
set wait = $1
set psfile = $2
set script_name=$3
set work_dir=$4
########################################################################
# Start SM
#
cd $work_dir
echo $work_dir
sm -S << FIN
define wait $wait
# do not alter the code above this line
# The following line outputs checks if the user has enetered a 'wait' time >1
# if the 'wait' time from the command line is 0, then the output is sent to a file
#
# to output to gif rather than postscript, uncomment the second line and
# comment out the first line
if (\$wait) {device x11} else {device postencap $psfile} # output to .ps
# if (\$wait) {device x11} else {device gif $psfile} # output to .gif
echo $psfile
# Type your sm commands here. An example is provided
inp $script_name
# do not alter below this line
!sleep \$wait
FIN
exit
---------------------------------------------------------------------------------------------------
sm script file:
#device postencap myfile.ps
set x = 1, 1000, 1
set y = (x**2.5 * lg(x)) + 1.
window 2 1 1 1 # define a window which is 2 in x, 1 in y
# and plot current data into bottom left window
ticksize 0 0 0 0
ctype black
expand 1.0
ptype 10 3 # use an open polygon with 10 sides
lweight 1.
limits x y
box 1 1 1 0 # label all of the axis except the RHS y-axis
expand 0.5 # make the points half the size
points x y # draw the points in x and y
expand 1.0 # set the rest of the plot back to normal size
relocate 50 8.5e7
label This is a more complex plot
ltype 4
relocate 50 8.2e7
draw 800 8.2e7
expand 1.3 # label the x and y axis
xlabel x-axis
ylabel y-axis
expand 1.0
# now make same plot, but in log scale in y
window 2 1 2 1 # now change windows and make a different plot
ticksize 0 0 -1 0 # plot the y-axis in log scale
set z = lg(y) # take the log of y
limits x z # reset the limits for the new data
box 1 0 1 1
ptype 4 3
ctype blue
points x z
ctype black
expand 0.8
xlabel x-axis
