edk_nc issueshttps://git.ufz.de/chs/progs/edk_nc/-/issues2023-03-10T11:34:39+01:00https://git.ufz.de/chs/progs/edk_nc/-/issues/15add range_maxdist_factor to namelist2023-03-10T11:34:39+01:00Sebastian Mülleradd range_maxdist_factor to namelistA `range_maxdist_factor` could be used to calculate maxdist from the variogram range to prevent artifacts in the output like reported in #14 .
Properties
- `range_maxdist_factor = -1.0` by default
- if `range_maxdist_factor > 0` set `ma...A `range_maxdist_factor` could be used to calculate maxdist from the variogram range to prevent artifacts in the output like reported in #14 .
Properties
- `range_maxdist_factor = -1.0` by default
- if `range_maxdist_factor > 0` set `maxDist = range * range_maxdist_factor`
- for spherical model `range_maxdist_factor = 1.0` should work perfectly
- for exponential or gaussian model one should use at least `range_maxdist_factor = 3.0`https://git.ufz.de/chs/progs/edk_nc/-/issues/14still artefacts in temperature fields2023-03-10T23:44:30+01:00Friedrich Boeingstill artefacts in temperature fieldsRounded artefacts still occur in temperature fields after solving bug in a4a19087 (which solved artefacts in precipitation fields).
I checked the station list, no doubled stations exist.
See example tavg field for 1 km and 4 km versio...Rounded artefacts still occur in temperature fields after solving bug in a4a19087 (which solved artefacts in precipitation fields).
I checked the station list, no doubled stations exist.
See example tavg field for 1 km and 4 km version:
![edk_tavg_artefacts](/uploads/5fce0ac19904b762277b45282fceb9af/edk_tavg_artefacts.png)https://git.ufz.de/chs/progs/edk_nc/-/issues/6fitting new variograms using the up-to-date dwd weather data (here for precip...2022-08-10T14:15:08+02:00Friedrich Boeingfitting new variograms using the up-to-date dwd weather data (here for precipitation)DWD weather data is downloaded from https://opendata.dwd.de/climate_environment/CDC/observations_germany/climate/daily/more_precip/historical/
A comparison to the precipitation station data stored in /data/stohyd/data/processed/Germany/...DWD weather data is downloaded from https://opendata.dwd.de/climate_environment/CDC/observations_germany/climate/daily/more_precip/historical/
A comparison to the precipitation station data stored in /data/stohyd/data/processed/Germany/DWD/pre/pre_1900-2015/ to the new station data showed that some stations changed:
```
Nr Stations| Old LUT | New LUT | Same | Abandoned | New Stations
pre |5394 |5616 |5380 |14 |236
tavg |1059 |1213 |1055 |4 |158
tmax |1059 |1213 |1055 |4 |158
tmin |1059 |1213 |1055 |4 |158
```
**However it is unclear which data was used for calculation of the variograms used in the current drought monitor version (status: 2019-11-28).**
Calculation of the variogram resulted in the following theoretical variogram parameters (using exponential variogram type)
```
nugget sill range Type Easting Northing BIAS RMSE r
0.3199E+00 0.9526E+00 0.6429E+05 2 0.0 0.0 0.0098 0.0009 0.9946
```
The Old variogram parameters were:
```
nugget sill range varTyp
0.2000E+00 0.9014E+00 0.3990E+05 2
```
In the graph the new variogram is plotted using new theoretical and experimental variogram data)
![variogram_new_pre](/uploads/ab8b0864211b9f5b95f0ef9b98674a21/variogram_new_pre.png)
Interestingly, the old nugget seems to fit the data better. The range is remarkable bigger in the new data. @thober any possible explanations?Friedrich BoeingFriedrich Boeing