|
|
|
## FoO of DOD and emissions calculation / wind erosion threshold retrieving / final output netCDF
|
|
|
|
|
|
|
|
This script loads the .RData created in the previous steps and does the following to the selected season (you can select the season by changing the variable "**season**" [1=DJF/2=MAM/3=JJA/4=SON/5=Annual]) 12UTC-data :
|
|
|
|
|
|
|
|
* Computes the FoO of DOD overpassing a threshold (0.2 is the threshold explicitly chosen in the code in the variable "**dodthr**")
|
|
|
|
|
|
|
|
* Computes the FoO of emissions overpassing a threshold (0 is the threshold explicitly chosen in the code in the variable "**emthr**")
|
|
|
|
|
|
|
|
* Computes the <img src="http://latex.codecogs.com/gif.latex?u^{*}_t" border="0"/> by applying the FoO method.
|
|
|
|
|
|
|
|
* Computes the <img src="http://latex.codecogs.com/gif.latex?u^{*}_{t_{sd}}" border="0"/> by applying the FoO method.
|
|
|
|
|
|
|
|
* Calculates and plots in screen the bias, RMSE and correlation coefficient between <img src="http://latex.codecogs.com/gif.latex?u^{*}_t" border="0"/> and <img src="http://latex.codecogs.com/gif.latex?u^{*}_{t_{sd}}" border="0"/> with its corresponding values used by the model.
|
|
|
|
|
|
|
|
* Computes the relative error of <img src="http://latex.codecogs.com/gif.latex?u^{*}_t" border="0"/> retrieved through the FoO method with respect to the used by the model.
|
|
|
|
|
|
|
|
* Plots an histogram (number of grid cells .vs. relative error) separating 0-emission and emission sources ["histo_northafrica_relative_error_DJF.png"]
|
|
|
|
|
|
|
|
* Plots two histograms used to have an example of how the method works ["histo_FoO_method.png"/"histo_FoO_method_em.png"]
|
|
|
|
|
|
|
|
* Adapts the results to finally save them in netCDF format.
|
|
|
|
|
|
|
|
Variables included in the final netCDF:
|
|
|
|
|
|
|
|
| Variable Name| Variable Description | Units |
|
|
|
|
| ------ | ------ | ------ |
|
|
|
|
| FoO_DOD | Frequency of occurrence of DOD overpassing a threshold (0.2 predetermined)|[%]|
|
|
|
|
| FoO_em | Frequency of occurrence of emissions overpassing a threshold (0.0 predetermined)|[%]|
|
|
|
|
| mc | Moisture correction parameter |[]|
|
|
|
|
| mean_em | Mean emissions in the period studied |[kg/m²]|
|
|
|
|
| rc | Roughness correction parameter |[]|
|
|
|
|
| rel_error | Relative error of u*t_FoO_method compared to the u*t_model|[%]|
|
|
|
|
| u*t_FoO_method | u*t retrieved by using the FoO method |[m/s]|
|
|
|
|
| u*t_model | mean u*t used by the model in this period |[m/s]|
|
|
|
|
| u*tsd_FoO_method_a |u*t for smooth and dry surface obtained by dividing u*t_FoO_method by mc and rc|[m/s]|
|
|
|
|
| u*tsd_FoO_method_b |u*t for smooth and dry surface obtained by trandforming u* to u*sd and then applying the FoO method|[m/s]|
|
|
|
|
| u*tsd_model |mean u*t for smooth and dry surface used by the model|[m/s]| |
|
|
|
\ No newline at end of file |