From af4a6889c36f9759a7637b0cd57050eed341a28a Mon Sep 17 00:00:00 2001 From: Sam Rabin Date: Wed, 25 Oct 2023 12:02:56 -0600 Subject: [PATCH] Added *C_TO_SEED_PERHARV hist fields. --- src/biogeochem/CNVegCarbonFluxType.F90 | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/src/biogeochem/CNVegCarbonFluxType.F90 b/src/biogeochem/CNVegCarbonFluxType.F90 index d6f5375b5d..c7aa3469e2 100644 --- a/src/biogeochem/CNVegCarbonFluxType.F90 +++ b/src/biogeochem/CNVegCarbonFluxType.F90 @@ -981,6 +981,7 @@ subroutine InitHistory(this, bounds, carbon_type) ptr_patch=this%livestemc_to_biofuelc_patch) this%repr_grainc_to_seed_patch(begp:endp,:) = spval + this%repr_grainc_to_seed_perharv_patch(begp:endp,:,:) = spval this%repr_grainc_to_seed_thisyr_patch(begp:endp,:) = spval do k = repr_grain_min, repr_grain_max data1dptr => this%repr_grainc_to_seed_patch(:,k) @@ -991,6 +992,16 @@ subroutine InitHistory(this, bounds, carbon_type) avgflag='A', & long_name=get_repr_longname(k)//' C to seed', & ptr_patch=data1dptr) + data2dptr => this%repr_grainc_to_seed_perharv_patch(:,:,k) + call hist_addfld2d ( & + ! e.g., GRAINC_TO_SEED_PERHARV + fname=get_repr_hist_fname(k)//'C_TO_SEED_PERHARV', & + units='gC/m^2', & + type2d='mxharvests', & + avgflag='I', & + long_name=get_repr_longname(k)//' C to seed per harvest; should only be output annually', & + ptr_patch=data2dptr, & + default='inactive') data1dptr => this%repr_grainc_to_seed_thisyr_patch(:,k) call hist_addfld1d ( & ! e.g., GRAINC_TO_SEED_ANN