diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index bd8edc3afb7322e8d8089be74f6c85729da0a835..ff7308644a7357484b6c98b7c96a235008e4d777 100644 --- a/modules/amv/intercompare.py +++ b/modules/amv/intercompare.py @@ -747,7 +747,7 @@ def direction_difference(dir_a, dir_b): return diff -def get_amv_winds_match(dist_threshold=150, cqi=50): +def get_amv_winds_match(dist_threshold=150, cqi=50, lat_range=[-60, 60]): # match all to EUM amvs_eum = get_amv_nd('/Users/rink/data/amv_intercompare/EUM321_csv.csv') amvs_brz = get_amv_nd('/Users/rink/data/amv_intercompare/BRZCPTECfin_121_csv.csv') @@ -756,12 +756,12 @@ def get_amv_winds_match(dist_threshold=150, cqi=50): amvs_noa = get_amv_nd('/Users/rink/data/amv_intercompare/NOA621_csv.csv') amvs_nwc = get_amv_nd('/Users/rink/data/amv_intercompare/NWC721_csv.csv') - amvs_eum = filter_amvs(amvs_eum, cqi=cqi, lat_range=[-60, 60]) - amvs_brz = filter_amvs(amvs_brz, cqi=cqi, lat_range=[-60, 60]) - amvs_jma = filter_amvs(amvs_jma, cqi=cqi, lat_range=[-60, 60]) - amvs_kma = filter_amvs(amvs_kma, cqi=cqi, lat_range=[-60, 60]) - amvs_nwc = filter_amvs(amvs_nwc, cqi=cqi, lat_range=[-60, 60]) - amvs_noa = filter_amvs(amvs_noa, cqi=cqi, lat_range=[-60, 60]) + amvs_eum = filter_amvs(amvs_eum, cqi=cqi, lat_range=lat_range) + amvs_brz = filter_amvs(amvs_brz, cqi=cqi, lat_range=lat_range) + amvs_jma = filter_amvs(amvs_jma, cqi=cqi, lat_range=lat_range) + amvs_kma = filter_amvs(amvs_kma, cqi=cqi, lat_range=lat_range) + amvs_nwc = filter_amvs(amvs_nwc, cqi=cqi, lat_range=lat_range) + amvs_noa = filter_amvs(amvs_noa, cqi=cqi, lat_range=lat_range) lat_idx = amv_hdr_list.index('lat') lon_idx = amv_hdr_list.index('lon')