PCR_GLOBWB_2

Runs PCR-GLOBWB 2 global hydrological simulations from forcing preparation through discharge validation.

155|6|Updated Mar 31, 2026
One-click install
npx skills add https://github.com/lzwei196/KISS-Knowledge-Infrastructure-for-Scientific-Simulation --skill pcr-globwb-2-lzwei196
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: PCR_GLOBWB_2
Source: https://github.com/lzwei196/KISS-Knowledge-Infrastructure-for-Scientific-Simulation/tree/main/models/PCR_GLOBWB_2
Command: npx skills add https://github.com/lzwei196/KISS-Knowledge-Infrastructure-for-Scientific-Simulation --skill pcr-globwb-2-lzwei196

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires numpy, netCDF4, pcraster, python-dateutil, six, and includes scripts (resource) and references (resource) components.

What problem does it solve? Running the PCR-GLOBWB 2 global water balance model requires operational expertise scattered across source code, docs, and examples: preparing NetCDF forcing in exact units and variable names, building clone maps from the LDD, configuring multi-section .ini files, and diagnosing silent failures like 1000x unit errors. This Skill packages that expertise so an agent can execute, check, and recover full model runs. ## Core Features & Use Cases - End-to-end pipeline tools: Six validated scripts cover clone/landmask creation, OPeNDAP input fetching, forcing conversion (CMFD/MSWX with FAO-56 Penman-Monteith refET), soil parameter conversion, model execution, and output parsing to CSV. - Diagnostic recovery: 32 symptom-to-remedy diagnostic triplets plus a preflight check catch known failure modes such as mm-vs-m precipitation traps, Kelvin temperature forcing, and missing mapattr on PATH. - Use Case: Simulate daily discharge at a new gauged basin (e.g., Elbe at Neu Darchau): trace the catchment on the model's own LDD, fetch a 30-arcmin input subset, build forcing on the clone grid, run the model, and score discharge against GRDC observations with NSE/KGE. ## Quick Start Run python preflight_check.py in this directory, then ask the agent to simulate daily discharge for your gauge coordinates using the PCR-GLOBWB 2 pipeline and validate it against observed discharge.

Frequently Asked Questions about PCR_GLOBWB_2

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I run a PCR-GLOBWB 2 simulation for a new basin?

Run preflight_check.py first, then use make_clone_map.py with the gauge coordinates and drainage area to build clone and landmask maps, fetch inputs with fetch_pcrglobwb_inputs.py, convert forcing, write the .ini, and execute run_pcrglobwb.py. Extract discharge at the snapped gauge cell with parse_pcrglobwb_output.py.

What forcing data and units does PCR-GLOBWB 2 require?

It needs daily NetCDF files named exactly precipitation, temperature, and referencePotET or evapotranspiration. All depth quantities must be in m/day and temperature in degrees Celsius; mm/day precipitation or Kelvin temperature cause silent 1000x errors.

Why is my simulated discharge 1000 times too high?

Precipitation was almost certainly supplied in mm/day instead of the required m/day, the most common silent failure (dt_001). Divide precipitation by 1000; the same trap applies to reference ET and soil storage depths.

Why does the model read zero precipitation and temperature?

The NetCDF variable names do not match the expected names exactly, so the model silently reads zeros (dt_009). Rename forcing variables to precipitation, temperature, and evapotranspiration or referencePotET.

Can PCR-GLOBWB 2 run at sub-daily timesteps or with lateral groundwater flow?

No. The model runs on a daily timestep only, and groundwater is a single vertical linear reservoir per cell with no lateral inter-cell flow. It also excludes glacier dynamics and backwater effects in routing.

How is simulated discharge validated against observations?

Extract daily discharge at the snapped gauge cell with --aggregation dailyTot and score against observed series using NSE and PBIAS conventions (satisfactory NSE 0.5, good 0.7 per Moriasi). The default global run is uncalibrated, with median gauge correlation around 0.60.