From a5470e32d7814719ed6700c3774d0581d66f98b4 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 3 Jun 2021 13:31:03 -0500
Subject: [PATCH] add new match file creator
---
modules/aeolus/aeolus_amv.py | 22 ++++++++++++----------
1 file changed, 12 insertions(+), 10 deletions(-)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index f4aaa3cf..c50578fc 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -1455,7 +1455,8 @@ def match_aeolus_to_clavrx(aeolus_dict, clvrx_files):
cc, ll = nav.earth_to_lc(lon, lat)
if cc is None or ll is None:
continue
-
+ if cc - 4 < 0 or ll - 4 < 0 or cc + 5 >= num_elems or ll + 5 >= num_lines:
+ continue
# c_rng, l_rng = get_search_box(nav, lon, lat)
# if c_rng is None:
# continue
@@ -1545,15 +1546,18 @@ def create_file_2(match_dct, filename, clvrx_params=['cld_height_acha', 'cld_pre
var.units = u
nc4_vars.append(var)
- # Write data to file...
+ # Write data to file ---------------------
+ prf_lon[:] = alons
+ prf_lat[:] = alats
+ prf_elem[:] = elems
+ prf_line[:] = lines
+ prf_time[:] = atimes
+
for idx, key in enumerate(keys):
tup = match_dct.get(key)
prof = tup[4]
param_nd = tup[5]
- for pidx, param in enumerate(clvrx_params):
- nc4_vars[pidx][idx, :, :] = param_nd[pidx, :, :]
-
nlevs = prof.shape[0]
for k in range(nlevs):
prf_spd[idx,k] = prof[k,6]
@@ -1561,11 +1565,9 @@ def create_file_2(match_dct, filename, clvrx_params=['cld_height_acha', 'cld_pre
prf_hht[idx,k] = prof[k,3]
prf_hhb[idx,k] = prof[k,4]
- prf_lon[:] = alons
- prf_lat[:] = alats
- prf_elem[:] = elems
- prf_line[:] = lines
- prf_time[:] = atimes
+ for pidx, param in enumerate(clvrx_params):
+ nc4_vars[pidx][idx, :, :] = param_nd[pidx, :, :]
+
rootgrp.close()
--
GitLab