Skip to content

Commit edd8e1f

Browse files
author
José Gómez-Dans
committed
Starting to build a package
1 parent 1fa838d commit edd8e1f

11 files changed

+1104
-0
lines changed

‎enwofost/__init__.py

Whitespace-only changes.

‎enwofost/data/CN1.008

+393
Large diffs are not rendered by default.

‎enwofost/data/CN1.009

+392
Large diffs are not rendered by default.

‎enwofost/data/Hengshui.soil

+52
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,52 @@
1+
** $Id: winter wheat soil data generated on 2013/2/28
2+
**
3+
** SOIL DATA FILE for use with WOFOST Version 6.0
4+
**
5+
** Some of data is from EC2-medium
6+
7+
SOLNAM='Hengshui.soil'
8+
9+
** physical soil characteristics
10+
11+
** soil water retention
12+
SMTAB = -1.000, 0.390, ! vol. soil moisture content
13+
1.000, 0.361, ! as function of pF [log (cm); cm3 cm-3]
14+
1.300, 0.360,
15+
1.491, 0.346,
16+
2.000, 0.309,
17+
2.400, 0.260,
18+
2.700, 0.216,
19+
3.400, 0.134,
20+
4.204, 0.099,
21+
6.000, 0.090
22+
SMW = 0.0915 ! soil moisture content at wilting point [cm3/cm3]
23+
SMFCF = 0.3 ! soil moisture content at field capacity [cm3/cm3]
24+
SM0 = 0.472 ! soil moisture content at saturation [cm3/cm3]
25+
CRAIRC = 0.075 ! critical soil air content for aeration [cm3/cm3]
26+
** hydraulic conductivity
27+
CONTAB = 0.000, 1.380, ! 10-log hydraulic conductivity
28+
1.000, 0.312, ! as function of pF [log (cm); log (cm/day)]
29+
1.300, 0.000,
30+
1.491, -0.268,
31+
1.700, -0.538,
32+
2.000, -1.114,
33+
2.400, -1.699,
34+
2.700, -2.658,
35+
3.000, -3.092,
36+
3.400, -3.620,
37+
3.700, -4.009,
38+
4.000, -4.377,
39+
4.204, -4.638
40+
41+
42+
K0 = 23.988 ! hydraulic conductivity of saturated soil [cm day-1]
43+
SOPE = 1.32 ! maximum percolation rate root zone[cm day-1]
44+
KSUB = 2.05 ! maximum percolation rate subsoil [cm day-1]
45+
RDMSOL = 80 ! maximum soil rootable depth [cm]
46+
47+
** soil workability parameters
48+
SPADS = 0.300 ! 1st topsoil seepage parameter deep seedbed
49+
SPODS = 0.035 ! 2nd topsoil seepage parameter deep seedbed
50+
SPASS = 0.400 ! 1st topsoil seepage parameter shallow seedbed
51+
SPOSS = 0.060 ! 2nd topsoil seepage parameter shallow seedbed
52+
DEFLIM = -0.200 ! required moisture deficit deep seedbed

‎enwofost/data/WWH108.CAB

