NCAR / pyDARTdiags

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

Bug: plot_profile may store the pressure level incorrectly #13

Open icastorm opened 3 months ago

icastorm commented 3 months ago

Issue

Users in the example code provided are encouraged to give pressure levels that define the bins in Pa. The plot_profile function intends to store the vertical pressure level in hPa, but then obtains these values by dividing by 1000 at line 101 of plots.py. The resulting incorrect vertical pressure levels can the be seen in the plots.

Proposed solutions

1) Divide by 100 instead of 1000 to correctly convert from Pa to hPa 2) Encourage users to give the bin ranges in hPa and convert to Pa when binning observations