module CNNStateUpdate3Mod !----------------------------------------------------------------------- ! !DESCRIPTION: ! Module for nitrogen state variable update, mortality fluxes. ! Also, sminn leaching flux. ! ! !USES: use shr_kind_mod , only: r8 => shr_kind_r8 use clm_varpar , only: nlevdecomp, ndecomp_pools use clm_time_manager , only : get_step_size_real use clm_varctl , only : iulog, use_nitrif_denitrif,use_matrixcn,use_soil_matrixcn use clm_varpar , only : i_cwd, i_met_lit, i_cel_lit, i_lig_lit use CNVegNitrogenStateType , only : cnveg_nitrogenstate_type use CNVegNitrogenFluxType , only : cnveg_nitrogenflux_type use SoilBiogeochemNitrogenStateType , only : soilbiogeochem_nitrogenstate_type use SoilBiogeochemNitrogenFluxType , only : soilbiogeochem_nitrogenflux_type ! implicit none private ! ! !PUBLIC MEMBER FUNCTIONS: public:: NStateUpdate3 !----------------------------------------------------------------------- contains !----------------------------------------------------------------------- subroutine NStateUpdate3(num_soilc, filter_soilc, num_soilp, filter_soilp, & cnveg_nitrogenflux_inst, cnveg_nitrogenstate_inst, & soilbiogeochem_nitrogenflux_inst, soilbiogeochem_nitrogenstate_inst) ! ! !DESCRIPTION: ! On the radiation time step, update all the prognostic nitrogen state ! variables affected by gap-phase mortality fluxes. Also the Sminn leaching flux. ! NOTE - associate statements have been removed where there are ! no science equations. This increases readability and maintainability. ! ! !ARGUMENTS: integer , intent(in) :: num_soilc ! number of soil columns in filter integer , intent(in) :: filter_soilc(:) ! filter for soil columns integer , intent(in) :: num_soilp ! number of soil patches in filter integer , intent(in) :: filter_soilp(:) ! filter for soil patches type(cnveg_nitrogenflux_type) , intent(in) :: cnveg_nitrogenflux_inst type(cnveg_nitrogenstate_type) , intent(inout) :: cnveg_nitrogenstate_inst type(soilbiogeochem_nitrogenflux_type) , intent(inout) :: soilbiogeochem_nitrogenflux_inst type(soilbiogeochem_nitrogenstate_type) , intent(inout) :: soilbiogeochem_nitrogenstate_inst ! ! !LOCAL VARIABLES: integer :: c,p,j,l,k ! indices integer :: fp,fc ! lake filter indices real(r8):: dt ! radiation time step (seconds) !----------------------------------------------------------------------- associate( & nf_veg => cnveg_nitrogenflux_inst , & ! Input ns_veg => cnveg_nitrogenstate_inst , & ! Output nf_soil => soilbiogeochem_nitrogenflux_inst , & ! Input ns_soil => soilbiogeochem_nitrogenstate_inst & ! Output ) ! set time steps dt = get_step_size_real() do j = 1, nlevdecomp ! column loop do fc = 1,num_soilc c = filter_soilc(fc) if (.not. use_nitrif_denitrif) then ! mineral N loss due to leaching ns_soil%sminn_vr_col(c,j) = ns_soil%sminn_vr_col(c,j) - nf_soil%sminn_leached_vr_col(c,j) * dt else ! mineral N loss due to leaching and runoff ns_soil%smin_no3_vr_col(c,j) = max( ns_soil%smin_no3_vr_col(c,j) - & ( nf_soil%smin_no3_leached_vr_col(c,j) + nf_soil%smin_no3_runoff_vr_col(c,j) ) * dt, 0._r8) ns_soil%sminn_vr_col(c,j) = ns_soil%smin_no3_vr_col(c,j) + ns_soil%smin_nh4_vr_col(c,j) end if ! column level nitrogen fluxes from fire ! patch-level wood to column-level CWD (uncombusted wood) if (.not. use_soil_matrixcn)then ns_soil%decomp_npools_vr_col(c,j,i_cwd) = ns_soil%decomp_npools_vr_col(c,j,i_cwd) + & nf_veg%fire_mortality_n_to_cwdn_col(c,j) * dt ! patch-level wood to column-level litter (uncombusted wood) ns_soil%decomp_npools_vr_col(c,j,i_met_lit) = ns_soil%decomp_npools_vr_col(c,j,i_met_lit) + & nf_veg%m_n_to_litr_met_fire_col(c,j)* dt ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) = ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) + & nf_veg%m_n_to_litr_cel_fire_col(c,j)* dt ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) = ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) + & nf_veg%m_n_to_litr_lig_fire_col(c,j)* dt else nf_soil%matrix_Ninput%V(c,j+(i_cwd-1)*nlevdecomp) = nf_soil%matrix_Ninput%V(c,j+(i_cwd-1)*nlevdecomp) + & nf_veg%fire_mortality_n_to_cwdn_col(c,j) * dt ! patch-level wood to column-level litter (uncombusted wood) nf_soil%matrix_Ninput%V(c,j+(i_met_lit-1)*nlevdecomp) = nf_soil%matrix_Ninput%V(c,j+(i_met_lit-1)*nlevdecomp) + & nf_veg%m_n_to_litr_met_fire_col(c,j)* dt nf_soil%matrix_Ninput%V(c,j+(i_cel_lit-1)*nlevdecomp) = nf_soil%matrix_Ninput%V(c,j+(i_cel_lit-1)*nlevdecomp) + & nf_veg%m_n_to_litr_cel_fire_col(c,j)* dt nf_soil%matrix_Ninput%V(c,j+(i_lig_lit-1)*nlevdecomp) = nf_soil%matrix_Ninput%V(c,j+(i_lig_lit-1)*nlevdecomp) + & nf_veg%m_n_to_litr_lig_fire_col(c,j)* dt end if ! not use_soil_matrix end do ! end of column loop end do ! litter and CWD losses to fire if(.not. use_soil_matrixcn)then do l = 1, ndecomp_pools do j = 1, nlevdecomp ! column loop do fc = 1,num_soilc c = filter_soilc(fc) ns_soil%decomp_npools_vr_col(c,j,l) = ns_soil%decomp_npools_vr_col(c,j,l) - & nf_veg%m_decomp_npools_to_fire_vr_col(c,j,l) * dt end do end do end do end if ! not use_soil_matrixcn ! patch-level nitrogen fluxes do fp = 1,num_soilp p = filter_soilp(fp) if(.not. use_matrixcn)then !from fire displayed pools ns_veg%leafn_patch(p) = ns_veg%leafn_patch(p) - & nf_veg%m_leafn_to_fire_patch(p) * dt ns_veg%frootn_patch(p) = ns_veg%frootn_patch(p) - & nf_veg%m_frootn_to_fire_patch(p) * dt ns_veg%livestemn_patch(p) = ns_veg%livestemn_patch(p) - & nf_veg%m_livestemn_to_fire_patch(p) * dt ns_veg%deadstemn_patch(p) = ns_veg%deadstemn_patch(p) - & nf_veg%m_deadstemn_to_fire_patch(p) * dt ns_veg%livecrootn_patch(p) = ns_veg%livecrootn_patch(p) - & nf_veg%m_livecrootn_to_fire_patch(p) * dt ns_veg%deadcrootn_patch(p) = ns_veg%deadcrootn_patch(p) - & nf_veg%m_deadcrootn_to_fire_patch(p) * dt ns_veg%leafn_patch(p) = ns_veg%leafn_patch(p) - & nf_veg%m_leafn_to_litter_fire_patch(p) * dt ns_veg%frootn_patch(p) = ns_veg%frootn_patch(p) - & nf_veg%m_frootn_to_litter_fire_patch(p) * dt ns_veg%livestemn_patch(p) = ns_veg%livestemn_patch(p) - & nf_veg%m_livestemn_to_litter_fire_patch(p) * dt - & nf_veg%m_livestemn_to_deadstemn_fire_patch(p) * dt ns_veg%deadstemn_patch(p) = ns_veg%deadstemn_patch(p) - & nf_veg%m_deadstemn_to_litter_fire_patch(p) * dt + & nf_veg%m_livestemn_to_deadstemn_fire_patch(p) * dt ns_veg%livecrootn_patch(p) = ns_veg%livecrootn_patch(p) - & nf_veg%m_livecrootn_to_litter_fire_patch(p) * dt - & nf_veg%m_livecrootn_to_deadcrootn_fire_patch(p) * dt ns_veg%deadcrootn_patch(p) = ns_veg%deadcrootn_patch(p) - & nf_veg%m_deadcrootn_to_litter_fire_patch(p) * dt + & nf_veg%m_livecrootn_to_deadcrootn_fire_patch(p) * dt ! storage pools ns_veg%leafn_storage_patch(p) = ns_veg%leafn_storage_patch(p) - & nf_veg%m_leafn_storage_to_fire_patch(p) * dt ns_veg%frootn_storage_patch(p) = ns_veg%frootn_storage_patch(p) - & nf_veg%m_frootn_storage_to_fire_patch(p) * dt ns_veg%livestemn_storage_patch(p) = ns_veg%livestemn_storage_patch(p) - & nf_veg%m_livestemn_storage_to_fire_patch(p) * dt ns_veg%deadstemn_storage_patch(p) = ns_veg%deadstemn_storage_patch(p) - & nf_veg%m_deadstemn_storage_to_fire_patch(p) * dt ns_veg%livecrootn_storage_patch(p) = ns_veg%livecrootn_storage_patch(p) - & nf_veg%m_livecrootn_storage_to_fire_patch(p) * dt ns_veg%deadcrootn_storage_patch(p) = ns_veg%deadcrootn_storage_patch(p) - & nf_veg%m_deadcrootn_storage_to_fire_patch(p) * dt ns_veg%leafn_storage_patch(p) = ns_veg%leafn_storage_patch(p) - & nf_veg%m_leafn_storage_to_litter_fire_patch(p) * dt ns_veg%frootn_storage_patch(p) = ns_veg%frootn_storage_patch(p) - & nf_veg%m_frootn_storage_to_litter_fire_patch(p) * dt ns_veg%livestemn_storage_patch(p) = ns_veg%livestemn_storage_patch(p) - & nf_veg%m_livestemn_storage_to_litter_fire_patch(p) * dt ns_veg%deadstemn_storage_patch(p) = ns_veg%deadstemn_storage_patch(p) - & nf_veg%m_deadstemn_storage_to_litter_fire_patch(p) * dt ns_veg%livecrootn_storage_patch(p) = ns_veg%livecrootn_storage_patch(p) - & nf_veg%m_livecrootn_storage_to_litter_fire_patch(p) * dt ns_veg%deadcrootn_storage_patch(p) = ns_veg%deadcrootn_storage_patch(p) - & nf_veg%m_deadcrootn_storage_to_litter_fire_patch(p) * dt ! transfer pools ns_veg%leafn_xfer_patch(p) = ns_veg%leafn_xfer_patch(p) - & nf_veg%m_leafn_xfer_to_fire_patch(p) * dt ns_veg%frootn_xfer_patch(p) = ns_veg%frootn_xfer_patch(p) - & nf_veg%m_frootn_xfer_to_fire_patch(p) * dt ns_veg%livestemn_xfer_patch(p) = ns_veg%livestemn_xfer_patch(p) - & nf_veg%m_livestemn_xfer_to_fire_patch(p) * dt ns_veg%deadstemn_xfer_patch(p) = ns_veg%deadstemn_xfer_patch(p) - & nf_veg%m_deadstemn_xfer_to_fire_patch(p) * dt ns_veg%livecrootn_xfer_patch(p) = ns_veg%livecrootn_xfer_patch(p) - & nf_veg%m_livecrootn_xfer_to_fire_patch(p) * dt ns_veg%deadcrootn_xfer_patch(p) = ns_veg%deadcrootn_xfer_patch(p) - & nf_veg%m_deadcrootn_xfer_to_fire_patch(p) * dt ns_veg%leafn_xfer_patch(p) = ns_veg%leafn_xfer_patch(p) - & nf_veg%m_leafn_xfer_to_litter_fire_patch(p) * dt ns_veg%frootn_xfer_patch(p) = ns_veg%frootn_xfer_patch(p) - & nf_veg%m_frootn_xfer_to_litter_fire_patch(p) * dt ns_veg%livestemn_xfer_patch(p) = ns_veg%livestemn_xfer_patch(p) - & nf_veg%m_livestemn_xfer_to_litter_fire_patch(p) * dt ns_veg%deadstemn_xfer_patch(p) = ns_veg%deadstemn_xfer_patch(p) - & nf_veg%m_deadstemn_xfer_to_litter_fire_patch(p) * dt ns_veg%livecrootn_xfer_patch(p) = ns_veg%livecrootn_xfer_patch(p) - & nf_veg%m_livecrootn_xfer_to_litter_fire_patch(p) * dt ns_veg%deadcrootn_xfer_patch(p) = ns_veg%deadcrootn_xfer_patch(p) - & nf_veg%m_deadcrootn_xfer_to_litter_fire_patch(p) * dt ! retranslocated N pool ns_veg%retransn_patch(p) = ns_veg%retransn_patch(p) - & nf_veg%m_retransn_to_fire_patch(p) * dt ns_veg%retransn_patch(p) = ns_veg%retransn_patch(p) - & nf_veg%m_retransn_to_litter_fire_patch(p) * dt else ! NOTE: The equivalent changes for matrix code are in CNFireBase and CNFireLi2014 codes EBK (11/26/2019) end if !.not. use_matrixcn end do end associate end subroutine NStateUpdate3 end module CNNStateUpdate3Mod