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

snapshot...

parent ac333735
Branches
Tags
No related merge requests found
...@@ -54,10 +54,15 @@ def infer_cloud_products(clvrx_path, output_dir, full_disk=True, satellite='GOES ...@@ -54,10 +54,15 @@ def infer_cloud_products(clvrx_path, output_dir, full_disk=True, satellite='GOES
cld_opd_1[clr_cat] = opd[clr_cat] cld_opd_1[clr_cat] = opd[clr_cat]
cld_opd_1[cld_cat] = opd[cld_cat] cld_opd_1[cld_cat] = opd[cld_cat]
cld_opd_2 = cld_opd_1.copy()
cld_opd_2 = np.where(cld_frac == 1, cld_opd_2 / 0.25, cld_opd_2)
cld_opd_2 = np.where(cld_frac == 2, cld_opd_2 / 0.50, cld_opd_2)
cld_opd_2 = np.where(cld_frac == 3, cld_opd_2 / 0.75, cld_opd_2)
cld_opd_1[np.invert(np.logical_and(cld_opd_1 < 0.0, cld_opd_1 > 160.0))] = -1.0 cld_opd_1[np.invert(np.logical_and(cld_opd_1 < 0.0, cld_opd_1 > 160.0))] = -1.0
cld_opd_1[np.isnan(opd)] = -1.0 cld_opd_1[np.isnan(opd)] = -1.0
write_cld_prods_file_nc4(clvrx_str_time, out_file, cld_frac, cld_opd_1, x_rad, y_rad, None, None, write_cld_prods_file_nc4(clvrx_str_time, out_file, cld_frac, cld_opd_1, cld_opd_2, x_rad, y_rad, None, None,
satellite=satellite, domain=domain, has_time=True) satellite=satellite, domain=domain, has_time=True)
t1 = time.time() t1 = time.time()
print('total time: ', (t1 - t0)) print('total time: ', (t1 - t0))
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment