dtcenter / MET

Model Evaluation Tools
https://dtcenter.org/community-code/model-evaluation-tools-met
Apache License 2.0
77 stars 24 forks source link

Bugfix #2902 main_v11.1 pb2nc_pbl #2903

Closed JohnHalleyGotway closed 4 months ago

JohnHalleyGotway commented 4 months ago

Expected Differences

Pull Request Testing

In addition, tested using main_v11.1 from the nightly build with -v 10. For station 71722 I see the following log messages:

DEBUG 8: compute_pbl() input to calpbl_ (buffer): index, P, Q, T, Z, U, V
DEBUG 8: compute_pbl()   0   15607511481187303424.00000 0.00000 -3152171068026237978910130176.00000 -0.00000 6.50000 -1.70000
DEBUG 8: compute_pbl()   1   93960.00000 0.00659 286.96384 610.00000 6.50000 -1.70000
DEBUG 8: compute_pbl()   2   97400.00000 0.00840 289.29169 305.00000 2.00000 0.70000
DEBUG 8: compute_pbl()   3   98700.00000 0.00906 290.14999 192.00000 -9999.00000 -9999.00000
DEBUG 7:  hpbl: 610, obs_arr[4]: 418   lat: 46.29999924, lon: -76, elv: 192 20190615_000000 ADPUPA 71722

Note the obvious problem with the pressure and temperature values for record index 0. This is the result of the bug.

Testing with this bugfix branch, I get the following corrected log messages which are an accurate reversal of the input pressure level values:

DEBUG 8: compute_pbl() input to calpbl_ (buffer): index, P, Q, T, Z, U, V
DEBUG 8: compute_pbl()   0   92500.00000 0.00580 285.95001 735.00000 6.50000 -1.70000
DEBUG 8: compute_pbl()   1   93960.00000 0.00659 286.96384 610.00000 6.50000 -1.70000
DEBUG 8: compute_pbl()   2   97400.00000 0.00840 289.29169 305.00000 2.00000 0.70000
DEBUG 8: compute_pbl()   3   98700.00000 0.00906 290.14999 192.00000 -9999.00000 -9999.00000
DEBUG 7:  hpbl: 735, obs_arr[4]: 543   lat: 46.29999924, lon: -76, elv: 192 20190615_000000 ADPUPA 71722

Pull Request Checklist

See the METplus Workflow for details.