bwinkel / pycraf

pycraf is a package that provides functions and procedures for various tasks in spectrum-management compatibility studies.
35 stars 15 forks source link

atten_slant_annex1 function for different humidity values #5

Closed YuankangLiu closed 3 years ago

YuankangLiu commented 6 years ago

Hi! It seems that I can't type in humidity values in the atten_slant_annex functions. I can only use the humidity value calculated from the profile function. In real observations, these humidity values are observed in the observatory. I'd like to test different humidity values.

How can I do that?

bwinkel commented 6 years ago

Hi @YuankangLiu,

your finding is correct - the atten_slant_annex1 function, which is the more accurate version of the two, can only be used with a height profile. This is because a (measured) ground humidity value is linked with the physical conditions in the layers of the atmosphere above the ground layer. P.676 unfortunately gives no advice how one could handle such a case. If you know how to construct a suitable profile function, please let me know and I might implement it.

For the annex-2 version, you can provide a wet scale height, which depends on the humidity. Not sure, if you could use that instead. (What's the accuracy you need?)

By-the-way, I also wrote a Python wrapper of the C++ AATM package. It is not yet on GitHub but I could upload it, if you want (will take a couple days).

YuankangLiu commented 6 years ago

Thanks for your reply. I want to build a model which describes the influence of humidity on opacity. So in the end, I decided to use the am model written by Harvard. Thanks anyway.