Skip to content

Rename hist fields to track them down more easily #2106

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

Merged
slevis-lmwg merged 5 commits into from
Aug 21, 2023
Merged

Rename hist fields to track them down more easily #2106

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
394f157
Reverse hist field strings such as MET_LIT to LIT_MET
slevis-lmwg Aug 15, 2023
2fb7513
Change C or N to _C and _N after hist field strings like LIT_MET
slevis-lmwg Aug 15, 2023
3e12e28
Update hist field names in user_nl_clm files
slevis-lmwg Aug 15, 2023
5650e9e
Merge remote-tracking branch 'escomp/master' into rename_mimics_hist_…
slevis-lmwg Aug 21, 2023
e0fd4f2
Update ChangeLog and ChangeSum
slevis-lmwg Aug 21, 2023
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
4 changes: 2 additions & 2 deletions cime_config/usermods_dirs/newton_krylov_spinup/user_nl_clm
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
hist_dov2xy = .true.,.false.
hist_nhtfrq = 0,-175200
hist_mfilt = 1,1
hist_fincl2 = 'FPI_vr', 'K_PAS_SOM', 'K_SLO_SOM', 'K_ACT_SOM',
'K_CWD', 'K_CEL_LIT', 'K_LIG_LIT', 'K_MET_LIT',
hist_fincl2 = 'FPI_vr', 'K_SOM_PAS', 'K_SOM_SLO', 'K_SOM_ACT',
'K_CWD', 'K_LIT_CEL', 'K_LIT_LIG', 'K_LIT_MET',
'CWD_PATHFRAC_L2_vr', 'CWD_RESP_FRAC_L2_vr',
'CWD_PATHFRAC_L3_vr', 'CWD_RESP_FRAC_L3_vr',
'L1_PATHFRAC_S1_vr', 'L1_RESP_FRAC_S1_vr',
Expand Down
6 changes: 3 additions & 3 deletions cime_config/usermods_dirs/output_bgc/user_nl_clm
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,8 +3,8 @@
!----------------------------------------------------------------------------------

! h0 stream (monthly average, gridcell-level)
hist_fexcl1 += 'ACT_SOMC_vr', 'ACT_SOMN_vr', 'SLO_SOMC_vr', 'SLO_SOMN_vr', 'PAS_SOMC_vr', 'PAS_SOMN_vr', 'SOILC_vr','SOILN_vr', 'CWDC_vr', 'MET_LITC_vr', 'CEL_LITC_vr', 'LIG_LITC_vr', 'MET_LITN_vr', 'CEL_LITN_vr', 'LIG_LITN_vr', 'CWDN_vr', 'SMIN_NO3_vr', 'CONC_O2_UNSAT', 'CONC_O2_SAT','SMIN_NH4_vr','SMINN_vr'
hist_fincl1 += 'LEAFC_TO_LITTER', 'FROOTC_TO_LITTER','MET_LITC_TO_ACT_SOMC','MET_LITN_TO_ACT_SOMN','CEL_LITC_TO_ACT_SOMC', 'CEL_LITN_TO_ACT_SOMN','LIG_LITC_TO_SLO_SOMC','LIG_LITN_TO_SLO_SOMN','DWT_WOOD_PRODUCTC_GAIN_PATCH'
hist_fexcl1 += 'SOM_ACT_C_vr', 'SOM_ACT_N_vr', 'SOM_SLO_C_vr', 'SOM_SLO_N_vr', 'SOM_PAS_C_vr', 'SOM_PAS_N_vr', 'SOILC_vr','SOILN_vr', 'CWD_C_vr', 'LIT_MET_C_vr', 'LIT_CEL_C_vr', 'LIT_LIG_C_vr', 'LIT_MET_N_vr', 'LIT_CEL_N_vr', 'LIT_LIG_N_vr', 'CWD_N_vr', 'SMIN_NO3_vr', 'CONC_O2_UNSAT', 'CONC_O2_SAT','SMIN_NH4_vr','SMINN_vr'
hist_fincl1 += 'LEAFC_TO_LITTER', 'FROOTC_TO_LITTER','LIT_MET_C_TO_SOM_ACT_C','LIT_MET_N_TO_SOM_ACT_N','LIT_CEL_C_TO_SOM_ACT_C', 'LIT_CEL_N_TO_SOM_ACT_N','LIT_LIG_C_TO_SOM_SLO_C','LIT_LIG_N_TO_SOM_SLO_N','DWT_WOOD_PRODUCTC_GAIN_PATCH'

