Commit 608fd9ca by Fabian Wachsmann

Enhanced description for percentile based ceis

parent 2bde0ded
cdo-cei.ipynb 100644 → 100755
 ... ... @@ -244,22 +244,22 @@ ## 5. Percentile based indices %% Cell type:markdown id: tags: The percentile based indices require a special calculation method for years that lie inside the base period. Then, bootstrapping must be applied where the base period is modified: The index year is taken from the base period and is replaced by another year. Then, the percentile as well as the index is calculated for the new 30 year base period. This is done 29 times so that each year from the base period will be accounted twice. In the end, the mean of 29 indices is taken.
The temperature related percentile based indices (*tx90p, tx10p, tn90p, tn10p*) require a special percentile calculation method for years that lie inside the base period. For that years, bootstrapping must be applied where the base period is modified: The index year is taken from the base period and is replaced by another year. Then, the percentile as well as the index is calculated for the new 30 year base period. This is done 29 times so that each year from the base period will be accounted twice. In the end, the mean of 29 indices is taken.
Therefore, the operators need input arguments: 1. The window days 2. The start year of the bootstrapping interval 3. The end year of the bootstrapping interval 4. The output frequency %% Cell type:code id: tags: ``` python !export CDO_PCTL_NBINS=302 tn10p_values = cdo.etccdi_tn10p(5,1961,1990,"freq=year", tn10p_values = cdo.etccdi_tn10p(5,1961,1990, input=tasminHamburg+" tasmin_runmin.nc tasmin_runmax.nc", output="tn10p_hamburg.nc", returnCdf=True).variables["tn10pETCCDI"][:] tn10p_values = tn10p_values.flatten() ``` ... ... @@ -272,11 +272,5 @@ plt.grid() plt.xlabel("Year") plt.ylabel("Number of days with tmin < tmin90") plt.show() ``` %% Cell type:code id: tags: ``` python ``` ... ...
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!