35
35
use multiprocess h5py reader from io utilities module
36
36
use updated ATL11 readers within icesat2_toolkit
37
37
use pathlib to define and operate on paths
38
- use raster file to define valid grid cells in mask
38
+ use raster file to define valid grid cells in mask
39
39
include cell_area in output HDF5 file
40
40
Written 03/2021
41
41
"""
@@ -378,7 +378,7 @@ def tidal_histogram(tile_file,
378
378
# histogram bin
379
379
attributes ['bins' ] = {}
380
380
attributes ['bins' ]['long_name' ] = 'Histogram bins'
381
- attributes ['bins' ]['units' ] = 'meters '
381
+ attributes ['bins' ]['units' ] = '1 '
382
382
attributes ['bins' ]['description' ] = \
383
383
'Center of each height difference histogram bin'
384
384
fill_value ['bins' ] = None
@@ -390,10 +390,11 @@ def tidal_histogram(tile_file,
390
390
attributes ['cell_area' ]['units' ] = 'm^2'
391
391
attributes ['cell_area' ]['coordinates' ] = 'y x'
392
392
attributes ['cell_area' ]['grid_mapping' ] = 'crs'
393
- fill_value ['cell_area' ] = 0
393
+ fill_value ['cell_area' ] = 0
394
394
# height difference histogram
395
395
attributes ['dh_hist' ] = {}
396
396
attributes ['dh_hist' ]['long_name' ] = 'Height difference histogram'
397
+ attributes ['dh_hist' ]['description' ] = 'Histogram of height differences'
397
398
attributes ['dh_hist' ]['units' ] = 'meters'
398
399
attributes ['dh_hist' ]['coordinates' ] = 'y x'
399
400
attributes ['dh_hist' ]['grid_mapping' ] = 'crs'
@@ -613,7 +614,7 @@ def main():
613
614
loglevels = [logging .CRITICAL , logging .INFO , logging .DEBUG ]
614
615
logging .basicConfig (level = loglevels [args .verbose ])
615
616
616
- # run program
617
+ # run program
617
618
try :
618
619
info (args )
619
620
tidal_histogram (args .infile ,
0 commit comments