search

lanl / palm_lanl

LANL Contributions to PArallelized Large-eddy simulation Model (PALM)
2 stars 5 forks source link

Wall at bottom boundary #35

Closed cbegeman closed 5 years ago

cbegeman commented 5 years ago

@vanroekel @lconlon @qingli411 @xylar

As far as I can tell, a wall at the bottom boundary is hard-coded here: https://github.com/xylar/palm_les_lanl/blob/42c59e36d9fffe82b0f4b3159be04f1860408ef7/trunk/SOURCE/init_grid.f90#L1911-L1912 because the nzb+1 index is used instead of nzb. topo is then used to define wall_flags_0.

I propose a new namelist parameter that designates whether the bottom or the top boundary is a wall, rather than hard-coding a wall at the bottom boundary or implementing a wall only if the boundary is associated with a constant flux layer.

Thoughts?

cbegeman commented 5 years ago

The main effect seems to be decreasing the mixing length near the bottom boundary of the domain, which @vanroekel thinks is appropriate. Therefore, no changes to the code along these lines will be pursued.