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

minor

parent ff4c0f9f
No related merge requests found
......@@ -502,8 +502,8 @@ exmp_file_fd = '/Users/tomrink/data/OR_ABI-L1b-RadF-M6C16_G16_s20212521800223_e2
# # UR from Taiwan
# lon, lat = 135.0, 35.0
# elem_ur, line_ur = (2499, 995)
taiwan_x0 = 1079
taiwan_y0 = 995
taiwan_i0 = 1079
taiwan_j0 = 995
taiwan_lenx = 1420
taiwan_leny = 1020
# geos.transform_point(135.0, 35.0, ccrs.PlateCarree(), False)
......@@ -522,8 +522,8 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain
if satellite == 'H08':
xlen = taiwan_lenx
ylen = taiwan_leny
i_0 = taiwan_x0
j_0 = taiwan_y0
i_0 = taiwan_i0
j_0 = taiwan_j0
grd_dct = {name: None for name in name_list}
......@@ -552,16 +552,6 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain
i_ul = i * w_x
grd_dct_n[ds_name].append(grd_dct[ds_name][j_ul:j_ul+w_y, i_ul:i_ul+w_x])
# TODO: replace this with numpy meshgrid and flatten
# ll = []
# cc = []
# for j in range(n_y-1):
# j_y = j_0 + j*w_y
# for i in range(n_x-1):
# i_x = i_0 + i*w_x
# ll.append(j_y)
# cc.append(i_x)
grd_dct = {name: None for name in name_list}
for didx, ds_name in enumerate(name_list):
grd_dct[ds_name] = np.stack(grd_dct_n[ds_name])
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment