diff --git a/modules/util/util.py b/modules/util/util.py index e474ad8d440c424aee2fc004859eabd4e4fe2c38..f3bb1a22fa95bc9fd74f79fd5eae63caaaa72c20 100644 --- a/modules/util/util.py +++ b/modules/util/util.py @@ -517,6 +517,8 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain w_y = 16 i_0 = 0 j_0 = 0 + s_x = w_x + s_y = w_y geos, xlen, xmin, xmax, ylen, ymin, ymax = get_cartopy_crs(satellite, domain) if satellite == 'H08': @@ -539,17 +541,17 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain grd_dct_n = {name: [] for name in name_list} - n_x = int(xlen/w_x) - n_y = int(ylen/w_y) + n_x = int(xlen/s_x) + n_y = int(ylen/s_y) - ll = [j_0 + j*w_y for j in range(n_y-1)] - cc = [i_0 + i*w_x for i in range(n_x-1)] + ll = [j_0 + j*s_y for j in range(n_y-1)] + cc = [i_0 + i*s_x for i in range(n_x-1)] for ds_name in name_list: for j in range(n_y-1): - j_ul = j * w_y + j_ul = j * s_y for i in range(n_x-1): - i_ul = i * w_x + i_ul = i * s_x grd_dct_n[ds_name].append(grd_dct[ds_name][j_ul:j_ul+w_y, i_ul:i_ul+w_x]) grd_dct = {name: None for name in name_list}