Skip to content
Snippets Groups Projects
Commit 7e6c59e5 authored by Sumbel Shangareeva's avatar Sumbel Shangareeva
Browse files

samara station

parent c08c6618
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 1831 additions and 59 deletions
initial_value/carbon_intake/Samara_FYM_1.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
\ No newline at end of file
initial_value/carbon_intake/Samara_FYM_2.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
\ No newline at end of file
initial_value/carbon_intake/Samara_FYM_3.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
\ No newline at end of file
initial_value/carbon_intake/Samara_FYM_4.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
\ No newline at end of file
initial_value/carbon_intake/Samara_FYM_5.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
\ No newline at end of file
initial_value/carbon_intake/Samara_org_1.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
0.03458333333333335
initial_value/carbon_intake/Samara_org_2.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
0.08625
initial_value/carbon_intake/Samara_org_3.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
0.0808333333333335
initial_value/carbon_intake/Samara_org_4.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
0.0520833333333335
initial_value/carbon_intake/Samara_org_5.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
0.05875
initial_value/climate_intake/Samara_MOI.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0.29307555
0.2850379975
0.342754365
0.3745193475
0.24490356
0.2302856425
0.170772555
0.1622810375
0.1596031225
0.201217655
0.2790336625
0.274822225
0.26963425
0.258052835
0.32000732
0.3727264425
0.31160354
0.1765213
0.19037628
0.16722107
0.1698951725
0.2528533925
0.3651237475
0.3703422525
0.3481597825
0.342750555
0.3899383525
0.3823089625
0.31943893
0.1966438325
0.165931705
0.14996338
0.2668342575
0.2948303225
0.3369483975
0.334098825
0.310440065
0.3119239825
0.3584213275
0.3973312425
0.2971611025
0.352161405
0.2482605025
0.1590538025
0.153736115
0.200489045
0.33250045
0.310752875
0.30353546
0.2889518775
0.28258897
0.38412476
0.268013
0.183578495
0.1700859075
0.153987885
0.160037995
0.1998519875
0.252338405
0.24238205
0.23479843
0.230030055
0.32348633
0.3578529375
0.2633667
0.1675415
0.165218355
0.198852535
0.196079255
0.26552582
0.29351425
0.2844696025
0.28234482
0.306079865
0.39237214
0.349906925
0.258747095
0.2083091725
0.1520194975
0.17358398
0.180683135
0.184360505
0.2247009275
0.2219543425
0.21771622
0.2754364025
0.2775535625
0.35427093
0.2111282325
0.190109255
0.1544914225
0.1468009925
0.2078208925
0.2145423875
0.25331497
0.3212127675
0.291404725
0.273983005
0.27889252
0.3834877075
0.3301315325
0.29838562
0.1797943125
0.14901352
0.18980789
0.283710475
0.3589096075
0.3366470375
0.3459053075
0.325630185
0.3828277575
0.29605483
0.198360445
0.1755027775
0.17208481
0.1529502875
0.177223205
0.2443313575
0.36258698
0.3443756125
0.33743668
0.3154411425
0.32928848
0.36785888
0.267692565
0.19408798
0.1643257125
0.278987885
0.18255234
0.181537625
0.286800385
0.277633665
\ No newline at end of file
initial_value/climate_intake/Samara_Temp.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
-11.97117
-14.00886
-2.75069999999999
4.17532
18.29843
19.87832
21.35003
22.94216
14.56597
4.69723000000005
-2.11959999999999
-7.03396999999995
-12.05625
-8.29336999999998
-5.38369999999998
5.38778000000002
15.32467
24.40414
20.93765
19.272
17.51364
2.89193
-1.3768
-3.94632999999999
-12.27598
-3.45612
-1.61221999999998
8.89443
15.23672
20.58428
23.9778
26.97213
12.62875
4.28137000000004
-3.86478
-12.47938
-10.50376
-9.76864999999998
-3.35343
4.71664000000004
12.6664
16.4407
21.14285
22.82266
15.76333
5.404
0.499170000000049
-7.06881999999996
-11.72336
-12.16715
-10.8817
3.58785
15.49606
19.33697
25.0676
22.303
16.92498
7.46328
-3.29226999999997
-8.95956999999999
-11.77836
-8.83612999999997
-2.0317
6.61834000000005
16.37554
21.70846
21.83544
18.66116
12.21304
8.53976
-1.9074
-5.67779999999999
-4.02222
-5.21457999999996
0.503780000000006
6.65470000000005
15.3123000000001
19.31198
26.52523
20.20464
14.04537
7.28060000000005
-2.95456999999999
-13.00004
-10.23722
-14.45063
-6.52401999999995
7.44400000000002
20.79904
24.06576
25.91244
27.5048
12.44256
6.42935
-0.980799999999988
-7.72489999999999
-9.73439999999999
-4.38415999999995
-5.69265999999999
8.2441
10.76272
18.1514
23.49978
25.63683
14.60797
6.90800000000002
-1.27942999999999
-8.45056999999997
-16.01447
-7.84213999999997
1.44677000000001
9.57546000000002
18.12997
19.2197
23.81207
23.0865
15.6824
6.82089999999999
0.538540000000012
-8.58807999999999
-13.98026
-13.22697
-4.89274
11.5964
12.1121
22.76174
23.85494
18.88174
16.0881
7.14880000000005
-0.533639999999991
-5.86703999999997
\ No newline at end of file
initial_value/climate_intake/Samara_Veg.txt 0 → 100644
+ 132
− 0
View file @ 7e6c59e5
0
0
0
0
0
0
0
0
1
1
1
1
1
1
1
1
1
1
1
0
1
1
1
1
1
1
1
1
1
1
1
0
0
0
0
0
0
0
1
1
1
1
0
0
1
1
1
1
1
1
1
1
1
1
1
0
0
0
0
0
0
0
1
1
1
1
0
0
0
0
0
0
0
0
0
0
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
0
0
0
0
0
0
0
0
1
1
1
1
1
1
0
0
0
0
0
0
0
0
0
0
0
1
1
1
1
1
1
1
1
1
1
1
0
0
0
0
0
\ No newline at end of file
main_spotpy.py
+ 1
− 1
View file @ 7e6c59e5
......@@ -15,7 +15,7 @@ if not os.path.exists(res_path):
n_it = 100
setup = spot_setup(setup_path=config_nml_path, params_path=parameters_nml_path)
sampler = spotpy.algorithms.rope(setup, dbname=f'{res_path}/arkha_rope_{n_it}',
sampler = spotpy.algorithms.rope(setup, dbname=f'{res_path}/samara_rope_{n_it}',
dbformat='csv',
save_sim=False
)
......
parameters_config.nml
+ 5
− 4
View file @ 7e6c59e5
&parameters_config_namelist
kd_in(1:6,1) = 0.03, 0.03, 0.0402, 0.08581933, 0.012, 0.07115579
kd_in(1:6,2) = 0.03, 0.03, 0.05067, 0.08581933, 0.007, 0.007
kd_in(1:6,3) = 0.035, 0.035, 0.04185, 0.08581933, 0.007, 0.007
kd_in(1:6,4) = 0.03, 0.03, 0.0, 0.0, 0.0, 0.0
kd_in(1:7,1) = 0.03, 0.03, 0.0402, 0.08581933, 0.012, 0.07115579, 0.020638477
kd_in(1:7,2) = 0.03, 0.03, 0.05067, 0.08581933, 0.007, 0.007, 0.0
kd_in(1:7,3) = 0.035, 0.035, 0.04185, 0.08581933, 0.007, 0.007, 0.0
kd_in(1:7,4) = 0.03, 0.03, 0.0, 0.0, 0.0, 0.0, 0.0
kd_in(1:7,5) = 0.03, 0.03, 0.0, 0.0, 0.0, 0.0, 0.0
/
plot.py
+ 8
− 0
View file @ 7e6c59e5
......@@ -34,6 +34,14 @@ if (station == 'Arkha'):
f'results/{carbon_model_type}/{station}_1.txt'
]
if (station == 'Samara'):
model_data_paths = [
f'results/{carbon_model_type}/{station}_1.txt',
# f'results/{carbon_model_type}/{station}_2.txt',
# f'results/{carbon_model_type}/{station}_3.txt',
# f'results/{carbon_model_type}/{station}_4.txt',
# f'results/{carbon_model_type}/{station}_5.txt'
]
if (station == 'DAO3') or (station == 'DAO4'):
model_data_paths = [
......
source/carbon_pools.f90
+ 63
− 33
View file @ 7e6c59e5
......@@ -225,18 +225,7 @@ contains
Csoil_iv = Csoil_st_in(3,3)
Csoilb_iv = Csoilb_st_in(3,3)
case('4')
CDPM_iv = CDPM_st_in(3,4)
CRPM_iv = CRPM_st_in(3,4)
CBIO_iv = CBIO_st_in(3,4)
CHUM_iv = CHUM_st_in(3,4)
CIOM_iv = CIOM_st_in(3,4)
C_s1_iv = C_s1_st_in(3,4)
C_s2_iv = C_s2_st_in(3,4)
Csoil_iv = Csoil_st_in(3,4)
Csoilb_iv = Csoilb_st_in(3,4)
end select
case('VLDMR')
......@@ -277,18 +266,7 @@ contains
Csoil_iv = Csoil_st_in(4,3)
Csoilb_iv = Csoilb_st_in(4,3)
case('4')
CDPM_iv = CDPM_st_in(4,4)
CRPM_iv = CRPM_st_in(4,4)
CBIO_iv = CBIO_st_in(4,4)
CHUM_iv = CHUM_st_in(4,4)
CIOM_iv = CIOM_st_in(4,4)
C_s1_iv = C_s1_st_in(4,4)
C_s2_iv = C_s2_st_in(4,4)
Csoil_iv = Csoil_st_in(4,4)
Csoilb_iv = Csoilb_st_in(4,4)
end select
case('Rostov')
select case(station_opt)
......@@ -328,18 +306,7 @@ contains
Csoil_iv = Csoil_st_in(5,3)
Csoilb_iv = Csoilb_st_in(5,3)
case('4')
CDPM_iv = CDPM_st_in(5,4)
CRPM_iv = CRPM_st_in(5,4)
CBIO_iv = CBIO_st_in(5,4)
CHUM_iv = CHUM_st_in(5,4)
CIOM_iv = CIOM_st_in(5,4)
C_s1_iv = C_s1_st_in(5,4)
C_s2_iv = C_s2_st_in(5,4)
Csoil_iv = Csoil_st_in(5,4)
Csoilb_iv = Csoilb_st_in(5,4)
end select
case('Arkha')
CDPM_iv = CDPM_st_in(6,1)
......@@ -353,6 +320,69 @@ contains
Csoil_iv = Csoil_st_in(6,1)
Csoilb_iv = Csoilb_st_in(6,1)
case('Samara')
select case(station_opt)
case('1')
CDPM_iv = CDPM_st_in(7,1)
CRPM_iv = CRPM_st_in(7,1)
CBIO_iv = CBIO_st_in(7,1)
CHUM_iv = CHUM_st_in(7,1)
CIOM_iv = CIOM_st_in(7,1)
C_s1_iv = C_s1_st_in(7,1)
C_s2_iv = C_s2_st_in(7,1)
Csoil_iv = Csoil_st_in(7,1)
Csoilb_iv = Csoilb_st_in(7,1)
case('2')
CDPM_iv = CDPM_st_in(7,2)
CRPM_iv = CRPM_st_in(7,2)
CBIO_iv = CBIO_st_in(7,2)
CHUM_iv = CHUM_st_in(7,2)
CIOM_iv = CIOM_st_in(7,2)
C_s1_iv = C_s1_st_in(7,2)
C_s2_iv = C_s2_st_in(7,2)
Csoil_iv = Csoil_st_in(7,2)
Csoilb_iv = Csoilb_st_in(7,2)
case('3')
CDPM_iv = CDPM_st_in(7,3)
CRPM_iv = CRPM_st_in(7,3)
CBIO_iv = CBIO_st_in(7,3)
CHUM_iv = CHUM_st_in(7,3)
CIOM_iv = CIOM_st_in(7,3)
C_s1_iv = C_s1_st_in(7,3)
C_s2_iv = C_s2_st_in(7,3)
Csoil_iv = Csoil_st_in(7,3)
Csoilb_iv = Csoilb_st_in(7,3)
case('4')
CDPM_iv = CDPM_st_in(7,4)
CRPM_iv = CRPM_st_in(7,4)
CBIO_iv = CBIO_st_in(7,4)
CHUM_iv = CHUM_st_in(7,4)
CIOM_iv = CIOM_st_in(7,4)
C_s1_iv = C_s1_st_in(7,4)
C_s2_iv = C_s2_st_in(7,4)
Csoil_iv = Csoil_st_in(7,4)
Csoilb_iv = Csoilb_st_in(7,4)
case('5')
CDPM_iv = CDPM_st_in(7,5)
CRPM_iv = CRPM_st_in(7,5)
CBIO_iv = CBIO_st_in(7,5)
CHUM_iv = CHUM_st_in(7,5)
CIOM_iv = CIOM_st_in(7,5)
C_s1_iv = C_s1_st_in(7,5)
C_s2_iv = C_s2_st_in(7,5)
Csoil_iv = Csoil_st_in(7,5)
Csoilb_iv = Csoilb_st_in(7,5)
end select
end select
end if
......
source/config.f90
+ 22
− 12
View file @ 7e6c59e5
......@@ -59,6 +59,8 @@ module config
character(len_default) :: datetime_last_5
character(len_default) :: datetime_init_6
character(len_default) :: datetime_last_6
character(len_default) :: datetime_init_7
character(len_default) :: datetime_last_7
real :: CDPM_in
real :: CRPM_in
......@@ -70,15 +72,15 @@ module config
real :: C_s1_in
real :: C_s2_in
real :: CDPM_st_in(6,4)
real :: CRPM_st_in(6,4)
real :: CBIO_st_in(6,4)
real :: CHUM_st_in(6,4)
real :: CIOM_st_in(6,4)
real :: Csoil_st_in(6,4)
real :: Csoilb_st_in(6,4)
real :: C_s1_st_in(6,4)
real :: C_s2_st_in(6,4)
real :: CDPM_st_in(7,5)
real :: CRPM_st_in(7,5)
real :: CBIO_st_in(7,5)
real :: CHUM_st_in(7,5)
real :: CIOM_st_in(7,5)
real :: Csoil_st_in(7,5)
real :: Csoilb_st_in(7,5)
real :: C_s1_st_in(7,5)
real :: C_s2_st_in(7,5)
namelist /carbon_pools_config_namelist/ &
& CDPM_in, &
......@@ -135,6 +137,8 @@ module config
& datetime_last_5, &
& datetime_init_6, &
& datetime_last_6, &
& datetime_init_7, &
& datetime_last_7, &
& CDPM_st_in, &
& CRPM_st_in, &
& CBIO_st_in, &
......@@ -242,13 +246,13 @@ module config
end select
select case(station_name)
case('Rostov', 'DAO3', 'DAO4', 'VLDMR', 'TRGK', 'Arkha')
case('Rostov', 'DAO3', 'DAO4', 'VLDMR', 'TRGK', 'Arkha', 'Samara')
case default
stop "check failed : unknown station name"
end select
select case(station_opt)
case('1', '2','3', '4')
case('1', '2','3', '4', '5')
case default
stop "check failed : unknown station_opt"
end select
......@@ -304,6 +308,12 @@ module config
stop "check failed : this station name can't work with this station_opt"
endif
if ((station_name == 'Rostov' .or. station_name == 'VLDMR' .or. station_name == 'TRGK' .or. station_name == 'DAO3' &
& .or. station_name == 'DAO4' .or. station_name == 'Arkha') .and. station_opt == '5') then
stop "check failed : this station name can't work with this station_opt"
endif
if (station_name == 'Arkha' .and. (station_opt == '2' .or. station_opt == '3' .or. station_opt == '4')) then
stop "check failed : this station name can't work with this station_opt"
endif
......
source/environment/environment_data_station.f90
+ 11
− 9
View file @ 7e6c59e5
......@@ -20,13 +20,13 @@ module environment_data_station
! ------- Serve value -------
integer :: mnc !< Количество месяцев в расчете
integer :: mncX = 0 !< Номер месяца с начала работы программы
integer, parameter :: station_max = 6 !< Максимальное число станций
integer, parameter :: opt_max = 4 !< Максимальное количество способов подачи удобрения
integer, parameter :: station_max = 7 !< Максимальное число станций
integer, parameter :: opt_max = 5 !< Максимальное количество способов подачи удобрения
integer :: opt_n !< Номер способа подачи удобрения
integer :: station_n ! DAO4 DAO3 TRGK VLDR ROST !< Номер станции наблюдения за климатом
! 1 2 3 4 5 6
integer :: mnclot_fst(station_max) = (/1935,1937,1956,1968,1975,2014/) !< Даты начала сбора данных по климату, в рамках станции
integer :: mnclot_lst(station_max) = (/2012,2012,2018,2018,2018,2024/) !< Даты конца сбора данных по климату, в рамках станции
! 1 2 3 4 5 6 7
integer :: mnclot_fst(station_max) = (/1935,1937,1956,1968,1975,2014,2014/) !< Даты начала сбора данных по климату, в рамках станции
integer :: mnclot_lst(station_max) = (/2012,2012,2018,2018,2018,2024,2024/) !< Даты конца сбора данных по климату, в рамках станции
! ---------------------------------------------------------------------------------------------------------------------
! ------ Determined externally ------
real(kind = 16), allocatable :: in_temp (:) !< Поступление извне данных cредней температуры воздуха в месяц, [Celsius]
......@@ -39,13 +39,13 @@ module environment_data_station
real, allocatable :: in_rainfall(:) !< Поступление извне данных влажности почвы, [mm]
real, allocatable :: in_Evpar(:) !< Поступление извне данных влажности почвы, [mm]
! station = DAO4 DAO3 TRGK VLDR ROST Arkha
! station = DAO4 DAO3 TRGK VLDR ROST Arkha Samara
! ---- For INMCM ----
real(kind = 16) :: rhodry_in(station_max) = (/1.39, 0.8, 1.4, 1.35, 2.2, 1.05/)
real(kind = 16) :: rhodry_in(station_max) = (/1.39, 0.8, 1.4, 1.35, 2.2, 1.05, 1.3146 /)
!1.29 !1.15
! ---- For Rothc ----
real(kind = 16) :: sw_in (station_max) = (/0.120, 0.120, 0.120, 0.120, 0.120, 0.120/) !0.146
real(kind = 16) :: bettar_in(station_max) = (/20., 26., 7., 8., 44., 21.86/)
real(kind = 16) :: sw_in (station_max) = (/0.120, 0.120, 0.120, 0.120, 0.120, 0.120, 0.120/) !0.146
real(kind = 16) :: bettar_in(station_max) = (/20., 26., 7., 8., 44., 21.86, 49.86/)
contains
......@@ -85,6 +85,8 @@ contains
station_n = 5
case('Arkha')
station_n = 6
case('Samara')
station_n = 7
end select
! External carbon intake
......
source/grid.f90
+ 5
− 0
View file @ 7e6c59e5
......@@ -386,6 +386,8 @@ contains
& datetime_last_5, &
& datetime_init_6, &
& datetime_last_6, &
& datetime_init_7, &
& datetime_last_7, &
& station_name
integer :: ncid, dimid, varid
......@@ -427,6 +429,9 @@ contains
case('Arkha')
DI = datetime_init_6
DL = datetime_last_6
case('Samara')
DI = datetime_init_7
DL = datetime_last_7
end select
case default
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment