Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Rayleigh damping near lateral boundaries (hrdamp = 1) bug

prashn

New member
This concerns the parameter hrdamp in the namelist. Setting this parameter to 1 enables horizontal damping for the u, v, and w velocity fields according to the README.namelist file. However, if you look at the source code in param.F (line 8706), only the variable tauf is initialized when hrdamp = 1 which acts on the w field and not u and v in solve1.F (line 403).

Either the README.namelist needs to be updated stating that lateral damping only acts on w or the variable tauh must be initialized in the same manner as tauf in param.F.
 
Top