Skip to content
Snippets Groups Projects
Commit 3fdb2119 authored by tomrink's avatar tomrink
Browse files

snapshot...

parent f61e662c
No related branches found
No related tags found
No related merge requests found
......@@ -111,24 +111,29 @@ def create_file(filename, data_dct, ds_list, ds_types, lon_c, lat_c, time_s, fl_
lon_ds = h5f.create_dataset('longitude', data=lon_c, dtype='f4')
lon_ds.dims[0].label = 'time'
lon_ds.attrs.create('units', data='degrees_east')
lon_ds.attrs.create('long_name', data='PIREP longitude')
lat_ds = h5f.create_dataset('latitude', data=lat_c, dtype='f4')
lat_ds.dims[0].label = 'time'
lat_ds.attrs.create('units', data='degrees_north')
lat_ds.attrs.create('long_name', data='PIREP latitude')
time_ds = h5f.create_dataset('time', data=time_s)
time_ds.dims[0].label = 'time'
time_ds.attrs.create('units', data='seconds since 1970-1-1 00:00:00')
time_ds.attrs.create('long_name', data='PIREP time')
ice_alt_ds = h5f.create_dataset('icing_altitude', data=fl_alt_s, dtype='f4')
ice_alt_ds.dims[0].label = 'time'
ice_alt_ds.attrs.create('units', data='m')
ice_alt_ds.attrs.create('long_name', data='PIREP altitude')
if icing_intensity is not None:
icing_int_ds = h5f.create_dataset('icing_intensity', data=icing_intensity, dtype='i4')
icing_int_ds.attrs.create('long_name', data='0:No intensity report, 1:Trace, 2:Light, 3:Light Moderate, 4:Moderate, 5:Moderate Severe, 6:Severe')
icing_int_ds.attrs.create('long_name', data='From PIREP. 0:No intensity report, 1:Trace, 2:Light, 3:Light Moderate, 4:Moderate, 5:Moderate Severe, 6:Severe')
unq_ids_ds = h5f.create_dataset('unique_id', data=unq_ids, dtype='i4')
unq_ids_ds.attrs.create('long_name', data='ID mapping to PIREP icing dictionary: see pireps.py')
# copy relevant attributes
for ds_name in ds_list:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment