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)




示例#3
0
        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 + '_')