Skip to content
New issue

Have a question about this project? # for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “#”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? # to your account

Add notes about power law interpolation vs linear #186

Merged
merged 1 commit into from
Oct 18, 2024

Conversation

ppinchuk
Copy link
Collaborator

@ppinchuk ppinchuk commented Sep 5, 2024

@grantbuster Did I miss anything in the main text that should be added?

I tried to place the note in a couple different places so it would be easier to find, but let me know if you think I missed a good place to put it.

@ppinchuk ppinchuk added the documentation Improvements or additions to documentation label Sep 5, 2024
@ppinchuk ppinchuk requested a review from grantbuster September 5, 2024 16:27
@ppinchuk ppinchuk self-assigned this Sep 5, 2024
@grantbuster
Copy link
Member

@bnb32 for your awareness: the original WTK code does linear interpolation between WRF vertical levels. I was confused because the WTK report says power law interp and we had previously discussed the benefits of such an approach. I reached out to Michael Rossol who originally implemented WTK processing and he came back with this explanation why we want to use linear. This might change how we think about processing ERA data. Power law might be okay because the model might learn to produce negative wind shear, but linear interp might better encourage this.

@bnb32
Copy link
Collaborator

bnb32 commented Sep 9, 2024

@bnb32 for your awareness: the original WTK code does linear interpolation between WRF vertical levels. I was confused because the WTK report says power law interp and we had previously discussed the benefits of such an approach. I reached out to Michael Rossol who originally implemented WTK processing and he came back with this explanation why we want to use linear. This might change how we think about processing ERA data. Power law might be okay because the model might learn to produce negative wind shear, but linear interp might better encourage this.

You can still get negative wind shear if you use a negative exponent though. I was using log interp which does allow for negative wind shear in a similar way, but I switched back to lin interp a since I didn't see much difference. Interesting that the wtk paper says power law interp though.

@grantbuster
Copy link
Member

Notes from our discussion: difference between log/power law fit and log/power law interpolation between WRF levels. If we fit a power law exponent between two WRF levels we can have a continuous power law interpolation, but it might behave poorly with sharp gradients at local maximum windspeeds. Based on this, i think we can proceed with linear interpolation under the argument that it is guaranteed to be a well-behaved interpolation method under negative wind shear conditions.

@grantbuster grantbuster merged commit 5861f85 into main Oct 18, 2024
15 of 16 checks passed
@grantbuster grantbuster deleted the pp/interpolation_notes branch October 18, 2024 19:53
github-actions bot pushed a commit that referenced this pull request Oct 18, 2024
Add notes about power law interpolation vs linear
# for free to join this conversation on GitHub. Already have an account? # to comment
Labels
documentation Improvements or additions to documentation
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants