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