+175
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,175 @@
1+
** $Id: wwh108.cab 1.0 2017/11/2 14:07:04 $
2+
** File WWH108.CAB
3+
** CROP DATA FILE for use with WOFOST Version 7.1.7, June 1992
4+
**
5+
** WHEAT, WINTER 108
6+
** Regions : Shenzhou Agro-meteo station, Hebei Province, China
7+
** start date 1 January
8+
** mean date of flowering: 3-14 May, mature 10-15 June
9+
10+
** Calculated from observations from Shenzhou Agro-meteo station.
11+
** Calibrated for use in WOFOST model at Hebei for regional yield simulation.
12+
13+
14+
** Purpose of application: Crop model uncertainty measurement.
15+
16+
** Developed in CAU and UCL.
17+
18+
CRPNAM='Winter wheat, Shenzhou, China'
19+
20+
** emergence
21+
TBASEM = 0.0 ! lower threshold temp. for emergence [cel]
22+
TEFFMX = 30.0 ! max. eff. temp. for emergence [cel]
23+
TSUMEM = 141.0
24+
25+
** phenology
26+
IDSL = 0 ! indicates whether pre-anthesis development depends
27+
! on temp. (=0), daylength (=1) , or both (=2)
28+
DLO = 14.0 ! optimum daylength for development [hr]
29+
DLC = 8.0 ! critical daylength (lower threshold) [hr]
30+
TSUM1 = 1300.3
31+
TSUM2 = 785.0 ! temperature sum from anthesis to maturity [cel d]
32+
DTSMTB = 0.00, 0.00, ! daily increase in temp. sum
33+
8.00, 8.00,
34+
30.00, 30.00, ! as function of av. temp. [cel; cel d]
35+
45.00, 0.00
36+
DVSI = 0. ! initial DVS
37+
DVSEND = 2.00 ! development stage at harvest (= 2.0 at maturity [-])
38+
39+
** initial
40+
TDWI = 210.00 ! initial total crop dry weight [kg ha-1]
41+
LAIEM = 0.1365 ! leaf area index at emergence [ha ha-1]
42+
RGRLAI = 0.00867 ! maximum relative increase in LAI [ha ha-1 d-1]
43+
44+
** green area
45+
SLATB = 0.00, 0.00202,
46+
0.20, 0.00181,
47+
0.40, 0.00155,
48+
0.70, 0.00261,
49+
1.00, 0.00351,
50+
1.60, 0.00285,
51+
2.00, 0.00174 ! specific leaf area as a function of DVS [-; ha kg-1]
52+
SPA = 0.000 ! specific pod area [ha kg-1]
53+
SSATB = 0.0, 0.0, ! specific stem area [ha kg-1]
54+
2.0, 0.0 ! as function of DVS
55+
SPAN = 30.0 ! life span of leaves growing at 35 Celsius [d]
56+
TBASE = 0.0 ! lower threshold temp. for ageing of leaves [cel]
57+
58+
** assimilation
59+
KDIFTB = 0.0, 0.61, ! extinction coefficient for diffuse visible light [-]
60+
2.0, 0.60 ! as function of DVS
61+
EFFTB = 0.0, 0.50, ! light-use effic. single leaf [kg ha-1 hr-1 j-1 m2 s]
62+
10.0, 0.50, ! as function of daily mean temp.
63+
25.0, 0.50,
64+
30.0, 0.45,
65+
40.0, 0.35
66+
AMAXTB = 0.00, 35.83,
67+
1.00, 45.83,
68+
1.30, 38.83,
69+
1.50, 35.83,
70+
2.00, 4.48
71+
TMPFTB = 0.00, 0.01, ! reduction factor of AMAX
72+
10.00, 0.60, ! as function of av. temp. [cel; -]
73+
15.00, 1.00,
74+
25.00, 1.00,
75+
35.00, 0.01
76+
TMNFTB = -5.00, 0.00, ! red. factor of gross assim. rate
77+
0.50, 1.00, ! as function of low min. temp. [cel; -]
78+
3.00, 1.00
79+
80+
** conversion of assimilates into biomass
81+
CVL = 0.700 ! efficiency of conversion into leaves [kg kg-1]
82+
CVO = 0.80 ! efficiency of conversion into storage org. [kg kg-1]
83+
CVR = 0.694 ! efficiency of conversion into roots [kg kg-1]
84+
CVS = 0.662 ! efficiency of conversion into stems [kg kg-1]
85+
86+
** maintenance respiration
87+
Q10 = 2.0 ! rel. incr. in resp. rate per 10 Cel temp. incr. [-]
88+
RML = 0.0250 ! rel. maint. resp. rate leaves [kg CH2O kg-1 d-1]
89+
RMO = 0.0090 ! rel. maint. resp. rate stor.org. [kg CH2O kg-1 d-1]
90+
RMR = 0.0100 ! rel. maint. resp. rate roots [kg CH2O kg-1 d-1]
91+
RMS = 0.0150 ! rel. maint. resp. rate stems [kg CH2O kg-1 d-1]
92+
RFSETB = 0.00, 1.00, ! red. factor for senescence
93+
2.00, 1.00 ! as function of DVS [-; -]
94+
95+
** partitioning
96+
FRTB = 0.00, 0.50, ! fraction of total dry matter to roots
97+
0.10, 0.50, ! as a function of DVS [-; kg kg-1]
98+
0.20, 0.40,
99+
0.35, 0.22,
100+
0.40, 0.17,
101+
0.50, 0.13,
102+
0.70, 0.07,
103+
0.90, 0.03,
104+
1.20, 0.00,
105+
2.00, 0.00
106+
FLTB = 0.00, 0.682,
107+
0.15, 0.560,
108+
0.25, 0.550,
109+
0.50, 0.310,
110+
0.64, 0.220,
111+
0.95, 0.00,
112+
2.00, 0.00
113+
FSTB = 0.00, 0.318,
114+
0.15, 0.440,
115+
0.25, 0.450,
116+
0.50, 0.690,
117+
0.64, 0.780,
118+
0.95, 0.30,
119+
2.00, 0.00
120+
FOTB = 0.00, 0.00,
121+
0.15, 0.00,
122+
0.25, 0.00,
123+
0.50, 0.00,
124+
0.64, 0.00,
125+
0.95, 0.70,
126+
2.00, 1.00
127+
128+
** death rates
129+
PERDL = 0.03 ! max. rel. death rate of leaves due to water stress
130+
RDRRTB = 0.00, 0.000, ! rel. death rate of stems
131+
1.50, 0.000, ! as a function of DVS [-; kg kg-1 d-1]
132+
1.5001, 0.020,
133+
2.00, 0.020
134+
RDRSTB = 0.00, 0.000, ! rel. death rate of roots
135+
1.50, 0.000, ! as a function of DVS [-; kg kg-1 d-1]
136+
1.5001, 0.020,
137+
2.00, 0.020
138+
139+
** water use
140+
CFET = 1.00 ! correction factor transpiration rate [-]
141+
DEPNR = 4.5 ! crop group number for soil water depletion [-]
142+
IAIRDU = 0 ! air ducts in roots present (=1) or not (=0)
143+
IOX = 0 ! oxygen stress disabled
144+
145+
** rooting
146+
RDI = 10. ! initial rooting depth [cm]
147+
RRI = 1.2 ! maximum daily increase in rooting depth [cm d-1]
148+
RDMCR = 125. ! maximum rooting depth [cm]
149+
150+
** nutrients
151+
** maximum and minimum concentrations of N, P, and K
152+
** in storage organs in vegetative organs [kg kg-1]
153+
NMINSO = 0.0110 ; NMINVE = 0.0030
154+
NMAXSO = 0.0310 ; NMAXVE = 0.0105
155+
PMINSO = 0.0016 ; PMINVE = 0.0004
156+
PMAXSO = 0.0060 ; PMAXVE = 0.0020
157+
KMINSO = 0.0030 ; KMINVE = 0.0070
158+
KMAXSO = 0.0080 ; KMAXVE = 0.0280
159+
YZERO = 200. ! max. amount veg. organs at zero yield [kg ha-1]
160+
NFIX = 0.00 ! fraction of N-uptake from biol. fixation [kg kg-1]
161+
162+
LAI_OBS = -99.
163+
PAR1 = -99.
164+
* PAR2 = -99.
165+
166+
IPARTB = 1
167+
NCTB1 = 'SLATB'
168+
ICTB1 = 2
169+
IICTB1 = 2
170+
* NCTB2 = 'AMAXTB'
171+
* ICTB2 = 2
172+
* IICTB2 = 2
173+
174+
175+

