elevation_angle = 0.
observations['dr'] = dr
observations['r1_slice'] = slice(0., 15.e3 + 1., dr * di)
observations['phi_slice'] = slice(0., 2. * np.pi, 2. * np.pi / 72.)
phi, r1 = np.mgrid[observations['phi_slice'], observations['r1_slice']]
observations['azel_r1_m'] = deepcopy(np.ndarray.flatten(r1))
observations['azel_alpha_rad'] = deepcopy(np.ndarray.flatten(phi))
del r1, phi
observations['n'] = len(observations['azel_r1_m'])
observations['enu_radar_location_x'] = np.zeros(observations['n'])
observations['enu_radar_location_y'] = np.zeros(observations['n'])
observations['enu_radar_location_z'] = np.zeros(observations['n']) + 500.
observations['enu_radar_time_t'] = np.zeros(observations['n'])
observations['azel_r2_m'] = observations['azel_r1_m'] + observations['dr']
observations['azel_gamma_rad'] = np.zeros(observations['n']) + np.deg2rad(elevation_angle)
observations['beam_FWHM0_rad'] = np.zeros(observations['n']) + np.deg2rad(2.1)
observations['beam_FWHM1_rad'] = np.zeros(observations['n']) + np.deg2rad(2.1)
observations['dt'] = np.ones(observations['n'])
ao = radarfilter.calc_radar_meas(observations)
else:
ao = additional_output.ao_dct(observations['additional_output_filename'])
fun_plot_scanning.plot_scanning(observations['additional_output_filename'], 'scanning_plots/scanning_'+myway+'_')
#if not os.path.exists(additional_output_filename):
if True:
observations = {}
#general configuration files
observations['cfg_filename'] = "../../../input_files/general/standard_output.cfg;"
observations['cfg_filename'] += "../../../input_files/general/water_refractive_index_segelstein.cfg;"
observations['cfg_filename'] += "../../../input_files/general/white1999_integral.cfg;"
observations['cfg_filename'] += "../../../input_files/general/atmosphere/US1976.cfg;"
observations['cfg_filename'] += "../../../input_files/general/instruments/tara.cfg;"
observations['cfg_filename'] += "../../../input_files/retrieval/algorithm/windvectors_fdvar_horizontal_hdir_solution.cfg;"
observations['additional_output_filename'] = additional_output_filename
observations['measurements_filename'] = measurements_filename
ao = retrieval.retrieval(observations)
else:
ao = additional_output.ao_dct(additional_output_filename)
opts = {
'Doppler_velocity_ms_min': -10.,
'Doppler_velocity_ms_max': 10.,
}
fun_plot_scanning.plot_scanning(additional_output_filename, 'fdvar_plots_scanning/scanning_'+myway+'_', opts)
observations['phi_slice'] = slice(0., 2. * np.pi, 2. * np.pi / 72.)
phi, r1 = np.mgrid[observations['phi_slice'], observations['r1_slice']]
observations['azel_r1_m'] = deepcopy(np.ndarray.flatten(r1))
observations['azel_alpha_rad'] = deepcopy(np.ndarray.flatten(phi))
del r1, phi
observations['n'] = len(observations['azel_r1_m'])
observations['enu_radar_location_x'] = np.zeros(observations['n'])
observations['enu_radar_location_y'] = np.zeros(observations['n'])
observations['enu_radar_location_z'] = np.zeros(
observations['n']) + 500.
observations['enu_radar_time_t'] = np.zeros(observations['n'])
observations[
'azel_r2_m'] = observations['azel_r1_m'] + observations['dr']
observations['azel_gamma_rad'] = np.zeros(
observations['n']) + np.deg2rad(elevation_angle)
observations['beam_FWHM0_rad'] = np.zeros(
observations['n']) + np.deg2rad(2.1)
observations['beam_FWHM1_rad'] = np.zeros(
observations['n']) + np.deg2rad(2.1)
observations['dt'] = np.ones(observations['n'])
ao = radarfilter.calc_radar_meas(observations)
else:
ao = additional_output.ao_dct(
observations['additional_output_filename'])
fun_plot_scanning.plot_scanning(observations['additional_output_filename'],
'scanning_plots/scanning_' + myway + '_')
