Skip to content
Snippets Groups Projects
Commit 31ba12ad authored by tomrink's avatar tomrink
Browse files

snaphot...

parent 3b45c83f
Branches
No related tags found
No related merge requests found
...@@ -352,7 +352,7 @@ def analyze(filename): ...@@ -352,7 +352,7 @@ def analyze(filename):
return cc_ctp, acha_ctp, amv_press_avg, lats return cc_ctp, acha_ctp, amv_press_avg, lats
def compute_and_add_geo_hgt(path_to_amvs, path_to_gfs, amv_source, band='14'): def compute_and_add_geo_hgt(path_to_amvs, path_to_gfs, amv_source, band='14', out_file=None):
gfs_ds = get_datasource(path_to_gfs, 'GFS') gfs_ds = get_datasource(path_to_gfs, 'GFS')
amv_ds = get_datasource(path_to_amvs, amv_source, file_time_span=20, band=band) amv_ds = get_datasource(path_to_amvs, amv_source, file_time_span=20, band=band)
amv_filenames = amv_ds.flist amv_filenames = amv_ds.flist
...@@ -390,4 +390,8 @@ def compute_and_add_geo_hgt(path_to_amvs, path_to_gfs, amv_source, band='14'): ...@@ -390,4 +390,8 @@ def compute_and_add_geo_hgt(path_to_amvs, path_to_gfs, amv_source, band='14'):
prs_s = np.array(prs_s) prs_s = np.array(prs_s)
print('Done...') print('Done...')
if out_file is not None:
np.save(out_file, (prs_s, alt_s))
else:
return prs_s, alt_s return prs_s, alt_s
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment