Skip to content

Commit

Permalink
Merge pull request #25 from UU-Hydro/bug_fix/uly_sm_vmc
Browse files Browse the repository at this point in the history 
Bug fix/uly sm vmc
  • Loading branch information
@edwinkost
edwinkost authored Oct 31, 2023
2 parents 711ef30 + b9a71a0 commit 0511485
Show file tree
Hide file tree
Showing 6 changed files with 359 additions and 27 deletions.
11 changes: 5 additions & 6 deletions README.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -8,11 +8,11 @@ Please also see the file README.txt.

Main reference/paper: Sutanudjaja, E. H., van Beek, R., Wanders, N., Wada, Y., Bosmans, J. H. C., Drost, N., van der Ent, R. J., de Graaf, I. E. M., Hoch, J. M., de Jong, K., Karssenberg, D., López López, P., Peßenteiner, S., Schmitz, O., Straatsma, M. W., Vannametee, E., Wisser, D., and Bierkens, M. F. P.: PCR-GLOBWB 2: a 5 arcmin global hydrological and water resources model, Geosci. Model Dev., 11, 2429-2453, https://doi.org/10.5194/gmd-11-2429-2018, 2018.

## Input and output files (including OPeNDAP-based access)
## Input and output files (including OPeNDAP-based access: https://opendap.4tu.nl/thredds/catalog/data2/pcrglobwb/catalog.html)

PCR-GLOBWB input and output files for the runs made in Sutanudjaja et al. (2018, https://doi.org/10.5194/gmd-11-2429-2018) are available on https://geo.data.uu.nl/research-pcrglobwb/pcr-globwb_gmd_paper_sutanudjaja_et_al_2018/. For requesting access, please send an e-mail to [email protected].
The input files for the runs made in the aformentioned paper (Sutanudjaja et al., 2018) are available on the OPeNDAP server: https://opendap.4tu.nl/thredds/catalog/data2/pcrglobwb/version_2019_11_beta/pcrglobwb2_input/catalog.html. The OPeNDAP protocol (https://www.opendap.org) allow users to access PCR-GLOBWB input files from the remote server and perform PCR-GLOBWB runs **without** the need to download the input files (with total size ~250 GB for the global extent).

The input files (and some output files) are also available on the OPeNDAP server: https://opendap.4tu.nl/thredds/catalog/data2/pcrglobwb/catalog.html. The OPeNDAP protocol (https://www.opendap.org) allow users to access PCR-GLOBWB input files from the remote server and perform PCR-GLOBWB runs **without** the need to download the input files (with total size ~250 GB for the global extent).
Some output files are also provided: https://opendap.4tu.nl/thredds/catalog/data2/pcrglobwb/version_2019_11_beta/example_output/global_05min_gmd_paper_output/catalog.html. More output files are available on https://geo.data.uu.nl/research-pcrglobwb/pcr-globwb_gmd_paper_sutanudjaja_et_al_2018/ (for requesting access, please send an e-mail to [email protected]).

## How to install

Expand Down Expand Up @@ -58,8 +58,7 @@ Please make sure that the correct conda environment in a command prompt:

`conda activate pcrglobwb_python3`

Go to to the PCR-GLOBWB *model* directory. You can start a PCR-GLOBWB run using the following command:

`python deterministic_runner.py <ini_configuration_file>`
## Exercises/cooking recipes

where <ini_configuration_file> is the configuration file of PCR-GLOBWB.
We included some exercise/cooking recipes for running PCR-GLOBWB. You can find these documents in the folder [exercise](exercise) within this repository. While these exercises were generally designed for our own computing facilities (e.g. velocity and eejit servers), they should be adaptable for use on other computing machines.
37 changes: 24 additions & 13 deletions config/setup_05min.ini
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4,27 +4,37 @@
# - Please make sure that you have access to it.
#~ outputDir = /scratch/depfg/sutan101/pcrglobwb2_output/05min/
#~ outputDir = /scratch/sutan101/pcrglobwb2_output_test_on_velocity/05min/
outputDir = /scratch/sutan101/pcrglobwb2_output_test_on_velocity/05min_rhine-meuse/
#~ outputDir = /scratch/sutan101/pcrglobwb2_output_test_on_velocity/05min_rhine-meuse/
outputDir = /scratch-shared/edwinoxy/pcrglobwb2_output_test/05min_rhine-meuse/


# Set the input directory map in an absolute path. The input forcing and parameter directories and files will be relative to this.
# - The following is an example using files from the opendap server.
inputDir = https://opendap.4tu.nl/thredds/dodsC/data2/pcrglobwb/version_2019_11_beta/pcrglobwb2_input/
#~ # - The followings are examples using input files stored locally in your computing machine.
#~ inputDir = /projects/0/dfguu/data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/
#~ inputDir = /data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/


# Please set the clone map file (cloneMap), which defines the spatial resolution and extent of your study area.
# - Please make sure that the file is stored locally in your computing machine.
# - The file must be in the pcraster format.
cloneMap = /data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/global_05min/cloneMaps/RhineMeuse05min.clone.map
#~ cloneMap = /data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/global_05min/cloneMaps/RhineMeuse05min.clone.map
#~ cloneMap = /data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/global_05min/cloneMaps/clone_global_05min.map
cloneMap = /home/edwinoxy/tmp/RhineMeuse05min.clone.map


# Set the input directory map in an absolute path. The input forcing and parameter directories and files will be relative to this.
# - The following is an example using files from the opendap server.
inputDir = https://opendap.4tu.nl/thredds/dodsC/data2/pcrglobwb/version_2019_11_beta/pcrglobwb2_input/
#~ # - The following is an example using input files stored locally in your computing machine.
#~ inputDir = /data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/

#~ (pcrglobwb_python3_2023-10-31) [email protected]:/home/edwinoxy/tmp$ ls -lah
#~ total 2.0K
#~ drwxr-x--- 2 edwinoxy edwinoxy 4.0K Oct 31 13:48 .
#~ drwxr-x--- 8 edwinoxy edwinoxy 4.0K Oct 31 13:46 ..
#~ -rwxr-xr-x 1 edwinoxy edwinoxy 11K Dec 5 2019 RhineMeuse05min.clone.map
#~ -rwxr-xr-x 1 edwinoxy edwinoxy 11K Dec 5 2019 RhineMeuse05min.landmask.map

# The area/landmask of interest:
# If None, area/landmask is limited for cells with ldd value.
#~ landmask = None
landmask = /data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/global_05min/cloneMaps/RhineMeuse05min.landmask.map
#~ landmask = /data/hydroworld/pcrglobwb2_input_release/version_2019_11_beta/pcrglobwb2_input/global_05min/cloneMaps/RhineMeuse05min.landmask.map
landmask = /home/edwinoxy/tmp/RhineMeuse05min.landmask.map


# netcdf attributes for output files:
Expand Down Expand Up @@ -446,14 +456,15 @@ timestepsToAvgDischargeIni = global_05min/initialConditions/non-natural/1999

# output files that will be written in the disk in netcdf files:
# - daily resolution
outDailyTotNC = discharge,totalRunoff,gwRecharge,totalGroundwaterAbstraction,surfaceWaterStorage
#~ outDailyTotNC = discharge,totalRunoff,gwRecharge,totalGroundwaterAbstraction,surfaceWaterStorage
outDailyTotNC = discharge,totalRunoff,gwRecharge,totalGroundwaterAbstraction,surfaceWaterStorage,ulyssesSMUpp,ulyssesSMLow,ulyssesSM
# - monthly resolution
outMonthTotNC = actualET,irrPaddyWaterWithdrawal,irrNonPaddyWaterWithdrawal,domesticWaterWithdrawal,industryWaterWithdrawal,livestockWaterWithdrawal,runoff,totalRunoff,baseflow,directRunoff,interflowTotal,totalGroundwaterAbstraction,desalinationAbstraction,surfaceWaterAbstraction,nonFossilGroundwaterAbstraction,fossilGroundwaterAbstraction,irrGrossDemand,nonIrrGrossDemand,totalGrossDemand,nonIrrWaterConsumption,nonIrrReturnFlow,precipitation,gwRecharge,surfaceWaterInf,referencePotET,totalEvaporation,totalPotentialEvaporation,totLandSurfaceActuaET,totalLandSurfacePotET,waterBodyActEvaporation,waterBodyPotEvaporation
outMonthAvgNC = discharge,temperature,dynamicFracWat,surfaceWaterStorage,interceptStor,snowFreeWater,snowCoverSWE,topWaterLayer,storUppTotal,storLowTotal,storGroundwater,storGroundwaterFossil,totalActiveStorageThickness,totalWaterStorageThickness,satDegUpp,satDegLow,channelStorage,waterBodyStorage
outMonthAvgNC = discharge,temperature,dynamicFracWat,surfaceWaterStorage,interceptStor,snowFreeWater,snowCoverSWE,topWaterLayer,storUppTotal,storLowTotal,storGroundwater,storGroundwaterFossil,totalActiveStorageThickness,totalWaterStorageThickness,satDegUpp,satDegLow,channelStorage,waterBodyStorage,ulyssesSMUpp,ulyssesSMLow,ulyssesSM
outMonthEndNC = storGroundwater,storGroundwaterFossil,waterBodyStorage,channelStorage,totalWaterStorageThickness,totalActiveStorageThickness
# - annual resolution
outAnnuaTotNC = totalEvaporation,precipitation,gwRecharge,totalRunoff,baseflow,desalinationAbstraction,surfaceWaterAbstraction,nonFossilGroundwaterAbstraction,fossilGroundwaterAbstraction,totalGroundwaterAbstraction,totalAbstraction,irrGrossDemand,nonIrrGrossDemand,totalGrossDemand,nonIrrWaterConsumption,nonIrrReturnFlow,runoff,actualET,irrPaddyWaterWithdrawal,irrNonPaddyWaterWithdrawal,irrigationWaterWithdrawal,domesticWaterWithdrawal,industryWaterWithdrawal,livestockWaterWithdrawal,precipitation_at_irrigation,netLqWaterToSoil_at_irrigation,evaporation_from_irrigation,transpiration_from_irrigation,referencePotET
outAnnuaAvgNC = temperature,discharge,surfaceWaterStorage,waterBodyStorage,interceptStor,snowFreeWater,snowCoverSWE,topWaterLayer,storUppTotal,storLowTotal,storGroundwater,storGroundwaterFossil,totalWaterStorageThickness,satDegUpp,satDegLow,channelStorage,waterBodyStorage,fractionWaterBodyEvaporation,fractionTotalEvaporation,fracSurfaceWaterAllocation,fracDesalinatedWaterAllocation,gwRecharge
outAnnuaAvgNC = temperature,discharge,surfaceWaterStorage,waterBodyStorage,interceptStor,snowFreeWater,snowCoverSWE,topWaterLayer,storUppTotal,storLowTotal,storGroundwater,storGroundwaterFossil,totalWaterStorageThickness,satDegUpp,satDegLow,channelStorage,waterBodyStorage,fractionWaterBodyEvaporation,fractionTotalEvaporation,fracSurfaceWaterAllocation,fracDesalinatedWaterAllocation,gwRecharge,ulyssesSMUpp,ulyssesSMLow,ulyssesSM
outAnnuaEndNC = surfaceWaterStorage,interceptStor,snowFreeWater,snowCoverSWE,topWaterLayer,storUppTotal,storLowTotal,storGroundwater,storGroundwaterFossil,totalWaterStorageThickness
# - monthly and annual maxima
outMonthMaxNC = channelStorage,dynamicFracWat,floodVolume,floodDepth,surfaceWaterLevel,discharge,totalRunoff
Expand Down
Loading

0 comments on commit 0511485

Please sign in to comment.