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

snapshot...

parent 11f16a60
Branches
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment