Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Missing barometer height in BUFR files #218

Closed
ladsmund opened this issue Jan 2, 2024 · 2 comments
Closed

Missing barometer height in BUFR files #218

ladsmund opened this issue Jan 2, 2024 · 2 comments
Assignees

Comments

@ladsmund
Copy link
Contributor

ladsmund commented Jan 2, 2024

The BUFR files are missing heightOfBarometerAboveMeanSeaLevel when z_boom_u_smooth I nan.

pypromice/src/pypromice/postprocess/csv2bufr.py:

#Set measurement heights
if math.isnan(row['z_boom_u_smooth']) is False:
    codes_set(ibufr,
              '#1#heightOfSensorAboveLocalGroundOrDeckOfMarinePlatform',
              row['z_boom_u_smooth']-0.1) # For air temp and RH
    codes_set(ibufr,
              '#7#heightOfSensorAboveLocalGroundOrDeckOfMarinePlatform',
              row['z_boom_u_smooth']+0.4) # For wind speed
    if math.isnan(row['gps_alt_fit']) is False:
        codes_set(ibufr, 'heightOfBarometerAboveMeanSeaLevel',
                  row['gps_alt_fit']+row['z_boom_u_smooth']) # For pressure

I don't understand why z_boom_u_smooth is used to compute heightOfBarometerAboveMeanSeaLevel. z_boom_u_smooth is relative to the snow height and not the barometer height above sea level.
I expected to see a constant offset defining the barometer height relative height between the gps receiver:

barometer_height = row['gps_alt_fit'] + barometer_height_relative_to_gps_receiver
@ladsmund ladsmund self-assigned this Jan 2, 2024
@ladsmund
Copy link
Contributor Author

ladsmund commented Jan 2, 2024

@PennyHow @patrickjwright Do you have some inputs to this?

@PennyHow
Copy link
Member

PennyHow commented Jan 8, 2024

Looks like this is all resolved now with the merged PR. Great catch @ladsmund!

@PennyHow PennyHow closed this as completed Jan 8, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants