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

Correct version of new station Arkha

parent 943eb2e7
No related branches found
No related tags found
No related merge requests found
date,1
2014-01-01,4.75
2024-01-01,4.75
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
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
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
0.04925
\ No newline at end of file
0.3461875875
0.339229575
0.3659286525
0.3923759475
0.36354065
0.3393821725
0.3203887925
0.281555175
0.3216056825
0.3275375375
0.3486785925
0.32446671
0.31612396
0.3125381475
0.33184052
0.3681602525
0.38012313
0.3573417625
0.3455085775
0.3688049325
0.3767852775
0.38108063
0.371978755
0.3903923025
0.35346222
0.34226227
0.3509140075
0.40509033
0.3183326725
0.3044891375
0.296989435
0.33869934
0.368946075
0.3546180725
0.3441467275
0.33811188
0.3270492575
0.32263565
0.339263915
0.38507843
0.4077033975
0.3770942675
0.370227805
0.3640785175
0.3773078925
0.3861923175
0.3792114175
0.347515115
0.336383825
0.325767515
0.3180084325
0.365818025
0.38244629
0.329792025
0.315387735
0.271030435
0.3328552275
0.390029905
0.380268095
0.3490753225
0.33719635
0.33073043
0.32922744
0.3878746025
0.3603820825
0.3047904925
0.34820556
0.3729248075
0.363758085
0.4001121525
0.37821579
0.3751297025
0.36209488
0.355361945
0.3701629575
0.39277267
0.39996338
0.28543855
0.33232498
0.3464698775
0.3870697
0.38690186
0.3791275
0.351459505
0.3406028725
0.3336448675
0.330230715
0.4027099625
0.37491226
0.279994965
0.272499075
0.3793602
0.3768882725
0.3718337975
0.3790702775
0.34532547
0.3346328725
0.326454165
0.321262355
0.381240855
0.398002625
0.3615112325
0.313655855
0.27645874
0.35876846
0.36804962
0.375492095
0.3489723175
0.339355465
0.3336105325
0.3570709225
0.3972969025
0.3349876425
0.271015165
0.3349075275
0.313903805
0.325611115
0.3832321175
0.3794441175
0.3455734225
0.335918425
0.3293304475
0.3588943475
0.4143753075
0.3531532325
0.3122100875
0.24771118
0.2462806675
0.2103118925
0.28814697
0.3482742375
0.3442497275
\ No newline at end of file
-12.65085
-6.25689999999997
-2.57882999999998
0.751500000000021
10.17928
12.87832
16.31487
16.57302
10.0523
-1.86135999999999
-2.68405999999999
-6.58866
-11.44297
-5.34024999999997
-1.52238
1.23934000000003
12.12936
14.57406
13.67397
12.89504
10.94723
1.61850000000004
-3.52329999999995
-4.25686999999999
-15.96543
-3.78229999999996
-3.59269999999998
3.78310000000005
12.38516
14.2464000000001
19.65164
16.42135
8.91586000000001
2.02746000000002
-7.48782999999997
-7.69421999999997
-11.1268
-9.19247999999999
-1.52139999999997
-0.293129999999962
3.39297000000005
10.84304
16.762
15.8539
8.54653000000002
2.26336000000003
-2.16684999999995
-4.26609999999999
-6.57687999999996
-13.2824
-10.2313
0.994100000000003
9.9941
13.21978
18.46603
16.16628
10.79803
3.16445000000004
-1.74149999999997
-8.82675999999998
-10.06156
-6.28729999999996
-3.06099999999998
1.68280000000004
10.2193
14.7729
13.70654
11.34475
8.91030000000001
2.10420000000005
-3.72573999999997
-2.80084999999997
-4.94409999999999
-4.31809999999996
-1.90832999999998
0.0785200000000259
7.12286
14.60885
16.9529700000001
13.4312
10.2524
4.23373000000004
-0.520969999999977
-7.49612999999999
-12.81925
-19.71625
-5.66854999999998
3.29165
10.68185
19.01302
18.3948600000001
14.74502
6.83904000000001
3.90005000000002
-2.1859
-13.07255
-10.1679999999999
-5.75525999999996
-5.34889999999996
0.214800000000025
6.16116
14.62795
19.0662
18.43566
7.62945000000002
4.46270000000004
-4.56849999999997
-8.72985999999997
-9.34262999999999
-8.07846999999998
-4.27392999999995
2.07742000000002
11.13778
13.00094
16.09918
16.25643
12.83083
1.77987000000002
-6.29935
-11.36738
-14.00175
-10.09025
-1.68959999999998
2.38937000000004
6.43434999999999
17.13674
18.83737
15.98525
14.6584
3.97106000000002
-0.85979999999995
-5.78695999999997
\ No newline at end of file
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
\ No newline at end of file
......@@ -43,7 +43,7 @@ module carbon_model_rothc_aux
! ------- Coefficients ------------
real :: ks(ntiles) = (/ 3.22*(1.E-7), 9.65*(1.E-9), 2.12*(1.E-8), 6.43*(1.E-10)/)
!< Скорость дыхания единицы массы каждого пула в стандартных условиях, [1-s]
real, parameter :: alphadr = 1.44 !< Определяет соотношение поступления опада между DPM и RPM, [dim],
real, parameter :: alphadr = 0.25 !< Определяет соотношение поступления опада между DPM и RPM, [dim],
! Для деревьев, для кустарников, для натуральной травы, для СХ (/ 0.25, 0.33, 0.67, 1.44 /)
real, allocatable :: SA(:,:) !< Определяет распределение выдыхаемой части пула по остальным пулам
......@@ -181,9 +181,8 @@ contains
use grid, only : date_c
!open(1, file='results/'//trim(carbon_model_type)//'/'//trim(station_name)//'_'//trim(station_opt)//'.txt', status='unknown')
!write(1,*) date_c%timestamp,';',CDPM(:,:,1),';',CRPM(:,:,1),';',CBIO(:,:,1),';',CHUM(:,:,1),';', CIOM(:,:,1)
!close(1)
open(1, file='results/'//trim(carbon_model_type)//'/'//trim(station_name)//'_'//trim(station_opt)//'.txt', status='unknown')
write(1,*) date_c%timestamp,';',CDPM(:,:,1),';',CRPM(:,:,1),';',CBIO(:,:,1),';',CHUM(:,:,1),';', CIOM(:,:,1)
! ---- Initialization ------
!print*, 'finish thing 1, ', CDPM(:,:,1)/cs*100,'%'
......
......@@ -83,7 +83,9 @@ contains
if (environment_data_type == 'station') then
select case(station_name)
case('Arkha')
case('DAO4')
select case(station_opt)
case('1')
CDPM_iv = CDPM_st_in(1,1)
CRPM_iv = CRPM_st_in(1,1)
CBIO_iv = CBIO_st_in(1,1)
......@@ -95,6 +97,43 @@ contains
Csoil_iv = Csoil_st_in(1,1)
Csoilb_iv = Csoilb_st_in(1,1)
case('2')
CDPM_iv = CDPM_st_in(1,2)
CRPM_iv = CRPM_st_in(1,2)
CBIO_iv = CBIO_st_in(1,2)
CHUM_iv = CHUM_st_in(1,2)
CIOM_iv = CIOM_st_in(1,2)
C_s1_iv = C_s1_st_in(1,2)
C_s2_iv = C_s2_st_in(1,2)
Csoil_iv = Csoil_st_in(1,2)
Csoilb_iv = Csoilb_st_in(1,2)
case('3')
CDPM_iv = CDPM_st_in(1,3)
CRPM_iv = CRPM_st_in(1,3)
CBIO_iv = CBIO_st_in(1,3)
CHUM_iv = CHUM_st_in(1,3)
CIOM_iv = CIOM_st_in(1,3)
C_s1_iv = C_s1_st_in(1,3)
C_s2_iv = C_s2_st_in(1,3)
Csoil_iv = Csoil_st_in(1,3)
Csoilb_iv = Csoilb_st_in(1,3)
case('4')
CDPM_iv = CDPM_st_in(1,4)
CRPM_iv = CRPM_st_in(1,4)
CBIO_iv = CBIO_st_in(1,4)
CHUM_iv = CHUM_st_in(1,4)
CIOM_iv = CIOM_st_in(1,4)
C_s1_iv = C_s1_st_in(1,4)
C_s2_iv = C_s2_st_in(1,4)
Csoil_iv = Csoil_st_in(1,4)
Csoilb_iv = Csoilb_st_in(1,4)
end select
case('DAO3')
select case(station_opt)
......@@ -302,6 +341,18 @@ contains
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)
CRPM_iv = CRPM_st_in(6,1)
CBIO_iv = CBIO_st_in(6,1)
CHUM_iv = CHUM_st_in(6,1)
CIOM_iv = CIOM_st_in(6,1)
C_s1_iv = C_s1_st_in(6,1)
C_s2_iv = C_s2_st_in(6,1)
Csoil_iv = Csoil_st_in(6,1)
Csoilb_iv = Csoilb_st_in(6,1)
end select
end if
......
......@@ -56,6 +56,8 @@ module config
character(len_default) :: datetime_last_4
character(len_default) :: datetime_init_5
character(len_default) :: datetime_last_5
character(len_default) :: datetime_init_6
character(len_default) :: datetime_last_6
real :: CDPM_in
real :: CRPM_in
......@@ -67,15 +69,15 @@ module config
real :: C_s1_in
real :: C_s2_in
real :: CDPM_st_in(5,4)
real :: CRPM_st_in(5,4)
real :: CBIO_st_in(5,4)
real :: CHUM_st_in(5,4)
real :: CIOM_st_in(5,4)
real :: Csoil_st_in(5,4)
real :: Csoilb_st_in(5,4)
real :: C_s1_st_in(5,4)
real :: C_s2_st_in(5,4)
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)
namelist /carbon_pools_config_namelist/ &
& CDPM_in, &
......@@ -130,6 +132,8 @@ module config
& datetime_last_4, &
& datetime_init_5, &
& datetime_last_5, &
& datetime_init_6, &
& datetime_last_6, &
& CDPM_st_in, &
& CRPM_st_in, &
& CBIO_st_in, &
......@@ -230,7 +234,7 @@ module config
end select
select case(station_name)
case('Rostov', 'DAO3','Arkha', 'VLDMR', 'TRGK')
case('Rostov', 'DAO3', 'DAO4', 'VLDMR', 'TRGK', 'Arkha')
case default
stop "check failed : unknown station name"
end select
......
......@@ -19,13 +19,13 @@ module environment_data_station
! ------- Serve value -------
integer :: mnc !< Количество месяцев в расчете
integer :: mncX = 0 !< Номер месяца с начала работы программы
integer, parameter :: station_max = 5 !< Максимальное число станций
integer, parameter :: station_max = 6 !< Максимальное число станций
integer, parameter :: opt_max = 4 !< Максимальное количество способов подачи удобрения
integer :: opt_n !< Номер способа подачи удобрения
integer :: station_n ! Arkha DAO3 TRGK VLDR ROST !< Номер станции наблюдения за климатом
! 1 2 3 4 5
integer :: mnclot_fst(station_max) = (/2014,1937,1956,1968,1975/) !< Даты начала сбора данных по климату, в рамках станции
integer :: mnclot_lst(station_max) = (/2025,2012,2018,2018,2018/) !< Даты конца сбора данных по климату, в рамках станции
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/) !< Даты конца сбора данных по климату, в рамках станции
! ---------------------------------------------------------------------------------------------------------------------
! ------ Determined externally ------
real(kind = 16), allocatable :: in_temp (:) !< Поступление извне данных cредней температуры воздуха в месяц, [Celsius]
......@@ -38,13 +38,13 @@ module environment_data_station
real, allocatable :: in_rainfall(:) !< Поступление извне данных влажности почвы, [mm]
real, allocatable :: in_Evpar(:) !< Поступление извне данных влажности почвы, [mm]
! station = Arkha DAO3 TRGK VLDR ROST
! station = DAO4 DAO3 TRGK VLDR ROST Arkha
! ---- For INMCM ----
real(kind = 16) :: rhodry_in(station_max) = (/1.05, 0.8, 1.4, 1.35, 2.2/)
real(kind = 16) :: rhodry_in(station_max) = (/1.39, 0.8, 1.4, 1.35, 2.2, 1.05/)
!1.29 !1.15
! ---- For Rothc ----
real(kind = 16) :: sw_in (station_max) = (/0.120, 0.120, 0.120, 0.120, 0.120/) !0.146
real(kind = 16) :: bettar_in(station_max) = (/21.86, 26., 7., 8., 44./)
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/)
! ---- For SOCs ----
real :: kd_in (station_max,opt_max) !< Коэф. перехода углерода из C2 в C1, десорбация, разрушение агрегатов
......@@ -54,6 +54,7 @@ module environment_data_station
data kd_in(03,:) /0.0402, 0.05067, 0.04185, 0./ !0.015, 0.008, 0.008, 0.
data kd_in(04,:) /0.08581933, 0.08581933, 0.08581933, 0./ !0.015, 0.008, 0.008, 0.
data kd_in(05,:) /0.012, 0.007, 0.007, 0./
data kd_in(06,:) /0.012, 0.007, 0.007, 0./
contains
......@@ -81,7 +82,7 @@ contains
read(station_opt,*) opt_n
select case (station_name)
case('Arkha')
case('DAO4')
station_n = 1
case('DAO3')
station_n = 2
......@@ -91,6 +92,8 @@ contains
station_n = 4
case('Rostov')
station_n = 5
case('Arkha')
station_n = 6
end select
! External carbon intake
......
......@@ -384,6 +384,8 @@ contains
& datetime_last_4, &
& datetime_init_5, &
& datetime_last_5, &
& datetime_init_6, &
& datetime_last_6, &
& station_name
integer :: ncid, dimid, varid
......@@ -422,6 +424,9 @@ contains
case('Rostov')
DI = datetime_init_5
DL = datetime_last_5
case('Arkha')
DI = datetime_init_6
DL = datetime_last_6
end select
case default
......
......@@ -11,23 +11,23 @@
!> station_name = 'DAO4':
!> Date
datetime_init_1 = '2014-01-01 00:00:00'
datetime_last_1 = '2025-01-01 00:00:00'
datetime_init_1 = '1935-01-01 00:00:00'
datetime_last_1 = '2012-01-01 00:00:00'
!> Initial value
! RothC
! opt 1
CDPM_st_in(1,:) = 0.017849986
CRPM_st_in(1,:) = 0.570733963
CBIO_st_in(1,:) = 0.09814079
CHUM_st_in(1,:) = 3.77426532
CIOM_st_in(1,:) = 0.289035
! opt 1 2 3 4
CDPM_st_in(1,:) = 0.01944, 0.02037, 0.02037, 0.0192
CRPM_st_in(1,:) = 1.73544, 1.75515, 1.75515, 1.73087
CBIO_st_in(1,:) = 0.06188, 0.06464, 0.06464, 0.0612
CHUM_st_in(1,:) = 1.0416 , 1.14567, 1.14567, 1.01558
CIOM_st_in(1,:) = 0.245 , 0.257 , 0.257 , 0.24
! INMCM
Csoil_st_in(1,:) = 0.095
Csoilb_st_in(1,:) = 4.655
Csoil_st_in(1,:) = 0.061994, 0.064774, 0.064496, 0.060882
Csoilb_st_in(1,:) = 3.037706, 3.173926, 3.160304, 2.983218
! SOCS
C_s1_st_in(1,:) = 4.655
C_s2_st_in(1,:) = 0.095
C_s1_st_in(1,:) = 3.037706, 3.173926, 3.160304, 2.983218
C_s2_st_in(1,:) = 0.061994, 0.064774, 0.064496, 0.060882
!> station_name = 'DAO3':
datetime_init_2 = '1937-01-01 00:00:00'
......@@ -105,5 +105,26 @@
C_s1_st_in(5,:) = 8.730624, 8.730624, 8.730624
C_s2_st_in(5,:) = 0.178176, 0.178176, 0.178176
!> station_name = 'Arkha':
!> Date
datetime_init_6 = '2014-01-01 00:00:00'
datetime_last_6 = '2025-01-01 00:00:00'
!> Initial value
! RothC
! opt 1
CDPM_st_in(6,:) = 0.017849986
CRPM_st_in(6,:) = 0.570733963
CBIO_st_in(6,:) = 0.09814079
CHUM_st_in(6,:) = 3.77426532
CIOM_st_in(6,:) = 0.289035
! INMCM
Csoil_st_in(6,:) = 4.655
Csoilb_st_in(6,:) = 0.095
! SOCS
C_s1_st_in(6,:) = 0.095
C_s2_st_in(6,:) = 4.655
! -----------------------------------------------------------------
/
......@@ -19,7 +19,7 @@
! 'calc' - расчёт согласно настройкам всё как обычно
!> Углеродная модель:
carbon_model_type = 'rothc'
carbon_model_type = 'inmcm'
!
! 'inmcm' - модель inmc [1,2]
! 'socs' - модель SOCS [3]
......@@ -41,11 +41,12 @@
! дополнительно для environment_data_type = 'station':
!> Имя станции:
station_name = 'Arkha'
! 'Arkha' - Долгопрудная опытная станция ДАОС 4
! 'DAO4' - Долгопрудная опытная станция ДАОС 4
! 'DAO3' - Долгопрудная опытная станция ДАОС 3
! 'TRGK' - Федеральный научный центр любяных культур "Торжок"
! 'VLDMR' - Верхневолжский ФАНЦ "Владимир"
! 'Rostov' - Донской Федеральный аграрный научный центр (ФАНЦ) "Ростов"
! 'Arkha' - Комбинация данных на основе [4], [5] и ERA5
!
!> Тип подачи удобрения:
station_opt = '1'
......@@ -107,7 +108,7 @@
! '2629800' - 1 месяц
!> Дата и время на момент старта:
datetime_init_mode = 'manual'
datetime_init_mode = 'auto'
!
! 'auto' - получить из входного файла (для lsm_offline)
! 'auto' - указать автоматически дату начала наблюдения станции
......@@ -117,16 +118,16 @@
ntime_mode = 'datetime_last'
!
! 'auto' - до конца входного файла (для lsm_offline)
! 'auto' - указать автоматически дату конца наблюдения станции
! 'auto' - автоматически указать дату конца наблюдения станции
! 'ntime' - указанное число шагов
ntime = 1000
! 'datetime_last' - до достижения указанной даты
datetime_last = '2175-01-01 00:00:00'
!Номер станции наблюдения за климатом
!Arkha DAO3 TRGK VLDR ROST
!1935,1937,1956,1968,1975
!2012,2012,2018,2018,2018
!DAO4 DAO3 TRGK VLDR ROST Arkha
!1935,1937,1956,1968,1975,2014
!2012,2012,2018,2018,2018,2024
! продвинутые настройки
......@@ -142,7 +143,7 @@
testing_log_mode = 'read' ! none, write, read
nv_singlecolumn = 1 ! Для модели INMCM нужно поставить 12, для RothC и SOCS 1
nv_singlecolumn = 12 ! Для модели INMCM нужно поставить 12, для RothC и SOCS 1
environment_model_type = 'inmcm' ! inmcm
......@@ -159,3 +160,9 @@
!
! [3] Рыжова И.М. Анализ динамики органического вещества почв на основе минимальных моделей круговорота углерода // Почвы-стратегический ресурс России: тезисы докладов VIII cъезда Общества почвоведов им. В.В. Докучаева и Школы молодых ученых по морфологии и классификации почв (Сыктывкар, 2020-2022 гг.). Часть 2. — Москва–Сыктывкар: ИБ ФИЦ Коми НЦ УрО РАН, 2021. — с.130–131.
! url : https://ib.komisc.ru/en/component/phocadownload/category/35-knigi?download=1313:soils-strategic-resource-of-russia-part-2
!
! [4] Почвенная карта России
! url: https://soil-db.ru/map
!
! [5] Продуктивность экосистем
! url: https://soil-db.ru/map?name=ecosystem-productivity&zoom=7
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment