3d1bee4532
Also update existing sea level pressure values stored in hdf5 files. (cherry picked from commit b733a89e8a6a2804bab9faaeee72496d8a65fba1 [formerly 73606213c72ba86a2a0f96232b1d1471fb6bd1c4]) Change-Id: I6924acf4fbe3fba7057f69306d582fc5c368812e Former-commit-id: a9f91757e4a9cc994b24c025d20aa63433e7c2f7
34 lines
904 B
Python
Executable file
34 lines
904 B
Python
Executable file
#!/usr/bin/env python2
|
|
|
|
# #5345
|
|
# Convert sea level pressure from hPa to Pa in an hdf5 file
|
|
# Do nothing if all values are already in Pa
|
|
|
|
import sys
|
|
import h5py
|
|
|
|
def main():
|
|
if len(sys.argv) != 2:
|
|
print "usage: {} filename.h5".format(sys.argv[0])
|
|
sys.exit(1)
|
|
|
|
didstuff = False
|
|
try:
|
|
with h5py.File(sys.argv[1], 'r+') as f:
|
|
if 'seaLevelPress' in f:
|
|
for i, data in enumerate(f['seaLevelPress']):
|
|
if data > 0 and data < 1999:
|
|
f['seaLevelPress'][i] = data * 100.0
|
|
didstuff = True
|
|
except Exception as e:
|
|
print "ERROR: " + str(sys.exc_info()[0]) + ": " + str(e)
|
|
sys.exit(1)
|
|
|
|
if didstuff:
|
|
print "INFO: {}: updated".format(sys.argv[1])
|
|
else:
|
|
print "INFO: {}: no update needed".format(sys.argv[1])
|
|
|
|
|
|
if __name__ == "__main__":
|
|
main()
|