NCAR / pyDARTdiags

Python observation space diagnostics for the Data Assimilation Research Testbed (DART)
Apache License 2.0
3 stars 1 forks source link

Feat request: Plot profiles of RMSE and bias of obs with the height vertical coordinate (and general vertical coordinate handling) #14

Open icastorm opened 2 months ago

icastorm commented 2 months ago

Description

While most observations in DART are likely to use the pressure vertical coordinate, the height vertical coordinate is also often used. In my previous experiments, for example, I used height exclusively. Some sets of obs may even use both, depending on the types of observations being used.

Currently, the profile plot can only handle observations with the pressure vertical coordinate. I'd like to see the plotting function able to handle obs_seq files with the height vertical coordinate as well.

Solution(s)

icastorm commented 2 months ago

This is a fairly easy add, so I can work on a branch with this feature and test it some time this week!

icastorm commented 2 months ago

I have a version of this working, will look at getting the plotting tests up to speed next. Currently, the new version of plot_profile defaults to using pressure coordinates, but can take in any vert_unit value desired. It then slices the given data down to just the rows that use the selected coordinate and returns an error if none are found.