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

Some little fix

parent 1de1b2ad
No related branches found
No related tags found
No related merge requests found
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
...@@ -18,42 +18,32 @@ implicit none ...@@ -18,42 +18,32 @@ implicit none
! ------------------------------------------------------------------------------------------------------------------- ! -------------------------------------------------------------------------------------------------------------------
! Main part of SOCS ! Main part of SOCS
! пулы ! ------- Pools -------
real, dimension(:,:,:), pointer :: Catm !< атмосфера real, dimension(:,:,:), pointer :: Catm !< атмосфера
real, dimension(:,:,:), pointer :: Cveg !< растительность real, dimension(:,:,:), pointer :: Cveg !< растительность
real, dimension(:,:,:), pointer :: C1 !< огранический углерод в почве / soil organic carbon (SOC) real, dimension(:,:,:), pointer :: C1 !< огранический углерод в почве / soil organic carbon (SOC)
real, dimension(:,:,:), pointer :: C2 !< минерализованный углерод в почве / mineral organic carbon (MOC) real, dimension(:,:,:), pointer :: C2 !< минерализованный углерод в почве / mineral organic carbon (MOC)
! потоки ! ------- Flows -------
real, dimension(:,:,:), pointer :: F_microbal_respiration !< Между почвой1 и атмосферой (микробное дыхание) real, dimension(:,:,:), pointer :: F_microbal_respiration !< Между почвой1 и атмосферой (микробное дыхание)
real, dimension(:,:,:), pointer :: F_litterfall !< Между растением и почвой1 (отмирание и опад) real, dimension(:,:,:), pointer :: F_litterfall !< Между растением и почвой1 (отмирание и опад)
real, dimension(:,:,:), pointer :: F_mineralization !< Между почвой1 и почвой2 (минерализация) real, dimension(:,:,:), pointer :: F_mineralization !< Между почвой1 и почвой2 (минерализация)
real, dimension(:,:,:), pointer :: F_destabilization !< Между почвой2 и почвой1 (дестабилизация?) real, dimension(:,:,:), pointer :: F_destabilization !< Между почвой2 и почвой1 (дестабилизация?)
! ------- Coefficients ------------
integer, parameter :: ntiles = 2
! 1 - целина
! 2 - пашня
real :: rare = 0.45 !< Доля углерода переходящего в пул защищенного C1 real :: rare = 0.45 !< Доля углерода переходящего в пул защищенного C1
real :: Cm = 12. !< Max кол-во органического углерода которое может быть защищено в почве real :: Cm = 12. !< Max кол-во органического углерода которое может быть защищено в почве
real :: kirk = 7.5/yrs !< Коэф. скорости разложения C1 !(ntiles) = (/ 7.5/yrs, 4./yrs /) real :: kirk = 7.5/yrs !< Коэф. скорости разложения C1 !(ntiles) = (/ 7.5/yrs, 4./yrs /)
real :: kd !< Коэф. перехода углерода из C2 в C1, десорбация, разрушение агрегатов real :: kd !< Коэф. перехода углерода из C2 в C1, десорбация, разрушение агрегатов
!real :: Irk(ntiles) = (/ 0.8/yrs, 0.2/yrs /) !< Поступление органического углерода в почву
!real :: kd(ntiles) = (/ 0.02/yrs, 0.008/yrs /) !< Коэф. перехода углерода из C2 в C1, десорбация, разрушение агрегатов
! ------ Litterfall -------- ! ------ Litterfall --------
real :: in_lambd(0:mnc) !< Поступление извне данных по поступлению углерода в почву, [kg/m**3 / year] real :: in_lambd(0:mnc) !< Поступление извне данных по поступлению углерода в почву, [kg/m**3 / year]
real :: lambd !< Поступление углерода в почву в интервал времени dt (lambdac) real :: lambd !< Поступление углерода в почву в интервал времени dt (lambdac)
contains contains
subroutine carbon_model_init() subroutine carbon_model_init()
! ---- Part of enviromental carbon_model_init() ---- ! ---- Part of enviromental carbon_model_init() ----
open (unit = 1, file = 'initial_value/'//trim(station)//'_'//trim(opt)//'.txt', status='unknown') open (unit = 1, file = 'initial_value/'//trim(station)//'_'//trim(opt)//'.txt', status='unknown')
read(1,*) in_lambd(1:mnc) read(1,*) in_lambd(1:mnc)
close (1) close (1)
!in_lambd(0:mnc) = 0.8/yrs
select case (station) select case (station)
case('Rostov') case('Rostov')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment