NCAR / geocat-comp

GeoCAT-comp provides implementations of computational functions for operating on geosciences data. Many of these functions originated in NCL and were translated into Python.
https://geocat-comp.readthedocs.io
Apache License 2.0
122 stars 55 forks source link

delta_pressure clarification question #635

Open kafitzgerald opened 1 month ago

kafitzgerald commented 1 month ago

Hello I am confuse, does computing thickness between two pressure level here based on hypsometric equation or its jjust simple difference between two pressure level. As in the shared result above , the 3rd value why it is 2500 if its the difference between 85000-80000, similarly for the 16th value why its 3750 and so on. Could you kindly please clarify this. I am completely confuse.

Originally posted by @bnepal22 in https://github.com/NCAR/geocat-comp/issues/538#issuecomment-2257575928

kafitzgerald commented 1 month ago

Hi @bnepal22, thanks for your question!

I went ahead and moved your comment to a new issue. For new clarification and usage questions about geocat-comp functions we recommend starting a discussion post or creating a new issue and linking to any relevant content rather than commenting on an existing one. For the most part you shouldn't need to worry about mentioning "@" people either unless there's a particular person you'd like to notify.

As for your question about delta_pressure, you can find a bit more detailed information in the NCL documentation of the function this replicates. It's computing the thickness (roughly a simple difference) with some exceptions for handling the top of the atmosphere and the surface. Let me know if you need further clarification on this.

This is probably something we could make a bit more clear in the documentation.

bnepal22 commented 1 month ago

I'm sorry for not getting back to you sooner. Further clarification would be appreciated.

kafitzgerald commented 1 month ago

Sure, let us know if there's anything in particular that would be helpful.

You can see here how delta_pressure is calculated in geocat-comp. It's just altering the calculations to use a specified top and surface pressure. https://github.com/NCAR/geocat-comp/blob/c3df44f40b436c4bbec59d08a9b9352dc2e6715d/geocat/comp/meteorology.py#L1382-L1448

We'll work on getting this updated in the documentation so it's more clear and I'll leave this issue open until that's been addressed. It may take a moment to get to this though.