Merge pull request #1 from ckoven/radtemp

added summary N litterfall flux variable
HPSCTerrSys · Feb 24, 2018 · 386dffd · 386dffd
2 parents 45f3d77 + 4947261
commit 386dffd
Showing 1 changed file with 51 additions and 0 deletions.
diff --git a/src/biogeochem/CNVegNitrogenFluxType.F90 b/src/biogeochem/CNVegNitrogenFluxType.F90
@@ -178,6 +178,7 @@ module CNVegNitrogenFluxType
      real(r8), pointer :: ndeploy_patch                             (:)     ! patch total N deployed to growth and storage (gN/m2/s)
      real(r8), pointer :: wood_harvestn_patch                       (:)     ! patch total N losses to wood product pools (gN/m2/s)
      real(r8), pointer :: wood_harvestn_col                         (:)     ! col total N losses to wood product pools (gN/m2/s) (p2c)
+     real(r8), pointer :: litfall_n_patch                           (:)     ! total N passed from plant to litter (gN/m2/s)
 
      ! phenology: litterfall and crop fluxes
      real(r8), pointer :: phenology_n_to_litr_met_n_col             (:,:)   ! col N fluxes associated with phenology (litterfall and crop) to litter metabolic pool (gN/m3/s)
@@ -441,6 +442,8 @@ subroutine InitAllocate(this, bounds)
     allocate(this%dwt_crop_productn_gain_patch (begp:endp))                   ; this%dwt_crop_productn_gain_patch (:)   = nan
     allocate(this%wood_harvestn_col            (begc:endc))                   ; this%wood_harvestn_col            (:)   = nan
 
+    allocate(this%litfall_n_patch (begp:endp))                                ; this%litfall_n_patch (:)                = nan
+
     allocate(this%dwt_frootn_to_litr_met_n_col (begc:endc,1:nlevdecomp_full)) ; this%dwt_frootn_to_litr_met_n_col (:,:) = nan
     allocate(this%dwt_frootn_to_litr_cel_n_col (begc:endc,1:nlevdecomp_full)) ; this%dwt_frootn_to_litr_cel_n_col (:,:) = nan
     allocate(this%dwt_frootn_to_litr_lig_n_col (begc:endc,1:nlevdecomp_full)) ; this%dwt_frootn_to_litr_lig_n_col (:,:) = nan
@@ -1034,6 +1037,11 @@ subroutine InitHistory(this, bounds)
          '(per-area-gridcell; only makes sense with dov2xy=.false.)', &
          ptr_patch=this%dwt_conv_nflux_patch, default='inactive')
 
+    this%litfall_n_patch(begp:endp) = spval
+    call hist_addfld1d (fname='LITFALL_N', units='gN/m^2/s', &
+         avgflag='A', long_name='N Flux from vegetation to soil', &
+         ptr_patch=this%litfall_n_patch)
+
     this%dwt_frootn_to_litr_met_n_col(begc:endc,:) = spval
     call hist_addfld_decomp (fname='DWT_FROOTN_TO_LITR_MET_N', units='gN/m^2/s',  type2d='levdcmp', &
          avgflag='A', long_name='fine root to litter due to landcover change', &
@@ -1820,6 +1828,49 @@ subroutine Summary_nitrogenflux(this, bounds, num_soilc, filter_soilc, num_soilp
             this%m_deadcrootn_xfer_to_fire_patch(p)     + &
             this%m_retransn_to_fire_patch(p)
 
+       ! total N fluxes from veg to litter
+       this%litfall_n_patch(p) = &
+            this%m_leafn_to_litter_patch(p)     + &
+            this%m_frootn_to_litter_patch(p)     + &
+            this%m_leafn_storage_to_litter_patch(p)     + &
+            this%m_frootn_storage_to_litter_patch(p)     + &
+            this%m_livestemn_storage_to_litter_patch(p)     + &
+            this%m_deadstemn_storage_to_litter_patch(p)     + &
+            this%m_livecrootn_storage_to_litter_patch(p)     + &
+            this%m_deadcrootn_storage_to_litter_patch(p)     + &
+            this%m_leafn_xfer_to_litter_patch(p)     + &
+            this%m_frootn_xfer_to_litter_patch(p)     + &
+            this%m_livestemn_xfer_to_litter_patch(p)     + &
+            this%m_deadstemn_xfer_to_litter_patch(p)     + &
+            this%m_livecrootn_xfer_to_litter_patch(p)     + &
+            this%m_deadcrootn_xfer_to_litter_patch(p)     + &
+            this%m_livestemn_to_litter_patch(p)     + &
+            this%m_deadstemn_to_litter_patch(p)     + &
+            this%m_livecrootn_to_litter_patch(p)     + &
+            this%m_deadcrootn_to_litter_patch(p)     + &
+            this%m_retransn_to_litter_patch(p)     + &
+            this%hrv_leafn_to_litter_patch(p)     + &
+            this%hrv_frootn_to_litter_patch(p)     + &
+            this%hrv_leafn_storage_to_litter_patch(p)     + &
+            this%hrv_frootn_storage_to_litter_patch(p)     + &
+            this%hrv_livestemn_storage_to_litter_patch(p)     + &
+            this%hrv_deadstemn_storage_to_litter_patch(p)     + &
+            this%hrv_livecrootn_storage_to_litter_patch(p)     + &
+            this%hrv_deadcrootn_storage_to_litter_patch(p)     + &
+            this%hrv_leafn_xfer_to_litter_patch(p)     + &
+            this%hrv_frootn_xfer_to_litter_patch(p)     + &
+            this%hrv_livestemn_xfer_to_litter_patch(p)     + &
+            this%hrv_deadstemn_xfer_to_litter_patch(p)     + &
+            this%hrv_livecrootn_xfer_to_litter_patch(p)     + &
+            this%hrv_deadcrootn_xfer_to_litter_patch(p)     + &
+            this%hrv_livestemn_to_litter_patch(p)     + &
+            this%hrv_livecrootn_to_litter_patch(p)     + &
+            this%hrv_deadcrootn_to_litter_patch(p)     + &
+            this%hrv_retransn_to_litter_patch(p)     + &
+            this%livestemn_to_litter_patch(p)     + &
+            this%leafn_to_litter_patch(p)     + &
+            this%frootn_to_litter_patch(p)
+
     end do
 
     call p2c(bounds, num_soilc, filter_soilc, &