Skip to content

ELM-side changes for outputting biophysics variables by land use type #7783

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

Draft
wants to merge 1 commit into
base: master
Choose a base branch
from
Draft

ELM-side changes for outputting biophysics variables by land use type #7783

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
2 changes: 1 addition & 1 deletion components/elm/src/main/elm_driver.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1291,7 +1291,7 @@ subroutine elm_drv(doalb, nextsw_cday, declinp1, declin, rstwr, nlend, rdate)
call alm_fates%WrapUpdateFatesRmean(nc)

! Update high-frequency history diagnostics for FATES
call alm_fates%wrap_update_hifrq_hist(bounds_clump)
call alm_fates%wrap_update_hifrq_hist(bounds_clump, solarabs_vars)
if ( is_beg_curr_day() ) then ! run ED at the start of each day
call alm_fates%dynamics_driv( bounds_clump, top_as, &
top_af, atm2lnd_vars, soilstate_vars, &
Expand Down
29 changes: 26 additions & 3 deletions components/elm/src/main/elmfates_interfaceMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -121,7 +121,7 @@ module ELMFatesInterfaceMod
use ColumnType , only : col_pp
use ColumnDataType , only : col_es, col_ws, col_wf, col_cs, col_cf
use ColumnDataType , only : col_nf, col_pf
use VegetationDataType, only : veg_es, veg_wf, veg_ws
use VegetationDataType, only : veg_es, veg_wf, veg_ws, veg_ef
use LandunitType , only : lun_pp

use landunit_varcon , only : istsoil
Expand Down Expand Up @@ -2994,24 +2994,31 @@ end subroutine WrapUpdateFatesSeedInOut

! ======================================================================================

subroutine wrap_update_hifrq_hist(this, bounds_clump )
subroutine wrap_update_hifrq_hist(this, bounds_clump, solarabs_inst)



! Arguments
class(hlm_fates_interface_type), intent(inout) :: this
type(bounds_type), intent(in) :: bounds_clump
type(solarabs_type) , intent(in) :: solarabs_inst

! locals
real(r8) :: dtime
integer :: nstep
logical :: is_beg_day
integer :: s,c,nc
integer :: ifp,p

associate(&
hr => col_cf%hr, & ! (gC/m2/s) total heterotrophic respiration
totsomc => col_cs%totsomc, & ! (gC/m2) total soil organic matter carbon
totlitc => col_cs%totlitc) ! (gC/m2) total litter carbon in BGC pools
totlitc => col_cs%totlitc, & ! (gC/m2) total litter carbon in BGC pools
eflx_lh_tot => veg_ef%eflx_lh_tot, & ! (W/m2) latent heat flux
eflx_sh_tot => veg_ef%eflx_sh_tot, & ! (W/m2) sensible heat flux
fsa_patch => solarabs_inst%fsa_patch, & ! (W/m2) absorbed solar flux
eflx_lwrad_net=> veg_ef%eflx_lwrad_net, & ! (W/m2) net longwave radiative flux
t_ref2m => veg_es%t_ref2m) ! (K) 2-m air temperature

nc = bounds_clump%clump_index
dtime = real(get_step_size(),r8)
Expand All @@ -3024,6 +3031,22 @@ subroutine wrap_update_hifrq_hist(this, bounds_clump )
this%fates(nc)%bc_in(s)%tot_litc = totlitc(c)
end do

! summarize biophysical variables that we want ot output on FATES dimensions

do s = 1, this%fates(nc)%nsites
c = this%f2hmap(nc)%fcolumn(s)
do ifp = 0, this%fates(nc)%sites(s)%youngest_patch%patchno !!!CDK was 1
p = ifp+col_pp%pfti(c)

this%fates(nc)%bc_in(s)%lhflux_pa(ifp) = eflx_lh_tot(p)
this%fates(nc)%bc_in(s)%shflux_pa(ifp) = eflx_sh_tot(p)
this%fates(nc)%bc_in(s)%swabs_pa(ifp) = fsa_patch(p)
this%fates(nc)%bc_in(s)%netlw_pa(ifp) = eflx_lwrad_net(p)
this%fates(nc)%bc_in(s)%t2m_pa(ifp) = t_ref2m(p)

end do
end do

! Update history variables that track these variables
call fates_hist%update_history_hifrq(nc, &
this%fates(nc)%nsites, &
Expand Down
Loading