From 5467852e814fb117b22b3a5ddaee5676910aaeec Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 7 Dec 2020 20:57:38 -0600
Subject: [PATCH] snapshot...
---
modules/aeolus/aeolus_amv.py | 23 +++++++++++++++++++++--
1 file changed, 21 insertions(+), 2 deletions(-)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 35c16db7..7589d962 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -83,6 +83,25 @@ class Framework(AMVFiles):
def __init__(self, files_path, file_time_span, band='14'):
super().__init__(files_path, file_time_span, '*WINDS_AMV_EN-'+band+'*.nc', band)
+ self.elem_name = 'Element'
+ self.line_name = 'Line'
+ self.lon_name = 'Longitude'
+ self.lat_name = 'Latitude'
+
+ self.out_params = ['Lon', 'Lat', 'Element', 'Line', 'pressure', 'wind_speed', 'wind_direction']
+ self.params = ['MedianPress', 'Wind_Speed', 'Wind_Dir']
+ self.meta_dict = {'Lon': ('degrees east', 'f4'), 'Lat': ('degrees north', 'f4'), 'Element': (None, 'i4'), 'Line': (None, 'i4'),
+ 'pressure': ('hPa', 'f4'), 'wind_speed': ('m s-1', 'f4'), 'wind_direction': ('degrees', 'f4')}
+
+ def get_parameters(self):
+ return self.params
+
+ def get_out_parameters(self):
+ return self.out_params
+
+ def get_meta_dict(self):
+ return self.meta_dict
+
def get_navigation(self):
GEOSNavigation(sub_lon=-75.0)
@@ -486,7 +505,7 @@ def analyze2(raob_to_amv_dct, raob_dct, gfs_filename=None):
mad = np.average(np.abs(diff))
bias = np.average(diff)
print('********************************************************')
- print('Number of good best fits: ', bf_p.shape[0])
+ print('Number of good best fits to RAOB: ', bf_p.shape[0])
print('press, MAD: {0:.2f}'.format(mad))
print('press, bias: {0:.2f}'.format(bias))
pd_std = np.std(diff)
@@ -556,7 +575,7 @@ def analyze2(raob_to_amv_dct, raob_dct, gfs_filename=None):
mad = np.average(np.abs(diff))
bias = np.average(diff)
print('********************************************************')
- print('Number of good best fits: ', bf_p.shape[0])
+ print('Number of good best fits to GFS: ', bf_p.shape[0])
print('press, MAD: {0:.2f}'.format(mad))
print('press, bias: {0:.2f}'.format(bias))
pd_std = np.std(diff)
--
GitLab