[cig-commits] [commit] master: Fixed example_geotherms (36053b4)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri Dec 12 17:39:39 PST 2014
Repository : https://github.com/geodynamics/burnman
On branch : master
Link : https://github.com/geodynamics/burnman/compare/1d2396e3c9d632ea6c4c7f8cbc61788c9e87ff80...a4ab63ec34d5a2c5c5f7373d058a34bea198b797
>---------------------------------------------------------------
commit 36053b4e24b7092a490dbc49d791e67be0bfce48
Author: I <kmanunterborn at gmail.com>
Date: Fri Dec 12 17:15:50 2014 -0800
Fixed example_geotherms
Fixed example_geotherms so as to not call B&S at too low of depth.
changed depth cutoff
Changed depth cutoff
>---------------------------------------------------------------
36053b4e24b7092a490dbc49d791e67be0bfce48
examples/example_geotherms.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/examples/example_geotherms.py b/examples/example_geotherms.py
index 5b028ac..11aa190 100644
--- a/examples/example_geotherms.py
+++ b/examples/example_geotherms.py
@@ -40,7 +40,7 @@ from burnman import minerals
if __name__ == "__main__":
# we want to evaluate several geotherms at these values
- pressures = np.arange(1e9,128e9,3e9)
+ pressures = np.arange(8.5e9,128e9,3e9)
#load two builtin geotherms and evaluate the temperatures at all pressures
temperature1 = burnman.geotherm.brown_shankland(pressures)
@@ -89,7 +89,7 @@ if __name__ == "__main__":
plt.plot(pressures/1e9,temperature5,'-m',label="Adiabat with pv (70%) and fp(30%)")
plt.legend(loc='lower right')
- plt.xlim([0, 130])
+ plt.xlim([8.5, 130])
plt.xlabel('Pressure/GPa')
plt.ylabel('Temperature')
plt.savefig("output_figures/example_geotherm.png")
More information about the CIG-COMMITS
mailing list