Skip to content
Snippets Groups Projects
Commit 7bef3404 authored by Georgiy Faikin's avatar Georgiy Faikin
Browse files

The rothc model rainfall and evaporation method for weather calculations

parent 68a99435
No related branches found
No related tags found
No related merge requests found
Showing
with 6413 additions and 73 deletions
-7.00000000
0.412552953
0,0
0,0
0,0
37,0
100,0
177,0
127,0
134,0
73,0
45,0
0,0
0,0
0,0
0,0
0,0
13,0
75,0
148,0
221,0
133,0
58,0
40,0
0,0
0,0
0,0
0,0
0,0
4,0
58,0
173,0
198,0
133,0
66,0
0,0
0,0
0,0
0,0
0,0
0,0
9,0
79,0
160,0
162,0
127,0
67,0
14,0
0,0
0,0
0,0
0,0
0,0
0,0
29,0
105,0
163,0
132,0
60,0
3,0
0,0
0,0
0,0
0,0
0,0
12,0
74,0
97,0
145,0
104,0
95,0
44,0
0,0
0,0
0,0
0,0
0,0
18,0
84,0
158,0
116,0
106,0
86,0
47,0
0,0
0,0
0,0
0,0
0,0
0,0
50,0
41,0
142,0
103,0
89,0
41,0
0,0
0,0
0,0
0,0
0,0
7,0
52,0
136,0
154,0
132,0
74,0
11,0
0,0
0,0
0,0
0,0
0,0
11,0
61,0
206,0
154,0
114,0
79,0
2,0
0,0
0,0
0,0
0,0
0,0
13,0
79,0
176,0
177,0
105,0
78,0
30,0
0,0
0,0
0,0
0,0
0,0
16,0
105,0
187,0
123,0
111,0
68,0
33,0
0,0
0,0
0,0
0,0
0,0
15,0
101,0
148,0
128,0
104,0
89,0
44,0
0,0
0,0
0,0
0,0
0,0
41,0
89,0
111,0
122,0
65,0
79,0
48,0
0,0
0,0
0,0
0,0
0,0
44,0
78,0
180,0
163,0
148,0
67,0
32,0
0,0
0,0
0,0
0,0
0,0
15,0
78,0
178,0
120,0
105,0
87,0
24,0
0,0
0,0
0,0
0,0
0,0
41,0
89,0
167,0
131,0
98,0
49,0
44,0
0,0
0,0
0,0
0,0
0,0
9,0
85,0
193,0
189,0
131,0
55,0
39,0
0,0
0,0
0,0
0,0
0,0
0,0
64,0
148,0
176,0
147,0
88,0
53,0
0,0
0,0
0,0
0,0
0,0
0,0
51,0
191,0
120,0
84,0
85,0
57,0
0,0
0,0
0,0
0,0
0,0
27,0
121,0
123,0
163,0
136,0
63,0
50,0
0,0
0,0
0,0
0,0
0,0
15,0
97,0
120,0
158,0
95,0
65,0
44,0
0,0
0,0
0,0
0,0
0,0
30,0
73,0
161,0
192,0
123,0
51,0
0,0
0,0
0,0
0,0
0,0
0,0
16,0
77,0
187,0
195,0
93,0
51,0
17,0
0,0
0,0
0,0
0,0
0,0
11,0
79,0
205,0
173,0
100,0
65,0
0,0
0,0
0,0
0,0
0,0
0,0
35,0
88,0
107,0
118,0
83,0
75,0
57,0
0,0
0,0
0,0
0,0
0,0
18,0
132,0
108,0
156,0
150,0
66,0
49,0
0,0
0,0
0,0
0,0
0,0
36,0
123,0
192,0
155,0
108,0
56,0
46,0
0,0
0,0
0,0
0,0
0,0
2,0
55,0
149,0
102,0
115,0
94,0
36,0
0,0
0,0
0,0
0,0
0,0
51,0
131,0
163,0
149,0
119,0
54,0
62,0
0,0
0,0
0,0
0,0
0,0
27,0
132,0
154,0
160,0
99,0
59,0
64,0
0,0
0,0
0,0
0,0
0,0
14,0
84,0
191,0
110,0
118,0
76,0
14,0
0,0
0,0
0,0
0,0
0,0
13,0
67,0
119,0
159,0
123,0
84,0
52,0
0,0
0,0
0,0
0,0
0,0
16,0
104,0
160,0
187,0
106,0
67,0
51,0
0,0
0,0
0,0
0,0
0,0
12,0
85,0
181,0
165,0
132,0
63,0
8,0
0,0
0,0
0,0
0,0
0,0
22,0
100,0
192,0
193,0
125,0
49,0
29,0
0,0
0,0
0,0
0,0
0,0
30,0
132,0
202,0
140,0
104,0
57,0
39,0
0,0
0,0
0,0
0,0
0,0
12,0
65,0
189,0
181,0
133,0
79,0
90,0
0,0
0,0
0,0
0,0
0,0
141,0
156,0
157,0
111,0
110,0
85,0
40,0
0,0
0,0
0,0
0,0
0,0
11,0
45,0
94,0
76,0
83,0
49,0
8,0
0,0
0,0
0,0
0,0
0,0
31,0
89,0
158,0
160,0
109,0
54,0
0,0
0,0
0,0
0,0
0,0
0,0
12,0
58,0
119,0
128,0
81,0
62,0
8,0
0,0
0,0
0,0
0,0
0,0
5,0
106,0
186,0
119,0
149,0
73,0
28,0
0,0
0,0
0,0
0,0
0,0
11,0
42,0
155,0
128,0
69,0
48,0
14,0
0,0
0,0
0,0
0,0
0,0
8,0
83,0
181,0
200,0
499,0
33,0
45,0
0,0
0,0
0,0
0,0
0,0
15,0
88,0
120,0
153,0
104,0
84,0
21,0
0,0
0,0
0,0
0,0
0,0
27,0
134,0
132,0
121,0
112,0
69,0
38,0
0,0
0,0
0,0
0,0
0,0
32,0
102,0
116,0
122,0
139,0
34,0
16,0
0,0
0,0
0,0
0,0
0,0
78,0
135,0
140,0
116,0
164,0
74,0
6,0
0,0
0,0
0,0
0,0
0,0
0,0
16,0
84,0
73,0
62,0
26,0
0,0
0,0
0,0
0,0
0,0
0,0
8,0
76,0
180,0
152,0
98,0
78,0
45,0
0,0
0,0
0,0
0,0
0,0
3,0
66,0
159,0
161,0
94,0
61,0
3,0
0,0
0,0
0,0
0,0
0,0
48,0
125,0
182,0
152,0
87,0
67,0
45,0
0,0
0,0
0,0
0,0
0,0
51,0
100,0
129,0
130,0
104,0
55,0
40,0
0,0
0,0
0,0
0,0
0,0
35,0
110,0
157,0
138,0
93,0
56,0
54,0
0,0
0,0
0,0
0,0
0,0
23,0
90,0
179,0
189,0
135,0
69,0
18,0
0,0
0,0
0,0
0,0
0,0
24,0
120,0
102,0
98,0
85,0
54,0
39,0
0,0
0,0
0,0
0,0
0,0
50,0
85,0
119,0
153,0
112,0
83,0
51,0
0,0
0,0
0,0
0,0
0,0
53,0
134,0
181,0
129,0
121,0
75,0
67,0
0,0
0,0
0,0
0,0
0,0
37,0
129,0
144,0
154,0
137,0
49,0
72,0
0,0
0,0
0,0
0,0
0,0
15,0
79,0
166,0
190,0
122,0
63,0
5,0
0,0
0,0
0,0
0,0
0,0
10,0
73,0
199,0
139,0
71,0
76,0
46,0
0,0
0,0
0,0
0,0
0,0
62,0
97,0
206,0
155,0
85,0
59,0
57,0
0,0
0,0
0,0
0,0
0,0
66,0
120,0
123,0
125,0
113,0
58,0
64,0
0,0
0,0
0,0
0,0
0,0
65,0
91,0
139,0
183,0
114,0
66,0
48,0
0,0
0,0
0,0
0,0
0,0
52,0
121,0
167,0
193,0
116,0
51,0
14,0
0,0
0,0
0,0
0,0
0,0
16,0
104,0
111,0
172,0
92,0
53,0
47,0
0,0
0,0
0,0
0,0
0,0
12,0
67,0
142,0
137,0
117,0
77,0
42,0
0,0
0,0
0,0
0,0
0,0
33,0
106,0
142,0
153,0
148,0
81,0
62,0
0,0
0,0
0,0
0,0
0,0
27,0
91,0
175,0
170,0
103,0
70,0
56,0
0,0
0,0
0,0
0,0
0,0
72,0
137,0
166,0
146,0
126,0
60,0
47,0
0,0
0,0
0,0
0,0
0,0
57,0
98,0
134,0
122,0
110,0
58,0
71,0
0,0
0,0
0,0
0,0
0,0
15,0
87,0
174,0
155,0
98,0
80,0
42,0
0,0
0,0
0,0
0,0
0,0
42,0
132,0
177,0
190,0
107,0
47,0
17,0
0,0
0,0
0,0
0,0
0,0
32,0
124,0
186,0
172,0
105,0
54,0
41,0
0,0
0,0
\ No newline at end of file
6.2
26
75.7
13.6
46.7
17
140.6
48
34
55.3
38.2
56.4
26
15.6
10.5
28.6
37.4
72.3
4.6
4.8
40.6
43.1
66.4
7.2
24.2
30.4
24.3
56.7
65.2
44.6
24.8
5.8
18.2
45.5
35.3
31.8
31.3
25.7
31.9
26.6
36.2
26.4
75.9
36.7
105.3
45.9
25
20.6
22.2
19.7
7.8
51.9
90.9
30.7
26.9
68.9
67.6
48.5
3.6
25.5
1.1
9.8
2.8
34.3
29.7
146
52
52.1
20.7
68.8
7.4
21
9.5
3.1
15.2
26.8
45.2
49
130.1
45.4
42.4
13.4
17
52.4
11.9
29
17.5
14.5
82.8
53.9
75.2
51.1
54.4
20.5
2.1
24.8
9.5
4
16.2
21.4
49.8
75.9
84.2
63.7
108.8
51.4
0.7
18
4.5
12.4
37.4
17.7
78.8
16.2
58
57.2
75.7
9.5
32.4
16.2
12.4
11.8
32.3
26
21
73.2
41.8
74
43.8
21.4
81.5
19.8
14.8
8.2
3.7
10.7
41
37.8
86.5
74.8
48.1
50
28.1
10.9
11.7
9.5
42.5
14
23.5
104.8
98.4
68.8
3.5
26.3
30.4
39.8
4.3
22.2
11
15.5
76.9
110.8
66.2
137
73.2
42
63.5
16.1
15.8
6
64.5
18.8
89.1
9.4
98.6
38.9
23.9
5.2
17.1
25.6
29.7
33.7
30.4
16.5
39
37
174.9
59.8
42.4
145.5
66.5
33.8
23.5
5.7
12
25.6
70.3
80.5
152.5
136.1
93.3
64.2
15
8
14.7
3.8
30
23.8
33.3
54.3
48.8
30.7
66.9
67.7
25.4
25.3
27.5
29.2
47.9
66.2
32.6
51.5
28.9
14.8
34.9
36.5
12.1
20.6
33.5
11.5
14.4
39
98.1
89.3
78.9
88.8
25
46.6
19.4
28.2
29.3
64.7
30.4
21.9
29.7
95.7
61.2
52.2
69.9
25.4
28.8
35.2
54.8
45.7
33.1
46.5
35.5
133.3
44.4
106.3
59.8
72.4
25.3
35.4
49.9
27.7
30.8
36.9
94.5
69.5
44.1
69.8
46.3
49.5
22.3
49.2
79.2
53.7
14
6.1
71.7
41.3
47.3
107.6
69
104.1
44.6
48.6
29.5
25.1
48.7
50.1
41
31.3
71.7
92.8
27.8
5
29.5
68.3
38.1
21.4
26
27.2
104.7
85.1
103.4
130.1
45
84.1
40.3
36.5
21.3
36.9
24.8
22.8
17.5
115.2
83.5
30.6
65
26.3
26.3
12.6
14.9
16.6
11.9
14
59.8
4.6
27
47.2
31
21.2
47.8
54.1
34.4
23.6
12.5
41.4
58.7
76.8
157.5
45.5
88.7
40.4
49.2
99.4
51.9
97.2
73.7
28.2
63.9
55.2
104.9
49
81
31.4
49.1
82.4
43.1
19.7
54.2
33.6
26
42.1
40.4
166.6
21.3
39.7
51.5
66.8
50.4
37.5
28.5
43
51.9
14
130.5
71.5
51
96.5
26.7
36.9
25.1
35
10.4
56.2
63.7
103.2
48.7
24.9
43.1
47.4
93.9
30.3
98.2
73.6
27.3
117.1
10.6
61.9
37.6
26.9
53.2
114.8
10.9
32.7
37.5
15.5
45.9
48.4
30.5
25.9
58.9
26.8
81.6
128.3
66.7
35.6
2.8
5.9
31.2
32.1
60.1
43.2
18.3
30.4
79.9
53.5
42.9
20.3
11.1
52.6
17.8
70.3
39.7
15.7
99.6
89.9
49.8
79.1
58.3
45.1
16.5
23.2
27.5
18
80.2
64.2
92.3
39.5
11.3
70.4
73.4
42.8
24.1
35.4
7.5
25.1
23.5
91.7
92.7
84.1
22.5
25.1
37.7
52.1
41.5
11.8
31
61
121.4
119.3
108.6
67.5
18.5
70.4
32
81.8
22.8
46.2
16.5
40.9
90
41.5
69.5
41.1
40
27.9
175.4
96
18.7
56.2
31.7
35.8
64.5
80.5
73
89.7
77.3
106.1
57.7
54.3
28.2
16
50.4
6
11.4
20.8
138.6
17.7
13.9
123.4
38.8
50.1
31.9
18
13
30.3
82.7
90.9
88.4
154.3
51.5
11.7
39.2
40.5
62.1
17.2
30.5
0.9
32.7
59
45
131
158.2
46.7
40.5
87.5
62.1
3
5
36
31.9
68.9
103.1
70.1
22.8
27
26.6
20.1
43.4
28
19.5
81.1
14.5
85.6
78.2
10
41
75.4
80.6
29
17.1
2
20.4
9.4
77
90.9
100.7
36.1
52.1
48.3
29.4
34.5
99.5
33
11.4
20.3
52.5
85.4
87.8
15.7
81.5
58.3
70
20.9
105.9
14
3.9
64.1
9.5
103.9
87
106
56.5
57.5
24.5
16.5
31.9
5
12
18.2
87.4
80
76.5
59
48.5
0
14.5
28.8
19.5
73.3
36.5
27.9
38.3
95.7
74.1
107.3
37.9
28.5
55.8
81.3
44.5
34.9
48.9
47.1
33.2
71.2
37.4
137.9
22.3
129.8
55.4
75.5
56.8
67.3
73
25.6
48.1
60.7
114.6
94.3
116.9
84.8
85.7
17
52.3
33.6
23.3
50.7
47.7
121.6
101.2
147.1
73.5
113.7
43.9
42.4
54.1
38.2
41.1
35.5
37.7
31.3
18.1
49.6
51.3
143.6
58.9
21.7
79.8
47.3
26.1
34.5
16.1
139.8
155.9
114.4
112.1
49.6
7
86
72.6
17.8
52.4
11.3
57.3
94.4
48.3
79.3
49.3
71.2
62.4
45.5
62.7
56
13.6
47.5
26.5
68.4
56.4
58.5
39.4
20.6
78.7
36.3
18.5
30.8
17.2
16.8
48.1
89.9
72.8
25.4
124.2
29.2
54.2
27.7
49.2
29.7
46.8
21.7
37.7
107.5
8.9
54.9
59.6
153.2
60.2
40.3
43
38.4
60.7
60.1
92.3
60.8
134.4
144.8
49.9
70.9
54.6
44
66.8
53.6
34.4
27.5
33.2
9.5
68.4
85.9
47.7
33.2
39.9
54.8
37.7
100.6
45
26.8
36.5
115.8
167.5
83.4
50.9
31.5
26.4
81.7
38
90.9
47.9
29
111.6
72.1
89.2
56.6
37.5
70.4
66.1
51.2
47.5
48.8
29.5
13.5
17.8
49.7
23.8
21.3
79.6
131.7
40.5
36.8
44.6
9.3
21
46.4
40.9
80.1
66.3
144.6
100.3
52.9
52.4
42.4
88.2
50.5
30.3
37.2
58
76.4
152.5
80.1
88.4
83.5
63.2
56.9
190.194677419355
42.7
40.1
49.1
80.9
88.9
86.7
28.7
12.2
44.5
31
78.7
26.7
36.2
46.9
44.4
55.6
54.6
34.5
129.3
60.4
52.5
51.4
32
69.8
46.5
39.3
24.4
34
21.2
81.2
86.5
56.6
109.2
65.1
14.3
32.5
59
56.1
51.2
68.3
76.3
181
115.1
69.6
62.5
99.6
52.4
41.7
51.4
35.7
25.1
58.8
45.8
90.9
86.4
42.7
135.5
69.8
55.8
18.2
63.7
22.1
30.7
59.1
62.3
11.9
68
70.3
42
74.9
78.7
37.2
38.1
21.2
44.1
24.9
38.4
70.7
61.5
74.9
62.5
45.2
75.5
\ No newline at end of file
0,0
0,0
0,0
11,0
69,0
171,0
100,0
81,0
78,0
53,0
0,0
0,0
0,0
0,0
0,0
17,0
95,0
203,0
193,0
110,0
46,0
5,0
0,0
0,0
0,0
0,0
0,0
37,0
100,0
177,0
127,0
134,0
73,0
45,0
0,0
0,0
0,0
0,0
0,0
13,0
75,0
148,0
221,0
133,0
58,0
40,0
0,0
0,0
0,0
0,0
0,0
4,0
58,0
173,0
198,0
133,0
66,0
0,0
0,0
0,0
0,0
0,0
0,0
9,0
79,0
160,0
162,0
127,0
67,0
14,0
0,0
0,0
0,0
0,0
0,0
0,0
29,0
105,0
163,0
132,0
60,0
3,0
0,0
0,0
0,0
0,0
0,0
12,0
74,0
97,0
145,0
104,0
95,0
44,0
0,0
0,0
0,0
0,0
0,0
18,0
84,0
158,0
116,0
106,0
86,0
47,0
0,0
0,0
0,0
0,0
0,0
0,0
50,0
41,0
142,0
103,0
89,0
41,0
0,0
0,0
0,0
0,0
0,0
7,0
52,0
136,0
154,0
132,0
74,0
11,0
0,0
0,0
0,0
0,0
0,0
11,0
61,0
206,0
154,0
114,0
79,0
2,0
0,0
0,0
0,0
0,0
0,0
13,0
79,0
176,0
177,0
105,0
78,0
30,0
0,0
0,0
0,0
0,0
0,0
16,0
105,0
187,0
123,0
111,0
68,0
33,0
0,0
0,0
0,0
0,0
0,0
15,0
101,0
148,0
128,0
104,0
89,0
44,0
0,0
0,0
0,0
0,0
0,0
41,0
89,0
111,0
122,0
65,0
79,0
48,0
0,0
0,0
0,0
0,0
0,0
44,0
78,0
180,0
163,0
148,0
67,0
32,0
0,0
0,0
0,0
0,0
0,0
15,0
78,0
178,0
120,0
105,0
87,0
24,0
0,0
0,0
0,0
0,0
0,0
41,0
89,0
167,0
131,0
98,0
49,0
44,0
0,0
0,0
0,0
0,0
0,0
9,0
85,0
193,0
189,0
131,0
55,0
39,0
0,0
0,0
0,0
0,0
0,0
0,0
64,0
148,0
176,0
147,0
88,0
53,0
0,0
0,0
0,0
0,0
0,0
0,0
51,0
191,0
120,0
84,0
85,0
57,0
0,0
0,0
0,0
0,0
0,0
27,0
121,0
123,0
163,0
136,0
63,0
50,0
0,0
0,0
0,0
0,0
0,0
15,0
97,0
120,0
158,0
95,0
65,0
44,0
0,0
0,0
0,0
0,0
0,0
30,0
73,0
161,0
192,0
123,0
51,0
0,0
0,0
0,0
0,0
0,0
0,0
16,0
77,0
187,0
195,0
93,0
51,0
17,0
0,0
0,0
0,0
0,0
0,0
11,0
79,0
205,0
173,0
100,0
65,0
0,0
0,0
0,0
0,0
0,0
0,0
35,0
88,0
107,0
118,0
83,0
75,0
57,0
0,0
0,0
0,0
0,0
0,0
18,0
132,0
108,0
156,0
150,0
66,0
49,0
0,0
0,0
0,0
0,0
0,0
36,0
123,0
192,0
155,0
108,0
56,0
46,0
0,0
0,0
0,0
0,0
0,0
2,0
55,0
149,0
102,0
115,0
94,0
36,0
0,0
0,0
0,0
0,0
0,0
51,0
131,0
163,0
149,0
119,0
54,0
62,0
0,0
0,0
0,0
0,0
0,0
27,0
132,0
154,0
160,0
99,0
59,0
64,0
0,0
0,0
0,0
0,0
0,0
14,0
84,0
191,0
110,0
118,0
76,0
14,0
0,0
0,0
0,0
0,0
0,0
13,0
67,0
119,0
159,0
123,0
84,0
52,0
0,0
0,0
0,0
0,0
0,0
16,0
104,0
160,0
187,0
106,0
67,0
51,0
0,0
0,0
0,0
0,0
0,0
12,0
85,0
181,0
165,0
132,0
63,0
8,0
0,0
0,0
0,0
0,0
0,0
22,0
100,0
192,0
193,0
125,0
49,0
29,0
0,0
0,0
0,0
0,0
0,0
30,0
132,0
202,0
140,0
104,0
57,0
39,0
0,0
0,0
0,0
0,0
0,0
12,0
65,0
189,0
181,0
133,0
79,0
90,0
0,0
0,0
0,0
0,0
0,0
141,0
156,0
157,0
111,0
110,0
85,0
40,0
0,0
0,0
0,0
0,0
0,0
11,0
45,0
94,0
76,0
83,0
49,0
8,0
0,0
0,0
0,0
0,0
0,0
31,0
89,0
158,0
160,0
109,0
54,0
0,0
0,0
0,0
0,0
0,0
0,0
12,0
58,0
119,0
128,0
81,0
62,0
8,0
0,0
0,0
0,0
0,0
0,0
5,0
106,0
186,0
119,0
149,0
73,0
28,0
0,0
0,0
0,0
0,0
0,0
11,0
42,0
155,0
128,0
69,0
48,0
14,0
0,0
0,0
0,0
0,0
0,0
8,0
83,0
181,0
200,0
499,0
33,0
45,0
0,0
0,0
0,0
0,0
0,0
15,0
88,0
120,0
153,0
104,0
84,0
21,0
0,0
0,0
0,0
0,0
0,0
27,0
134,0
132,0
121,0
112,0
69,0
38,0
0,0
0,0
0,0
0,0
0,0
32,0
102,0
116,0
122,0
139,0
34,0
16,0
0,0
0,0
0,0
0,0
0,0
78,0
135,0
140,0
116,0
164,0
74,0
6,0
0,0
0,0
0,0
0,0
0,0
0,0
16,0
84,0
73,0
62,0
26,0
0,0
0,0
0,0
0,0
0,0
0,0
8,0
76,0
180,0
152,0
98,0
78,0
45,0
0,0
0,0
0,0
0,0
0,0
3,0
66,0
159,0
161,0
94,0
61,0
3,0
0,0
0,0
0,0
0,0
0,0
48,0
125,0
182,0
152,0
87,0
67,0
45,0
0,0
0,0
0,0
0,0
0,0
51,0
100,0
129,0
130,0
104,0
55,0
40,0
0,0
0,0
0,0
0,0
0,0
35,0
110,0
157,0
138,0
93,0
56,0
54,0
0,0
0,0
0,0
0,0
0,0
23,0
90,0
179,0
189,0
135,0
69,0
18,0
0,0
0,0
0,0
0,0
0,0
24,0
120,0
102,0
98,0
85,0
54,0
39,0
0,0
0,0
0,0
0,0
0,0
50,0
85,0
119,0
153,0
112,0
83,0
51,0
0,0
0,0
0,0
0,0
0,0
53,0
134,0
181,0
129,0
121,0
75,0
67,0
0,0
0,0
0,0
0,0
0,0
37,0
129,0
144,0
154,0
137,0
49,0
72,0
0,0
0,0
0,0
0,0
0,0
15,0
79,0
166,0
190,0
122,0
63,0
5,0
0,0
0,0
0,0
0,0
0,0
10,0
73,0
199,0
139,0
71,0
76,0
46,0
0,0
0,0
0,0
0,0
0,0
62,0
97,0
206,0
155,0
85,0
59,0
57,0
0,0
0,0
0,0
0,0
0,0
66,0
120,0
123,0
125,0
113,0
58,0
64,0
0,0
0,0
0,0
0,0
0,0
65,0
91,0
139,0
183,0
114,0
66,0
48,0
0,0
0,0
0,0
0,0
0,0
52,0
121,0
167,0
193,0
116,0
51,0
14,0
0,0
0,0
0,0
0,0
0,0
16,0
104,0
111,0
172,0
92,0
53,0
47,0
0,0
0,0
0,0
0,0
0,0
12,0
67,0
142,0
137,0
117,0
77,0
42,0
0,0
0,0
0,0
0,0
0,0
33,0
106,0
142,0
153,0
148,0
81,0
62,0
0,0
0,0
0,0
0,0
0,0
27,0
91,0
175,0
170,0
103,0
70,0
56,0
0,0
0,0
0,0
0,0
0,0
72,0
137,0
166,0
146,0
126,0
60,0
47,0
0,0
0,0
0,0
0,0
0,0
57,0
98,0
134,0
122,0
110,0
58,0
71,0
0,0
0,0
0,0
0,0
0,0
15,0
87,0
174,0
155,0
98,0
80,0
42,0
0,0
0,0
0
0
0
42
132
177
190
107
47
17
0
0
0
0
0
32
124
186
172
105
54
41
0
0
\ No newline at end of file
9.5
19.7
5.3
50.6
32
58.7
153.4
131.6
92.7
96.3
17.2
53.8
10.3
24
11
26
27.1
33.2
40.5
49.8
63.3
43.7
38.9
32.5
6.2
26
75.7
13.6
46.7
17
140.6
48
34
55.3
38.2
56.4
26
15.6
10.5
28.6
37.4
72.3
4.6
4.8
40.6
43.1
66.4
7.2
24.2
30.4
24.3
56.7
65.2
44.6
24.8
5.8
18.2
45.5
35.3
31.8
31.3
25.7
31.9
26.6
36.2
26.4
75.9
36.7
105.3
45.9
25
20.6
22.2
19.7
7.8
51.9
90.9
30.7
26.9
68.9
67.6
48.5
3.6
25.5
1.1
9.8
2.8
34.3
29.7
146
52
52.1
20.7
68.8
7.4
21
9.5
3.1
15.2
26.8
45.2
49
130.1
45.4
42.4
13.4
17
52.4
11.9
29
17.5
14.5
82.8
53.9
75.2
51.1
54.4
20.5
2.1
24.8
9.5
4
16.2
21.4
49.8
75.9
84.2
63.7
108.8
51.4
0.7
18
4.5
12.4
37.4
17.7
78.8
16.2
58
57.2
75.7
9.5
32.4
16.2
12.4
11.8
32.3
26
21
73.2
41.8
74
43.8
21.4
81.5
19.8
14.8
8.2
3.7
10.7
41
37.8
86.5
74.8
48.1
50
28.1
10.9
11.7
9.5
42.5
14
23.5
104.8
98.4
68.8
3.5
26.3
30.4
39.8
4.3
22.2
11
15.5
76.9
110.8
66.2
137
73.2
42
63.5
16.1
15.8
6
64.5
18.8
89.1
9.4
98.6
38.9
23.9
5.2
17.1
25.6
29.7
33.7
30.4
16.5
39
37
174.9
59.8
42.4
145.5
66.5
33.8
23.5
5.7
12
25.6
70.3
80.5
152.5
136.1
93.3
64.2
15
8
14.7
3.8
30
23.8
33.3
54.3
48.8
30.7
66.9
67.7
25.4
25.3
27.5
29.2
47.9
66.2
32.6
51.5
28.9
14.8
34.9
36.5
12.1
20.6
33.5
11.5
14.4
39
98.1
89.3
78.9
88.8
25
46.6
19.4
28.2
29.3
64.7
30.4
21.9
29.7
95.7
61.2
52.2
69.9
25.4
28.8
35.2
54.8
45.7
33.1
46.5
35.5
133.3
44.4
106.3
59.8
72.4
25.3
35.4
49.9
27.7
30.8
36.9
94.5
69.5
44.1
69.8
46.3
49.5
22.3
49.2
79.2
53.7
14
6.1
71.7
41.3
47.3
107.6
69
104.1
44.6
48.6
29.5
25.1
48.7
50.1
41
31.3
71.7
92.8
27.8
5
29.5
68.3
38.1
21.4
26
27.2
104.7
85.1
103.4
130.1
45
84.1
40.3
36.5
21.3
36.9
24.8
22.8
17.5
115.2
83.5
30.6
65
26.3
26.3
12.6
14.9
16.6
11.9
14
59.8
4.6
27
47.2
31
21.2
47.8
54.1
34.4
23.6
12.5
41.4
58.7
76.8
157.5
45.5
88.7
40.4
49.2
99.4
51.9
97.2
73.7
28.2
63.9
55.2
104.9
49
81
31.4
49.1
82.4
43.1
19.7
54.2
33.6
26
42.1
40.4
166.6
21.3
39.7
51.5
66.8
50.4
37.5
28.5
43
51.9
14
130.5
71.5
51
96.5
26.7
36.9
25.1
35
10.4
56.2
63.7
103.2
48.7
24.9
43.1
47.4
93.9
30.3
98.2
73.6
27.3
117.1
10.6
61.9
37.6
26.9
53.2
114.8
10.9
32.7
37.5
15.5
45.9
48.4
30.5
25.9
58.9
26.8
81.6
128.3
66.7
35.6
2.8
5.9
31.2
32.1
60.1
43.2
18.3
30.4
79.9
53.5
42.9
20.3
11.1
52.6
17.8
70.3
39.7
15.7
99.6
89.9
49.8
79.1
58.3
45.1
16.5
23.2
27.5
18
80.2
64.2
92.3
39.5
11.3
70.4
73.4
42.8
24.1
35.4
7.5
25.1
23.5
91.7
92.7
84.1
22.5
25.1
37.7
52.1
41.5
11.8
31
61
121.4
119.3
108.6
67.5
18.5
70.4
32
81.8
22.8
46.2
16.5
40.9
90
41.5
69.5
41.1
40
27.9
175.4
96
18.7
56.2
31.7
35.8
64.5
80.5
73
89.7
77.3
106.1
57.7
54.3
28.2
16
50.4
6
11.4
20.8
138.6
17.7
13.9
123.4
38.8
50.1
31.9
18
13
30.3
82.7
90.9
88.4
154.3
51.5
11.7
39.2
40.5
62.1
17.2
30.5
0.9
32.7
59
45
131
158.2
46.7
40.5
87.5
62.1
3
5
36
31.9
68.9
103.1
70.1
22.8
27
26.6
20.1
43.4
28
19.5
81.1
14.5
85.6
78.2
10
41
75.4
80.6
29
17.1
2
20.4
9.4
77
90.9
100.7
36.1
52.1
48.3
29.4
34.5
99.5
33
11.4
20.3
52.5
85.4
87.8
15.7
81.5
58.3
70
20.9
105.9
14
3.9
64.1
9.5
103.9
87
106
56.5
57.5
24.5
16.5
31.9
5
12
18.2
87.4
80
76.5
59
48.5
0
14.5
28.8
19.5
73.3
36.5
27.9
38.3
95.7
74.1
107.3
37.9
28.5
55.8
81.3
44.5
34.9
48.9
47.1
33.2
71.2
37.4
137.9
22.3
129.8
55.4
75.5
56.8
67.3
73
25.6
48.1
60.7
114.6
94.3
116.9
84.8
85.7
17
52.3
33.6
23.3
50.7
47.7
121.6
101.2
147.1
73.5
113.7
43.9
42.4
54.1
38.2
41.1
35.5
37.7
31.3
18.1
49.6
51.3
143.6
58.9
21.7
79.8
47.3
26.1
34.5
16.1
139.8
155.9
114.4
112.1
49.6
7
86
72.6
17.8
52.4
11.3
57.3
94.4
48.3
79.3
49.3
71.2
62.4
45.5
62.7
56
13.6
47.5
26.5
68.4
56.4
58.5
39.4
20.6
78.7
36.3
18.5
30.8
17.2
16.8
48.1
89.9
72.8
25.4
124.2
29.2
54.2
27.7
49.2
29.7
46.8
21.7
37.7
107.5
8.9
54.9
59.6
153.2
60.2
40.3
43
38.4
60.7
60.1
92.3
60.8
134.4
144.8
49.9
70.9
54.6
44
66.8
53.6
34.4
27.5
33.2
9.5
68.4
85.9
47.7
33.2
39.9
54.8
37.7
100.6
45
26.8
36.5
115.8
167.5
83.4
50.9
31.5
26.4
81.7
38
90.9
47.9
29
111.6
72.1
89.2
56.6
37.5
70.4
66.1
51.2
47.5
48.8
29.5
13.5
17.8
49.7
23.8
21.3
79.6
131.7
40.5
36.8
44.6
9.3
21
46.4
40.9
80.1
66.3
144.6
100.3
52.9
52.4
42.4
88.2
50.5
30.3
37.2
58
76.4
152.5
80.1
88.4
83.5
63.2
56.9
190.194677419355
42.7
40.1
49.1
80.9
88.9
86.7
28.7
12.2
44.5
31
78.7
26.7
36.2
46.9
44.4
55.6
54.6
34.5
129.3
60.4
52.5
51.4
32
69.8
46.5
39.3
24.4
34
21.2
81.2
86.5
56.6
109.2
65.1
14.3
32.5
59
56.1
51.2
68.3
76.3
181
115.1
69.6
62.5
99.6
52.4
41.7
51.4
35.7
25.1
58.8
45.8
90.9
86.4
42.7
135.5
69.8
55.8
18.2
63.7
22.1
30.7
59.1
62.3
11.9
68
70.3
42
74.9
78.7
37.2
38.1
21.2
44.1
24.9
38.4
70.7
61.5
74.9
62.5
45.2
75.5
\ No newline at end of file
File moved
File moved
-7
-2.6
2.2
7.8
15.5
20.5
22
21.4
17.7
14
4
0.7
-1.1 -1.1
-2.9 -2.9
3 3
......
...@@ -10,6 +10,18 @@ ...@@ -10,6 +10,18 @@
1 1
1 1
1 1
0
0
0
0
0
0
0
0
1
1
1
1
1 1
1 1
1 1
......
0,0
0,0
0,0
0,0
47,0
171,0
87,0
63,0
69,0
32,0
0,0
0,0
0,0
0,0
16,0
74,0
122,0
137,0
106,0
78,0
39,0
0,0
0,0
0,0
0,0
0,0
7,0
64,0
113,0
108,0
71,0
72,0
54,0
0,0
0,0
0,0
0,0
0,0
13,0
74,0
154,0
178,0
117,0
62,0
0,0
0,0
0,0
0,0
0,0
0,0
11,0
71,0
168,0
185,0
67,0
59,0
16,0
0,0
0,0
0,0
0,0
0,0
5,0
55,0
180,0
152,0
80,0
79,0
40,0
0,0
0,0
0,0
0,0
0,0
22,0
67,0
98,0
76,0
70,0
81,0
68,0
0,0
0,0
0,0
0,0
0,0
16,0
107,0
108,0
143,0
122,0
71,0
45,0
0,0
0,0
0,0
0,0
0,0
8,0
50,0
175,0
175,0
87,0
82,0
46,0
0,0
0,0
0,0
0,0
0,0
0,0
36,0
139,0
91,0
79,0
89,0
55,0
0,0
0,0
0,0
0,0
0,0
20,0
94,0
146,0
127,0
110,0
67,0
43,0
0,0
0,0
0,0
0,0
0,0
13,0
101,0
152,0
122,0
83,0
65,0
58,0
0,0
0,0
0,0
0,0
0,0
15,0
59,0
175,0
113,0
108,0
88,0
11,0
0,0
0,0
0,0
0,0
0,0
15,0
66,0
106,0
141,0
77,0
73,0
24,0
0,0
0,0
0,0
0,0
0,0
10,0
74,0
146,0
155,0
111,0
69,0
31,0
0,0
0,0
0,0
0,0
0,0
6,0
76,0
151,0
123,0
103,0
70,0
11,0
0,0
0,0
0,0
0,0
0,0
11,0
63,0
176,0
200,0
140,0
53,0
27,0
0,0
0,0
0,0
0,0
0,0
0,0
59,3
83,0
84,6
75,8
73,6
56,3
24,6
15,1
0,0
0,0
0,0
0,0
37,6
140,2
69,0
106,2
102,7
62,1
20,6
0,0
0,0
0,0
0,0
0,0
33,1
75,6
131,8
129,2
95,4
54,9
26,8
0,0
0,0
0,0
0,0
0,0
46,8
63,7
103,0
82,2
91,1
53,7
21,7
0,0
0,0
0,0
0,0
0,0
48,3
112,6
112,2
107,0
82,4
50,2
31,7
0,0
0,0
0,0
0,0
0,0
45,8
128,0
110,0
123,8
101,4
53,7
42,4
17,3
0,0
0,0
0,0
0,0
56,2
96,7
137,5
74,2
91,3
49,8
21,7
0,0
0,0
0,0
0,0
0,0
43,3
83,7
94,0
94,6
93,6
52,2
29,2
12,3
0,0
0,0
0,0
0,0
47,5
95,5
99,5
114,5
90,3
48,1
21,3
0,0
0,0
0,0
0,0
0,0
47,5
109,2
94,7
89,4
94,1
47,4
20,0
0,0
0,0
0,0
0,0
0,0
48,4
91,7
121,8
148,4
164,3
59,9
27,2
0,0
0,0
0,0
0,0
0,0
49,6
102,0
130,2
115,1
78,2
41,7
25,9
9,8
0,0
0,0
0,0
0,0
49,9
58,7
108,9
85,5
74,9
74,4
42,1
16,0
10,5
0,0
0,0
12,1
78,6
117,6
112,8
106,7
73,7
67,1
31,0
0,0
0,0
0,0
0,0
0,0
54,7
70,5
66,9
76,7
63,4
55,2
0,0
0,0
0,0
0,0
0,0
0,0
51,4
87,4
99,2
100,2
76,4
39,4
21,8
15,0
0,0
0,0
0,0
0,0
50,9
90,1
75,3
80,0
67,6
38,6
24,8
20,6
0,0
0,0
0,0
0,0
38,4
181,9
118,0
81,4
93,3
51,6
24,2
0,0
0,0
0,0
0,0
0,0
46,3
57,0
114,9
78,3
63,4
39,7
29,1
0,0
0,0
0,0
0,0
0,0
44,2
101,8
119,7
150,6
94,6
37,3
32,7
17,4
0,0
0,0
0,0
0,0
46,9
79,9
81,1
101,6
80,5
55,5
13,0
21,5
0,0
0,0
0,0
0,0
53,2
100,7
77,0
85,5
91,4
67,0
28,9
0,0
0,0
0,0
0,0
0,0
65,5
126,3
89,4
78,7
60,3
43,0
20,2
0,0
0,0
0,0
0,0
0,0
39,2
87,4
59,9
65,9
96,4
42,5
22,9
0,0
0,0
0,0
0,0
0,0
44,6
100,8
97,4
79,7
68,1
32,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
81,6
76,7
74,1
56,9
25,2
32,7
0,0
0,0
0,0
0,0
0,0
54,7
103,3
95,1
118,0
64,5
46,6
25,4
0
0
0
0
30,8
51,9
83,1
97,8
87,9
67,3
51,6
21,7
0
0
0
14,7
28,1
65,6
76,2
80,4
68,6
61,8
27,9
25
8,7
0
0
0
0
54,4
74,2
78,7
87,1
61,4
37,3
15,7
12,4
0
0
0
25,4
42,8
87,4
106,6
107,7
94,6
58,6
14,2
0
0
0
0
8
56
133
115
107
75
34
0
0
0
0
0
28
76
116
152
132
81
46
0
0
0
0
0
18
71
151
140
77
75
47
0
0
0
0
0
69
115
151
151
133
58
44
0
0
0
0
0
43
77
124
129
103
64
59
0
0
0
0
0
12
54
141
151
98
87
36
0
0
0
0
0
29
110
146
205
103
48
18
0
0
0
0
0
25
86
157
172
74
49
48
0
0
0
0
0
23,77777778
85,77777778
156,7777778
169,1111111
79,44444444
50,77777778
42,66666667
0
0
0
0
0
22,55555556
85,55555556
156,5555556
166,2222222
84,88888889
52,55555556
37,33333333
0
0
0
0
0
21,33333333
85,33333333
156,3333333
163,3333333
90,33333333
54,33333333
32
0
0
0
0
0
20,11111111
85,11111111
156,1111111
160,4444444
95,77777778
56,11111111
26,66666667
0
0
0
0
0
18,88888889
84,88888889
155,8888889
157,5555556
101,2222222
57,88888889
21,33333333
0
0
0
0
0
17,66666667
84,66666667
155,6666667
154,6666667
106,6666667
59,66666667
16
0
0
0,0
\ No newline at end of file
26.5
9.5
6.5
36.6
57.7
70.6
101.3
100
33.8
33.5
21.4
19
38.9
27.9
23.7
73.6
84.8
64
45.1
55.2
37.5
18.5
18.7
35
35.2
24.7
45
47
81
98.2
104.3
48.5
71
24.5
37.7
36.6
21.4
16.4
28.9
23.8
70.2
29.1
27.5
60.8
54.2
21.9
14.2
42.7
30.7
18.2
9
48.9
39.4
41.8
138
83.6
43.7
37.7
54
33.4
7.5
42.6
20.5
77.3
15.3
60.3
81.3
30.8
4.6
25.6
57.2
15
19.6
32.3
25.6
89.9
92.3
146.4
60.6
42.5
55.1
24.8
27.1
20.6
21.3
16.4
20.2
31.2
70.1
53.3
54.5
31.4
48.9
35.5
20.2
10.9
48.2
24.4
25.3
96.7
25.5
38
48.6
40
19.3
51.3
78.8
53.6
27.9
35.6
21.2
72.5
37.7
98.7
93.6
30.9
36.3
27.5
84.5
54.5
70.5
61.6
47.8
43.8
81
100.1
42
54.6
16.8
27.5
78.4
42.9
15.1
38.1
56.5
27.9
45.9
79.3
146.4
22.3
76.7
36.4
44.5
51.4
25.4
26
60.2
92.3
12.3
82.4
46.2
30.3
48.9
25.1
40.2
14.9
17.8
15.3
23.8
39.3
114.7
40.7
95.1
69.5
42.8
89.5
42
144.8
51.4
27.1
87
33.5
58.4
39.8
30.1
75.9
76.4
38.4
44.9
32.2
31.9
36.3
53.9
19.1
64.1
62.4
37.6
94.9
59.9
54
49.8
3.1
4.6
18.6
42.7
73.1
37.5
29.4
14.7
85.4
63.7
63.4
26.6
17.5
52.7
43
39.9
52.1
43.3
45.1
38.3
72.8
54.6
64.5
33.1
19.9
34.7
11.9
15.8
80.6
65.9
117
25.7
33.6
46.8
63.9
56.9
36.2
41.5
28
42.6
28.9
30.7
56.6
118.7
23.8
22
33.5
46.4
20.9
10.8
45.1
49.9
109.6
88.1
100.1
68.2
10.3
52.8
29.6
65.6
15
58.3
16.7
37
44.1
73.9
55.7
65.3
47
40.5
97.8
29.2
11.1
17.4
24.9
16.5
84.9
93.7
51.9
100
89
40.9
48
25.6
49.4
14.7
40.9
26
29.2
19.2
100.8
34.1
65.4
29.9
34.7
28
32.2
22.5
11.9
23.5
80.2
26.7
199.3
143.8
47.2
79
16.5
31.3
30.8
18.2
35.8
16.8
28
121.1
16.7
118.1
113.7
106.2
64.1
50.1
22.6
4.6
12.4
23.7
42.9
113
98.2
57.3
36.7
22.6
38.7
40.2
38.6
25.6
25
71.4
46.6
101.7
180
32.8
70.2
33
54.9
29.6
19.6
5.6
8
22.1
108.6
130.3
101.4
45.9
108.6
41
20.7
22.3
51.1
22.5
13.1
32.5
55.7
118.9
86.1
53.2
80.7
42.4
43.2
33.4
60.7
13.5
8.3
74.1
22.5
84.4
59.5
47.1
90.7
22
16.1
39.2
14.3
9.1
19.2
14.2
74
147.2
69.5
79.4
85.2
0.4
49.4
27.8
7.7
47.9
50.5
37.7
22
91.1
32.8
108.9
28.2
38.4
54.9
44.5
29.2
18.3
27.8
29.7
43
111.9
81.7
122.6
21.1
62.6
48.7
37.9
39.7
35.9
43.9
27.6
34
55.1
143.1
54.1
108.6
68.7
54.1
11.8
22.8
41.2
28.4
20.4
125.1
81.5
117.9
54.2
45.6
27
16.5
15.3
31.8
28.5
31.5
35.6
33.9
43.2
39.7
49.7
70.4
104.5
21.9
32
36.9
19.7
24
28.5
20.6
47.2
114.1
88.4
106.3
28.3
1.8
58
44.5
10.2
34.9
20.1
87.9
55.6
43.1
48.6
52.4
43.9
49.8
26.5
53.2
46.7
19.2
52.4
107.4
83.4
32.7
57.8
51.4
28
30.3
28.3
13.1
27.1
10.9
20.2
54.1
80.3
106.2
8.2
74.9
38.3
47.4
20.7
38.6
32.9
33.1
54.1
76
62.7
27.2
96.9
81.2
151.2
53.4
37.1
38
29.7
30.1
37.6
49.6
68.4
145.5
128.8
32.2
45.8
38.8
52.5
55.6
42.5
43.9
18.9
28.2
1.1
39.8
174.7
46.9
37.6
29.7
59.3
35.8
62.9
29.9
69.4
61.4
127.3
203.2
35.3
28.5
30.5
46
50.8
33.1
65.3
42.9
28.1
32.3
90.5
57.9
45.9
45.6
57.6
60
34.3
27.9
53.8
36.2
9.4
10.8
68.6
12.6
14.6
22.7
87.4
33.7
38.4
30
8.4
14.5
53.2
83.3
111.7
81.4
142.6
26.1
89
62
44.2
85
68
29
16
69
57
144
81
98
86
35
52
49
32
34
10
102
109
71
28
26
50
45
61
26
22
31
27
63
106
67
153
45
73
37
30
53
24
21
24
39
44
98
69
40
68
52
26
30
37
48
43
80
76
102
101
46
54
67
56
55
50
47
27
132
84
59
70
32
97
86
55
13
66
24
42
69
108
12
122
52
61
80
74
51
35
13
27
64
90
92
195
64
33
32
73
52.6111111111111
37.9888888888889
17.5666666666667
29.0888888888889
63.4333333333333
89.1
91.3
181.966666666667
64.7222222222222
36.2777777777778
35.6555555555556
72.2777777777778
54.2222222222222
40.9777777777778
22.1333333333333
31.1777777777778
62.8666666666667
88.2
90.6
168.933333333333
65.4444444444445
39.5555555555556
39.3111111111111
71.5555555555555
55.8333333333333
43.9666666666667
26.7
33.2666666666667
62.3
87.3
89.9
155.9
66.1666666666667
42.8333333333333
42.9666666666667
70.8333333333333
57.4444444444445
46.9555555555555
31.2666666666667
35.3555555555556
61.7333333333333
86.4
89.2
142.866666666667
66.8888888888889
46.1111111111111
46.6222222222222
70.1111111111111
59.0555555555556
49.9444444444444
35.8333333333333
37.4444444444444
61.1666666666667
85.5
88.5
129.833333333333
67.6111111111111
49.3888888888889
50.2777777777778
69.3888888888889
60.6666666666667
52.9333333333333
40.4
39.5333333333333
60.6
84.6
87.8
116.8
68.3333333333334
52.6666666666667
53.9333333333334
68.6666666666666
62.2777777777778
\ No newline at end of file
0,0
0,0
0,0
3,0
74,0
182,0
110,0
103,0
85,0
20,0
0,0
0,0
0,0
0,0
0,0
26,0
76,0
131,0
151,0
103,0
78,0
32,0
0,0
0,0
0,0
0,0
0,0
11,0
81,0
146,0
187,0
112,0
77,0
34,0
0,0
0,0
0,0
0,0
0,0
7,0
78,0
165,0
135,0
117,0
85,0
9,0
0,0
0,0
0,0
0,0
0,0
14,0
95,0
181,0
199,0
138,0
51,0
32,0
0,0
0,0
0,0
0,0
0,0
38,0
116,0
184,0
155,0
93,0
52,0
7,0
0,0
0,0
0,0
0,0
0,0
4,0
52,0
159,0
158,0
110,0
96,0
60,0
0,0
0,0
0,0
0,0
0,0
54,0
131,0
163,0
152,0
103,0
69,0
40,0
0,0
0,0
0,0
0,0
0,0
11,0
58,0
150,0
167,0
98,0
84,0
35,0
0,0
0,0
0,0
0,0
0,0
32,0
113,0
157,0
155,0
125,0
53,0
0,0
0,0
0,0
0,0
0,0
0,0
6,0
53,0
122,0
95,0
107,0
85,0
0,0
0,0
0,0
0,0
0,0
0,0
11,0
121,0
164,0
106,0
125,0
65,0
21,0
0,0
0,0
0,0
0,0
0,0
14,0
46,0
168,0
137,0
58,0
80,0
51,0
0,0
0,0
0,0
0,0
0,0
9,0
85,0
211,0
195,0
98,0
47,0
43,0
0,0
0,0
0,0
0,0
0,0
14,0
82,0
118,0
157,0
109,0
78,0
20,0
0,0
0,0
0,0
0,0
0,0
52,0
129,0
103,0
121,0
118,0
72,0
50,0
0,0
0,0
0,0
0,0
0,0
38,0
141,0
144,0
151,0
89,0
64,0
51,0
0,0
0,0
0,0
0,0
0,0
11,0
78,0
121,0
120,0
142,0
70,0
38,0
0,0
0,0
0,0
0,0
0,0
45,0
115,0
169,0
143,0
68,0
47,0
0,0
0,0
0,0
0,0
0,0
0,0
7,0
70,0
160,0
133,0
77,0
74,0
25,0
0,0
0,0
0,0
0,0
0,0
17,0
93,0
199,0
177,0
100,0
52,0
39,0
0,0
0,0
0,0
0,0
0,0
36,0
101,0
185,0
149,0
92,0
62,0
38,0
0,0
0,0
0,0
0,0
0,0
37,0
77,0
129,0
135,0
107,0
72,0
26,0
0,0
0,0
0,0
0,0
0,0
33,0
109,0
166,0
142,0
109,0
58,0
52,0
0,0
0,0
0,0
0,0
0,0
12,0
94,0
173,0
172,0
128,0
69,0
23,0
0,0
0,0
0,0
0,0
0,0
17,0
111,0
126,0
112,0
109,0
68,0
34,0
0,0
0,0
0,0
0,0
0,0
36,0
77,0
126,0
151,0
91,0
94,0
48,0
0,0
0,0
0,0
0,0
0,0
44,0
122,0
172,0
136,0
121,0
71,0
57,0
0,0
0,0
0,0
0,0
0,0
17,0
97,0
130,0
138,0
112,0
71,0
14,0
0,0
0,0
0,0
0,0
0,0
8,0
67,0
146,0
182,0
114,0
72,0
3,0
0,0
0,0
0,0
0,0
0,0
3,0
65,0
185,0
168,0
68,0
82,0
36,0
0,0
0,0
0,0
0,0
0,0
48,0
76,0
214,0
182,0
88,0
58,0
50,0
0,0
0,0
0,0
0,0
0,0
58,0
113,0
135,0
154,0
102,0
57,0
50,0
0,0
0,0
0,0
0,0
0,0
56,0
98,0
123,0
201,0
124,0
60,0
29,0
0,0
0,0
0,0
0,0
0,0
33,0
100,0
165,0
175,0
111,0
57,0
10,0
0,0
0,0
0,0
0,0
0,0
16,0
99,0
107,0
175,0
102,0
68,0
50,0
0,0
0,0
0,0
0,0
0,0
6,0
63,0
107,0
169,0
125,0
89,0
44,0
0,0
0,0
0,0
0,0
0,0
17,0
108,0
123,0
155,0
140,0
82,0
51,0
0,0
0,0
0,0
0,0
0,0
14,0
83,0
197,0
159,0
95,0
84,0
43,0
0,0
0,0
0,0
0,0
0,0
58,0
124,0
157,0
159,0
135,0
50,0
42,0
0,0
0,0
0,0
0,0
0,0
50,0
81,0
127,0
127,0
105,0
66,0
63,0
0,0
0,0
0,0
0,0
0,0
15,0
92,0
154,0
156,0
137,0
86,0
56,0
0,0
0,0
0,0
0,0
0,0
38,0
114,0
186,0
195,0
119,0
60,0
14,0
0,0
0,0
0,0
0,0
0,0
15,0
97,0
183,0
204,0
134,0
56,0
49,0
0
0
0
0
0
34
22
141
161
114
79
69
0
0
0
0
0
18
97
180
146
152
51
55
0
0
0
0
0
24
109
124
157
155
88
35
0
0
0
0
0
13
89
150
136
147
100
60
0
0
0
0
0
32
126
167
149
118
70
48
0
0
0
0
0
12
57
114
117
168
79
53
0
0
\ No newline at end of file
56.3
38.7
31.2
53.2
43.5
24.4
133.6
40.8
76.7
84.2
45.7
49.5
4.9
12.5
21.1
41.7
55.9
62.4
58.5
58.2
50.3
28.2
86.5
21
96.3
51.3
62.9
97.2
20.7
97.8
31.3
34.1
65.4
130.1
49.9
25.7
46.9
15.6
9
68.1
24
57.2
113.8
16.7
62.3
135.7
127.6
32.7
2.6
16.5
36.8
52.2
20.2
70.5
30.1
21.1
47.4
74.9
54.4
40.4
29.4
65
39
45
48.1
10
53.6
100.6
70.7
83.2
70.2
64.8
12.9
41.6
22.9
37.7
85
43
87.9
55
19.8
35.7
49.4
56.8
33.7
35
37.5
41.8
57.6
60.6
31.9
63
27.5
41
33.3
36.3
38.6
18.9
45.1
46.9
98.5
63.1
41
78.5
13.7
33.2
41.6
45.7
17.8
75.7
18.3
44.9
51.9
67.1
76.8
37
66.3
60.8
102.7
54.9
26.8
19.9
42.7
12
74.6
72.1
134.1
41.5
45
77.6
50.4
25.6
74.5
20.9
64.4
32.2
9.5
23.4
178
19.5
99.4
50.1
37.6
85.2
28.9
34.3
16.9
39.5
115.3
50.9
88.3
156.9
60
21.9
48.1
68.7
53.3
40.3
51.3
27.7
29.9
22.6
7.5
87.2
75.3
74.5
46
83.6
56.5
8.4
9.7
44.1
54.8
105.3
64.6
58.6
55.2
92.9
48.6
67.3
55.7
52.7
38.3
48.1
43.3
113.5
102.3
52.9
26.4
91.3
76.2
42.5
36
0.6
19.5
16.2
16.4
76.6
66.8
102.2
73
50
35.3
51.4
96
35
6.2
79
38
127
80
26
88
52
81
64
51
26
23
45
14
61
65
170
92
62
27
58
9
19
16
47
75
117
64
105
95
0.2
48
43
25
30
13
7
29
56
75
62
76
54
35
38
33
25
42
36.1
46.4
62
65.9
100.8
13.4
132
63
46
42.8
39.1
54.9
35.3
59.3
76
86.2
42.9
89.4
121.9
67.6
17.1
30.4
40.6
13
38.5
52.7
94.6
68.6
55.2
81.7
77.1
26.6
53.6
52.9
36.3
31.7
60.7
31.4
18.9
39.6
29.2
82.8
182.7
97.1
19.2
46.3
17.2
36.3
30.3
9.9
80.8
122.2
57.4
75.2
60
5.9
-16.1
60.4
19.2
48.7
13.8
40.7
88.8
32.1
70.5
26
54.1
78.8
43.2
53.8
41.2
21.4
79.8
50.8
56.2
49.2
47.8
38.6
18.2
116.7
36
19.4
19.7
9.4
16.2
53.1
130.2
55.2
48.7
68.9
18.8
32.5
19.4
40
27.7
32.6
32.6
32.9
144.1
2.3
29.7
65.6
166.7
45.7
29.1
24.5
36.7
41.3
42.3
81.7
47.7
66.3
120.7
45.6
87.2
45.4
50.7
68
41
15
39
56
5.6
34
76
28.6
41.9
18.6
46.5
32
33.9
44.3
24
20.1
85.5
101.1
82.3
75.8
23.6
26.5
90
41.8
74.7
39.7
24.9
72.6
114.7
17.7
21.1
42.3
63.8
53.1
42.8
42
54.1
29.7
30.1
19.7
52
79.6
19.9
76.2
92.5
50.7
31.6
35
9.6
33
45.2
30.8
89.8
19
94.8
41.8
41
37
31
50
18
45
31
57
161.8
42
45
38
54
52
37
52.4
38
39
72
38
141
46
17
12.4
31
18
84
24
29
39
40
40
23
58
92
52
50
70
22
58
43
24
23
26
34
71
69
100
55
59
20
29
40
70
32.3
93
65
161
101
46
45
38
42
36.2
25.1
26.5
18
30
88.1
70.4
32.1
19.2
130.5
51.4
39.7
17.5
37.9
16.6
26.2
96.6
46.4
13.1
49.3
48.2
39
73.6
82.5
40.8
21.5
24.8
24.8
30.6
30.6
53.8
31.4
81.6
40.8
46.6
50.4
84.4
24.8
40.8
63.9
40.1
121.1
28.2
130.2
51
112.4
38.5
41.2
39.5
28.6
68.2
25.8
95.6
67.1
81.7
44.9
202.2
35.4
47.7
46
40.5
17.5
17.4
21.6
50.5
132.6
23.1
69.6
28.7
51.2
12.5
55.6
49.1
30.8
3.7
44.9
62.4
127.1
54.9
36.6
22.8
23.4
44.6
42.7
85.9
55.8
27.1
59.9
35
48.5
74.9
135.9
56.7
25.1
82.5
24.9
25.2
30.3
52.8
35.2
60
91.8
113.5
22.1
56.9
83.4
67.7
76.2
\ No newline at end of file
...@@ -26,10 +26,10 @@ module carbon_model_rothc ...@@ -26,10 +26,10 @@ module carbon_model_rothc
call set_pool(pid = n_CHUM, name = 'CHUM', initial_value = 8.11102631424 , alias = CHUM) call set_pool(pid = n_CHUM, name = 'CHUM', initial_value = 8.11102631424 , alias = CHUM)
! litterfall ! litterfall
call set_flux(fid = n_F, pid_out = n_Cveg, pid_in = n_CDPM, name = 'fdpm*lambdac') call set_flux(fid = n_F, pid_out = n_Cveg, pid_in = n_CDPM, name = 'fdpm*lambd')
call set_mult(n_F, 'lin', x_ij = lambd, k = fdpm) call set_mult(n_F, 'lin', x_ij = lambd, k = fdpm)
call set_flux(fid = n_F, pid_out = n_Cveg, pid_in = n_CRPM, name = '(1-fdpm)*lambdac') call set_flux(fid = n_F, pid_out = n_Cveg, pid_in = n_CRPM, name = '(1-fdpm)*lambd')
call set_mult(n_F, 'lin', x_ij = lambd, k = (1 - fdpm)) call set_mult(n_F, 'lin', x_ij = lambd, k = (1 - fdpm))
! to atmosphere ! to atmosphere
...@@ -46,11 +46,11 @@ module carbon_model_rothc ...@@ -46,11 +46,11 @@ module carbon_model_rothc
call set_mult(n_F, 'lin', x = RHUM) call set_mult(n_F, 'lin', x = RHUM)
! atmosphere fall ! atmosphere fall
call set_flux(fid = n_F, pid_out = n_Catm, pid_in = n_CBIO, name = '(0.46*bettar)*Rs') call set_flux(fid = n_F, pid_out = n_Catm, pid_in = n_CBIO, name = 'BIO_prop*Rs')
call set_mult(n_F, 'lin', x = Rs, k_ij = BIO_prop) ! (0.46*bettar/100) call set_mult(n_F, 'lin', x = Rs, k_ij = BIO_prop)
call set_flux(fid = n_F, pid_out = n_Catm, pid_in = n_CHUM, name = '(0.54*bettar)*Rs') call set_flux(fid = n_F, pid_out = n_Catm, pid_in = n_CHUM, name = 'HUM_prop*Rs')
call set_mult(n_F, 'lin', x = Rs, k_ij = HUM_prop) ! (0.54*bettar/100) call set_mult(n_F, 'lin', x = Rs, k_ij = HUM_prop)
end subroutine carbon_model_assembly end subroutine carbon_model_assembly
......
...@@ -2,7 +2,7 @@ module carbon_model_rothc_aux ...@@ -2,7 +2,7 @@ module carbon_model_rothc_aux
! интерфейс ! интерфейс
! ----------------------------------------------- Use pack ------------------------------------------------------------ ! ----------------------------------------------- Use pack ------------------------------------------------------------
use environment_core, only: Tsoil, Wsoil, lambd, veg, sw, bettar use environment_core, only: Tsoil, Wsoil, lambd, veg, sw, bettar, Evpar, pr
use config, only: station_name, station_opt, carbon_model_type use config, only: station_name, station_opt, carbon_model_type
use const, only: yrs, pi, nmonth use const, only: yrs, pi, nmonth
use grid, only: date_c, i0, i1, j0, j1, ml, dt use grid, only: date_c, i0, i1, j0, j1, ml, dt
...@@ -38,16 +38,14 @@ module carbon_model_rothc_aux ...@@ -38,16 +38,14 @@ module carbon_model_rothc_aux
real :: fdpm !< Функция определяющая разлагаемую часть опада, [dim] real :: fdpm !< Функция определяющая разлагаемую часть опада, [dim]
! ------- Coefficients ------------ ! ------- Coefficients ------------
real :: ks(ntiles) = (/ 3.22*(1.E-7), 9.65*(1.E-9), 2.12*(1.E-8), 6.43*(1.E-10)/) real :: ks(ntiles) = (/ 3.22*(1.E-7), 9.65*(1.E-9), 2.12*(1.E-8), 6.43*(1.E-10)/)
!real :: ks(ntiles) = (/ 3.22*(1.E-6), 9.65*(1.E-9), 2.12*(1.E-8), 6.43*(1.E-9)/)
!real :: ks(ntiles) = (/ 10., 0.3, 0.66, 0.02/)
!< Скорость дыхания единицы массы каждого пула в стандартных условиях, [1-s] !< Скорость дыхания единицы массы каждого пула в стандартных условиях, [1-s]
real, parameter :: alphadr = 1.44 !< Определяет соотношение поступления опада между DPM и RPM, [dim], real, parameter :: alphadr = 1.44 !< Определяет соотношение поступления опада между DPM и RPM, [dim],
! Для деревьев, для кустарников, для натуральной травы, для СХ (/ 0.25, 0.33, 0.67, 1.44 /) ! Для деревьев, для кустарников, для натуральной травы, для СХ (/ 0.25, 0.33, 0.67, 1.44 /)
real, allocatable :: SA(:,:) ! Определяет распределение выдыхаемой части пула по остальным пулам real, allocatable :: SA(:,:) !< Определяет распределение выдыхаемой части пула по остальным пулам
real, allocatable :: CO2_prop(:,:) ! Определяет распределение выдыхаемой части пула по остальным пулам real, allocatable :: CO2_prop(:,:) !< Определяет распределение выдыхаемой части пула в CO2
real, allocatable :: BIO_prop(:,:) ! Определяет распределение выдыхаемой части пула по остальным пулам real, allocatable :: BIO_prop(:,:) !< Определяет распределение выдыхаемой части пула в CBIO
real, allocatable :: HUM_prop(:,:) ! Определяет распределение выдыхаемой части пула по остальным пулам real, allocatable :: HUM_prop(:,:) !< Определяет распределение выдыхаемой части пула в CHUM
! ------ Climate variables -------- ! ------ Climate variables --------
! ------ Determined within ! ------ Determined within
...@@ -57,23 +55,15 @@ module carbon_model_rothc_aux ...@@ -57,23 +55,15 @@ module carbon_model_rothc_aux
!real :: bettar(i0:i1,j0:j1) !< Fraction of soil respiration, [%] from 0 to 1 !real :: bettar(i0:i1,j0:j1) !< Fraction of soil respiration, [%] from 0 to 1
!real :: sw !< Влажность увядания !real :: sw !< Влажность увядания
! ------ Initialization ------ ! ------ Initialization ------
real :: cs ! Сумма пулов real :: cs !< Сумма пулов
! ----- The topsoil moisture deficit ------- ! ----- The topsoil moisture deficit -------
real, allocatable :: TSMD(:,:) ! the topsoil moisture deficit real :: deep = 30. ! глубина слоя (@todo сопоставить с переменной Ml)
real :: deep = 30. ! глубина слоя real, allocatable :: TSMD (:,:) !< the topsoil moisture deficit
real, allocatable :: SMD(:,:) ! TSMD на актуальной глубине real, allocatable :: SMD (:,:) !< TSMD на актуальной глубине
real, allocatable :: bSMD(:,:) ! Случай TSMD для голой почвы (bare soil) real, allocatable :: bSMD (:,:) !< Случай TSMD для голой почвы (bare soil)
real, allocatable :: moi_bal(:,:) ! rainfall minus evaporation real, allocatable :: moi_bal (:,:) !< pr minus evaporation
real, allocatable :: aTSMD(:,:) ! accumulated topsoil mositure deficit real, allocatable :: aTSMD (:,:) !< accumulated topsoil mositure deficit
real, allocatable :: evpar_part(:,:) !< 75% from evaporation
real, allocatable :: in_rainfall(:)
real, allocatable :: in_Evpar(:)
real, allocatable :: rainfall(:,:) ! осадки (из внешних данных)
real, allocatable :: Evpar(:,:) ! Evaporation
real, allocatable :: evpar_part(:,:) ! 75% from evaporation
integer :: mncX
integer :: mnc = 516 ! Случай для Ростова
contains contains
...@@ -95,26 +85,12 @@ contains ...@@ -95,26 +85,12 @@ contains
allocate(bSMD(i0:i1,j0:j1)) allocate(bSMD(i0:i1,j0:j1))
allocate(moi_bal(i0:i1,j0:j1)) allocate(moi_bal(i0:i1,j0:j1))
allocate(aTSMD(i0:i1,j0:j1)) allocate(aTSMD(i0:i1,j0:j1))
allocate(in_rainfall(0:mnc))
allocate(in_Evpar(0:mnc))
allocate(rainfall(i0:i1,j0:j1))
allocate(Evpar(i0:i1,j0:j1))
allocate(evpar_part(i0:i1,j0:j1)) allocate(evpar_part(i0:i1,j0:j1))
TSMD(:,:) = -(20. + 1.3*bettar(:,:) - 0.01*bettar(:,:)**2) TSMD(:,:) = -(20. + 1.3*bettar(:,:) - 0.01*bettar(:,:)**2)
SMD(:,:) = TSMD(:,:)/23.*deep SMD(:,:) = TSMD(:,:)/23.*deep
bSMD(:,:) = SMD(:,:)/1.8 bSMD(:,:) = SMD(:,:)/1.8
open (unit = 1, file = 'initial_value/rainfall.txt', status='unknown')
read(1,*) in_rainfall(1:mnc)
close (1)
open (unit = 1, file = 'initial_value/Evpar.txt', status='unknown')
read(1,*) in_Evpar(1:mnc)
close (1)
mncX = 1
! main init part of Rothc ! main init part of Rothc
do i = 1, ml do i = 1, ml
s0 (:,:,i) = 0.5*(1. + sw(:,:)) s0 (:,:,i) = 0.5*(1. + sw(:,:))
...@@ -149,11 +125,9 @@ contains ...@@ -149,11 +125,9 @@ contains
integer, intent(in) :: ii, jj integer, intent(in) :: ii, jj
! ---- TPMD_calc_at_cell ---- ! ---- TPMD_calc_at_cell ----
rainfall(ii,jj) = in_rainfall(mncX)
Evpar(ii,jj) = in_Evpar(mncX)
evpar_part(ii,jj) = 0.75*Evpar(ii,jj) evpar_part(ii,jj) = 0.75*Evpar(ii,jj)
moi_bal(ii,jj) = rainfall(ii,jj) - evpar_part(ii,jj) moi_bal(ii,jj) = pr(ii,jj) - evpar_part(ii,jj)
if (veg(ii,jj) == 1) then if (veg(ii,jj) == 1) then
aTSMD(ii,jj) = max(SMD(ii,jj), min(0.,aTSMD(ii,jj) + moi_bal(ii,jj))) aTSMD(ii,jj) = max(SMD(ii,jj), min(0.,aTSMD(ii,jj) + moi_bal(ii,jj)))
...@@ -161,8 +135,6 @@ contains ...@@ -161,8 +135,6 @@ contains
aTSMD(ii,jj) = max(min(bSMD(ii,jj),aTSMD(ii,jj)),min(0.,aTSMD(ii,jj) + moi_bal(ii,jj))) aTSMD(ii,jj) = max(min(bSMD(ii,jj),aTSMD(ii,jj)),min(0.,aTSMD(ii,jj) + moi_bal(ii,jj)))
end if end if
mncX = mncX + 1
! --------------
! main part of rothc_calc_at_cell ! main part of rothc_calc_at_cell
if (Wsoil(ii,jj,1) >= s0(ii,jj,1)) then if (Wsoil(ii,jj,1) >= s0(ii,jj,1)) then
...@@ -187,7 +159,6 @@ contains ...@@ -187,7 +159,6 @@ contains
Ft = 0 Ft = 0
end if end if
Fv = 0.6 + 0.4*(1. - veg(ii,jj)) Fv = 0.6 + 0.4*(1. - veg(ii,jj))
RDPM = CDPM(ii,jj,1)*(1. - exp(-Ft*Fv*Fs*ks(1)*dt))/dt RDPM = CDPM(ii,jj,1)*(1. - exp(-Ft*Fv*Fs*ks(1)*dt))/dt
...@@ -195,16 +166,6 @@ contains ...@@ -195,16 +166,6 @@ contains
RBIO = CBIO(ii,jj,1)*(1. - exp(-Ft*Fv*Fs*ks(3)*dt))/dt RBIO = CBIO(ii,jj,1)*(1. - exp(-Ft*Fv*Fs*ks(3)*dt))/dt
RHUM = CHUM(ii,jj,1)*(1. - exp(-Ft*Fv*Fs*ks(4)*dt))/dt RHUM = CHUM(ii,jj,1)*(1. - exp(-Ft*Fv*Fs*ks(4)*dt))/dt
!RDPM = CDPM(ii,jj,1)*Ft*Fv*Fs*ks(1)
!RRPM = CRPM(ii,jj,1)*Ft*Fv*Fs*ks(2)
!RBIO = CBIO(ii,jj,1)*Ft*Fv*Fs*ks(3)
!RHUM = CHUM(ii,jj,1)*Ft*Fv*Fs*ks(4)
print*, 'RDPM:', RDPM
print*, 'RRPM:', RRPM
print*, 'RBIO:', RBIO
print*, 'RHUM:', RHUM
Rs = RDPM + RRPM + RBIO + RHUM Rs = RDPM + RRPM + RBIO + RHUM
! ---- Initialization ------ ! ---- Initialization ------
......
...@@ -8,7 +8,7 @@ module environment_core ...@@ -8,7 +8,7 @@ module environment_core
private private
public :: Temp, e, Rswd, Rlwd, p, pr, Wind, Tsrf, Tgr, snow, ra, Tsoil, Wsoil, Isoil public :: Temp, e, Rswd, Rlwd, p, pr, Wind, Tsrf, Tgr, snow, ra, Tsoil, Wsoil, Isoil
public :: sw, bettar, veg, kd, lambd, rhodry, g2gw_to_cmw public :: sw, bettar, veg, kd, lambd, rhodry, g2gw_to_cmw, Evpar
public :: environment_core_init public :: environment_core_init
! переменные ! переменные
...@@ -37,6 +37,8 @@ module environment_core ...@@ -37,6 +37,8 @@ module environment_core
real, allocatable, dimension(:,:,:), target :: rhodry !< плотность сухой почвы, г/см3 real, allocatable, dimension(:,:,:), target :: rhodry !< плотность сухой почвы, г/см3
real, allocatable, dimension(:,:,:), target :: g2gw_to_cmw !< г/г воды в см воды в слоях почвы real, allocatable, dimension(:,:,:), target :: g2gw_to_cmw !< г/г воды в см воды в слоях почвы
real, allocatable, dimension(:,:), target :: Evpar !< Evaporation
contains contains
...@@ -70,6 +72,7 @@ contains ...@@ -70,6 +72,7 @@ contains
allocate(lambd(i0:i1,j0:j1)) allocate(lambd(i0:i1,j0:j1))
allocate(rhodry(i0:i1,j0:j1,ml)) allocate(rhodry(i0:i1,j0:j1,ml))
allocate(g2gw_to_cmw(i0:i1,j0:j1,ml)) allocate(g2gw_to_cmw(i0:i1,j0:j1,ml))
allocate(Evpar(i0:i1,j0:j1))
g2gw_to_cmw = miss_v g2gw_to_cmw = miss_v
rhodry = miss_v rhodry = miss_v
...@@ -93,7 +96,7 @@ contains ...@@ -93,7 +96,7 @@ contains
Tsoil = miss_v Tsoil = miss_v
Wsoil = miss_v Wsoil = miss_v
Isoil = miss_v Isoil = miss_v
Evpar = miss_v
end subroutine end subroutine
......
...@@ -2,8 +2,8 @@ module environment_data_station ...@@ -2,8 +2,8 @@ module environment_data_station
! интерфейс ! интерфейс
! ----------------------------------------------- Use pack ------------------------------------------------------------ ! ----------------------------------------------- Use pack ------------------------------------------------------------
use environment_core, only: Temp, e, Rlwd, Rswd, p, pr, Wind, ra, Tsrf, Tsoil, Wsoil, Isoil, snow, Tgr, & use environment_core!, only: Temp, e, Rlwd, Rswd, p, pr, Wind, ra, Tsrf, Tsoil, Wsoil, Isoil, snow, Tgr, &
& sw, bettar, kd, veg, lambd, rhodry, g2gw_to_cmw !& sw, bettar, kd, veg, lambd, rhodry, g2gw_to_cmw
use grid, only: date_c, date_fst, date_lst, dt, i0, i1, j0, j1, ml, ms,dz use grid, only: date_c, date_fst, date_lst, dt, i0, i1, j0, j1, ml, ms,dz
use config, only: station_name, station_opt use config, only: station_name, station_opt
use const, only: yrs, nmonth, pi use const, only: yrs, nmonth, pi
...@@ -32,6 +32,9 @@ module environment_data_station ...@@ -32,6 +32,9 @@ module environment_data_station
real(kind = 16), allocatable :: in_veg (:) !< Поступление извне данных средней концентрации растений в месяц, [%] real(kind = 16), allocatable :: in_veg (:) !< Поступление извне данных средней концентрации растений в месяц, [%]
real(kind = 16), allocatable :: in_wsoil(:) !< Поступление извне данных влажности почвы, [dim] real(kind = 16), allocatable :: in_wsoil(:) !< Поступление извне данных влажности почвы, [dim]
real(kind = 16), allocatable :: in_lambd(:) !< Поступление извне данных по поступлению углерода в почву, [kg/m**3 / year] real(kind = 16), allocatable :: in_lambd(:) !< Поступление извне данных по поступлению углерода в почву, [kg/m**3 / year]
real, allocatable :: in_rainfall(:) !< Поступление извне данных влажности почвы, [mm]
real, allocatable :: in_Evpar(:) !< Поступление извне данных влажности почвы, [mm]
! station = DAO4 DAO3 TRGK VLDR ROST ! station = DAO4 DAO3 TRGK VLDR ROST
! ---- For INMCM ---- ! ---- For INMCM ----
real(kind = 16) :: rhodry_in(station_max) = (/1.26, 1.26, 1.23, 1.30, 1.28/) real(kind = 16) :: rhodry_in(station_max) = (/1.26, 1.26, 1.23, 1.30, 1.28/)
...@@ -39,6 +42,7 @@ module environment_data_station ...@@ -39,6 +42,7 @@ module environment_data_station
! ---- For Rothc ---- ! ---- For Rothc ----
real(kind = 16) :: sw_in (station_max) = (/0.120, 0.120, 0.120, 0.120, 0.146/) real(kind = 16) :: sw_in (station_max) = (/0.120, 0.120, 0.120, 0.120, 0.146/)
real(kind = 16) :: bettar_in(station_max) = (/25., 25., 7., 8., 44./) real(kind = 16) :: bettar_in(station_max) = (/25., 25., 7., 8., 44./)
! ---- For SOCs ---- ! ---- For SOCs ----
real :: kd_in (station_max,opt_max) real :: kd_in (station_max,opt_max)
! opt = 1 2 3 4 ! opt = 1 2 3 4
...@@ -64,6 +68,8 @@ contains ...@@ -64,6 +68,8 @@ contains
allocate(in_temp (0:mnc)) allocate(in_temp (0:mnc))
allocate(in_veg (0:mnc)) allocate(in_veg (0:mnc))
allocate(in_wsoil(0:mnc)) allocate(in_wsoil(0:mnc))
allocate(in_Evpar(0:mnc))
allocate(in_rainfall(0:mnc))
kd_in = kd_in/yrs kd_in = kd_in/yrs
...@@ -98,6 +104,14 @@ contains ...@@ -98,6 +104,14 @@ contains
read(1,*) in_veg(1:mnc) read(1,*) in_veg(1:mnc)
close (1) close (1)
open (unit = 1, file = 'initial_value/'//trim(station_name)//'_Rainfall.txt', status='unknown')
read(1,*) in_rainfall(1:mnc)
close (1)
open (unit = 1, file = 'initial_value/'//trim(station_name)//'_Evpar.txt', status='unknown')
read(1,*) in_Evpar(1:mnc)
close (1)
in_lambd = in_lambd/dt in_lambd = in_lambd/dt
rhodry(:,:,:) = rhodry_in(station_n) rhodry(:,:,:) = rhodry_in(station_n)
...@@ -117,6 +131,8 @@ contains ...@@ -117,6 +131,8 @@ contains
in_veg (0) = in_veg (1) in_veg (0) = in_veg (1)
in_wsoil(0) = in_wsoil(1) in_wsoil(0) = in_wsoil(1)
in_lambd(0) = in_lambd(1) in_lambd(0) = in_lambd(1)
in_Evpar(0) = in_Evpar(1)
in_rainfall(0) = in_rainfall(1)
else ! Для шага по времени месяц: else ! Для шага по времени месяц:
mncX = (date_fst%y - mnclot_fst(station_n))*nmonth + date_lst%m mncX = (date_fst%y - mnclot_fst(station_n))*nmonth + date_lst%m
end if end if
...@@ -124,7 +140,6 @@ contains ...@@ -124,7 +140,6 @@ contains
! Набор переменных которые должны быть, но не данны через качественные внешние данные ! Набор переменных которые должны быть, но не данны через качественные внешние данные
p(:,:) = 980. !< pressure p(:,:) = 980. !< pressure
Rlwd(:,:) = 350. !< longwave radiation Rlwd(:,:) = 350. !< longwave radiation
pr(:,:) = 0. !< precipitation
ra(:,:) = 20. !< aerodynamical resistance ra(:,:) = 20. !< aerodynamical resistance
snow(:,:) = 0. !< snow cover snow(:,:) = 0. !< snow cover
Isoil(:,:,:) = 0. !< part of ice in the soil Isoil(:,:,:) = 0. !< part of ice in the soil
...@@ -151,13 +166,16 @@ contains ...@@ -151,13 +166,16 @@ contains
Temp (:,:) = in_temp(mncX)*sin(pi*(nod-1)/(N(nom)-1)) + in_temp(mncX+1)*(1-sin(pi*(nod-1)/(N(nom)-1))) Temp (:,:) = in_temp(mncX)*sin(pi*(nod-1)/(N(nom)-1)) + in_temp(mncX+1)*(1-sin(pi*(nod-1)/(N(nom)-1)))
Wsoil(:,:,:) = in_wsoil(mncX)*sin(pi*(nod-1)/(N(nom)-1)) + in_wsoil(mncX+1)*(1-sin(pi*(nod-1)/(N(nom)-1))) Wsoil(:,:,:) = in_wsoil(mncX)*sin(pi*(nod-1)/(N(nom)-1)) + in_wsoil(mncX+1)*(1-sin(pi*(nod-1)/(N(nom)-1)))
veg (:,:) = in_veg(mncX) veg (:,:) = in_veg(mncX)
Evpar (:,:) = in_Evpar(mncX)/N(nom)
pr (:,:) = in_rainfall(mncX)/N(nom)
if (nod == (N(nom)/2)) mncX = mncX + 1 if (nod == (N(nom)/2)) mncX = mncX + 1
else ! Для шага по времени месяц: else ! Для шага по времени месяц:
lambd (:,:) = in_lambd(mncX) lambd (:,:) = in_lambd(mncX)
Temp (:,:) = in_temp(mncX) Temp (:,:) = in_temp(mncX)
Wsoil(:,:,:) = in_wsoil(mncX) Wsoil(:,:,:) = in_wsoil(mncX)
veg (:,:) = in_veg(mncX) veg (:,:) = in_veg(mncX)
Evpar (:,:) = in_Evpar(mncX)
pr (:,:) = in_rainfall(mncX)
mncX = mncX + 1 mncX = mncX + 1
end if end if
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment