Bug: interp_sigma_to_hybrid and WRF_Chem output

[anissa111]

Hi @erogluorhan , sorry for the late response. It took me some time before I could test interp_sigma_to_hybrid. However, I am facing two difficulties in using this file for WRF_Chem output

1. With the current implementation, it seems not possible to interpolate concentration fields directly. By default, the concentration fields have units of ppm which are not recognised by metpy. We need to convert the units into something that metpy recognizes (e.g. molecules cm-3). See the discussion here.
2. The current implementation is very slow as compared to ncl. For example, I am interpolating my WRF-Chem absolute humidity fields on TM5 grid. While with ncl, on my machine it took me less than 1 second, with interp_sigma_to_hybrid it took me more than 200s. For your reference, I have attached my ncl script, python scripts and a sample data file here.

Best regards and thanks a lot,
Vinod

Originally posted by @vinodkatmos in #142 (comment)

[anissa111]

@vinodkatmos, just wanted to give you an update on this issue you submitted.

1. I'm not entirely sure where the best location for this unit change would occur. If possible, could you point me to metpy's implementation or a usage example so I could better understand?
2. We will definitely be looking into optimizing this (and other!) functionality in the near future. This issue should provide links to that progress when it occurs.

Thank you again for taking the time to provide feedback to our project!

[anissa111]

I'm going to close this issue, as the performance issues will be included in a different issue shortly.

[erogluorhan]

Hey @anissa111 it was closed on Feb 21 but is now assigned to me. I am happy to look into it but just wanted to understand correctly if the intent is actually to work on this now.

[anissa111]

Oh, whoops! I think this was a mis-click while trying to archive it in our PL meeting yesterday my bad!