1	#! /bin/csh -f
    2	#
    3	# study a selected orbit from demo1, this script is 'tkrun' enabled
    4	#
    5	
    6	#> ENTRY run=run1 	
    7	#> SCALE orbit=0 	0:1000:1
    8	#> SCALE omega=-1.3	-2:2:0.01
    9	#> SCALE range=2        0:10:0.1
   10	#> SCALE tstop=10       0:100:0.1
   11	#> SCALE dt=0.01        0.0001:0.1:0.0001
   12	
   13	set run=run1
   14	set orbit=0
   15	set omega=1.3
   16	set range=2
   17	set tstop=10
   18	set dt=0.01
   19	
   20	foreach a ($*)
   21	  set $a
   22	end
   23	
   24	#	set some useful shorthands
   25	set xyrange=-${range}:${range}
   26	set nsteps=`nemoinp $tstop/$dt format=%d`
   27	
   28	#       no printouts to for tkrun  to run quiet
   29	setenv DEBUG -1
   30	
   31	snapmask $run.13 - $orbit |\
   32	 snapscale - - rscale=1,1,0 vscale=1,1,0 |\
   33	 stoo - - |\
   34	 orbint - $run.tmp nsteps=$nsteps dt=$dt potname=ccd potfile=$run.112 potpars=$omega
   35	orbplot $run.tmp xrange=$xyrange yrange=$xyrange "xlabel=x orbit $orbit" yapp=1/xs
   36	orbsos $run.tmp y | tabplot - 3 4 point=2,0.1 color=2 yapp=2/xs
   37	
   38	# 	cleanup temporary file(s)
   39	rm -f $run.tmp