diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index c45a14a077a14352628d471d28de25bc38212809..32ad6a80b3e59f15947b5a37157234b0cddc2f75 100644 --- a/modules/amv/intercompare.py +++ b/modules/amv/intercompare.py @@ -834,7 +834,7 @@ def get_amv_winds_match(dist_threshold=150, cqi=50): return amvs_all -def best_fit(amv_spd, amv_dir, amv_prs, amv_lat, amv_lon, fcst_spd, fcst_dir, fcst_prs): +def best_fit(amv_spd, amv_dir, amv_prs, amv_lat, amv_lon, fcst_spd, fcst_dir, fcst_prs, verbose=False): """Finds the background model best fit pressure associated with the AMV. The model best-fit pressure is the height (in pressure units) where the @@ -869,9 +869,6 @@ def best_fit(amv_spd, amv_dir, amv_prs, amv_lat, amv_lon, fcst_spd, fcst_dir, fc fcst_num_levels = fcst_spd.shape[0] -# verbose = True - verbose = False - SatwindBestFitPress = undef SatwindBestFitU = undef SatwindBestFitV = undef @@ -1022,8 +1019,7 @@ def best_fit(amv_spd, amv_dir, amv_prs, amv_lat, amv_lon, fcst_spd, fcst_dir, fc def best_fit_altitude(amv_spd, amv_dir, amv_alt, amv_lat, amv_lon, fcst_spd, fcst_dir, fcst_alt, amv_uwind=None, amv_vwind=None, fcst_uwind=None, fcst_vwind=None, - alt_top=25000.0, alt_bot=0.0, bf_half_width=1000.0, constraint_half_width=1000.0): - verbose = False + alt_top=25000.0, alt_bot=0.0, bf_half_width=1000.0, constraint_half_width=1000.0, verbose=False): fcst_num_levels = fcst_alt.shape[0] flag = 3