[oorb] r282 committed - Gnuplot script for plotting topocentric ranges

[oo...@googlecode.com]

Revision: 282
Author: jxvir...@gmail.com
Date: Fri May 8 14:32:56 2015 UTC
Log: Gnuplot script for plotting topocentric ranges
https://code.google.com/p/oorb/source/detail?r=282

Added:
/trunk/gnuplot/sor_plot_range.gp

=======================================
--- /dev/null
+++ /trunk/gnuplot/sor_plot_range.gp Fri May 8 14:32:56 2015 UTC
@@ -0,0 +1,57 @@
+#====================================================================#
+# #
+# Copyright 2002,2003,2004,2005,2006,2007,2008,2009,2010,2011 #
+# Mikael Granvik, Jenni Virtanen, Karri Muinonen, Teemu Laakso, #
+# Dagmara Oszkiewicz #
+# #
+# This file is part of OpenOrb. #
+# #
+# OpenOrb is free software: you can redistribute it and/or modify it #
+# under the terms of the GNU General Public License as published by #
+# the Free Software Foundation, either version 3 of the License, or #
+# (at your option) any later version. #
+# #
+# OpenOrb is distributed in the hope that it will be useful, but #
+# WITHOUT ANY WARRANTY; without even the implied warranty of #
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU #
+# General Public License for more details. #
+# #
+# You should have received a copy of the GNU General Public License #
+# along with OpenOrb. If not, see <http://www.gnu.org/licenses/>. #
+# #
+#====================================================================#
+#
+# Gnuplot script for plotting the orbital-element probability-density
+# function in Keplerian phase space.
+#
+# Author: JV
+# Version: 2015-04-15
+#
+reset
+set terminal postscript eps enhanced linewidth 1.0 12.0
+set out 'sor_ranges.eps'
+#set logscale x
+unset key
+set pointsize 0.3
+set size 1.0,1.0
+set origin 0.0,0.0
+set multiplot
+set lmargin 11
+set rmargin 2
+set tmargin 1
+set bmargin 3
+set size 0.5,0.5
+set origin 0.0,0.5
+set xlabel 'rho1 [au]'
+set ylabel 'number density'
+set style histogram rowstacked gap 0
+set style fill solid 0.5 border lt -1
+plot "sor_histo.out" using 1:2 smooth freq with boxes
+set size 0.5,0.5
+set origin 0.0,0.0
+set ylabel 'pdf'
+plot "sor_histo.out" using 1:3 smooth freq with boxes
+#plot 'sor_ranges.out' using 1:2 pt 7
+#set size 0.5,0.33
+#set origin 0.5,0.66
+reset