Ryan May

I'm not shocked that updating the formula causes test failures; you are producing different values, and thus plots.

None that I'm aware of. You could try comparing curves to e.g. [those from the National Weather Service](https://www.weather.gov/source/zhu/ZHU_Training_Page/convective_parameters/skewt/skewtinfo.html), but it's usually really hard to match since you don't know what...

@robwarrenwx I appreciate the feedback. I think using a temperature-dependent Lv in `moist_lapse()` is definitely something we can consider as we continue to refine our thermo calcs, e.g. implement Romps...

Is there benefit to publishing these when the netcdf c library is not thread safe and calling netcdf functions from multiple threads is likely to lead to segfaults?

Seems reasonable. Just wanted to make sure expectations were in line with reality.

I'm happy to update the MetPy parser if there are general problems that occur in the original products we should account for. My concern going beyond that would be that...

I think there is an enhancement we can make here in terms of being more fault-tolerant, if not aggressively cleaning. Since we're iterating line-by-line, we should catch `ValueError`s and issue...

@josin77 The problem with your suggested approach is that to calculate the parcel profile directly as a function of height would require a version of `moist_lapse` that uses height instead...

Thanks @akrherz ! That's really helpful and we can definitely include a few more fixes here.

@bvmcode That makes sense since it looks like the issue is in MetPy's xarray handling combined with what `tke` is doing.