From 6b27f210b6ce599f673bafdda92daddecf63b7e8 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 2 Nov 2020 15:04:12 -0600
Subject: [PATCH] snapshot...
---
modules/amv/intercompare.py | 10 ----------
1 file changed, 10 deletions(-)
diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py
index a4f0a111..0f74928c 100644
--- a/modules/amv/intercompare.py
+++ b/modules/amv/intercompare.py
@@ -772,14 +772,6 @@ def get_amv_winds_match(dist_threshold=150, qitype=amv_cqi_idx, qival=50, lat_ra
noa_idxs = np.array(noa_idxs)
nwc_idxs = np.array(nwc_idxs)
- # amvs_all = []
- # amvs_all.append(amvs_eum[eum_idxs, :])
- # amvs_all.append(amvs_brz[brz_idxs, :])
- # amvs_all.append(amvs_kma[kma_idxs, :])
- # amvs_all.append(amvs_jma[jma_idxs, :])
- # amvs_all.append(amvs_noa[noa_idxs, :])
- # amvs_all.append(amvs_nwc[nwc_idxs, :])
-
amvs_dict = {}
amvs_dict['EUM'] = amvs_eum[eum_idxs, :]
amvs_dict['BRZ'] = amvs_brz[brz_idxs, :]
@@ -940,7 +932,6 @@ def best_fit(amv_spd, amv_dir, amv_prs, amv_lat, amv_lon, fcst_spd, fcst_dir, fc
SatwindGoodConstraint = 1
flag = 0
# Check for overall good agreement
- #vdiff = np.sqrt((amv_uwind - SatwindBestFitU) ** 2 + (amv_vwind - SatwindBestFitV) ** 2)
vdiff = np.sqrt((fcst_uwind - SatwindBestFitU) ** 2 + (fcst_vwind - SatwindBestFitV) ** 2)
min_vdiff = min(vdiff)
if min_vdiff > 4.0:
@@ -951,7 +942,6 @@ def best_fit(amv_spd, amv_dir, amv_prs, amv_lat, amv_lon, fcst_spd, fcst_dir, fc
if SatwindGoodConstraint == 1:
mm = np.where(fcst_prs > (SatwindBestFitPress + 100))[0]
nn = np.where(fcst_prs < (SatwindBestFitPress - 100))[0]
- #if (np.sum(VecDiff[mm] < (vdiff + 2.0)) + np.sum(VecDiff[nn] < (vdiff + 2.0))) > 0:
if (np.sum(vdiff[mm] < (min_vdiff + 2.0)) + np.sum(vdiff[nn] < (min_vdiff + 2.0))) > 0:
SatwindGoodConstraint = 0
flag = 1
--
GitLab