Skip to content

Commit

Permalink
prepare fluxes
Browse files Browse the repository at this point in the history
  • Loading branch information
@mvdebolskiy
mvdebolskiy committed Apr 21, 2024
1 parent 50d6784 commit 902f1a5
Show file tree
Hide file tree
Showing 2 changed files with 89 additions and 15 deletions.
93 changes: 81 additions & 12 deletions src/biogeochem/CNFireEmissionsMod.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -12,12 +12,14 @@ module CNFireEmissionsMod
use decompMod, only : bounds_type
use shr_fire_emis_mod, only : shr_fire_emis_comps_n, shr_fire_emis_comp_t, shr_fire_emis_linkedlist
use shr_fire_emis_mod, only : shr_fire_emis_mechcomps_n, shr_fire_emis_mechcomps
use clm_varctl, only : use_fates_bgc
!
implicit none
private
!
! !PUBLIC MEMBER FUNCTIONS:
public :: CNFireEmisUpdate
public :: FatesFireEmisUpdate
!
! !PRIVATE TYPES:
type, private :: emis_t
Expand Down Expand Up @@ -60,20 +62,26 @@ subroutine Init(this, bounds)

if ( shr_fire_emis_mechcomps_n < 1) return

call fire_emis_factors_init( shr_fire_emis_factors_file )

emis_cmp => shr_fire_emis_linkedlist
do while(associated(emis_cmp))
allocate(emis_cmp%emis_factors(maxveg))
call fire_emis_factors_get( trim(emis_cmp%name), factors, molec_wght )
emis_cmp%emis_factors = factors*1.e-3_r8 ! convert g/kg dry fuel to kg/kg
emis_cmp%molec_weight = molec_wght
emis_cmp => emis_cmp%next_emiscomp
enddo
if (.not. use_fates_bgc) then
if ( shr_fire_emis_factors_file == ' ') then
call endrun('ERROR:: fire_emis_factors_file must be set to use fire emissions with use_cn' )
endif
call fire_emis_factors_init( shr_fire_emis_factors_file )
emis_cmp => shr_fire_emis_linkedlist
do while(associated(emis_cmp))
allocate(emis_cmp%emis_factors(maxveg))
call fire_emis_factors_get( trim(emis_cmp%name), factors, molec_wght )
emis_cmp%emis_factors = factors*1.e-3_r8 ! convert g/kg dry fuel to kg/kg
emis_cmp%molec_weight = molec_wght
emis_cmp => emis_cmp%next_emiscomp
enddo
endif

call this%InitAllocate(bounds)
call this%InitHistory(bounds)

! FATES has it's own history!
if (.not. use_fates_bgc) then
call this%InitHistory(bounds)
endif
end subroutine Init

!-----------------------------------------------------------------------
Expand Down Expand Up @@ -315,6 +323,67 @@ subroutine CNFireEmisUpdate(bounds, num_bgc_vegp, filter_bgc_vegp, cnveg_cf_inst

end subroutine CNFireEmisUpdate


subroutine FatesFireEmisUpdate(bounds, num_bgc_vegp, filter_bgc_vegp, clm_fates, fireemis_inst)

use CLMFatesInterfaceMod, only : hlm_fates_interface_type

!ARGUMENTS:
type(bounds_type), intent(in) :: bounds
integer, intent(in) :: num_bgc_vegp ! number of bgc veg patches
integer, intent(in) :: filter_bgc_vegp(:) ! filter for bgc veg patches
type(hlm_fates_interface_type), intent(in) :: clm_fates
type(fireemis_type), intent(inout) :: fireemis_inst

!LOCAL VARIABLES:
type(shr_fire_emis_comp_t), pointer :: emis_cmp
real(r8) :: emis_flux(shr_fire_emis_comps_n)
integer :: fp,p,g,c ! indices
integer :: i, ii, icomp, imech, n_emis_comps, l, j

if ( shr_fire_emis_mechcomps_n < 1) return

associate( &
fire_emis => fireemis_inst%fireflx_patch, &
totfire => fireemis_inst%totfire, &
mech => fireemis_inst%mech, &
comp => fireemis_inst%comp, &
ztop => fireemis_inst%ztop_patch &
)

! initialize to zero ...
fire_emis(bounds%begp:bounds%endp,:) = 0._r8
totfire%emis(bounds%begp:bounds%endp) = 0._r8
ztop(bounds%begp:bounds%endp) = 0._r8

do i = 1, shr_fire_emis_mechcomps_n
mech(i)%emis(bounds%begp:bounds%endp) = 0._r8
enddo
do i = 1, shr_fire_emis_comps_n
comp(i)%emis(bounds%begp:bounds%endp) = 0._r8
enddo

do fp = 1,num_bgc_vegp
p = filter_bgc_vegp(fp)

! sum up the emissions compontent fluxes for the fluxes of chem mechanism compounds
do imech = 1,shr_fire_emis_mechcomps_n
n_emis_comps = shr_fire_emis_mechcomps(imech)%n_emis_comps
do icomp = 1,n_emis_comps ! loop over number of emission components that make up the nth mechanism compoud
ii = shr_fire_emis_mechcomps(imech)%emis_comps(icomp)%ptr%index
fire_emis(p,imech) = fire_emis(p,imech) + shr_fire_emis_mechcomps(imech)%coeffs(icomp) * fates_emis(p,ii)
enddo
mech(imech)%emis(p) = fire_emis(p,imech)
enddo

!ztop(p) = vert_dist_top( patch%itype(p) ) TODO: this needs to be rewritten
enddo ! fp

end associate

end subroutine FatesFireEmisUpdate


! Private methods
!-----------------------------------------------------------------------
!ztop compiled from Val Martin et al ACP 2010, Tosca et al. JGR 2011 and Jian et al., ACP 2013
Expand Down
11 changes: 8 additions & 3 deletions src/biogeochem/CNVegetationFacade.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -82,7 +82,7 @@ module CNVegetationFacade
use SoilBiogeochemCarbonFluxType , only : soilBiogeochem_carbonflux_type
use SoilBiogeochemNitrogenStateType , only : soilbiogeochem_nitrogenstate_type
use SoilBiogeochemNitrogenFluxType , only : soilbiogeochem_nitrogenflux_type
use CNFireEmissionsMod , only : fireemis_type, CNFireEmisUpdate
use CNFireEmissionsMod , only : fireemis_type, CNFireEmisUpdate, FatesFireEmisUpdate
use CNDriverMod , only : CNDriverInit
use CNDriverMod , only : CNDriverSummarizeStates, CNDriverSummarizeFluxes
use CNDriverMod , only : CNDriverNoLeaching, CNDriverLeaching
Expand Down Expand Up @@ -1033,8 +1033,13 @@ subroutine EcosystemDynamicsPreDrainage(this, bounds, &
nutrient_competition_method, this%cnfire_method, this%dribble_crophrv_xsmrpool_2atm)

! fire carbon emissions
call CNFireEmisUpdate(bounds, num_bgc_vegp, filter_bgc_vegp, &
this%cnveg_carbonflux_inst, this%cnveg_carbonstate_inst, fireemis_inst )
if (use_fates_bgc) then
call FatesFireEmisUpdate(bounds, num_bgc_vegp, filter_bgc_vegp, &
clm_fates, fireemis_inst)
else
call CNFireEmisUpdate(bounds, num_bgc_vegp, filter_bgc_vegp, &
this%cnveg_carbonflux_inst, this%cnveg_carbonstate_inst, fireemis_inst )
endif

call CNAnnualUpdate(bounds, &
num_bgc_soilc, filter_bgc_soilc, &
Expand Down

0 comments on commit 902f1a5

Please sign in to comment.