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

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

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

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

set terminal x11 2 title "p"
unset xlabel
set key default
plot p(sigma), sigma title "p", ap(sigma), sigma title "ap", bp(sigma) * ps, \
     sigma title "bp * ps"
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
7	pa = 5e4
8	ps = 101325.
9	bp(sigma) = exp(1. - 1. / sigma**2)
10	ap(sigma) = pa * (sigma - bp(sigma))
11	p(sigma) = ap(sigma) + bp(sigma) * ps
12
13	set terminal x11 0 title "bp"
14	set parametric
15	set xlabel "bp"
16	set ylabel "sigma"
17	set yrange [] reverse
18	set dummy sigma
19	set trange [0.01:1]
20	unset key
21	plot bp(sigma), sigma
22
23	set terminal x11 1 title "ap"
24	set xlabel "ap"
25	plot ap(sigma), sigma
26
27	set terminal x11 2 title "p"
28	unset xlabel
29	set key default
30	plot p(sigma), sigma title "p", ap(sigma), sigma title "ap", bp(sigma) * ps, \
31	sigma title "bp * ps"