‎enwofost/data/alwdgg.tif

17.8 MB
Binary file not shown.

‎enwofost/data/ec3.soil

+53
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,53 @@
1+
** $Id: ec3.new 1.2 1997/09/18 17:33:54 LEM release $
2+
**
3+
** SOIL DATA FILE for use with WOFOST Version 5.0, June 1990
4+
**
5+
** EC3-medium fine
6+
7+
SOLNAM='EC3-medium fine'
8+
9+
** physical soil characteristics
10+
11+
** soil water retention
12+
SMTAB = -1.000, 0.410, ! vol. soil moisture content
13+
1.000, 0.398, ! as function of pF [log (cm); cm3 cm-3]
14+
1.300, 0.389,
15+
1.491, 0.380,
16+
2.000, 0.340,
17+
2.400, 0.287,
18+
2.700, 0.241,
19+
3.400, 0.148,
20+
4.204, 0.104,
21+
6.000, 0.090
22+
SMW = 0.104 ! soil moisture content at wilting point [cm3/cm3]
23+
SMFCF = 0.300 ! soil moisture content at field capacity [cm3/cm3]
24+
SM0 = 0.410 ! soil moisture content at saturation [cm3/cm3]
25+
CRAIRC = 0.060 ! critical soil air content for aeration [cm3/cm3]
26+
27+
** hydraulic conductivity
28+
CONTAB = 0.000, 1.408, ! 10-log hydraulic conductivity
29+
1.000, 0.167, ! as function of pF [log (cm); log (cm/day)]
30+
1.300, -0.215,
31+
1.491, -0.638,
32+
1.700, -0.854,
33+
2.000, -1.155,
34+
2.400, -1.796,
35+
2.700, -2.260,
36+
3.000, -2.745,
37+
3.400, -3.357,
38+
3.700, -3.824,
39+
4.000, -4.276,
40+
4.204, -4.678
41+
42+
43+
K0 = 25.586 ! hydraulic conductivity of saturated soil [cm day-1]
44+
SOPE = 1.47 ! maximum percolation rate root zone[cm day-1]
45+
KSUB = 1.47 ! maximum percolation rate subsoil [cm day-1]
46+
RDMSOL = 80 ! maximum soil rootable depth [cm]
47+
48+
** soil workability parameters
49+
SPADS = 0.100 ! 1st topsoil seepage parameter deep seedbed
50+
SPODS = 0.030 ! 2nd topsoil seepage parameter deep seedbed
51+
SPASS = 0.200 ! 1st topsoil seepage parameter shallow seedbed
52+
SPOSS = 0.050 ! 2nd topsoil seepage parameter shallow seedbed
53+
DEFLIM = -0.300 ! required moisture deficit deep seedbed

