Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix for walkthrough files #1294

Merged
merged 1 commit into from
Oct 15, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion amr-wind/utilities/sampling/PlaneSampler.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -65,7 +65,7 @@ void PlaneSampler::check_bounds()
const amrex::Vector<amrex::Real> points = {
point, point + m_axis1[d], point + m_axis2[d]};
for (const auto& pt : points) {
if ((pt <= prob_lo[d]) || (pt >= prob_hi[d])) {
if ((pt < prob_lo[d]) || (pt >= prob_hi[d])) {
amrex::Abort(
"PlaneSampler: Point out of domain. Redefine your "
"planes so they are completely inside the domain.");
Expand Down
28 changes: 14 additions & 14 deletions docs/sphinx/walkthrough/precursor_inp.txt
Original file line number Diff line number Diff line change
Expand Up @@ -30,7 +30,7 @@ transport.turbulent_prandtl = 0.3333
geometry.prob_lo = 0. 0. 0. # Lo corner coordinates
geometry.prob_hi = 2560. 2560. 1280. # Hi corner coordinates
amr.n_cell = 128 128 64 # Grid cells at coarsest AMRlevel
amr.max_level = 0 # Max AMR level in hierarchy
amr.max_level = 0 # Max AMR level in hierarchy
geometry.is_periodic = 1 1 0 # Periodicity x y z (0/1)

zlo.type = wall_model
Expand Down Expand Up @@ -74,21 +74,21 @@ sampling.output_frequency = 96
sampling.fields = velocity temperature

#---- sample defs ----
sampling.labels = xy-domain xz-domain
sampling.labels = xy-domain xz-domain

sampling.xy-domain.type = PlaneSampler
sampling.xy-domain.num_points = 256 256
sampling.xy-domain.origin = 0.0 0.0 86.5
sampling.xy-domain.axis1 = 2550.0 0.0 0.0
sampling.xy-domain.axis2 = 0.0 2550.0 0.0
sampling.xy-domain.offset_vector = 0.0 0.0 1.0
sampling.xy-domain.offsets = -63.45 0.0 63.45
sampling.xy-domain.type = PlaneSampler
sampling.xy-domain.num_points = 256 256
sampling.xy-domain.origin = 0.0 0.0 86.5
sampling.xy-domain.axis1 = 2550.0 0.0 0.0
sampling.xy-domain.axis2 = 0.0 2550.0 0.0
sampling.xy-domain.offset_vector = 0.0 0.0 1.0
sampling.xy-domain.offsets = -63.45 0.0 63.45

sampling.xz-domain.type = PlaneSampler
sampling.xz-domain.num_points = 256 128
sampling.xz-domain.origin = 0.0 1280.0 0.0
sampling.xz-domain.axis1 = 2550.0 0.0 0.0
sampling.xz-domain.axis2 = 0.0 0.0 1270.0
sampling.xz-domain.type = PlaneSampler
sampling.xz-domain.num_points = 256 128
sampling.xz-domain.origin = 0.0 1280.0 0.0
sampling.xz-domain.axis1 = 2550.0 0.0 0.0
sampling.xz-domain.axis2 = 0.0 0.0 1270.0

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# AVERAGING #
Expand Down
32 changes: 16 additions & 16 deletions docs/sphinx/walkthrough/spinup_inp.txt
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@ time.checkpoint_interval = 1800 # Steps between checkpoi

incflo.physics = ABL
turbulence.model = OneEqKsgsM84 # For neutral ABL, use "Smagorinsky"
TKE.source_terms = KsgsM84Src
TKE.source_terms = KsgsM84Src
incflo.gravity = 0. 0. -9.81 # Gravitational force (3D)
incflo.density = 1.225 # Reference density; make sure this agrees with OpenFAST values
transport.viscosity = 1.0e-5 # Dynamic viscosity [N-s/m^2]
Expand All @@ -24,7 +24,7 @@ transport.turbulent_prandtl = 0.3333
geometry.prob_lo = 0. 0. 0. # Lo corner coordinates
geometry.prob_hi = 2560. 2560. 1280. # Hi corner coordinates
amr.n_cell = 128 128 64 # Grid cells at coarsest AMRlevel
amr.max_level = 0 # Max AMR level in hierarchy
amr.max_level = 0 # Max AMR level in hierarchy
geometry.is_periodic = 1 1 0 # Periodicity x y z (0/1)

zlo.type = wall_model
Expand Down Expand Up @@ -63,22 +63,22 @@ ABL.surface_temp_flux = 0.05 # Surface temperature flux [K-m
incflo.post_processing = sampling

# --- Sampling parameters ---
sampling.output_frequency = 900
sampling.output_frequency = 900
sampling.fields = velocity temperature

#---- sample defs ----
sampling.labels = xy-domain xz-domain
sampling.labels = xy-domain xz-domain

sampling.xy-domain.type = PlaneSampler
sampling.xy-domain.num_points = 256 256
sampling.xy-domain.origin = 0.0 0.0 86.5
sampling.xy-domain.axis1 = 2550.0 0.0 0.0
sampling.xy-domain.axis2 = 0.0 2550.0 0.0
sampling.xy-domain.offset_vector = 0.0 0.0 1.0
sampling.xy-domain.offsets = -63.45 0.0 63.45
sampling.xy-domain.type = PlaneSampler
sampling.xy-domain.num_points = 256 256
sampling.xy-domain.origin = 0.0 0.0 86.5
sampling.xy-domain.axis1 = 2550.0 0.0 0.0
sampling.xy-domain.axis2 = 0.0 2550.0 0.0
sampling.xy-domain.offset_vector = 0.0 0.0 1.0
sampling.xy-domain.offsets = -63.45 0.0 63.45

sampling.xz-domain.type = PlaneSampler
sampling.xz-domain.num_points = 256 128
sampling.xz-domain.origin = 0.0 1280.0 0.0
sampling.xz-domain.axis1 = 2550.0 0.0 0.0
sampling.xz-domain.axis2 = 0.0 0.0 1270.0
sampling.xz-domain.type = PlaneSampler
sampling.xz-domain.num_points = 256 128
sampling.xz-domain.origin = 0.0 1280.0 0.0
sampling.xz-domain.axis1 = 2550.0 0.0 0.0
sampling.xz-domain.axis2 = 0.0 0.0 1270.0
Loading