fix wind scaling problem in prsflg==2 #5
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
In mod_random_forcing, the prsflg==2 option scales the wind perturbation magnitude to that specified in the namelist. The scaling analysis use the pressure perturbation magnitude and the pressure wind relationship to compute the expected wind magnitude. But the grid spacing used in the code is minscpx which is not spatially varying. Grid spacing is varying (scpx and scpy) in the computation of pressure gradients dpresx and dpresy, therefore changed to using varying grid spacing in wprsfac computation, so that resulting wind perturbation magnitude won't be affected by varying dx.
|
the specified wind speed amplitude still doesn't match the output exactly. The wprsfac is computed by "pressure_amp / hcorr / fcor /rhoa" the pressure amplitude doesn't exactly drop to zero over the distance "hcorr", roughly half that distance instead. So adding a compensating factor should improve the result, instead of "rhdx" use "0.54rh*dx" as the scale distance. |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
In mod_random_forcing, the prsflg==2 option scales the wind perturbation magnitude to that specified in the namelist. The scaling analysis use the pressure perturbation magnitude and the pressure wind relationship to compute the expected wind magnitude. But the grid spacing used in the code is minscpx which is not spatially varying. Grid spacing is varying (scpx and scpy) in the computation of pressure gradients dpresx and dpresy, therefore changed to using varying grid spacing in wprsfac computation, so that resulting wind perturbation magnitude won't be affected by varying dx.