Last version documentation
On commit 45dc479
:
- a loop over the
optical module
and a loop over therange resolution
have been added : the overlap model changes if one of those variables changes over time, the code produces an overlap model for each specific case - some unused code blocks have been removed
- some code blocks for figure production have been removed
- the figure producing the overlap model itself is not coded yet, but here is an already existing (very basic) python code :
import xarray as xr
import matplotlib.pyplot as plt
from matplotlib import figure
model = xr.open_dataset('path_to_model').to_dataframe()
model.reset_index(drop=True, level=1, inplace=True)
model = model[~model.index.duplicated(keep='first')]
plt.figure(figsize=(8,10))
plt.title('Overlap model', fontsize=20)
_ = plt.plot(model.a[0:600], model.index[0:40], linewidth=5)
plt.ylabel('Range [magl]', fontsize=20)
plt.xlabel('Coefficient a', fontsize=20)
plt.yticks(fontsize=20)
plt.xticks(fontsize=20)
plt.legend(fontsize=18)
plt.savefig('output_figure.png')
To-do list for me this week :
-
the figure PR2_*.png, produced at the end of the code, shows what a 'test day' looks like after applying the overlap model; the 'test day' needs to be specific to each optical module
andrange resolution
: it should be a clear-sky day, selected by the code itself in the L1 folder, it could be one of the L1 used to create a daily function; but this still needs to be implemented. -
a python script to decide whether artefact_mode
should be 0 or 1
Edited by Melania Van Hove