Commit 1c55dfd9 authored by Pierre Augier's avatar Pierre Augier
Browse files

spatiotemporal: handle probe_region with floats

parent 286726b55f15
Pipeline #22421 passed with stage
in 16 minutes and 45 seconds
...@@ -50,6 +50,7 @@ params.init_fields.type = "noise" ...@@ -50,6 +50,7 @@ params.init_fields.type = "noise"
params.init_fields.noise.velo_max = 0.05 * Uf params.init_fields.noise.velo_max = 0.05 * Uf
params.oper.NO_SHEAR_MODES = True params.oper.NO_SHEAR_MODES = True
params.oper.NO_KY0 = True
# forcing params # forcing params
# key_forced = "rot_fft" # key_forced = "rot_fft"
......
...@@ -178,11 +178,11 @@ class SpatioTemporalSpectra3D(SpecificOutput): ...@@ -178,11 +178,11 @@ class SpatioTemporalSpectra3D(SpecificOutput):
if self.nb_dim == 3: if self.nb_dim == 3:
ikzmax = 4 ikzmax = 4
ikxmax = min(ikxmax, params.oper.nx // 2) ikxmax = min(int(round(ikxmax)), params.oper.nx // 2)
ikymax = min(ikymax, params.oper.ny // 2) ikymax = min(int(round(ikymax)), params.oper.ny // 2)
if self.nb_dim == 3: if self.nb_dim == 3:
ikzmax = min(ikzmax, params.oper.nz // 2) ikzmax = min(int(round(ikzmax)), params.oper.nz // 2)
self.probes_region = ikxmax, ikymax, ikzmax self.probes_region = ikxmax, ikymax, ikzmax
else: else:
self.probes_region = ikxmax, ikymax self.probes_region = ikxmax, ikymax
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment