%matplotlib inline
import yt
from yt.utilities.fits_image import FITSProjection
ds = yt.load("GasSloshing/sloshing_nomag2_hdf5_plt_cnt_0100")
yt : [WARNING ] 2014-05-05 11:06:12,691 integer runtime parameter checkpointfilenumber overwrites a simulation scalar of the same name yt : [WARNING ] 2014-05-05 11:06:12,691 integer runtime parameter plotfilenumber overwrites a simulation scalar of the same name yt : [INFO ] 2014-05-05 11:06:12,709 Parameters: current_time = 7.89058001997e+16 yt : [INFO ] 2014-05-05 11:06:12,709 Parameters: domain_dimensions = [16 16 16] yt : [INFO ] 2014-05-05 11:06:12,711 Parameters: domain_left_edge = [ -3.70272000e+24 -3.70272000e+24 -3.70272000e+24] yt : [INFO ] 2014-05-05 11:06:12,711 Parameters: domain_right_edge = [ 3.70272000e+24 3.70272000e+24 3.70272000e+24] yt : [INFO ] 2014-05-05 11:06:12,712 Parameters: cosmological_simulation = 0.0
prj = yt.ProjectionPlot(ds, "z", ["density"], width=(500,"kpc"))
prj.show()
yt : [INFO ] 2014-05-05 11:06:13,258 Loading field plugins. yt : [INFO ] 2014-05-05 11:06:13,259 Loaded angular_momentum (8 new fields) yt : [INFO ] 2014-05-05 11:06:13,259 Loaded astro (14 new fields) yt : [INFO ] 2014-05-05 11:06:13,260 Loaded cosmology (20 new fields) yt : [INFO ] 2014-05-05 11:06:13,261 Loaded fluid (54 new fields) yt : [INFO ] 2014-05-05 11:06:13,262 Loaded fluid_vector (86 new fields) yt : [INFO ] 2014-05-05 11:06:13,263 Loaded geometric (101 new fields) yt : [INFO ] 2014-05-05 11:06:13,263 Loaded local (101 new fields) yt : [INFO ] 2014-05-05 11:06:13,264 Loaded magnetic_field (107 new fields) yt : [INFO ] 2014-05-05 11:06:13,264 Loaded species (107 new fields) yt : [INFO ] 2014-05-05 11:06:15,851 Projection completed yt : [INFO ] 2014-05-05 11:06:15,862 xlim = -771419395240581203492864.000000 771419395240581203492864.000000 yt : [INFO ] 2014-05-05 11:06:15,863 ylim = -771419395240581203492864.000000 771419395240581203492864.000000 yt : [INFO ] 2014-05-05 11:06:15,864 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800 yt : [INFO ] 2014-05-05 11:06:15,878 xlim = -771419395240581203492864.000000 771419395240581203492864.000000 yt : [INFO ] 2014-05-05 11:06:15,878 ylim = -771419395240581203492864.000000 771419395240581203492864.000000 yt : [INFO ] 2014-05-05 11:06:15,879 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800 yt : [INFO ] 2014-05-05 11:06:15,894 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800
fits_proj = FITSProjection(ds, "z", ["density"])
fits_proj.writeto("proj.fits", clobber=True)
yt : [INFO ] 2014-05-05 11:06:17,898 Projection completed yt : [INFO ] 2014-05-05 11:06:17,918 Making a FITS image of field density yt : [INFO ] 2014-05-05 11:06:17,919 Making a fixed resolution buffer of (density) 1024 by 1024 WARNING: Overwriting existing file 'proj.fits'. [astropy.io.fits.file] astropy: [WARNING ] 2014-05-05 11:06:17,965 Overwriting existing file 'proj.fits'.
ds = yt.load("proj.fits")
yt : [WARNING ] 2014-05-05 11:06:18,145 Cannot find time yt : [INFO ] 2014-05-05 11:06:18,206 Found length units of kpc. yt : [INFO ] 2014-05-05 11:06:18,214 Parameters: current_time = 0.0 yt : [INFO ] 2014-05-05 11:06:18,214 Parameters: domain_dimensions = [1024 1024 1] yt : [INFO ] 2014-05-05 11:06:18,215 Parameters: domain_left_edge = [ 0.5 0.5 0.5] yt : [INFO ] 2014-05-05 11:06:18,217 Parameters: domain_right_edge = [ 1024.5 1024.5 1.5] yt : [INFO ] 2014-05-05 11:06:18,219 Parameters: cosmological_simulation = 0.0
sp = ds.sphere("c", (500,"kpc"))
prof = yt.create_profile(sp, "radius", ["density"], weight_field="ones", n_bins=50, logs={"radius":False})
pp = yt.ProfilePlot(sp, "radius", ["density"], profiles=[prof])
pp.set_log("radius", True)
pp.show()