Closed cbegeman closed 5 years ago
@xylar It'd be great to have you look over this PR. The modifications to coriolis terms are small for the slopes we're looking at. If you want to test this, you can use the same namelist files for the constant flux layer cases but set alpha_surface
equal to a nonzero slope in degrees.
When a surface slope is used (
alpha_surface
in degrees), the whole computational domain is rotated. In PR https://github.com/xylar/palm_les_lanl/pull/4, the gravity vector was rotated accordingly. In this PR, the coriolis parameters are also rotated. Stokes drift forces are not modified and are thus inaccurate for sloped domain cases.