data_provider.add(ice.a_priori[0])
data_provider.add(ice.a_priori[1])
data_provider.add(snow.a_priori[0])
data_provider.add(snow.a_priori[1])
data_provider.add(rain.a_priori[0])
data_provider.add(rain.a_priori[1])
data_provider.add(cloud_water_a_priori)
data_provider.add(rh_a_priori)
#data_provider.add(temperature_a_priori)
data_provider.add(ObservationError(sensors))
#
# Run the retrieval.
#
retrieval = CloudRetrieval(hydrometeors, sensors, data_provider)
retrieval.setup()
retrieval.run(658)
def plot_misfit():
ws = retrieval.simulation.workspace
y = ws.y.value
yf = ws.yf.value
y_hamp = y[:59]
yf_hamp = yf[:59]
plt.plot(y_hamp)
plt.plot(yf_hamp)
plt.figure()
y_hamp = y[59:]
bins=bins)
y = 10 * np.log10(np.maximum(10**-2.6, avgs / counts))
y[counts == 0] = -26.0
inds = np.where(cloud_sat.range_bins > 1e3)[0][0] - 1
y[:inds] = y[inds]
return y
cloud_sat.y_min = -26.0
sensors = [cloud_sat]
hydrometeors = [ice, rain]
data_provider = DardarProvider(filename, i_start, i_end)
data_provider.add(ice.a_priori[0])
data_provider.add(ice.a_priori[1])
data_provider.add(rain.a_priori[0])
data_provider.add(rain.a_priori[1])
data_provider.add(mcrf.liras.ObservationError(sensors))
#
# Run the retrieval.
#
retrieval = CloudRetrieval(hydrometeors, sensors, data_provider)
retrieval.setup(verbosity=0)
output_file = "retrieval_results.nc"
retrieval.simulation.initialize_output_file(
output_file, [("profile", i_end - i_start, i_start)],
full_retrieval_output=False)
retrieval.simulation.run_ranges(range(i_start, i_end))
#
data_provider.add(ice.a_priori[0])
data_provider.add(ice.a_priori[1])
data_provider.add(snow.a_priori[0])
data_provider.add(snow.a_priori[1])
data_provider.add(rain.a_priori[0])
data_provider.add(rain.a_priori[1])
data_provider.add(cloud_water_a_priori)
data_provider.add(rh_a_priori)
data_provider.add(mcrf.liras.ObservationError(sensors))
data_provider.add(observations)
#
# Define hydrometeors and sensors.
#
hydrometeors = [ice, rain]
sensors = [lcpr]
#
# Run the retrieval.
#
retrieval = CloudRetrieval(hydrometeors,
sensors,
data_provider,
include_cloud_water=False)
retrieval.setup()
retrieval.run(ip)