 |
Plotting H on PlotXY for paper AA518L8Fig3
Source code name: "real_life_AA518L8Fig3.py"
Programming language: Python
Topic: Physics/Cosmology
DMelt Version 1. Last modified: 05/09/2015. License: Pro
https://datamelt.org/code/cache/real_life_AA518L8Fig3_1694.py
To run this script using the DMelt IDE,
copy the above URL link to the menu [File]→[Read script from URL] of the DMelt IDE.
from java.awt import Color
from org.jplot2d.data import *
from org.jplot2d.element import *
from org.jplot2d.transform import *
from org.jplot2d.util import *
from jhplot import *
from array import *
c1=HPlotXY(False)
c1.visible()
ef=c1.getFactory()
p=c1.getPlot()
title=ef.createTitle("Title")
title.setFontScale(2)
p.addTitle(title)
p.getLegend().setVisible(False)
p.setFontSize(16)
xaxis = ef.createAxis();
xaxis.getTitle().setText("r (cm)");
xaxis.getTitle().setFontScale(1);
xaxis.getTickManager().getAxisTransform().setTransform(TransformType.LOGARITHMIC);
xaxis.getTickManager().getAxisTransform().setRange( Range.Double(5e13, 2.5e17));
xaxisTop = ef.createAxis();
xaxisTop.setPosition(AxisPosition.POSITIVE_SIDE);
xaxisTop.getTitle().setText('r(\u2033)'.decode('unicode-escape'))
xaxisTop.getTitle().setFontScale(1);
xaxisTop.getTickManager().getAxisTransform().setTransform(TransformType.LOGARITHMIC);
xaxisTop.getTickManager().getAxisTransform().setRange(Range.Double(2.8e-2, 2e2));
xaxisTop.getTickManager().setLabelFormat("%.0m");
yaxisLeft = ef.createAxis();
yaxisLeft.setColor(Color.RED);
yaxisLeft.getTitle().setText("X");
yaxisLeft.getTitle().setFontScale(1);
yaxisLeft.getTickManager().getAxisTransform().setTransform(TransformType.LOGARITHMIC);
yaxisLeft.getTickManager().getAxisTransform().setRange(Range.Double(3e-10, 7e-7));
yaxisRight = ef.createAxis();
yaxisRight.setColor(Color.BLUE);
yaxisRight.setPosition(AxisPosition.POSITIVE_SIDE);
yaxisRight.getTitle().setText("$\\mathrm{T_K}$");
yaxisRight.getTitle().setFontScale(1);
yaxisRight.getTickManager().getAxisTransform().setTransform(TransformType.LINEAR);
yaxisRight.getTickManager().getAxisTransform().setRange(Range.Double(0, 1200));
yaxisRight.getTickManager().setTickInterval(500);
p.addXAxis(xaxis);
p.addXAxis(xaxisTop);
p.addYAxis(yaxisLeft);
p.addYAxis(yaxisRight);
tkx = array('d',[ 1.1e14, 1e15, 1e16, 1e17 ])
tky = array('d',[1200, 360, 80, 0 ])
tkl = ef.createLayer(ef.createXYGraph(tkx, tky));
tkl.setColor(Color.BLUE);
tka = ef.createSymbolAnnotation(5e14, 600, "T$\\mathrm{_K}$");
tka.setColor(Color.BLUE);
tka.setFontScale(1.2);
tkl.addAnnotation(tka);
p.addLayer(tkl, xaxis, yaxisRight);
# SiC2
scx = array('d',[ 5e13, 1e16, 2e16, 7e16, 1e17, 2.1e17])
scy = array('d',[2e-7, 2e-7, 5e-7, 2e-7, 5e-8, 3e-10])
scl = ef.createLayer(ef.createXYGraph(scx, scy));
scl.setColor(Color.RED);
sca = ef.createSymbolAnnotation(2e16, 1e-7, "SiC$\\mathrm{_2}$");
sca.setColor(Color.RED);
sca.setFontScale(1.2);
scl.addAnnotation(sca);
p.addLayer(scl, xaxis, yaxisLeft);
#SiC2 LTE
ltex = array('d',[5e13, 6e13, 9e13, 1.3e14, 2e14, 2.1e14, 2.7e14])
ltey = array('d',[5e-8, 3e-7, 1.8e-7, 3e-7, 4e-8, 6e-8, 3e-10 ])
ltel = ef.createLayer(ef.createXYGraph(ltex, ltey));
ltel.setColor(Color.GREEN);
ltea = ef.createSymbolAnnotation(7e13, 3e-9, "SiC$\\mathrm{_2}$\nLTE");
ltea.setColor(Color.GREEN);
ltea.setFontScale(1.2);
ltel.addAnnotation(ltea);
p.addLayer(ltel, xaxis, yaxisLeft);
You see the box below because you did not login.