‎enwofost/data/par_prior.csv

+20
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
#PARAM_CODE,Variation,PARAM_XVALUE,PARAM_YVALUE,Min,Max,StdDev,Distribution
2+
AMAXTB00,YA,0,35.83,10,70,12.5,Gaussian
3+
AMAXTB100,YA,1,55.83,20,70,12.5,Gaussian
4+
AMAXTB130,YA,1.3,55.83,20,70,12.5,Gaussian
5+
AMAXTB200,YA,2,6.48,1,20,5.25,Gaussian
6+
CVL,S,-1,0.7,0.6,0.76,0.07,Gaussian
7+
CVO,S,-1,0.81,0.55,0.85,0.07,Gaussian
8+
KDIFTB100,YP,1,0.65,0.4,1,0.05,Gaussian
9+
RGRLAI,S,-1,0.00817,0.007,0.1,0.002,Gaussian
10+
SLATB00,YA,0,0.003,0.0007,0.0042,0.0003,Gaussian
11+
SLATB040,YA,0.4,0.0028,0.0007,0.0042,0.0003,Gaussian
12+
SLATB090,YA,0.9,0.0025,0.0007,0.0042,0.0003,Gaussian
13+
SLATB200,YA,2,0.0023,0.0007,0.0042,0.0003,Gaussian
14+
SPAN,S,-1,35,17,50,5,Gaussian
15+
TBASE,S,-1,0,-10,10,2,Gaussian
16+
TDWI,S,-1,250,50,400,50,Gaussian
17+
TMNFTB_min,XP,-3,0,-7,1.5,1.5,Gaussian
18+
TSUM1,S,-1,1300,1000,1500,30,Gaussian
19+
TSUM2,S,-1,785,600,900,30,Gaussian
20+
TMPFTB_max,XP,37,0,30,45,3,Gaussian

‎enwofost/data/shenzhou_wheat.amgt

+13
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
Version: 1.0
2+
AgroManagement:
3+
- 2008-10-12:
4+
CropCalendar:
5+
crop_name: 'winter-wheat'
6+
variety_name: 'Shenzhou_wheat'
7+
crop_start_date: 2008-10-12
8+
crop_start_type: sowing
9+
crop_end_date: 2009-06-30
10+
crop_end_type: harvest
11+
max_duration: 300
12+
TimedEvents: null
13+
StateEvents: null

‎enwofost.py ‎enwofost/enwofost.py

File renamed without changes.

‎requirements.txt

+6
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
pcse
2+
numpy
3+
gdal
4+
netCDF4
5+
cdsapi
6+
scipy

0 commit comments

Comments
 (0)