I have done a 3D accretion disk simulation in cartesian coordinate. The
disk is centered at the origin and the disk plane is x-y plane. I want to
know the disk structure at a fixed radius R, i.e. sqrt(x*x+y*y) equals a
constant. It is something like the surface of a cylinder. A slice plot of
cylinder surface is helpful. But sliceplot only applies to a flat plane. So
how to do a slice plot of a cylinder surface?
Thanks a lot!
Department of Astronomy&Astrophysics,UC Santa Cruz
I am trying to analyze a unigrid simulation, 512^3 zones. Each dump has a size of ~ 7 GB.
I use the following script to get the angular momenta:
I am running it on a cluster with 1 core and 16GB. It typically runs fine for the first 2 dumps, then crashes on the third one when asking for AngularMomentumZ. I assume that I must not erase variables correctly. Does someone have any suggestion about how to optimize the script?
Also, would that make sense to run the script in parallel or would the root grid be duplicated on all processors?
Sorry for the naive questions, and thanks a lot for the help!
I downloaded the example script of light cone projectoin(
to my laptop and changed the Enzo simulation parameter file and redshifts.
But when I ran it the following error occurred:
File "light_cone_projection.py", line 19, in <module>
line 158, in __init__
line 44, in __init__
File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/convenience.py", line
131, in simulation
line 81, in __init__
line 343, in __init__
line 273, in _parse_parameter_file
param, vals = (i.strip() for i in line.split("="))
ValueError: need more than 1 value to unpack
Does the parameters "near_redshift" and "far_redshift" in LightCone need to
be the same with "CosmologyFinalRedshift" and "CosmologyInitialRedshift" in
enzo parameter file?
Does anyone know how to solve this?
Thank you all!