diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py index 35c16db7158a2d255844fb660c44a0332965592e..7589d96294f7f77dcde608fb0c070e2e84027b18 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)