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

Update for Rothc, and SOCS model climate type "VLDMR" added

parent f1f6069b
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 75 additions and 74 deletions
data/obs_data_DAO3.csv
+ 8
− 8
View file @ c77d5cec
date,Control,NPK1.5+FYM,NPK3,NPK3+FYM
15-05-1939,2.838,2.8638,3.0444,3.0014
01-05-1962,2.881,3.3798,3.1046,3.2336
10-05-1966,2.7864,3.2336,3.01215,3.6249
12-05-1970,2.35425,2.9799,2.709,3.3411
08-05-1973,2.40585,3.1218,2.8155,3.32175
01-05-1976,2.451,3.14115,2.869089,3.326265
01-05-2006,2.734671,3.381864,3.3580137,3.400569
01-05-2009,2.701002,2.734671,3.32949,3.194814
1939-05-15,2.838,2.8638,3.0444,3.0014
1962-05-01,2.881,3.3798,3.1046,3.2336
1966-05-10,2.7864,3.2336,3.01215,3.6249
1970-05-12,2.35425,2.9799,2.709,3.3411
1973-05-08,2.40585,3.1218,2.8155,3.32175
1976-05-01,2.451,3.14115,2.869089,3.326265
2006-05-01,2.734671,3.381864,3.3580137,3.400569
2009-05-01,2.701002,2.734671,3.32949,3.194814
data/obs_data_DAO4.csv
+ 8
− 8
View file @ c77d5cec
date,Control,NPK1.5+FYM,NPK3,NPK3+FYM
10-05-1937,3.0997,3.2387,3.2248,3.0441
10-05-1949,2.65026,3.07653,2.5854,2.75
10-05-1952,2.53906,2.919,2.5298,2.43713
10-05-1958,2.53906,3.06726,2.63173,2.5298
10-05-1966,2.10353,2.8726,2.2796,2.34446
10-05-1972,2.09426,2.8356,2.2796,2.21473
10-05-1977,2.0294,2.7336,2.24253,2.15913
10-05-2009,2.74108,2.910382,2.90232,2.9399426
1937-05-10,3.0997,3.2387,3.2248,3.0441
1949-05-10,2.65026,3.07653,2.5854,2.75
1952-05-10,2.53906,2.919,2.5298,2.43713
1958-05-10,2.53906,3.06726,2.63173,2.5298
1966-05-10,2.10353,2.8726,2.2796,2.34446
1972-05-10,2.09426,2.8356,2.2796,2.21473
1977-05-10,2.0294,2.7336,2.24253,2.15913
2009-05-10,2.74108,2.910382,2.90232,2.9399426
data/obs_data_VLDMR.csv
+ 11
− 11
View file @ c77d5cec
date,1,10,16
16-05-1968,1.647,1.647,1.647
28-08-1971,1.539,1.674,1.728
18-09-1979,1.566,1.647,1.701
20-10-1983,1.485,1.485,1.728
10-05-1990,1.242,1.323,1.539
10-09-1993,1.323,1.485,1.809
15-09-1997,1.35,1.485,1.782
28-09-2001,1.404,1.593,1.944
05-10-2005,1.431,1.62,2.025
25-09-2009,1.458,1.566,1.944
28-08-2013,1.323,1.62,2.025
1968-05-16,1.647,1.647,1.647
1971-08-28,1.539,1.674,1.728
1979-09-18,1.566,1.647,1.701
1983-10-20,1.485,1.485,1.728
1990-05-10,1.242,1.323,1.539
1993-09-10,1.323,1.485,1.809
1997-09-15,1.35,1.485,1.782
2001-09-28,1.404,1.593,1.944
2005-10-05,1.431,1.62,2.025
2009-09-25,1.458,1.566,1.944
2013-08-28,1.323,1.62,2.025
plot.py
+ 2
− 1
View file @ c77d5cec
......@@ -51,7 +51,8 @@ if (carbon_model_type == 'rothc'):
pools = {'CDPM': 'Разлагаемый растительный материал',
'CRPM': 'Устойчивый растительный материал',
'CBIO': 'Пул микробной биомассы',
'CHUM': 'Долгоживущий гумусовый пул'
'CHUM': 'Долгоживущий гумусовый пул',
'CIOM': 'Пул инертного органического вещества'
}
elif (carbon_model_type == 'socs'):
pools = {'csoil1': 'Огранический углерод в почве',
......
source/carbon/carbon_models/carbon_model_ROTHC.f90
+ 23
− 23
View file @ c77d5cec
......@@ -13,35 +13,35 @@ module carbon_model_rothc
subroutine carbon_model_assembly()
integer :: n_Catm, n_Cveg, n_CDPM, n_CRPM, n_CBIO, n_CHUM !< id пулов
integer :: n_Catm, n_Cveg, n_CDPM, n_CRPM, n_CBIO, n_CHUM, n_CIOM !< id пулов
integer :: n_F !< id потока
call set_tiles(ntiles)
call set_pool(pid = n_Catm, name = 'catm')
call set_pool(pid = n_Cveg, name = 'cveg')
call set_pool(pid = n_CDPM, name = 'CDPM', initial_value = 0.025581370134274, alias = CDPM)
call set_pool(pid = n_CRPM, name = 'CRPM', initial_value = 0.4060806529106 , alias = CRPM)
call set_pool(pid = n_CBIO, name = 'CBIO', initial_value = 0.06760704463772 , alias = CBIO)
call set_pool(pid = n_CHUM, name = 'CHUM', initial_value = 2.51250900729 , alias = CHUM)
! station: opt: CDPM: CRPM: CBIO: CHUM:
! ROST 1 0.043915506788352 0.57587725621248 0.20176158117888 8.11102631424
! VLDR 1 0.0349378381755 0.572185167615 0.0250107214104 1.028654864073
! TRGK 1 0.0987695862776 1.239004653124 0.046724721204 1.812566071456
! 2 0.088876707141306 1.11490447467219 0.04204471760199 1.63101729974136
! 3 0.0935251834664 1.173216795436 0.044243762556 1.716323625184
! DAO3 1 0.02418868809258 0.383973110202 0.0639264318924 2.3757248493
! 2 0.024408585257058 0.3874637748402 0.06450758127324 2.39732234793
! 3 0.025947865408404 0.4118984273076 0.06857562693912 2.54850483834
! 4 0.025581370134274 0.4060806529106 0.06760704463772 2.51250900729
! DAO4 1 0.012403056043085 0.3965742789296 0.06819309095991 2.6225268288181
! 2 0.012959246897035 0.4143578788816 0.07125107710161 2.7401289287651
! 3 0.01290362781164 0.4125795188864 0.07094527848744 2.7283687187704
! 4 0.012180579701505 0.3894608389488 0.06696989650323 2.5754859888393
call set_pool(pid = n_CDPM, name = 'CDPM', initial_value = 0.0196, alias = CDPM)
call set_pool(pid = n_CRPM, name = 'CRPM', initial_value = 1.08615, alias = CRPM)
call set_pool(pid = n_CBIO, name = 'CBIO', initial_value = 0.05455, alias = CBIO)
call set_pool(pid = n_CHUM, name = 'CHUM', initial_value = 1.46279, alias = CHUM)
call set_pool(pid = n_CIOM, name = 'CIOM', initial_value = 0.224, alias = CIOM)
! station: opt: CDPM: CRPM: CBIO: CHUM: CIOM:
! DAO4 1 0.01944 1.73544 0.06188 1.0416 0.245
! 2 0.02037 1.75515 0.06464 1.14567 0.257
! 3 0.02037 1.75515 0.06464 1.14567 0.257
! 4 0.0192 1.73087 0.0612 1.01558 0.24
! DAO3 1 0.01792 1.10183 0.05705 1.45947 0.221
! 2 0.01792 1.10183 0.05705 1.45947 0.224
! 3 0.01804 1.10426 0.05741 1.47248 0.225
! 4 0.01792 1.10183 0.05705 1.45947 0.223
! TRGK 1 0.02132 1.12986 0.05946 1.64844 0.238
! 2 0.0196 1.08615 0.05455 1.46279 0.224
! 3 0.02132 1.12986 0.05946 1.64844 0.238
! VLDR 1 0.00777 1.38303 0.03275 0.14961 0.119
! ROST 1 0.04936 0.799 0.11873 7.07025 0.814
! litterfall
call set_flux(fid = n_F, pid_out = n_Cveg, pid_in = n_CDPM, name = 'fdpm*organic')
......
source/carbon/carbon_models/carbon_model_ROTHC_aux.f90
+ 2
− 1
View file @ c77d5cec
......@@ -16,6 +16,7 @@ module carbon_model_rothc_aux
real, dimension(:,:,:), pointer :: CRPM !< устойчивый растительный материал / resistant plant material
real, dimension(:,:,:), pointer :: CBIO !< пул микробной биомассы / microbial biomass pool
real, dimension(:,:,:), pointer :: CHUM !< долгоживущий гумусовый пул / long lived humified pool
real, dimension(:,:,:), pointer :: CIOM !< инертное органическое вещество / inert organic matter
! ------- Flows -------
real, dimension(:,:,:), pointer :: litterfall1 !< Между растениями и почвой DPM
real, dimension(:,:,:), pointer :: litterfall2 !< Между растениеми и почвой RPM
......@@ -203,7 +204,7 @@ contains
use grid, only : date_c
open(unit=500, file='results/'//trim(carbon_model_type)//'/'//trim(station_name)//'_'//trim(station_opt)//'.txt', status='unknown')
write(500,*) date_c%timestamp,';',CDPM(:,:,1),';',CRPM(:,:,1),';',CBIO(:,:,1),';',CHUM(:,:,1)
write(500,*) date_c%timestamp,';',CDPM(:,:,1),';',CRPM(:,:,1),';',CBIO(:,:,1),';',CHUM(:,:,1),';', CIOM(:,:,1)
! ---- Initialization ------
!print*, 'finish thing 1, ', CDPM(:,:,1)/cs*100,'%'
......
source/carbon/carbon_models/carbon_model_SOCS.f90
+ 11
− 12
View file @ c77d5cec
......@@ -25,24 +25,23 @@ module carbon_model_socs
call set_pool(pid = n_Catm, name = 'catm')
call set_pool(pid = n_Cveg, name = 'cveg')
call set_pool(pid = n_Csoil1, name = 'csoil1',initial_value = 0.178176, alias = C1)
call set_pool(pid = n_Csoil2, name = 'csoil2',initial_value = 8.730624, alias = C2)
call set_pool(pid = n_Csoil1, name = 'csoil1',initial_value = 0.060028, alias = C1)
call set_pool(pid = n_Csoil2, name = 'csoil2',initial_value = 2.941372, alias = C2)
! station: opt: С1: С2:
! ROST 1 0.178176 8.730624
! VLDR 1 0.0338432 1.6583168
! TRGK 1 0.0619416 3.0351384
! 2 0.0569418 2.7901482
! 3 0.0619416 3.0351384
! DAO3 1 0.05676 2.78124
! 2 0.057276 2.806524
! 3 0.060888 2.983512
! 4 0.060028 2.941372
! DAO4 1 0.061994 3.037706
! 2 0.064774 3.173926
! 3 0.064496 3.160304
! 4 0.060882 2.983218
! DAO3 1 0.05676 2.78124
! 2 0.057276 2.806524
! 3 0.060888 2.983512
! 4 0.060028 2.941372
! TRGK 1 0.0619416 3.0351384
! 2 0.0569418 2.7901482
! 3 0.0619416 3.0351384
! VLDR 1 0.0338432 1.6583168
! ROST 1 0.178176 8.730624
call set_flux(fid = n_F, pid_out = n_Cveg, pid_in = n_Csoil1, name = 'litterfall')
call set_mult(n_F, 'lin', x_ij = lambd)
......
source/carbon/carbon_models/carbon_model_SOCS_aux.f90
+ 1
− 1
View file @ c77d5cec
......@@ -90,7 +90,7 @@ contains
end if
!kirk = 7.5/yrs
kirk = Ft*Fs*7.5/yrs
kirk = Ft*Fs*7.5/yrs !7.5
!print*, kirk
end subroutine
......
source/environment/environment_data_station.f90
+ 5
− 5
View file @ c77d5cec
......@@ -43,16 +43,16 @@ module environment_data_station
real(kind = 16) :: rhodry_in(station_max) = (/1.39, 1.29, 1.4, 1.35, 2.2/)
! 1.15
! ---- 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.120/) !0.146
real(kind = 16) :: bettar_in(station_max) = (/20., 26., 7., 8., 44./)
! ---- For SOCs ----
real :: kd_in (station_max,opt_max)
! opt = 1 2 3 4
data kd_in(01,:) /0.03, 0.03, 0.035, 0.03/
data kd_in(02,:) /0.03, 0.03, 0.035, 0.03/
data kd_in(02,:) /0.03, 0.03, 0.035, 0.03/ !0.08581933
data kd_in(03,:) /0.015, 0.008, 0.008, 0./ !0.015, 0.008, 0.008, 0.
data kd_in(04,:) /0.015, 0.008, 0.008, 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./
contains
......
ui1_config.nml
+ 4
− 4
View file @ c77d5cec
......@@ -12,7 +12,7 @@
! -----------------------------------------------------------------
!> Углеродная модель:
carbon_model_type = 'rothc'
carbon_model_type = 'socs'
!
! 'inmcm' - модель inmc [1,2]
! 'socs' - модель SOCS [3]
......@@ -34,11 +34,11 @@
! дополнительно для environment_data_type = 'station':
!> Имя станции:
station_name = 'DAO3'
! 'Rostov' - Станция ФАНЦ
! 'DAO3' - Станция в Долгопрудном 1
! 'DAO4' - Станция в Долгопрудном 2
! 'VLDMR' - Данные Владимир
! 'DAO3' - Станция в Долгопрудном 1
! 'TRGK' - Данные Торжок
! 'VLDMR' - Данные Владимир
! 'Rostov' - Станция ФАНЦ
!
!> Тип подачи удобрения:
station_opt = '4'
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment