diff --git a/book/_toc.yml b/book/_toc.yml index b62df30..25157f5 100644 --- a/book/_toc.yml +++ b/book/_toc.yml @@ -4,12 +4,10 @@ root: intro.md parts: - caption: Contents chapters: - - file: oneModel/oneModel.md - sections: - - file: oneModel/example_model_run_HBV.ipynb + - file: oneModel/example_model_run_HBV.ipynb - file: oneModel/additionalExperiments.md sections: - - file: oneModel/example_model_run_HBV_CMIP_forcing.ipynb - - file: oneModel/calibration_HBV/example_calibrate_HBV.ipynb - - file: oneModel/pcrglobwb/example_model_run_pcrglobwb.ipynb + - file: oneModel/additional/example_model_run_HBV_CMIP_forcing.ipynb + - file: oneModel/additional/calibration_HBV/example_calibrate_HBV.ipynb + - file: oneModel/additional/pcrglobwb/example_model_run_pcrglobwb.ipynb - file: references.md \ No newline at end of file diff --git a/book/oneModel/calibration_HBV/example_calibrate_HBV.ipynb b/book/oneModel/additional/calibration_HBV/example_calibrate_HBV.ipynb similarity index 100% rename from book/oneModel/calibration_HBV/example_calibrate_HBV.ipynb rename to book/oneModel/additional/calibration_HBV/example_calibrate_HBV.ipynb diff --git a/book/oneModel/example_model_run_HBV_CMIP_forcing.ipynb b/book/oneModel/additional/example_model_run_HBV_CMIP_forcing.ipynb similarity index 99% rename from book/oneModel/example_model_run_HBV_CMIP_forcing.ipynb rename to book/oneModel/additional/example_model_run_HBV_CMIP_forcing.ipynb index 173393b..d7bdaf2 100644 --- a/book/oneModel/example_model_run_HBV_CMIP_forcing.ipynb +++ b/book/oneModel/additional/example_model_run_HBV_CMIP_forcing.ipynb @@ -5,7 +5,7 @@ "id": "e9cd8d85-eb4a-45d6-925a-e3e770592043", "metadata": {}, "source": [ - "# Leven at Newby Bridge using HBV model in eWaterCycle\n", + "# Climate Change impact for River Leven at Newby Bridge using HBV model in eWaterCycle\n", "This notebook demonstrates how eWaterCycle is used to run a model. For the EuroCSDMS workshop we, of course, choose the beautiful region we are in as a case study. River discharge in river Leven is measured at the weirs at Newby Bridge, see the photo below. This observational data is available in the CamelsGB dataset (Coxon 20XX). The larger Caravan (collection of Camels...) dataset is available through eWaterCycle and will be used below.\n", "\n", "![image](https://upload.wikimedia.org/wikipedia/commons/7/76/Weirs_on_the_River_Leven_at_Newby_Bridge_-_geograph.org.uk_-_5455774.jpg)\n", diff --git a/book/oneModel/pcrglobwb/example_model_run_pcrglobwb.ipynb b/book/oneModel/additional/pcrglobwb/example_model_run_pcrglobwb.ipynb similarity index 100% rename from book/oneModel/pcrglobwb/example_model_run_pcrglobwb.ipynb rename to book/oneModel/additional/pcrglobwb/example_model_run_pcrglobwb.ipynb diff --git a/book/oneModel/pcrglobwb/generate_clonemaps.ipynb b/book/oneModel/additional/pcrglobwb/generate_clonemaps.ipynb similarity index 100% rename from book/oneModel/pcrglobwb/generate_clonemaps.ipynb rename to book/oneModel/additional/pcrglobwb/generate_clonemaps.ipynb diff --git a/book/oneModel/pcrglobwb/generate_forcing.ipynb b/book/oneModel/additional/pcrglobwb/generate_forcing.ipynb similarity index 100% rename from book/oneModel/pcrglobwb/generate_forcing.ipynb rename to book/oneModel/additional/pcrglobwb/generate_forcing.ipynb diff --git a/book/oneModel/additionalExperiments.md b/book/oneModel/additionalExperiments.md index 992a9e7..b7b8002 100644 --- a/book/oneModel/additionalExperiments.md +++ b/book/oneModel/additionalExperiments.md @@ -1,8 +1,8 @@ (intro)= # Additional Experiments -Next to the simple "run a model" tutorial we provide a number of notebooks to demonstrate more advanced uses of eWaterCycle. The purpose of all these notebooks is to demonstrate what little changes are needed to go from a +Next to the simple "run a model" tutorial we provide a number of notebooks to demonstrate more advanced uses of eWaterCycle. The purpose of all these notebooks is to demonstrate what little changes are needed to go from a simple "Hello World" of running a model, to more advanced use cases. -- Generating forcing from CMIP for climate change impact studies. In this notebook we make little changes to the tutorial notebook to not use ERA5 re-analysis data as forcing input, but rather use data from the CMIP models to do climate change impact studies. -- Calibrate a model using observational data. In this notebook we run the tutorial model a thousand times with different parameter sets to find the optimal parameter set linking model output to observations. -- Running a distributed, more complex model. In [this notebook](pcrglobwb/example_model_run_pcrglobwb.ipynb) we run the PCRGlobWB distributed model developed by Utrecht University for a lot of global hydrology studies, focussed on the UK. +- Generating forcing from CMIP for climate change impact studies. In [this notebook](additional/example_model_run_HBV_CMIP_forcing.ipynb) we make little changes to the tutorial notebook to not use ERA5 re-analysis data as forcing input, but rather use data from the CMIP models to do climate change impact studies. +- Calibrate a model using observational data. In [this notebook](additional/calibration_HBV/example_calibrate_HBV.ipynb) we run the tutorial model a thousand times with different parameter sets to find the optimal parameter set linking model output to observations. +- Running a distributed, more complex model. In [this notebook](additional/pcrglobwb/example_model_run_pcrglobwb.ipynb) we run the PCRGlobWB distributed model developed by Utrecht University for a lot of global hydrology studies, focussed on the UK. diff --git a/book/oneModel/example_model_run_HBV.ipynb b/book/oneModel/example_model_run_HBV.ipynb index bb08ba4..5199e8b 100644 --- a/book/oneModel/example_model_run_HBV.ipynb +++ b/book/oneModel/example_model_run_HBV.ipynb @@ -274,7 +274,7 @@ { "data": { "text/plain": [ - "[]" + "[]" ] }, "execution_count": 9, diff --git a/book/oneModel/oneModel.md b/book/oneModel/oneModel.md deleted file mode 100644 index 3b23454..0000000 --- a/book/oneModel/oneModel.md +++ /dev/null @@ -1,6 +0,0 @@ -(intro)= -# Running a Model in eWaterCycle - -_Write book intro here. This is the first page the student will see when opening the url._ - -Text as a test, to be added more later. -Rolf- \ No newline at end of file diff --git a/book/oneModel/pcrglobwb/pcrglobwb_uk_05min.ini b/book/oneModel/pcrglobwb/pcrglobwb_uk_05min.ini deleted file mode 100644 index 66700a1..0000000 --- a/book/oneModel/pcrglobwb/pcrglobwb_uk_05min.ini +++ /dev/null @@ -1,466 +0,0 @@ -[globalOptions] - -# Set the pcrglobwb output directory in an absolute path. -outputDir = /data/output -# Set the input directory map in an absolute path. -# - The input forcing and parameter directories and files will be relative to this. -inputDir = /data/input - -# Map of clone (must be provided in PCRaster maps) -# - Spatial resolution and coverage are based on this map: -cloneMap = global_05min/cloneMaps/uk_05min.map -# The area/landmask of interest: -# If None, area/landmask is limited for cells with ldd value. -landmask = global_05min/cloneMaps/uk_05min.map - - -# netcdf attributes for output files: -institution = Department of Physical Geography, Utrecht University -title = PCR-GLOBWB 2 output (not coupled to MODFLOW) -description = by Edwin H. Sutanudjaja (contact: e.h.sutanudjaja@uu.nl) - -startTime = 2002-01-01 -endTime = 2002-12-31 -# Format: YYYY-MM-DD ; The model runs on daily time step. - - -# spinning up options: -maxSpinUpsInYears = 0 -minConvForSoilSto = 0.0 -minConvForGwatSto = 0.0 -minConvForChanSto = 0.0 -minConvForTotlSto = 0.0 - - - -[meteoOptions] - -# Set the forcing temperature and precipitation files (relative to inputDir) -precipitationNC = global_05min/meteo/pcrglobwb_OBS6_ERA5_reanaly_1_day_pr_2002-2016_merrimack.nc -temperatureNC = global_05min/meteo/pcrglobwb_OBS6_ERA5_reanaly_1_day_tas_2002-2016_merrimack.nc - -# Method to calculate referencePotETP (reference potential evaporation+transpiration) -# options are "Hamon" and "Input" ; If "Input", the netcdf input file must be given: -referenceETPotMethod = Hamon -refETPotFileNC = 'test' - -# variable names in the forcing files (optional) -precipitationVariableName = pr -temperatureVariableName = tas -referenceEPotVariableName = potentialEvaporation - -# conversion constants and factors to correct forcing values (optional) so that the units are in m.day-1 and degree Celcius -precipitationConstant = 0.0 -precipitationFactor = 1.0 -temperatureConstant = -273.15 -temperatureFactor = 1.0 -referenceEPotConstant = 0.0 -referenceEPotFactor = 1.0 - - -[meteoDownscalingOptions] -# This section is for a 5 arcmin run, for downscaling meteorological forcing at 30 arcmin to 5 arcmin. - -downscalePrecipitation = True -downscaleTemperature = True -downscaleReferenceETPot = True - -# downscaling (based on the digital elevation model): -# The downscaling will be performed by providing the "cellIds" (meteoDownscaleIds) of lower resolution cells. -meteoDownscaleIds = global_05min/meteo/downscaling_from_30min/uniqueIds_30min.map -highResolutionDEM = global_05min/meteo/downscaling_from_30min/gtopo05min.map - -# lapse rates: -temperLapseRateNC = global_05min/meteo/downscaling_from_30min/temperature_slope.nc -precipLapseRateNC = global_05min/meteo/downscaling_from_30min/precipitation_slope.nc - -# downscaling criteria (TODO: remove these): -temperatCorrelNC = global_05min/meteo/downscaling_from_30min/temperature_correl.nc -precipitCorrelNC = global_05min/meteo/downscaling_from_30min/precipitation_correl.nc - -# windows length (unit: arc-degree) for smoothing/averaging forcing data (not recommended): -smoothingWindowsLength = 0 - - -[landSurfaceOptions] -debugWaterBalance = True - -numberOfUpperSoilLayers = 2 - -# soil and parameters -# - they are used for all land cover types, unless they are are defined in certain land cover type options -# (e.g. different/various soil types for agriculture areas) -topographyNC = global_05min/landSurface/topography/topography_parameters_5_arcmin_october_2015.nc -soilPropertiesNC = global_05min/landSurface/soil/soilProperties5ArcMin.nc - - -includeIrrigation = True - -# netcdf time series for historical expansion of irrigation areas (unit: hectares). -# Note: The resolution of this map must be consisten with the resolution of cellArea. -historicalIrrigationArea = global_05min/waterUse/irrigation/irrigated_areas/irrigationArea05ArcMin.nc - -# a pcraster map/value defining irrigation efficiency (dimensionless) - optional -irrigationEfficiency = global_30min/waterUse/irrigation/irrigation_efficiency/efficiency.map - - -includeDomesticWaterDemand = True -includeIndustryWaterDemand = True -includeLivestockWaterDemand = True - -# domestic, industrial and livestock water demand data (unit must be in m.day-1) -domesticWaterDemandFile = global_05min/waterUse/waterDemand/domestic/domestic_water_demand_version_april_2015.nc -industryWaterDemandFile = global_05min/waterUse/waterDemand/industry/industry_water_demand_version_april_2015.nc -livestockWaterDemandFile = global_05min/waterUse/waterDemand/livestock/livestock_water_demand_version_april_2015.nc - - -# desalination water supply (maximum/potential/capacity) -desalinationWater = global_05min/waterUse/desalination/desalination_water_version_april_2015.nc - - -# zone IDs (scale) at which allocations of groundwater and surface water (as well as desalinated water) are performed -allocationSegmentsForGroundSurfaceWater = global_05min/waterUse/abstraction_zones/abstraction_zones_60min_05min.map - -# pcraster maps defining the partitioning of groundwater - surface water source -# -# - predefined surface water - groundwater partitioning for irrigation demand (e.g. based on Siebert, Global Map of Irrigation Areas version 5) -irrigationSurfaceWaterAbstractionFractionData = global_05min/waterUse/source_partitioning/surface_water_fraction_for_irrigation/AEI_SWFRAC.map -# -- quality map -irrigationSurfaceWaterAbstractionFractionDataQuality = global_05min/waterUse/source_partitioning/surface_water_fraction_for_irrigation/AEI_QUAL.map -# -# - threshold values defining the preference for surface water source for irrigation purpose -# -- treshold to maximize surface water irrigation use (cells with irrSurfaceWaterAbstractionFraction above this will prioritize irrigation surface water use) -treshold_to_maximize_irrigation_surface_water = 0.50 -# -- treshold to minimize fossil water withdrawal for irrigation (cells with irrSurfaceWaterAbstractionFraction below this have no fossil withdrawal for irrigation) -treshold_to_minimize_fossil_groundwater_irrigation = 0.70 -# -# - predefined surface water - groundwater partitioning for non irrigation demand (e.g. based on McDonald, 2014) -maximumNonIrrigationSurfaceWaterAbstractionFractionData = global_30min/waterUse/source_partitioning/surface_water_fraction_for_non_irrigation/max_city_sw_fraction.map - - -[forestOptions] - -name = forest -debugWaterBalance = True - -# snow module properties -snowModuleType = Simple -freezingT = 0.0 -degreeDayFactor = 0.0025 -snowWaterHoldingCap = 0.1 -refreezingCoeff = 0.05 - -# other paramater values -minTopWaterLayer = 0.0 -minCropKC = 0.2 - -cropCoefficientNC = global_05min/landSurface/landCover/naturalTall/cropCoefficientForest.nc -interceptCapNC = global_05min/landSurface/landCover/naturalTall/interceptCapInputForest.nc -coverFractionNC = global_05min/landSurface/landCover/naturalTall/coverFractionInputForest.nc - -landCoverMapsNC = None -# If NC file is not provided, we have to provide the following pcraster maps: -fracVegCover = global_05min/landSurface/landCover/naturalTall/vegf_tall.map -minSoilDepthFrac = global_30min/landSurface/landCover/naturalTall/minf_tall_permafrost.map -maxSoilDepthFrac = global_30min/landSurface/landCover/naturalTall/maxf_tall.map -rootFraction1 = global_05min/landSurface/landCover/naturalTall/rfrac1_tall.map -rootFraction2 = global_05min/landSurface/landCover/naturalTall/rfrac2_tall.map -maxRootDepth = 1.0 -# Note: The maxRootDepth is not used for non irrigated land cover type. -# -# Parameters for the Arno's scheme: -arnoBeta = None -# If arnoBeta is defined, the soil water capacity distribution is based on this. -# If arnoBeta is NOT defined, maxSoilDepthFrac must be defined such that arnoBeta will be calculated based on maxSoilDepthFrac and minSoilDepthFrac. - -# initial conditions: -interceptStorIni = global_05min/initialConditions/interceptStor_forest_1999-12-31.map -snowCoverSWEIni = global_05min/initialConditions/snowCoverSWE_forest_1999-12-31.map -snowFreeWaterIni = global_05min/initialConditions/snowFreeWater_forest_1999-12-31.map -topWaterLayerIni = global_05min/initialConditions/topWaterLayer_forest_1999-12-31.map -storUppIni = global_05min/initialConditions/storUpp_forest_1999-12-31.map -storLowIni = global_05min/initialConditions/storLow_forest_1999-12-31.map -interflowIni = global_05min/initialConditions/interflow_forest_1999-12-31.map - - -[grasslandOptions] - -name = grassland -debugWaterBalance = True - -# snow module properties -snowModuleType = Simple -freezingT = 0.0 -degreeDayFactor = 0.0025 -snowWaterHoldingCap = 0.1 -refreezingCoeff = 0.05 - -# other paramater values -minTopWaterLayer = 0.0 -minCropKC = 0.2 - -cropCoefficientNC = global_05min/landSurface/landCover/naturalShort/cropCoefficientGrassland.nc -interceptCapNC = global_05min/landSurface/landCover/naturalShort/interceptCapInputGrassland.nc -coverFractionNC = global_05min/landSurface/landCover/naturalShort/coverFractionInputGrassland.nc - -landCoverMapsNC = None -# If NC file is not provided, we have to provide the following values: -fracVegCover = global_05min/landSurface/landCover/naturalShort/vegf_short.map -minSoilDepthFrac = global_30min/landSurface/landCover/naturalShort/minf_short_permafrost.map -maxSoilDepthFrac = global_30min/landSurface/landCover/naturalShort/maxf_short.map -rootFraction1 = global_05min/landSurface/landCover/naturalShort/rfrac1_short.map -rootFraction2 = global_05min/landSurface/landCover/naturalShort/rfrac2_short.map -maxRootDepth = 0.5 -# Note: The maxRootDepth is not used for non irrigated land cover type. -# -# Parameters for the Arno's scheme: -arnoBeta = None -# If arnoBeta is defined, the soil water capacity distribution is based on this. -# If arnoBeta is NOT defined, maxSoilDepthFrac must be defined such that arnoBeta will be calculated based on maxSoilDepthFrac and minSoilDepthFrac. - -# initial conditions: -interceptStorIni = global_05min/initialConditions/interceptStor_grassland_1999-12-31.map -snowCoverSWEIni = global_05min/initialConditions/snowCoverSWE_grassland_1999-12-31.map -snowFreeWaterIni = global_05min/initialConditions/snowFreeWater_grassland_1999-12-31.map -topWaterLayerIni = global_05min/initialConditions/topWaterLayer_grassland_1999-12-31.map -storUppIni = global_05min/initialConditions/storUpp_grassland_1999-12-31.map -storLowIni = global_05min/initialConditions/storLow_grassland_1999-12-31.map -interflowIni = global_05min/initialConditions/interflow_grassland_1999-12-31.map - - -[irrPaddyOptions] - -name = irrPaddy -debugWaterBalance = True - -# snow module properties -snowModuleType = Simple -freezingT = 0.0 -degreeDayFactor = 0.0025 -snowWaterHoldingCap = 0.1 -refreezingCoeff = 0.05 -# -landCoverMapsNC = None -# If NC file is not provided, we have to provide the following values: -fracVegCover = global_05min/landSurface/landCover/irrPaddy/fractionPaddy.map -minSoilDepthFrac = global_30min/landSurface/landCover/irrPaddy/minf_paddy_permafrost.map -maxSoilDepthFrac = global_30min/landSurface/landCover/irrPaddy/maxf_paddy.map -rootFraction1 = global_30min/landSurface/landCover/irrPaddy/rfrac1_paddy.map -rootFraction2 = global_30min/landSurface/landCover/irrPaddy/rfrac2_paddy.map -maxRootDepth = 0.5 -# -# Parameters for the Arno's scheme: -arnoBeta = None -# If arnoBeta is defined, the soil water capacity distribution is based on this. -# If arnoBeta is NOT defined, maxSoilDepthFrac must be defined such that arnoBeta will be calculated based on maxSoilDepthFrac and minSoilDepthFrac. -# -# other paramater values -minTopWaterLayer = 0.05 -minCropKC = 0.2 -cropDeplFactor = 0.2 -minInterceptCap = 0.0002 - -cropCoefficientNC = global_30min/landSurface/landCover/irrPaddy/Global_CropCoefficientKc-IrrPaddy_30min.nc - -# initial conditions: -interceptStorIni = global_05min/initialConditions/interceptStor_irrPaddy_1999-12-31.map -snowCoverSWEIni = global_05min/initialConditions/snowCoverSWE_irrPaddy_1999-12-31.map -snowFreeWaterIni = global_05min/initialConditions/snowFreeWater_irrPaddy_1999-12-31.map -topWaterLayerIni = global_05min/initialConditions/topWaterLayer_irrPaddy_1999-12-31.map -storUppIni = global_05min/initialConditions/storUpp_irrPaddy_1999-12-31.map -storLowIni = global_05min/initialConditions/storLow_irrPaddy_1999-12-31.map -interflowIni = global_05min/initialConditions/interflow_irrPaddy_1999-12-31.map - - -[irrNonPaddyOptions] - -name = irrNonPaddy -debugWaterBalance = True - -# snow module properties -snowModuleType = Simple -freezingT = 0.0 -degreeDayFactor = 0.0025 -snowWaterHoldingCap = 0.1 -refreezingCoeff = 0.05 -# -landCoverMapsNC = None -# If NC file is not provided, we have to provide the following values: -fracVegCover = global_05min/landSurface/landCover/irrNonPaddy/fractionNonPaddy.map -minSoilDepthFrac = global_30min/landSurface/landCover/irrNonPaddy/minf_nonpaddy_permafrost.map -maxSoilDepthFrac = global_30min/landSurface/landCover/irrNonPaddy/maxf_nonpaddy.map -rootFraction1 = global_30min/landSurface/landCover/irrNonPaddy/rfrac1_nonpaddy.map -rootFraction2 = global_30min/landSurface/landCover/irrNonPaddy/rfrac2_nonpaddy.map -maxRootDepth = 1.0 -# -# Parameters for the Arno's scheme: -arnoBeta = None -# If arnoBeta is defined, the soil water capacity distribution is based on this. -# If arnoBeta is NOT defined, maxSoilDepthFrac must be defined such that arnoBeta will be calculated based on maxSoilDepthFrac and minSoilDepthFrac. -# -# other paramater values -minTopWaterLayer = 0.0 -minCropKC = 0.2 -cropDeplFactor = 0.5 -minInterceptCap = 0.0002 - -cropCoefficientNC = global_30min/landSurface/landCover/irrNonPaddy/Global_CropCoefficientKc-IrrNonPaddy_30min.nc - -# initial conditions: -interceptStorIni = global_05min/initialConditions/interceptStor_irrNonPaddy_1999-12-31.map -snowCoverSWEIni = global_05min/initialConditions/snowCoverSWE_irrNonPaddy_1999-12-31.map -snowFreeWaterIni = global_05min/initialConditions/snowFreeWater_irrNonPaddy_1999-12-31.map -topWaterLayerIni = global_05min/initialConditions/topWaterLayer_irrNonPaddy_1999-12-31.map -storUppIni = global_05min/initialConditions/storUpp_irrNonPaddy_1999-12-31.map -storLowIni = global_05min/initialConditions/storLow_irrNonPaddy_1999-12-31.map -interflowIni = global_05min/initialConditions/interflow_irrNonPaddy_1999-12-31.map - - - - -[groundwaterOptions] - -debugWaterBalance = True - -groundwaterPropertiesNC = global_05min/groundwater/properties/groundwaterProperties5ArcMin.nc -# The file will containspecificYield (m3.m-3), kSatAquifer, recessionCoeff (day-1) -# -# - minimum value for groundwater recession coefficient (day-1) -minRecessionCoeff = 1.0e-4 - -# some options for constraining groundwater abstraction -limitFossilGroundWaterAbstraction = True -estimateOfRenewableGroundwaterCapacity = 0.0 -estimateOfTotalGroundwaterThickness = global_05min/groundwater/aquifer_thickness_estimate/thickness_05min.map -# minimum and maximum total groundwater thickness -minimumTotalGroundwaterThickness = 100. -maximumTotalGroundwaterThickness = None - -# annual pumping capacity for each region (unit: billion cubic meter per year), should be given in a netcdf file -pumpingCapacityNC = global_30min/waterUse/groundwater_pumping_capacity/regional_abstraction_limit.nc - -# initial conditions: -storGroundwaterIni = global_05min/initialConditions/storGroundwater_1999-12-31.map -storGroundwaterFossilIni = global_05min/initialConditions/storGroundwaterFossil_1999-12-31.map -# -# additional initial conditions for pumping behaviors -avgNonFossilGroundwaterAllocationLongIni = global_05min/initialConditions/avgNonFossilGroundwaterAllocationLong_1999-12-31.map -avgNonFossilGroundwaterAllocationShortIni = global_05min/initialConditions/avgNonFossilGroundwaterAllocationShort_1999-12-31.map -avgTotalGroundwaterAbstractionIni = global_05min/initialConditions/avgTotalGroundwaterAbstraction_1999-12-31.map -avgTotalGroundwaterAllocationLongIni = global_05min/initialConditions/avgTotalGroundwaterAllocationLong_1999-12-31.map -avgTotalGroundwaterAllocationShortIni = global_05min/initialConditions/avgTotalGroundwaterAllocationShort_1999-12-31.map -# -# additional initial conditions (needed only for MODFLOW run) -relativeGroundwaterHeadIni = global_05min/initialConditions/relativeGroundwaterHead_1999-12-31.map -baseflowIni = global_05min/initialConditions/baseflow_1999-12-31.map - -# zonal IDs (scale) at which zonal allocation of groundwater is performed -allocationSegmentsForGroundwater = global_05min/waterUse/abstraction_zones/abstraction_zones_30min_05min.map - -# assumption for the minimum transmissivity value (unit: m2/day) that can be extracted (via capillary rise and/or groundwater abstraction) - optional -minTransimissivityForProductiveAquifer = 50.0 -# - if None, abstraction and/or capillary rise can always occur everywhere and unlimited - - - -[routingOptions] - -debugWaterBalance = True - -# drainage direction map -lddMap = global_05min/routing/ldd_and_cell_area/lddsound_05min.map - -# cell area (unit: m2) -cellAreaMap = global_05min/routing/ldd_and_cell_area/cellsize05min.correct.map - -# routing method: #accuTravelTime or kinematicWave -routingMethod = accuTravelTime - -# manning coefficient -manningsN = 0.04 - -# Option for flood plain simulation -dynamicFloodPlain = False - -# manning coefficient for floodplain -floodplainManningsN = 0.07 - - -# channel gradient -gradient = global_05min/routing/channel_properties/channel_gradient.map - -# constant channel depth -constantChannelDepth = global_05min/routing/channel_properties/bankfull_depth.map - -# constant channel width (optional) -constantChannelWidth = global_05min/routing/channel_properties/bankfull_width.map - -# minimum channel width (optional) -minimumChannelWidth = global_05min/routing/channel_properties/bankfull_width.map - -# channel properties for flooding -bankfullCapacity = None -# - If None, it will be estimated from (bankfull) channel depth (m) and width (m) - - -# files for relative elevation (above minimum dem) -relativeElevationFiles = global_05min/routing/channel_properties/dzRel%04d.map -relativeElevationLevels = 0.0, 0.01, 0.05, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, 0.90, 1.00 - - -# composite crop factors for WaterBodies: -cropCoefficientWaterNC = global_30min/routing/kc_surface_water/cropCoefficientForOpenWater.nc -minCropWaterKC = 1.00 - - -# lake and reservoir parameters -waterBodyInputNC = global_05min/routing/surface_water_bodies/waterBodies5ArcMin.nc -onlyNaturalWaterBodies = False - - -# initial conditions: -waterBodyStorageIni = global_05min/initialConditions/waterBodyStorage_1999-12-31.map -channelStorageIni = global_05min/initialConditions/channelStorage_1999-12-31.map -readAvlChannelStorageIni = global_05min/initialConditions/readAvlChannelStorage_1999-12-31.map -avgDischargeLongIni = global_05min/initialConditions/avgDischargeLong_1999-12-31.map -avgDischargeShortIni = global_05min/initialConditions/avgDischargeShort_1999-12-31.map -m2tDischargeLongIni = global_05min/initialConditions/m2tDischargeLong_1999-12-31.map -avgBaseflowLongIni = global_05min/initialConditions/avgBaseflowLong_1999-12-31.map -riverbedExchangeIni = global_05min/initialConditions/riverbedExchange_1999-12-31.map -# -# initial condition of sub-time step discharge (needed for estimating number of time steps in kinematic wave methods) -subDischargeIni = global_05min/initialConditions/subDischarge_1999-12-31.map -# -avgLakeReservoirInflowShortIni = global_05min/initialConditions/avgLakeReservoirInflowShort_1999-12-31.map -avgLakeReservoirOutflowLongIni = global_05min/initialConditions/avgLakeReservoirOutflowLong_1999-12-31.map -# -# number of days (timesteps) that have been performed for spinning up initial conditions in the routing module (i.e. channelStorageIni, avgDischargeLongIni, avgDischargeShortIni, etc.) -timestepsToAvgDischargeIni = global_05min/initialConditions/timestepsToAvgDischarge_1999-12-31.map -# Note that: -# - maximum number of days (timesteps) to calculate long term average flow values (default: 5 years = 5 * 365 days = 1825) -# - maximum number of days (timesteps) to calculate short term average values (default: 1 month = 1 * 30 days = 30) - - - - -[reportingOptions] - -# output files that will be written in the disk in netcdf files: -# - daily resolution -#outDailyTotNC = temperature,precipitation,totalEvaporation,gwRecharge,totalRunoff,totalGroundwaterAbstraction,discharge,surfaceWaterStorage,interceptStor,snowFreeWater,snowCoverSWE,topWaterLayer,storUppTotal,storLowTotal,storGroundwater,totalWaterStorageThickness,channelStorage,surfaceWaterAbstraction,nonIrrGrossDemand,nonIrrWaterConsumption,irrPaddyWaterWithdrawal,domesticWaterWithdrawal,industryWaterWithdrawal,precipitation_at_irrigation,evaporation_from_irrigation,surfaceWaterStorage,dynamicFracWat,floodVolume,floodDepth,surfaceWaterLevel -outDailyTotNC = temperature,precipitation,totalEvaporation,gwRecharge,totalRunoff,totalGroundwaterAbstraction,discharge,surfaceWaterStorage,interceptStor,snowFreeWater,snowCoverSWE,topWaterLayer,storUppTotal,storLowTotal,storGroundwater,totalWaterStorageThickness,channelStorage,surfaceWaterAbstraction,nonIrrGrossDemand,nonIrrWaterConsumption,irrPaddyWaterWithdrawal,domesticWaterWithdrawal,industryWaterWithdrawal,precipitation_at_irrigation,evaporation_from_irrigation,surfaceWaterStorage,dynamicFracWat,surfaceWaterLevel -# - monthly resolution -outMonthTotNC = None -outMonthAvgNC = None -outMonthEndNC = None -# - annual resolution -outAnnuaTotNC = None -outAnnuaAvgNC = None -outAnnuaEndNC = None -# - monthly and annual maxima -outMonthMaxNC = None -outAnnuaMaxNC = None - -# netcdf format and zlib setup -formatNetCDF = NETCDF4 -zlib = True diff --git a/book/oneModel/pcrglobwb/uk_05min.map b/book/oneModel/pcrglobwb/uk_05min.map deleted file mode 100644 index b446f5f..0000000 Binary files a/book/oneModel/pcrglobwb/uk_05min.map and /dev/null differ