Manuel_LMDZE.texfol/Graphiques/disvert.gp

# This is a script for Gnuplot. This script plots the coefficients
# "ap" and "bp" used for the vertical grid in LMDZ.

reset
##set terminal x11 reset
set terminal epslatex color
set output "disvert.eps"

pa = 5e4
ps = 101325.
bp(s) = exp(1. - 1. / s**2)
ap(s) = pa * (s - bp(s))
p(s) = ap(s) + bp(s) * ps

##set terminal x11 0 title "bp"
set parametric
##set xlabel "bp"
set ylabel "$s$"
set yrange [] reverse
set dummy s
set trange [0.01:1]
##unset key
##plot bp(s), s

# # set terminal x11 1 title "ap"
# # set xlabel "ap"
# # plot ap(s), s

##set terminal x11 2 title "p"
##unset xlabel
##set key default
plot p(s), s title "$p$", ap(s), s title "$a$", \
     bp(s) * ps, s title "$b \\times p_s$"

set output
1	# This is a script for Gnuplot. This script plots the coefficients
2	# "ap" and "bp" used for the vertical grid in LMDZ.
3
4	reset
5	##set terminal x11 reset
6	set terminal epslatex color
7	set output "disvert.eps"
8
9	pa = 5e4
10	ps = 101325.
11	bp(s) = exp(1. - 1. / s**2)
12	ap(s) = pa * (s - bp(s))
13	p(s) = ap(s) + bp(s) * ps
14
15	##set terminal x11 0 title "bp"
16	set parametric
17	##set xlabel "bp"
18	set ylabel "$s$"
19	set yrange [] reverse
20	set dummy s
21	set trange [0.01:1]
22	##unset key
23	##plot bp(s), s
24
25	# # set terminal x11 1 title "ap"
26	# # set xlabel "ap"
27	# # plot ap(s), s
28
29	##set terminal x11 2 title "p"
30	##unset xlabel
31	##set key default
32	plot p(s), s title "$p$", ap(s), s title "$a$", \
33	bp(s) * ps, s title "$b \\times p_s$"
34
35	set output