result_o6_z = coldens.main(coords, hsmooth, mass*hydrogen*o6, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.o6.z.l%3.1f.png'%(snapname,lgrid),Vmin=12.7, Vmax=16.7,ion='O VI',mhalo=mhalo,extralabel='O VI',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_o7_z = coldens.main(coords, hsmooth, mass*hydrogen*o7, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.o7.z.l%3.1f.png'%(snapname,lgrid),Vmin=13.0, Vmax=17.0,ion='O VII',mhalo=mhalo,extralabel='O VII',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_o8_z = coldens.main(coords, hsmooth, mass*hydrogen*o8, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.o8.z.l%3.1f.png'%(snapname,lgrid),Vmin=13.0, Vmax=17.0,ion='O VIII',mhalo=mhalo,extralabel='O VIII',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_mg2_z = coldens.main(coords, hsmooth, mass*hydrogen*mg2, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.mg2.z.l%3.1f.png'%(snapname,lgrid),Vmin=11, Vmax=15,ion='Mg II',mhalo=mhalo,extralabel='Mg II',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_si2_z = coldens.main(coords, hsmooth, mass*hydrogen*si2, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.si2.z.l%3.1f.png'%(snapname,lgrid),Vmin=11.5, Vmax=15.5,ion='Si II',mhalo=mhalo,extralabel='Si II',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_si3_z = coldens.main(coords, hsmooth, mass*hydrogen*si3, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.si3.z.l%3.1f.png'%(snapname,lgrid),Vmin=11.5, Vmax=15.5,ion='Si III',mhalo=mhalo,extralabel='Si III',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_si4_z = coldens.main(coords, hsmooth, mass*hydrogen*si4, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.si4.z.l%3.1f.png'%(snapname,lgrid),Vmin=11.5, Vmax=15.5,ion='Si IV',mhalo=mhalo,extralabel='Si IV',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
if(snip=="0"):
result_ne8_z = coldens.main(coords, hsmooth, mass*hydrogen*ne8, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.ne8.z.l%3.1f.png'%(snapname,lgrid),Vmin=11, Vmax=15,ion='Ne VIII',mhalo=mhalo,extralabel='Ne VIII',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_mg10_z = coldens.main(coords, hsmooth, mass*hydrogen*mg10, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.mg10.z.l%3.1f.png'%(snapname,lgrid),Vmin=11, Vmax=15,ion='Mg X',mhalo=mhalo,extralabel='Mg X',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_si12_z = coldens.main(coords, hsmooth, mass*hydrogen*si12, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.si12.z.l%3.1f.png'%(snapname,lgrid),Vmin=11, Vmax=15,ion='Si XII',mhalo=mhalo,extralabel='Si XII',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_n5_z = coldens.main(coords, hsmooth, mass*hydrogen*n5, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.n5.z.l%3.1f.png'%(snapname,lgrid),Vmin=11, Vmax=15,ion='N V',mhalo=mhalo,extralabel='N V',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
else:
result_ne8_z = result_mg10_z = result_si12_z = result_n5_z = 0*result_o6_z
coldens.make_colourmap('coldens.%s.fo6.z.l%3.1f.png'%(snapname,lgrid), result_o6_z-result_o_z, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -4.0, -1.0, 'xOVI', ngrid, mhalo=mhalo, extralabel='OVI/O',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
coldens.make_colourmap('coldens.%s.fo7.z.l%3.1f.png'%(snapname,lgrid), result_o7_z-result_o_z, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -3.0, 0.0, 'xOVII', ngrid, mhalo=mhalo, extralabel='OVII/O',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
coldens.make_colourmap('coldens.%s.fo8.z.l%3.1f.png'%(snapname,lgrid), result_o8_z-result_o_z, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -3.0, 0.0, 'xOVIII', ngrid, mhalo=mhalo, extralabel='OVIII/O',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
f = file('colion_map.%s.z.l%3.1f.dat'%(snapname,lgrid), 'w')
for i in xrange(ngrid):
for j in xrange(ngrid):
f.write('%3d %3d % 5.2f % 5.2f % 5.2f % 5.2f % 5.2f % 5.2f % 5.2f % 5.2f %5.1f % 5.2f % 5.2f % 5.2f % 5.2f\n'%(i,j,result_h_z[i,j],result_o_z[i,j],result_h1_z[i,j],result_c4_z[i,j],result_o6_z[i,j],result_o7_z[i,j],result_o8_z[i,j],result_mg2_z[i,j],np.sqrt((ngrid/2.-0.5-j)**2+(ngrid/2.-0.5-i)**2)/ngrid*lgrid*1e+03,result_ne8_z[i,j],result_mg10_z[i,j],result_si12_z[i,j],result_n5_z[i,j]))
f.close()
write_radial_column('columnave.%s.z.%3.1f.dat'%(snapname,lgrid),result_h_z,result_o_z,result_h1_z,result_c4_z,result_o6_z,result_o7_z,result_o8_z,result_mg2_z,result_ne8_z,result_mg10_z,result_si12_z,result_n5_z,ngrid,lgrid,0)
write_radial_column('columnave_linsum.%s.z.%3.1f.dat'%(snapname,lgrid),10**result_h_z,10**result_o_z,10**result_h1_z,10**result_c4_z,10**result_o6_z,10**result_o7_z,10**result_o8_z,10**result_mg2_z,10**result_ne8_z,10**result_mg10_z,10**result_si12_z,10**result_n5_z,ngrid,lgrid,1)
result_h_y = coldens.main(coords, hsmooth, mass*hydrogen, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.hydrogen.y.l%3.1f.png'%(snapname,lgrid),Vmin=17.5, Vmax=21.5,ion='H',theta=90,mhalo=mhalo,extralabel='H',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
result_o_y = coldens.main(coords, hsmooth, mass*oxygen/16.0, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.oxygen.y.l%3.1f.png'%(snapname,lgrid),Vmin=14, Vmax=18,ion='O',theta=90,mhalo=mhalo,extralabel='O',haloinfostr=haloinfostr, docbar=docbar, R200=R200)
result_h1_y = coldens.main(coords, hsmooth, mass*hydrogen*h1, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.h1.y.l%3.1f.png'%(snapname,lgrid),Vmin=13, Vmax=21,ion='H I',theta=90,mhalo=mhalo,extralabel='H I',haloinfostr=haloinfostr,docbar=docbar, R200=R200)
(snapname, direction, lgrid),
Vmin=22.0,
Vmax=29.0,
ion='Tmass',
npix=ngrid,
theta=thetaangle,
phi=phiangle,
redshift=redshift,
extralabel='T')
coldens.make_colourmap('coldens_num.%s.T.%s.l%3.1f.png' %
(snapname, direction, lgrid),
result_tempmass - result_mass,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
4.0,
7.0,
"T",
ngrid,
redshift=redshift,
extralabel='T',
haloinfostr=haloinfostr)
print "result_tempmass= ", result_tempmass
result_nHmass = coldens.main(
coords_pos,
hsmooth_pos,
nH_mass,
center,
lgrid,
lgrid,
lgridz,
if(direction=="z"):
thetaangle = 0
phiangle = 0
if(direction=="y"):
thetaangle = 90
phiangle = 0
if(direction=="x"):
thetaangle = 0
phiangle = 90
result_den = coldens.main(coords[indexes], hsmooth[indexes], denominator, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.%s.%s.l%3.1f.png'%(snapname,ion,direction,lgrid),Vmin=colmin, Vmax=colmax,ion=ion,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
result_num_T = coldens.main(coords[indexes], hsmooth[indexes], numerator_T, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='num.png',Vmin=colmin, Vmax=colmax,ion=ion,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
coldens.make_colourmap('coldens.%s.%s_T.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), result_num_T-result_den, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, 3.0, 7.0, ion, ngrid, redshift=redshift, extralabel=ionname)
result_num_nH = coldens.main(coords[indexes], hsmooth[indexes], numerator_nH, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='num.png',Vmin=colmin, Vmax=colmax,ion=ion,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
coldens.make_colourmap('coldens.%s.%s_nH.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), result_num_nH-result_den, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -6.0, -1.0, ion, ngrid, redshift=redshift, extralabel=ionname)
result_num_nHWD = coldens.main(coords[indexes], hsmooth[indexes], numerator_nHWD, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='num.png',Vmin=colmin, Vmax=colmax,ion=ion,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
coldens.make_colourmap('coldens.%s.%s_nHWD.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), result_num_nHWD-result_den, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -6.0, -1.0, ion, ngrid, redshift=redshift, extralabel=ionname)
result_num_P = coldens.main(coords[indexes], hsmooth[indexes], numerator_P, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='num.png',Vmin=colmin, Vmax=colmax,ion=ion,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
coldens.make_colourmap('coldens.%s.%s_P.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), result_num_P-result_den, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, 0, 5, ion, ngrid, redshift=redshift, extralabel=ionname)
result_num_PWD = coldens.main(coords[indexes], hsmooth[indexes], numerator_PWD, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='num.png',Vmin=colmin, Vmax=colmax,ion=ion,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
coldens.make_colourmap('coldens.%s.%s_PWD.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), result_num_PWD-result_den, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, 0, 5, ion, ngrid, redshift=redshift, extralabel=ionname)
result_num_rad = coldens.main(coords[indexes], hsmooth[indexes], numerator_rad, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='num.png',Vmin=colmin, Vmax=colmax,ion=ion,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
coldens.make_colourmap('coldens.%s.%s_rad.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), 10**(result_num_rad-result_den), -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, 0, 1000.0, ion, ngrid, redshift=redshift, extralabel=ionname)
fig_name='num.png',
Vmin=colmin,
Vmax=colmax,
ion=ion,
theta=thetaangle,
phi=phiangle,
redshift=redshift,
extralabel=ionname)
#make_colourmap(fig_name, ResultW, Xmin, Xmax, Ymin, Ymax, Vmin, Vmax, ion, npix, redshift)
coldens.make_colourmap('coldens.%s.%s.%s.l%3.1f.png' %
(snapname, ion, direction, lgrid),
result_num - result_den,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
colmin,
colmax,
ion,
ngrid,
redshift=redshift,
extralabel=ionname,
kpc='True')
print "result_num= ", result_num
print "result_den= ", result_den
N10lgrid = 0.0
nN10lgrid = 0
N05lgrid = 0.0
nN05lgrid = 0
N025lgrid = 0.0
nN025lgrid = 0
if (i==1):
direction="y"
thetaangle = 90
phiangle = 0
if (i==2):
direction="z"
thetaangle = 0
phiangle = 0
lgridz = lgrid*2 # *2 #Now 2 times as of 12/30/14 # Four times longer than Lgrid in zdirection.
result_Nh = coldens.main(coords, hsmooth, mass*hydrogen, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens_num.%s.hydrogen.%s.l%3.1f.png'%(snapname,direction,lgrid),Vmin=18.0+np.log10(lgridz/2./lgrid), Vmax=22.0+np.log10(lgridz/2./lgrid),ion='H',npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel='H')
print "result_Nh= ", result_Nh
result_mass = coldens.main(coords, hsmooth, mass, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens_num.%s.mass.%s.l%3.1f.png'%(snapname,direction,lgrid),Vmin=18.0+np.log10(lgridz/2./lgrid), Vmax=22.0+np.log10(lgridz/2./lgrid),ion='mass',npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel='mass')
result_tempmass = coldens.main(coords, hsmooth, temp_mass, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens_num.%s.Tmass.%s.l%3.1f.png'%(snapname,direction,lgrid),Vmin=22.0, Vmax=29.0,ion='Tmass',npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel='T')
coldens.make_colourmap('coldens_num.%s.T.%s.l%3.1f.png'%(snapname,direction,lgrid), result_tempmass-result_mass, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, 4.0,7.0, "T", ngrid, redshift=redshift, extralabel='T',haloinfostr=haloinfostr)
print "result_tempmass= ", result_tempmass
result_nhmass = coldens.main(coords, hsmooth, nh_mass, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens_num.%s.nHmass.%s.l%3.1f.png'%(snapname,direction,lgrid),Vmin=14.0, Vmax=22.0,ion='nhmass',npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel='')
coldens.make_colourmap('coldens_num.%s.nh.%s.l%3.1f.png'%(snapname,direction,lgrid), result_nhmass-result_mass, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -7.0,-1.0, "nH", ngrid, redshift=redshift, extralabel='n$_{\mathrm{H}}$',haloinfostr=haloinfostr)
print "result_nhmass= ", result_nhmass
result_vdotroverrmass = coldens.main(coords, hsmooth, vdotroverr_mass, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens_num.%s.vdotroverrmass.%s.l%3.1f.png'%(snapname,direction,lgrid),Vmin=-vlimit, Vmax=vlimit,ion='vdotroverrmass',npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel='')
print "result_vdotroverrmass= ", result_vdotroverrmass
print "10**result_mass= ", 10**result_mass
coldens.make_colourmap('coldens_num.%s.vdotroverr.%s.l%3.1f.png'%(snapname,direction,lgrid), result_vdotroverrmass/10**result_mass, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -vlimit,vlimit, "vdotroverr", ngrid, redshift=redshift, extralabel='v',haloinfostr=haloinfostr)
print "result_vdotroverrmass/10**result_mass= ", result_vdotroverrmass/10**result_mass
result_vxmass = coldens.main(coords, hsmooth, vx_mass, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens_num.%s.vxmass.%s.l%3.1f.png'%(snapname,direction,lgrid),Vmin=-vlimit, Vmax=vlimit,ion='vxmass',npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel='')
print "result_vxmass= ", result_vxmass
print "10**result_mass= ", 10**result_mass
denominator = mass*hydrogen*o7_neg
ionname = 'OVII'
colmin = 13.5
colmax = 16.5
if(ion=="o8"):
numerator = mass*hydrogen*o8_pos
denominator = mass*hydrogen*o8_neg
ionname = 'OVIII'
colmin = 13.5
colmax = 16.5
if(direction=="z"):
thetaangle = 0
phiangle = 0
if(direction=="y"):
thetaangle = 90
phiangle = 0
if(direction=="x"):
thetaangle = 0
phiangle = 90
result_den = coldens.main(coords, hsmooth, denominator, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens_diff.%s.%s.%s.l%3.1f.png'%(snapname,ion,direction,lgrid),Vmin=colmin, Vmax=colmax,ion=ionname,npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
result_num = coldens.main(coords, hsmooth, numerator, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='num.png',Vmin=colmin, Vmax=colmax,ion=ionname,npix=ngrid,theta=thetaangle,phi=phiangle,redshift=redshift,extralabel=ionname)
coldens.make_colourmap('coldens_diff.%s.%s.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), result_num-result_den, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -0.5, 0.5, ionname, ngrid, redshift=redshift, extralabel=ionname'$_{\mathrm{NEQ-Equ}}')
print "result_num= ", result_num
print "result_den= ", result_den
Vmin=colmin,
Vmax=colmax,
ion=ion,
theta=thetaangle,
phi=phiangle,
redshift=redshift,
extralabel=ionname)
#make_colourmap(fig_name, ResultW, Xmin, Xmax, Ymin, Ymax, Vmin, Vmax, ion, npix, redshift)
###coldens.make_colourmap('coldens.%s.%s.%s.l%3.1f.png'%(snapname,ion,direction,lgrid), result_num-result_den, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, colmin, colmax, ion, ngrid, redshift=redshift, extralabel=ionname)
coldens.make_colourmap('coldens.%s.%s.%s.l%3.1f.png' %
(snapname, ion, direction, lgrid),
result_num - result_den,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
colmin,
colmax,
ion,
ngrid,
topbottomlabel=modelname,
extralabel=ionname)
print "result_num= ", result_num
print "result_den= ", result_den
N10lgrid = 0.0
nN10lgrid = 0
N05lgrid = 0.0
nN05lgrid = 0
N025lgrid = 0.0
nN025lgrid = 0
fig_name='coldens.%s.xraydens.x.l%3.1f.png' %
(snapname, lgrid),
Vmin=-30,
Vmax=30,
ion='dens',
phi=90,
extralabel='dens',
haloinfostr=haloinfostr,
docbar=docbar)
coldens.make_colourmap('coldens.%s.clump.x.l%3.1f.png' % (snapname, lgrid),
result_xraydens2_x - result_dens_x - result_dens_x,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
-1.0 * lgrid / 2.0,
lgrid / 2.0,
-5,
5,
'Xray',
ngrid,
extralabel='Xray',
haloinfostr=haloinfostr,
docbar=docbar)
result_hot_x = coldens.main(coords[indexes_xray],
hsmooth[indexes_xray],
mass[indexes_xray],
center,
lgrid,
lgrid,
lgridz,
ngrid,
o6 = chem[:,28]
o7 = chem[:,29]
o8 = chem[:,30]
mg2 = chem[:,45]
ne8 = chem[:,40]
mg10 = chem[:,53]
si12 = chem[:,68]
n5 = chem[:,19]
result_o6_z = coldens.main(coords, hsmooth, mass*hydrogen*o6, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.o6.z.l%3.1f.png'%(snapname,lgrid),Vmin=11, Vmax=15,ion='OVI',redshift=redshift,extralabel='OVI',haloinfostr=haloinfostr,docbar=docbar)
result_o7_z = coldens.main(coords, hsmooth, mass*hydrogen*o7, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.o7.z.l%3.1f.png'%(snapname,lgrid),Vmin=12.5, Vmax=16.5,ion='OVII',redshift=redshift,extralabel='OVII',haloinfostr=haloinfostr,docbar=docbar)
result_o8_z = coldens.main(coords, hsmooth, mass*hydrogen*o8, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.o8.z.l%3.1f.png'%(snapname,lgrid),Vmin=12.5, Vmax=16.5,ion='OVIII',redshift=redshift,extralabel='OVIII',haloinfostr=haloinfostr,docbar=docbar)
result_h_z = coldens.main(coords, hsmooth, mass*hydrogen, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.hydrogen.z.l%3.1f.png'%(snapname,lgrid),Vmin=17, Vmax=21,ion='H',redshift=redshift,extralabel='H',haloinfostr=haloinfostr,docbar=docbar)
result_o_z = coldens.main(coords, hsmooth, mass*oxygen/16.0, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.oxygen.z.l%3.1f.png'%(snapname,lgrid),Vmin=14, Vmax=18,ion='O',redshift=redshift,extralabel='O',haloinfostr=haloinfostr,docbar=docbar)
coldens.make_colourmap('coldens.%s.fo6.z.l%3.1f.png'%(snapname,lgrid), result_o6_z-result_o_z, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -4.0, -1.0, 'fOVI', ngrid, redshift=redshift, extralabel='f$_\mathrm{OVI}$',haloinfostr=haloinfostr,docbar=docbar)
coldens.make_colourmap('coldens.%s.fo7.z.l%3.1f.png'%(snapname,lgrid), result_o7_z-result_o_z, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -3.0, 0.0, 'fOVII', ngrid, redshift=redshift, extralabel='f$_\mathrm{OVII}$',haloinfostr=haloinfostr,docbar=docbar)
coldens.make_colourmap('coldens.%s.fo8.z.l%3.1f.png'%(snapname,lgrid), result_o8_z-result_o_z, -1.0*lgrid/2.0, lgrid/2.0, -1.0*lgrid/2.0, lgrid/2.0, -3.0, 0.0, 'fOVIII', ngrid, redshift=redshift, extralabel='f$_\mathrm{OVIII}$',haloinfostr=haloinfostr,docbar=docbar)
f = file('colion_map.%s.z.l%3.1f.dat'%(snapname,lgrid), 'w')
for i in xrange(ngrid):
for j in xrange(ngrid):
f.write('%3d %3d % 5.2f % 5.2f % 5.2f % 5.2f % 5.2f % 5.2f % 5.2f\n'%(i,j,result_h_z[i,j],result_o_z[i,j],result_h_z[i,j],result_o_z[i,j],result_o6_z[i,j],result_o7_z[i,j],result_o8_z[i,j]))
f.close()
write_radial_column('columnave.%s.z.%3.1f.dat'%(snapname,lgrid),result_h_z,result_o_z,result_o6_z,result_o7_z,result_o8_z,ngrid,lgrid,0)
write_radial_column('columnave_linsum.%s.z.%3.1f.dat'%(snapname,lgrid),10**result_h_z,10**result_o_z,10**result_o6_z,10**result_o7_z,10**result_o8_z,ngrid,lgrid,1)
result_o6_y = coldens.main(coords, hsmooth, mass*hydrogen*o6, center, lgrid, lgrid, lgridz, ngrid, ngrid, 58, boxsize,fig_name='coldens.%s.o6.y.l%3.1f.png'%(snapname,lgrid),Vmin=11, Vmax=15,ion='OVI',redshift=redshift,extralabel='OVI',haloinfostr=haloinfostr,docbar=docbar,theta=90)