! h1 stream (monthly average, finest sub-grid)
hist_fincl2 += 'GPP', 'NPP', 'AGNPP', 'TOTVEGC', 'NPP_NUPTAKE', 'AR', 'HR', 'HTOP'
Expand All @@ -14,7 +14,7 @@ hist_fincl2 += 'GPP', 'NPP', 'AGNPP', 'TOTVEGC', 'NPP_NUPTAKE', 'AR', 'HR', 'HTO
hist_fincl3 += 'GPP', 'NPP', 'AR', 'HR', 'DWT_CONV_CFLUX_PATCH', 'WOOD_HARVESTC', 'DWT_WOOD_PRODUCTC_GAIN_PATCH', 'SLASH_HARVESTC', 'COL_FIRE_CLOSS', 'FROOTC:I', 'HTOP'

! h3 stream (yearly average, gridcell-level)
hist_fincl4 += 'SOILC_vr', 'SOILN_vr', 'CWDC_vr', 'MET_LITC_vr', 'CEL_LITC_vr', 'LIG_LITC_vr', 'MET_LITN_vr', 'CEL_LITN_vr', 'LIG_LITN_vr','CWDN_vr', 'TOTLITC:I', 'TOT_WOODPRODC:I', 'TOTSOMC:I','TOTVEGC:I'
hist_fincl4 += 'SOILC_vr', 'SOILN_vr', 'CWD_C_vr', 'LIT_MET_C_vr', 'LIT_CEL_C_vr', 'LIT_LIG_C_vr', 'LIT_MET_N_vr', 'LIT_CEL_N_vr', 'LIT_LIG_N_vr','CWD_N_vr', 'TOTLITC:I', 'TOT_WOODPRODC:I', 'TOTSOMC:I','TOTVEGC:I'

! h4 stream (yearly average, landunit-level)
hist_fincl5 += 'TOTSOMC:I', 'TOTSOMC_1m:I', 'TOTECOSYSC:I', 'TOTVEGC:I', 'WOODC:I', 'TOTLITC:I', 'LIVECROOTC:I', 'DEADCROOTC:I', 'FROOTC:I'
Expand Down
66 changes: 66 additions & 0 deletions doc/ChangeLog
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,70 @@
===============================================================
Tag name: ctsm5.1.dev135
Originator(s): slevis (Samuel Levis,UCAR/TSS,303-665-1310)
Date: Mon Aug 21 15:06:35 MDT 2023
One-line Summary: Rename hist fields to track them down more easily

Purpose and description of changes
----------------------------------

Renaming history fields to make easier to find in lists, e.g. when
using ncview. For example, litter fields like MET_LIT and STR_LIT
will be LIT_MET and LIT_STR.


Significant changes to scientifically-supported configurations
--------------------------------------------------------------

Does this tag change answers significantly for any of the following physics configurations?
(Details of any changes will be given in the "Answer changes" section below.)

[Put an [X] in the box for any configuration with significant answer changes.]

[ ] clm5_1

[ ] clm5_0

[ ] ctsm5_0-nwp

[ ] clm4_5


Bugs fixed or introduced
------------------------
CTSM issues fixed (include CTSM Issue #):
Fixes #2095


Testing summary:
----------------
[Remove any lines that don't apply.]

[PASS means all tests PASS; OK means tests PASS other than expected fails.]

regular tests (aux_clm: https://github.com/ESCOMP/CTSM/wiki/System-Testing-Guide#pre-merge-system-testing):

cheyenne ---- OK
izumi ------- OK

any other testing (give details below):
Sam L. ran the LMWG diag. pkg and found only one plot affected by this
PR's changes. In particular, set 6 CWD_C, which was CWDC


Answer changes
--------------

Changes answers relative to baseline:
No. Field lists differ. In some tests, the namelists differ.


Other details
-------------
Pull Requests that document the changes (include PR ids):
https://github.com/ESCOMP/ctsm/pull/2106

===============================================================
===============================================================
Tag name: ctsm5.1.dev134
Originator(s): rgknox (Ryan Knox,LBNL EESA), erik (Erik Kluzek,UCAR/TSS,303-497-1326)
Date: Wed Aug 16 17:20:27 MDT 2023
Expand Down
1 change: 1 addition & 0 deletions doc/ChangeSum
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,6 @@
Tag Who Date Summary
============================================================================================================================
ctsm5.1.dev135 slevis 08/21/2023 Rename hist fields to track them down more easily
ctsm5.1.dev134 rgknox 08/16/2023 Migration of FATES to share normal soil BGC call sequence and functionality
ctsm5.1.dev133 glemieux 08/09/2023 FATES API update to facilitate fates refactor
ctsm5.1.dev132 slevis 08/04/2023 Add parameterization to allow excess ice in soil and subsidence
Expand Down
12 changes: 6 additions & 6 deletions src/biogeochem/CNVegCarbonFluxType.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2967,15 +2967,15 @@ subroutine InitHistory(this, bounds, carbon_type)
do k = 1, ndecomp_pools
if ( decomp_cascade_con%is_litter(k) .or. decomp_cascade_con%is_cwd(k) ) then
data1dptr => this%m_decomp_cpools_to_fire_col(:,k)
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TO_FIRE'
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TO_FIRE'
longname = trim(decomp_cascade_con%decomp_pool_name_long(k))//' C fire loss'
call hist_addfld1d (fname=fieldname, units='gC/m^2/s', &
avgflag='A', long_name=longname, &
ptr_col=data1dptr, default='inactive')

if ( nlevdecomp_full > 1 ) then
data2dptr => this%m_decomp_cpools_to_fire_vr_col(:,:,k)
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TO_FIRE'//trim(vr_suffix)
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TO_FIRE'//trim(vr_suffix)
longname = trim(decomp_cascade_con%decomp_pool_name_long(k))//' C fire loss'
call hist_addfld_decomp (fname=fieldname, units='gC/m^3/s', type2d='levdcmp', &
avgflag='A', long_name=longname, &
Expand Down Expand Up @@ -3176,15 +3176,15 @@ subroutine InitHistory(this, bounds, carbon_type)
do k = 1, ndecomp_pools
if ( decomp_cascade_con%is_litter(k) .or. decomp_cascade_con%is_cwd(k) ) then
data1dptr => this%m_decomp_cpools_to_fire_col(:,k)
fieldname = 'C13_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TO_FIRE'
fieldname = 'C13_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TO_FIRE'
longname = 'C13 '//trim(decomp_cascade_con%decomp_pool_name_long(k))//' C fire loss'
call hist_addfld1d (fname=fieldname, units='gC13/m^2', &
avgflag='A', long_name=longname, &
ptr_col=data1dptr, default='inactive')

if ( nlevdecomp_full > 1 ) then
data2dptr => this%m_decomp_cpools_to_fire_vr_col(:,:,k)
fieldname = 'C13_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TO_FIRE'//trim(vr_suffix)
fieldname = 'C13_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TO_FIRE'//trim(vr_suffix)
longname = 'C13 '//trim(decomp_cascade_con%decomp_pool_name_long(k))//' C fire loss'
call hist_addfld_decomp (fname=fieldname, units='gC13/m^3', type2d='levdcmp', &
avgflag='A', long_name=longname, &
Expand Down Expand Up @@ -3336,15 +3336,15 @@ subroutine InitHistory(this, bounds, carbon_type)
do k = 1, ndecomp_pools
if ( decomp_cascade_con%is_litter(k) .or. decomp_cascade_con%is_cwd(k) ) then
data1dptr => this%m_decomp_cpools_to_fire_col(:,k)
fieldname = 'C14_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TO_FIRE'
fieldname = 'C14_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TO_FIRE'
longname = 'C14 '//trim(decomp_cascade_con%decomp_pool_name_long(k))//' C fire loss'
call hist_addfld1d (fname=fieldname, units='gC14/m^2', &
avgflag='A', long_name=longname, &
ptr_col=data1dptr, default='inactive')

if ( nlevdecomp_full > 1 ) then
data2dptr => this%m_decomp_cpools_to_fire_vr_col(:,:,k)
fieldname = 'C14_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TO_FIRE'//trim(vr_suffix)
fieldname = 'C14_M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TO_FIRE'//trim(vr_suffix)
longname = 'C14 '//trim(decomp_cascade_con%decomp_pool_name_long(k))//' C fire loss'
call hist_addfld_decomp (fname=fieldname, units='gC14/m^3', type2d='levdcmp', &
avgflag='A', long_name=longname, &
Expand Down
4 changes: 2 additions & 2 deletions src/biogeochem/CNVegNitrogenFluxType.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1058,7 +1058,7 @@ subroutine InitHistory(this, bounds)
if ( decomp_cascade_con%is_litter(k) .or. decomp_cascade_con%is_cwd(k) ) then
this%m_decomp_npools_to_fire_col(begc:endc,k) = spval
data1dptr => this%m_decomp_npools_to_fire_col(:,k)
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'N_TO_FIRE'
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_N_TO_FIRE'
longname = trim(decomp_cascade_con%decomp_pool_name_long(k))//' N fire loss'
call hist_addfld1d (fname=fieldname, units='gN/m^2', &
avgflag='A', long_name=longname, &
Expand All @@ -1067,7 +1067,7 @@ subroutine InitHistory(this, bounds)
if ( nlevdecomp_full > 1 ) then
this%m_decomp_npools_to_fire_vr_col(begc:endc,:,k) = spval
data2dptr => this%m_decomp_npools_to_fire_vr_col(:,:,k)
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'N_TO_FIRE'//trim(vr_suffix)
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_N_TO_FIRE'//trim(vr_suffix)
longname = trim(decomp_cascade_con%decomp_pool_name_long(k))//' N fire loss'
call hist_addfld_decomp (fname=fieldname, units='gN/m^3', type2d='levdcmp', &
avgflag='A', long_name=longname, &
Expand Down
20 changes: 10 additions & 10 deletions src/soilbiogeochem/SoilBiogeochemCarbonFluxType.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -304,8 +304,8 @@ subroutine InitHistory(this, bounds, carbon_type)
if ( decomp_cascade_con%cascade_receiver_pool(l) /= 0 ) then
data1dptr => this%decomp_cascade_ctransfer_col(:,l)
fieldname = &
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_donor_pool(l)))//'C_TO_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_receiver_pool(l)))//'C'
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_donor_pool(l)))//'_C_TO_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_receiver_pool(l)))//'_C'
longname = 'decomp. of '//trim(decomp_cascade_con%decomp_pool_name_long(decomp_cascade_con%cascade_donor_pool(l)))//&
' C to '//trim(decomp_cascade_con%decomp_pool_name_long(decomp_cascade_con%cascade_receiver_pool(l)))//' C'
call hist_addfld1d (fname=fieldname, units='gC/m^2/s', &
Expand Down Expand Up @@ -342,9 +342,9 @@ subroutine InitHistory(this, bounds, carbon_type)
if ( decomp_cascade_con%cascade_receiver_pool(l) /= 0 ) then
data2dptr => this%decomp_cascade_ctransfer_vr_col(:,:,l)
fieldname = &
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_donor_pool(l)))//'C_TO_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_donor_pool(l)))//'_C_TO_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_receiver_pool(l)))&
//'C'//trim(vr_suffix)
//'_C'//trim(vr_suffix)
longname = 'decomp. of '//&
trim(decomp_cascade_con%decomp_pool_name_long(decomp_cascade_con%cascade_donor_pool(l)))//&
' C to '//&
Expand Down Expand Up @@ -412,14 +412,14 @@ subroutine InitHistory(this, bounds, carbon_type)
do k = 1, ndecomp_pools ! none from CWD
if ( .not. decomp_cascade_con%is_cwd(k) ) then
data1dptr => this%decomp_cpools_leached_col(:,k)
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TO_LEACHING'
fieldname = 'M_'//trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TO_LEACHING'
longname = trim(decomp_cascade_con%decomp_pool_name_long(k))//' C leaching loss'
call hist_addfld1d (fname=fieldname, units='gC/m^2/s', &
avgflag='A', long_name=longname, &
ptr_col=data1dptr, default='inactive')

data2dptr => this%decomp_cpools_transport_tendency_col(:,:,k)
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(k))//'C_TNDNCY_VERT_TRANSPORT'
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(k))//'_C_TNDNCY_VERT_TRANSPORT'
longname = trim(decomp_cascade_con%decomp_pool_name_long(k))//' C tendency due to vertical transport'
call hist_addfld_decomp (fname=fieldname, units='gC/m^3/s', type2d='levdcmp', &
avgflag='A', long_name=longname, &
Expand Down Expand Up @@ -502,9 +502,9 @@ subroutine InitHistory(this, bounds, carbon_type)
data2dptr => this%decomp_cascade_ctransfer_vr_col(:,:,l)
fieldname = 'C13_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_donor_pool(l)))&
//'C_TO_'//&
//'_C_TO_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_receiver_pool(l)))&
//'C'//trim(vr_suffix)
//'_C'//trim(vr_suffix)
longname = 'C13 decomp. of '&
//trim(decomp_cascade_con%decomp_pool_name_long(decomp_cascade_con%cascade_donor_pool(l)))&
//' C to '//&
Expand Down Expand Up @@ -586,9 +586,9 @@ subroutine InitHistory(this, bounds, carbon_type)

fieldname = 'C14_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_donor_pool(l)))&
//'C_TO_'//&
//'_C_TO_'//&
trim(decomp_cascade_con%decomp_pool_name_history(decomp_cascade_con%cascade_receiver_pool(l)))&
//'C'//trim(vr_suffix)
//'_C'//trim(vr_suffix)
longname = 'C14 decomp. of '&
//trim(decomp_cascade_con%decomp_pool_name_long(decomp_cascade_con%cascade_donor_pool(l)))//&
' C to '//trim(decomp_cascade_con%decomp_pool_name_long(decomp_cascade_con%cascade_receiver_pool(l)))//' C'
Expand Down
16 changes: 8 additions & 8 deletions src/soilbiogeochem/SoilBiogeochemCarbonStateType.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -184,23 +184,23 @@ subroutine InitHistory(this, bounds, carbon_type)
do l = 1, ndecomp_pools
if ( nlevdecomp_full > 1 ) then
data2dptr => this%decomp_cpools_vr_col(:,1:nlevsoi,l)
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(l))//'C_vr'
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C_vr'
longname = trim(decomp_cascade_con%decomp_pool_name_history(l))//' C (vertically resolved)'
call hist_addfld2d (fname=fieldname, units='gC/m^3', type2d='levsoi', &
avgflag='A', long_name=longname, &
ptr_col=data2dptr)
endif

data1dptr => this%decomp_cpools_col(:,l)
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(l))//'C'
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C'
longname = trim(decomp_cascade_con%decomp_pool_name_history(l))//' C'
call hist_addfld1d (fname=fieldname, units='gC/m^2', &
avgflag='A', long_name=longname, &
ptr_col=data1dptr)

if ( nlevdecomp_full > 1 ) then
data1dptr => this%decomp_cpools_1m_col(:,l)
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(l))//'C_1m'
fieldname = trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C_1m'
longname = trim(decomp_cascade_con%decomp_pool_name_history(l))//' C to 1 meter'
call hist_addfld1d (fname=fieldname, units='gC/m^2', &
avgflag='A', long_name=longname, &
Expand Down Expand Up @@ -293,15 +293,15 @@ subroutine InitHistory(this, bounds, carbon_type)
do l = 1, ndecomp_pools
if ( nlevdecomp_full > 1 ) then
data2dptr => this%decomp_cpools_vr_col(:,1:nlevsoi,l)
fieldname = 'C13_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'C_vr'
fieldname = 'C13_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C_vr'
longname = 'C13 '//trim(decomp_cascade_con%decomp_pool_name_history(l))//' C (vertically resolved)'
call hist_addfld2d (fname=fieldname, units='gC13/m^3', type2d='levsoi', &
avgflag='A', long_name=longname, &
ptr_col=data2dptr, default='inactive')
endif

data1dptr => this%decomp_cpools_col(:,l)
fieldname = 'C13_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'C'
fieldname = 'C13_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C'
longname = 'C13 '//trim(decomp_cascade_con%decomp_pool_name_history(l))//' C'
call hist_addfld1d (fname=fieldname, units='gC13/m^2', &
avgflag='A', long_name=longname, &
Expand Down Expand Up @@ -391,21 +391,21 @@ subroutine InitHistory(this, bounds, carbon_type)
do l = 1, ndecomp_pools
if ( nlevdecomp_full > 1 ) then
data2dptr => this%decomp_cpools_vr_col(:,1:nlevsoi,l)
fieldname = 'C14_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'C_vr'
fieldname = 'C14_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C_vr'
longname = 'C14 '//trim(decomp_cascade_con%decomp_pool_name_history(l))//' C (vertically resolved)'
call hist_addfld2d (fname=fieldname, units='gC14/m^3', type2d='levsoi', &
avgflag='A', long_name=longname, ptr_col=data2dptr, default='inactive')
endif

data1dptr => this%decomp_cpools_col(:,l)
fieldname = 'C14_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'C'
fieldname = 'C14_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C'
longname = 'C14 '//trim(decomp_cascade_con%decomp_pool_name_history(l))//' C'
call hist_addfld1d (fname=fieldname, units='gC14/m^2', &
avgflag='A', long_name=longname, ptr_col=data1dptr, default='inactive')

if ( nlevdecomp_full > 1 ) then
data1dptr => this%decomp_cpools_1m_col(:,l)
fieldname = 'C14_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'C_1m'
fieldname = 'C14_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//'_C_1m'
longname = 'C14_'//trim(decomp_cascade_con%decomp_pool_name_history(l))//' C to 1 meter'
call hist_addfld1d (fname=fieldname, units='gC/m^2', &
avgflag='A', long_name=longname, ptr_col=data1dptr, default='inactive')
Expand Down
Loading