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.
 
Back
Top