Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Making P code compatible with Century decomposition submodel #371

Merged
merged 2 commits into from
Nov 23, 2015
Merged

Making P code compatible with Century decomposition submodel #371

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
f2f6664
Making P code working with Century decomposition submodel
xyk-ornl Oct 9, 2015
0150e5d
cleanup comment
Oct 17, 2015
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
54 changes: 54 additions & 0 deletions models/lnd/clm/src/biogeochem/CNDecompCascadeBGCMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,6 +45,16 @@ module CNDecompCascadeBGCMod
real(r8):: cn_s2_bgc !C:N for SOM 2
real(r8):: cn_s3_bgc !C:N for SOM 3

real(r8):: np_s1_new_bgc !N:P for SOM 1
real(r8):: np_s2_new_bgc !N:P for SOM 2
real(r8):: np_s3_new_bgc !N:P for SOM 3

real(r8):: cp_s1_new_bgc !C:P for SOM 1
real(r8):: cp_s2_new_bgc !C:P for SOM 2
real(r8):: cp_s3_new_bgc !C:P for SOM 3



real(r8):: rf_l1s1_bgc !respiration fraction litter 1 -> SOM 1
real(r8):: rf_l2s1_bgc
real(r8):: rf_l3s2_bgc
Expand Down Expand Up @@ -148,6 +158,24 @@ subroutine readCNDecompBgcParams ( ncid )
if ( .not. readv ) call endrun(msg=trim(errCode)//trim(tString)//errMsg(__FILE__, __LINE__))
CNDecompBgcParamsInst%cn_s3_bgc=tempr

!!! read in phosphorus variables - note that these NP ratio parameters for BGC will have
!!! to be added in the parameter file

tString='np_s1_new'
call ncd_io(trim(tString),tempr, 'read', ncid, readvar=readv)
if ( .not. readv ) call endrun(msg=trim(errCode)//trim(tString)//errMsg(__FILE__, __LINE__))
CNDecompBgcParamsInst%np_s1_new_bgc=tempr

tString='np_s2_new'
call ncd_io(trim(tString),tempr, 'read', ncid, readvar=readv)
if ( .not. readv ) call endrun(msg=trim(errCode)//trim(tString)//errMsg(__FILE__, __LINE__))
CNDecompBgcParamsInst%np_s2_new_bgc=tempr

tString='np_s3_new'
call ncd_io(trim(tString),tempr, 'read', ncid, readvar=readv)
if ( .not. readv ) call endrun(msg=trim(errCode)//trim(tString)//errMsg(__FILE__, __LINE__))
CNDecompBgcParamsInst%np_s3_new_bgc=tempr

tString='rf_l1s1_bgc'
call ncd_io(trim(tString),tempr, 'read', ncid, readvar=readv)
if ( .not. readv ) call endrun(msg=trim(errCode)//trim(tString)//errMsg(__FILE__, __LINE__))
Expand Down Expand Up @@ -244,6 +272,11 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
real(r8) :: cn_s1
real(r8) :: cn_s2
real(r8) :: cn_s3
real(r8) :: np_s1_new
real(r8) :: np_s2_new
real(r8) :: np_s3_new


!real(r8) :: f_s1s2(bounds%begc:bounds%endc,1:nlevdecomp)
!real(r8) :: f_s1s3(bounds%begc:bounds%endc,1:nlevdecomp)
real(r8), allocatable :: f_s1s2(:,:)
Expand Down Expand Up @@ -281,6 +314,7 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
cascade_donor_pool => decomp_cascade_con%cascade_donor_pool , & ! Output: [integer (:) ] which pool is C taken from for a given decomposition step
cascade_receiver_pool => decomp_cascade_con%cascade_receiver_pool , & ! Output: [integer (:) ] which pool is C added to for a given decomposition step
floating_cn_ratio_decomp_pools => decomp_cascade_con%floating_cn_ratio_decomp_pools , & ! Output: [logical (:) ] TRUE => pool has fixed C:N ratio
floating_cp_ratio_decomp_pools => decomp_cascade_con%floating_cp_ratio_decomp_pools , & ! Output: [logical (:) ] TRUE => pool has fixed C:P ratio
decomp_pool_name_restart => decomp_cascade_con%decomp_pool_name_restart , & ! Output: [character(len=8) (:) ] name of pool for restart files
decomp_pool_name_history => decomp_cascade_con%decomp_pool_name_history , & ! Output: [character(len=8) (:) ] name of pool for history files
decomp_pool_name_long => decomp_cascade_con%decomp_pool_name_long , & ! Output: [character(len=20) (:) ] name of pool for netcdf long names
Expand All @@ -289,6 +323,7 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil => decomp_cascade_con%is_soil , & ! Output: [logical (:) ] TRUE => pool is a soil pool
is_cwd => decomp_cascade_con%is_cwd , & ! Output: [logical (:) ] TRUE => pool is a cwd pool
initial_cn_ratio => decomp_cascade_con%initial_cn_ratio , & ! Output: [real(r8) (:) ] c:n ratio for initialization of pools
initial_cp_ratio => decomp_cascade_con%initial_cp_ratio , & ! Output: [real(r8) (:) ] c:p ratio for initialization of pools
initial_stock => decomp_cascade_con%initial_stock , & ! Output: [real(r8) (:) ] initial concentration for seeding at spinup
is_metabolic => decomp_cascade_con%is_metabolic , & ! Output: [logical (:) ] TRUE => pool is metabolic material
is_cellulose => decomp_cascade_con%is_cellulose , & ! Output: [logical (:) ] TRUE => pool is cellulose
Expand All @@ -307,6 +342,11 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
cn_s2 = CNDecompBgcParamsInst%cn_s2_bgc
cn_s3 = CNDecompBgcParamsInst%cn_s3_bgc

! set soil organic matter C:P ratios -X. YANG
np_s1_new=CNDecompBgcParamsInst%np_s1_new_bgc
np_s2_new=CNDecompBgcParamsInst%np_s2_new_bgc
np_s3_new=CNDecompBgcParamsInst%np_s3_new_bgc

! set respiration fractions for fluxes between compartments
rf_l1s1 = CNDecompBgcParamsInst%rf_l1s1_bgc
rf_l2s1 = CNDecompBgcParamsInst%rf_l2s1_bgc
Expand Down Expand Up @@ -340,6 +380,7 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
!------------------- list of pools and their attributes ------------
i_litr1 = i_met_lit
floating_cn_ratio_decomp_pools(i_litr1) = .true.
floating_cp_ratio_decomp_pools(i_litr1) = .true.
decomp_pool_name_restart(i_litr1) = 'litr1'
decomp_pool_name_history(i_litr1) = 'LITR1'
decomp_pool_name_long(i_litr1) = 'litter 1'
Expand All @@ -348,13 +389,15 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil(i_litr1) = .false.
is_cwd(i_litr1) = .false.
initial_cn_ratio(i_litr1) = 90._r8
initial_cp_ratio(i_litr1) = 900._r8
initial_stock(i_litr1) = 0._r8
is_metabolic(i_litr1) = .true.
is_cellulose(i_litr1) = .false.
is_lignin(i_litr1) = .false.

i_litr2 = i_cel_lit
floating_cn_ratio_decomp_pools(i_litr2) = .true.
floating_cp_ratio_decomp_pools(i_litr2) = .true.
decomp_pool_name_restart(i_litr2) = 'litr2'
decomp_pool_name_history(i_litr2) = 'LITR2'
decomp_pool_name_long(i_litr2) = 'litter 2'
Expand All @@ -363,13 +406,15 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil(i_litr2) = .false.
is_cwd(i_litr2) = .false.
initial_cn_ratio(i_litr2) = 90._r8
initial_cp_ratio(i_litr2) = 900._r8
initial_stock(i_litr2) = 0._r8
is_metabolic(i_litr2) = .false.
is_cellulose(i_litr2) = .true.
is_lignin(i_litr2) = .false.

i_litr3 = i_lig_lit
floating_cn_ratio_decomp_pools(i_litr3) = .true.
floating_cp_ratio_decomp_pools(i_litr3) = .true.
decomp_pool_name_restart(i_litr3) = 'litr3'
decomp_pool_name_history(i_litr3) = 'LITR3'
decomp_pool_name_long(i_litr3) = 'litter 3'
Expand All @@ -378,13 +423,15 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil(i_litr3) = .false.
is_cwd(i_litr3) = .false.
initial_cn_ratio(i_litr3) = 90._r8
initial_cp_ratio(i_litr3) = 900._r8
initial_stock(i_litr3) = 0._r8
is_metabolic(i_litr3) = .false.
is_cellulose(i_litr3) = .false.
is_lignin(i_litr3) = .true.

! CWD
floating_cn_ratio_decomp_pools(i_cwd) = .true.
floating_cp_ratio_decomp_pools(i_cwd) = .true.
decomp_pool_name_restart(i_cwd) = 'cwd'
decomp_pool_name_history(i_cwd) = 'CWD'
decomp_pool_name_long(i_cwd) = 'coarse woody debris'
Expand All @@ -393,13 +440,15 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil(i_cwd) = .false.
is_cwd(i_cwd) = .true.
initial_cn_ratio(i_cwd) = 90._r8
initial_cp_ratio(i_cwd) = 900._r8
initial_stock(i_cwd) = 0._r8
is_metabolic(i_cwd) = .false.
is_cellulose(i_cwd) = .false.
is_lignin(i_cwd) = .false.

i_soil1 = 5
floating_cn_ratio_decomp_pools(i_soil1) = .false.
floating_cp_ratio_decomp_pools(i_soil1) = .true.
decomp_pool_name_restart(i_soil1) = 'soil1'
decomp_pool_name_history(i_soil1) = 'SOIL1'
decomp_pool_name_long(i_soil1) = 'soil 1'
Expand All @@ -408,13 +457,15 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil(i_soil1) = .true.
is_cwd(i_soil1) = .false.
initial_cn_ratio(i_soil1) = cn_s1
initial_cp_ratio(i_soil1) = cn_s1*np_s1_new
initial_stock(i_soil1) = 20._r8
is_metabolic(i_soil1) = .false.
is_cellulose(i_soil1) = .false.
is_lignin(i_soil1) = .false.

i_soil2 = 6
floating_cn_ratio_decomp_pools(i_soil2) = .false.
floating_cp_ratio_decomp_pools(i_soil2) = .true.
decomp_pool_name_restart(i_soil2) = 'soil2'
decomp_pool_name_history(i_soil2) = 'SOIL2'
decomp_pool_name_long(i_soil2) = 'soil 2'
Expand All @@ -423,13 +474,15 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil(i_soil2) = .true.
is_cwd(i_soil2) = .false.
initial_cn_ratio(i_soil2) = cn_s2
initial_cp_ratio(i_soil2) = cn_s2*np_s2_new
initial_stock(i_soil2) = 20._r8
is_metabolic(i_soil2) = .false.
is_cellulose(i_soil2) = .false.
is_lignin(i_soil2) = .false.

i_soil3 = 7
floating_cn_ratio_decomp_pools(i_soil3) = .false.
floating_cp_ratio_decomp_pools(i_soil3) = .true.
decomp_pool_name_restart(i_soil3) = 'soil3'
decomp_pool_name_history(i_soil3) = 'SOIL3'
decomp_pool_name_long(i_soil3) = 'soil 3'
Expand All @@ -438,6 +491,7 @@ subroutine init_decompcascade_bgc(bounds, cnstate_vars, soilstate_vars)
is_soil(i_soil3) = .true.
is_cwd(i_soil3) = .false.
initial_cn_ratio(i_soil3) = cn_s3
initial_cp_ratio(i_soil3) = cn_s3*np_s3_new
initial_stock(i_soil3) = 20._r8
is_metabolic(i_soil3) = .false.
is_cellulose(i_soil3) = .false.
Expand Down