2. Hydro: Load Selected Variables and Data Ranges
Simulation Overview
using Mera
info = getinfo(420, "../../testing/simulations/manu_sim_sf_L10");┌ Info: Precompiling Mera [02f895e8-fdb1-4346-8fe6-c721699f5126]
└ @ Base loading.jl:1273
*__ __ _______ ______ _______
| |_| | | _ | | _ |
| | ___| | || | |_| |
| | |___| |_||_| |
| | ___| __ | |
| ||_|| | |___| | | | _ |
|_| |_|_______|___| |_|__| |__|
[Mera]: 2020-02-08T13:41:33.944
Code: RAMSES
output [420] summary:
mtime: 2017-07-27T01:22:09
ctime: 2019-12-24T09:57:04.822
=======================================================
simulation time: 624.91 [Myr]
boxlen: 48.0 [kpc]
ncpu: 1024
ndim: 3
-------------------------------------------------------
amr: true
level(s): 6 - 10 --> cellsize(s): 750.0 [pc] - 46.88 [pc]
-------------------------------------------------------
hydro: true
hydro-variables: 6 --> (:rho, :vx, :vy, :vz, :p, :var6)
hydro-descriptor: (:density, :velocity_x, :velocity_y, :velocity_z, :thermal_pressure, :passive_scalar_1)
γ: 1.01
-------------------------------------------------------
gravity: true
gravity-variables: (:epot, :ax, :ay, :az)
-------------------------------------------------------
particles: true
particle variables: (:vx, :vy, :vz, :mass, :birth)
-------------------------------------------------------
clumps: true
clump-variables: (:index, :lev, :parent, :ncell, :peak_x, :peak_y, :peak_z, Symbol("rho-"), Symbol("rho+"), :rho_av, :mass_cl, :relevance)
-------------------------------------------------------
namelist-file: false
timer-file: false
compilation-file: true
makefile: true
patchfile: true
=======================================================Select Variables
Choose from the existing hydro variables listed in the simulation-info. Use the quoted Symbols: :varn1 or :cpu (=neg. one), :var1 or :rho, :var2 or :vx, :var3 or :vy, :var4 or :vz, :var5 or :p. Variables above 5 can be selected by :var6, :var7 etc. . No order is required. The selection of the variable's names from the descriptor files will be implemented in the future.
Read all variables (default)
gas = gethydro(info); [Mera]: Get hydro data: 2020-02-08T13:08:41.471
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
domain:
xmin::xmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
ymin::ymax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
zmin::zmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:48
Memory used for data table :85.94877052307129 MB
-------------------------------------------------------gas.dataTable with 1126532 rows, 10 columns:
Columns:
# colname type
────────────────────
1 level Int64
2 cx Int64
3 cy Int64
4 cz Int64
5 rho Float64
6 vx Float64
7 vy Float64
8 vz Float64
9 p Float64
10 var6 Float64Select several variables w/o a keyword
gas_a = gethydro(info, vars=[:rho, :p], smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:09:34.25
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 5) = (:rho, :p)
domain:
xmin::xmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
ymin::ymax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
zmin::zmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:40
Memory used for data table :51.5693416595459 MB
-------------------------------------------------------The same variables can be read by using the var-number:
gas_a = gethydro(info, vars=[:var1, :var5], smallr = 1e-4 / info.unit_d); [Mera]: Get hydro data: 2020-02-08T13:10:15.474
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 5) = (:rho, :p)
domain:
xmin::xmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
ymin::ymax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
zmin::zmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:50
Memory used for data table :51.5693416595459 MB
-------------------------------------------------------A keyword argument for the variables is not needed if the following order is preserved: InfoType-object, variables:
gas_a = gethydro(info, [:rho, :p], smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:11:05.814
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 5) = (:rho, :p)
domain:
xmin::xmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
ymin::ymax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
zmin::zmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:53
Memory used for data table :51.5693416595459 MB
-------------------------------------------------------gas_a.dataTable with 1126532 rows, 6 columns:
level cx cy cz rho p
────────────────────────────────────────────
6 1 1 1 0.0001 2.85781e-7
6 1 1 2 0.0001 2.85781e-7
6 1 1 3 0.0001 2.85781e-7
6 1 1 4 0.0001 2.85781e-7
6 1 1 5 0.0001 2.85781e-7
6 1 1 6 0.0001 2.85781e-7
6 1 1 7 0.0001 2.85781e-7
6 1 1 8 0.0001 2.85781e-7
6 1 1 9 0.0001 2.85781e-7
6 1 1 10 0.0001 2.85781e-7
6 1 1 11 0.0001 2.85781e-7
6 1 1 12 0.0001 2.85781e-7
⋮
10 822 507 516 0.0305045 0.000589997
10 822 508 511 0.0551132 0.00106596
10 822 508 512 0.0551132 0.00106596
10 822 508 513 0.0845289 0.0016349
10 822 508 514 0.0788161 0.00152441
10 822 508 515 0.0305045 0.000589997
10 822 508 516 0.0305045 0.000589997
10 822 509 513 0.0861783 0.0016668
10 822 509 514 0.0861783 0.0016668
10 822 510 513 0.0861783 0.0016668
10 822 510 514 0.0861783 0.0016668Select one variable
In this case, no array and keyword is necessary, but preserve the following order: InfoType-object, variable:
gas_c = gethydro(info, :vx ); [Mera]: Get hydro data: 2020-02-08T13:11:59.783
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(2,) = (:vx,)
domain:
xmin::xmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
ymin::ymax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
zmin::zmax: 0.0 :: 1.0 ==> 0.0 [kpc] :: 48.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:43
Memory used for data table :42.97448444366455 MB
-------------------------------------------------------gas_c.dataTable with 1126532 rows, 5 columns:
level cx cy cz vx
────────────────────────────────
6 1 1 1 0.0990548
6 1 1 2 0.110889
6 1 1 3 0.117122
6 1 1 4 0.120163
6 1 1 5 0.122688
6 1 1 6 0.124978
6 1 1 7 0.12718
6 1 1 8 0.129342
6 1 1 9 0.131503
6 1 1 10 0.133663
6 1 1 11 0.135799
6 1 1 12 0.137892
⋮
10 822 507 516 -0.0840177
10 822 508 511 -0.08318
10 822 508 512 -0.08318
10 822 508 513 -0.0948844
10 822 508 514 -0.0986642
10 822 508 515 -0.0840177
10 822 508 516 -0.0840177
10 822 509 513 -0.136379
10 822 509 514 -0.136379
10 822 510 513 -0.136379
10 822 510 514 -0.136379Select Spatial Ranges
Use RAMSES Standard Notation
Ranges correspond to the domain [0:1]^3 and are related to the box corner at [0., 0., 0.] by default. Here, we limit the loading of the data to a maximum level of 8:
gas = gethydro(info, lmax=8,
xrange=[0.2,0.8],
yrange=[0.2,0.8],
zrange=[0.4,0.6],
smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:12:44.551
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
domain:
xmin::xmax: 0.2 :: 0.8 ==> 9.6 [kpc] :: 38.4 [kpc]
ymin::ymax: 0.2 :: 0.8 ==> 9.6 [kpc] :: 38.4 [kpc]
zmin::zmax: 0.4 :: 0.6 ==> 19.2 [kpc] :: 28.8 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:35
Memory used for data table :23.961313247680664 MB
-------------------------------------------------------The loaded data ranges are assigned to the field ranges as an array in RAMSES standard notation (domain: [0:1]^3):
gas.ranges6-element Array{Float64,1}:
0.2
0.8
0.2
0.8
0.4
0.6Ranges relative to a given center:
gas = gethydro(info, lmax=8,
xrange=[-0.3, 0.3],
yrange=[-0.3, 0.3],
zrange=[-0.1, 0.1],
center=[0.5, 0.5, 0.5],
smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:13:20.861
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
center: [0.5, 0.5, 0.5] ==> [24.0 [kpc] :: 24.0 [kpc] :: 24.0 [kpc]]
domain:
xmin::xmax: 0.2 :: 0.8 ==> 9.6 [kpc] :: 38.4 [kpc]
ymin::ymax: 0.2 :: 0.8 ==> 9.6 [kpc] :: 38.4 [kpc]
zmin::zmax: 0.4 :: 0.6 ==> 19.2 [kpc] :: 28.8 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:36
Memory used for data table :23.961313247680664 MB
-------------------------------------------------------Use notation in physical units
In the following example the ranges are given in unit "kpc", relative to the box corner 0., 0., 0.:
gas = gethydro(info, lmax=8,
xrange=[2.,22.],
yrange=[2.,22.],
zrange=[22.,26.],
range_unit=:kpc,
smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:13:57.226
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
domain:
xmin::xmax: 0.0416667 :: 0.4583333 ==> 2.0 [kpc] :: 22.0 [kpc]
ymin::ymax: 0.0416667 :: 0.4583333 ==> 2.0 [kpc] :: 22.0 [kpc]
zmin::zmax: 0.4583333 :: 0.5416667 ==> 22.0 [kpc] :: 26.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:32
Memory used for data table :4.347360610961914 MB
-------------------------------------------------------The possible physical length units for the keyword range_unit are defined in the field scale :
viewfields(info.scale) # or e.g.: gas.info.scale [Mera]: Fields to scale from user/code units to selected units
=======================================================================
Mpc = 0.0010000000000006482
kpc = 1.0000000000006481
pc = 1000.0000000006482
mpc = 1.0000000000006482e6
ly = 3261.5637769461323
Au = 2.0626480623310105e23
km = 3.0856775812820004e16
m = 3.085677581282e19
cm = 3.085677581282e21
mm = 3.085677581282e22
μm = 3.085677581282e25
Msol_pc3 = 0.9997234790001649
g_cm3 = 6.76838218451376e-23
Msol_pc2 = 999.7234790008131
g_cm2 = 0.20885045168302602
Gyr = 0.014910986463557083
Myr = 14.910986463557084
yr = 1.4910986463557083e7
s = 4.70554946422349e14
ms = 4.70554946422349e17
Msol = 9.99723479002109e8
Mearth = 3.329677459032007e14
Mjupiter = 1.0476363431814971e12
g = 1.9885499720830952e42
km_s = 65.57528732282063
m_s = 65575.28732282063
cm_s = 6.557528732282063e6
nH = 30.987773856809987
erg = 8.551000140274429e55
g_cms2 = 2.9104844143584656e-9
T_mu = 517028.3199143136
Ba = 2.910484414358466e-9Ranges relative to a given center e.g. in unit "kpc":
gas = gethydro(info, lmax=8,
xrange=[-16.,16.],
yrange=[-16.,16.],
zrange=[-2.,2.],
center=[24.,24.,24.],
range_unit=:kpc,
smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:41:54.547
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
center: [0.5, 0.5, 0.5] ==> [24.0 [kpc] :: 24.0 [kpc] :: 24.0 [kpc]]
domain:
xmin::xmax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
ymin::ymax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
zmin::zmax: 0.4583333 :: 0.5416667 ==> 22.0 [kpc] :: 26.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:37
Memory used for data table :22.76700782775879 MB
-------------------------------------------------------Use the short notation for the box center :bc or :boxcenter for all dimensions (x,y,z):
gas = gethydro(info, lmax=8,
xrange=[-16., 16.],
yrange=[-16., 16.],
zrange=[-2., 2.],
center=[:boxcenter],
range_unit=:kpc,
smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:42:56.283
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
center: [0.5, 0.5, 0.5] ==> [24.0 [kpc] :: 24.0 [kpc] :: 24.0 [kpc]]
domain:
xmin::xmax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
ymin::ymax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
zmin::zmax: 0.4583333 :: 0.5416667 ==> 22.0 [kpc] :: 26.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:28
Memory used for data table :22.76700782775879 MB
-------------------------------------------------------gas = gethydro(info, lmax=8,
xrange=[-16., 16.],
yrange=[-16., 16.],
zrange=[-2., 2.],
center=[:bc],
range_unit=:kpc,
smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:15:01.751
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
center: [0.5, 0.5, 0.5] ==> [24.0 [kpc] :: 24.0 [kpc] :: 24.0 [kpc]]
domain:
xmin::xmax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
ymin::ymax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
zmin::zmax: 0.4583333 :: 0.5416667 ==> 22.0 [kpc] :: 26.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:31
Memory used for data table :22.76700782775879 MB
-------------------------------------------------------Use the box center notation for individual dimensions, here x,z:
gas = gethydro(info, lmax=8,
xrange=[-16., 16.],
yrange=[-16., 16.],
zrange=[-2., 2.],
center=[:bc, 24., :bc],
range_unit=:kpc,
smallr = 1e-4); [Mera]: Get hydro data: 2020-02-08T13:44:27.734
Key vars=(:level, :cx, :cy, :cz)
Using var(s)=(1, 2, 3, 4, 5, 6) = (:rho, :vx, :vy, :vz, :p, :var6)
center: [0.5, 0.5, 0.5] ==> [24.0 [kpc] :: 24.0 [kpc] :: 24.0 [kpc]]
domain:
xmin::xmax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
ymin::ymax: 0.1666667 :: 0.8333333 ==> 8.0 [kpc] :: 40.0 [kpc]
zmin::zmax: 0.4583333 :: 0.5416667 ==> 22.0 [kpc] :: 26.0 [kpc]
Reading data...
100%|███████████████████████████████████████████████████| Time: 0:00:27
Memory used for data table :22.76700782775879 MB
-------------------------------------------------------