Contour Plot (xfig example)

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");