forked from computationalgeography/pycatch
-
Notifications
You must be signed in to change notification settings - Fork 0
/
readme.txt
101 lines (71 loc) · 3.74 KB
/
readme.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
User manual
=============
Copyright of this software: Derek Karssenberg & Noemi Lana-Renault Monreal
Contact: [email protected] or [email protected]
Installation
================
Create a conda environment by running pcraster_pycatch.yaml
To run the model with 1 h timestep, run main.py
To run the model with 1 week timestep, run main_weekly.py (this model includes erosion, not documented et)
Configuring the main.py model
=============================
To remove all output, run clean.sh, this does NOT remove inputs, so there is no major risk, but check what is in clean.sh first
Reports are defined in the modules, in pcrasterPythonModules
As PCRaster maps:
- switch on or off reporting of whole modules, comment/uncomment in reportComponentsDynamic (in main.py)
- select maps reported by a module, edit in the module component def report
- note, in module component there is 'full' and 'filtering' by selecting one of these when instantiating
the module in main.py, you can also make selections in the reporting
As numpy arrays (stored as .txt files and then combined in a numpy array in the postmcloop)
- switch on or off reporting of whole modules, comment/uncomment in reportComponentsDynamic (in main.py)
- select maps reported as numpy by a module, edit in the module component, reportAsNumpy
en reportAsNumpyComponentsPostmcloop
Output file name conventions
----------------------------------
# first letter:
exchange vars starting with X
precipitation files starting with P
interception files starting with V (from vegetation)
surface store files starting with S
infiltration files staring with I
evapotranspiration files starting with E
subsurface store files starting with G (from groundwater)
runoff files starting with R
shading files starting with M (from shading due to Mountains )
budgets files starting with B
soil wash files starting with W
regolith files starting with A
bedrock weathering files starting with C
base level files starting with L
creep files starting with D
# second letter:
store (s, unit m)
actual flux in (i, unit m/h, for geomorphology m/year)
potential flux in (j, unit m/h, for geomorphology m/year)
actual flux out (o, unit m/h, for geomorphology m/year)
potential flux out (p, unit m/h, for geomorphology m/year)
actual flux in (in (positive) or out(negative)) (c, of change, unit m/h, for geomorphology m/year)
another flux (x, unit m/h, for geomorphology m/year)
lateral flux (q, cubic metre per hour)
two extra letters: other values, eg Ecl, cloud factor
see respective classes for details !! E.g. actual flux out is not always
all out fluxes but sometimes only one of two.
For particle filtering (not yet tested)
-------------------------------------------
To create observed discharge, run createObservedDischarge.sh, it creates maps in the folder 'observations'. This is only required
for particle filtering
todo and changes are in changes.txt
at the bottom of changes is also the names of files used
file name conventions
~~~~~~~~~~~~~~~~~~~~~~
# reports as numpy arrays
Got.npy self.totalActualAbstractionInUpstreamAreaCubicMetrePerHour, from subsurface
Gxt.npy self.totalUpwardSeepageInUpstreamAreaCubicMetrePerHour,from subsurface
Vot.npy self.totalActualAbstractionInUpstreamAreaCubicMetrePerHour, from canopy
Rq.npy discharge m3/h
Rqs.npy discharge m3/h, averaged over 2 hours
RPic.npy self.maximumInterceptionCapacityPerLAI,
RPks.npy self.ksat,
RPmm.npy self.multiplierMaxStomatalConductance
RPrt.npy self.regolithThicknessHomogeneous,
RPsc.npy self.saturatedConductivityMetrePerDay,