Contour Plot (xfig example)
Revision as of 11:05, 17 April 2018 by Obst (talk | contribs) (Obst moved page Example: Contour Plot to Contour Plot (xfig example) without leaving a redirect)
Contour Plot of XTE J1946+274
require("isisscripts"); variable qu1 = 6; variable qu2 = 3; variable vars = 14; variable colmap=png_get_colormap_names[7]; variable raw = String_Type[0]; variable i,j,k; k=0; for (i=1 ; i<=1; i=i+1) { for (j=i+1; j<=vars; j=j+1) { raw = [raw,sprintf("/userdata/data/mueller/1946/epochs/epochI/contoursCRSFIII/V%dV%d.fits",i,j)]; }} variable map2 = load_conf(raw[0]);map2.chisqr = log(map2.chisqr); variable p2 = xfig_plot_new (qu1,qu2); p2.world(map2.px.min,map2.px.max,map2.py.min,map2.py.max); p2.ylabel("$\\sigma_{\\text{G}}$ [keV]"); p2.plot_png(map2.chisqr; cmap=colmap); p2.axis(;color="white"); p2.y2axis(;ticlabels=0); p2=xfig_plot_contours(raw[0],p2;width=2,line=[0,1,2],col=["white","white","white"]); variable map3 = load_conf(raw[1]);map3.chisqr = log(map3.chisqr); variable p3 = xfig_plot_new (qu1,qu2); p3.world(map3.px.min,map3.px.max,map3.py.min,map3.py.max); p3.plot_png(map3.chisqr; cmap=colmap); p3.axis(;color="white"); p3.y2axis(;ticlabels=0); p3.ylabel("$\\Gamma$"); p3=xfig_plot_contours(raw[1],p3;width=2,line=[0,1,2],col=["white","white","white"]); variable map4 = load_conf(raw[2]);map4.chisqr = log(map4.chisqr); variable p4 = xfig_plot_new (qu1,qu2); p4.world(map4.px.min,map4.px.max,map4.py.min,map4.py.max); p4.plot_png(map4.chisqr; cmap=colmap); p4.axis(;color="white"); p4.y2axis(;ticlabels=0); p4.x2axis(;ticlabels=0); p4.ylabel("$N_{\\text H}\\ [10^{22}\\,\\text{cm}^{-2}]$"); p4=xfig_plot_contours(raw[2],p4;width=2,line=[0,1,2],col=["white","white","white"]); variable map5 = load_conf(raw[3]);map5.chisqr = log(map5.chisqr); variable p5 = xfig_plot_new (qu1,qu2); p5.world(map5.px.min,map5.px.max,map5.py.min,map5.py.max); p5.ylabel("$E_{\\text{cut}}$ [keV]"); p5.plot_png(map5.chisqr; cmap=colmap); p5.axis(;color="white"); p5.y2axis(;ticlabels=0); p5=xfig_plot_contours(raw[3],p5;width=2,line=[0,1,2],col=["white","white","white"]); variable map6 = load_conf(raw[4]);map6.chisqr = log(map6.chisqr); variable p6 = xfig_plot_new (qu1,qu2); p6.world(map6.px.min,map6.px.max,map6.py.min,map6.py.max); p6.ylabel("$E_{\\text{fold}}$ [keV]"); p6.xlabel("$E_{\\text{CRSF}}$ [keV]"); p6.plot_png(map6.chisqr; cmap=colmap); p6.axis(;color="white"); p6.y2axis(;ticlabels=0); p6=xfig_plot_contours(raw[4],p6;width=2,line=[0,1,2],col=["white","white","white"]); variable map7 = load_conf(raw[9]);map7.chisqr = log(map7.chisqr); variable p7 = xfig_plot_new (qu1,qu2); p7.world(map7.px.min,map7.px.max,map7.py.min,map7.py.max); p7.ylabel("$\\tau_{\\text{G}}$"); p7.plot_png(map7.chisqr; cmap=colmap); p7.axis(;color="white"); p7.y2axis(;ticlabels=0); p7=xfig_plot_contours(raw[9],p7;width=2,line=[0,1,2],col=["white","white","white"]); variable map8 = load_conf(raw[11]);map8.chisqr = log(map8.chisqr); variable p8 = xfig_plot_new (qu1,qu2); p8.world(map8.px.min,map8.px.max,map8.py.min,map8.py.max); p8.ylabel("$\\xi_{\\text{b}}$"); p8.xlabel("$E_{\\text{CRSF}}$ [keV]"); p8.axis(;color="white"); p8.y2axis(;ticlabels=0); p8.plot_png(map8.chisqr; cmap=colmap); p8=xfig_plot_contours(raw[11],p8;width=2,line=[0,1,2],col=["white","white","white"]); variable map9 = load_conf(raw[10]);map9.chisqr = log(map9.chisqr); variable p9 = xfig_plot_new (qu1,qu2); p9.world(map9.px.min,map9.px.max,map9.py.min,map9.py.max); p9.ylabel("$E_{\\text{G}}$ [keV]"); p9.axis(;color="white"); p9.x2axis(;ticlabels=0); p9.y2axis(;ticlabels=0); p9.plot_png(map9.chisqr; cmap=colmap); p9=xfig_plot_contours(raw[10],p9;width=2,line=[0,1,2],col=["white","white","white"]); variable finpl1 = xfig_multiplot(p9,p2,p7,p8); variable finpl2 = xfig_multiplot(p4,p3,p5,p6); variable finplot = xfig_new_hbox_compound(finpl1,finpl2,0.); finplot.render("CRSFcontours.eps"); ()=system("epstopdf CRSFcontours.eps"); ()=system("pdftops -eps -level3 CRSFcontours.pdf");