Skip to content
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

Changes the units of water flux between root and soil in ALM #1409

Merged
merged 1 commit into from
Apr 21, 2017
Merged

Changes the units of water flux between root and soil in ALM #1409

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/clm/src/betr/BetrBGCMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -717,7 +717,7 @@ subroutine do_tracer_advection(bounds, lbj, ubj, jtops, num_soilc, filter_soilc,
qflx_adv_local(c,jtops(c)-1) = safe_div(qflx_adv(c,jtops(c)-1),aqu2bulkcef_mobile_col(c,jtops(c),j),eps=loc_eps)
do l = jtops(c), ubj
qflx_adv_local(c,l) = safe_div(qflx_adv(c,l),aqu2bulkcef_mobile_col(c,l,j),eps=loc_eps)
qflx_rootsoi_local(c,l) = safe_div(qflx_rootsoi(c,l),aqu2bulkcef_mobile_col(c,l,j),eps=loc_eps)
qflx_rootsoi_local(c,l) = safe_div(qflx_rootsoi(c,l)*1.e-3_r8,aqu2bulkcef_mobile_col(c,l,j),eps=loc_eps)
enddo
enddo

Expand Down
4 changes: 2 additions & 2 deletions components/clm/src/betr/betr_core/TracerParamsMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1368,7 +1368,7 @@ subroutine diagnose_advect_water_flux(bounds, num_hydrologyc, filter_hydrologyc,
if(j==nlevsoi)then
qflx_adv(c,j) = qcharge(c) * 1.e-3_r8
else
qflx_adv(c,j) = 1.e-3_r8 * (h2osoi_liq(c,j+1)-h2osoi_liq_copy(c,j+1))/dtime + qflx_adv(c,j+1) + qflx_rootsoi(c,j+1)
qflx_adv(c,j) = 1.e-3_r8 * (h2osoi_liq(c,j+1)-h2osoi_liq_copy(c,j+1))/dtime + qflx_adv(c,j+1) + qflx_rootsoi(c,j+1)*1.e-3_r8
endif

enddo
Expand All @@ -1380,7 +1380,7 @@ subroutine diagnose_advect_water_flux(bounds, num_hydrologyc, filter_hydrologyc,
c = filter_hydrologyc(fc)

!obtain the corrected infiltration
qflx_infl(c) = (h2osoi_liq(c,1)-h2osoi_liq_copy(c,1))/dtime + (qflx_rootsoi(c,1)+qflx_adv(c,1))*1.e3_r8
qflx_infl(c) = (h2osoi_liq(c,1)-h2osoi_liq_copy(c,1))/dtime + qflx_rootsoi(c,1) + qflx_adv(c,1)*1.e3_r8
!the predicted net infiltration
infl_tmp=qflx_gross_infl_soil(c)-qflx_gross_evap_soil(c)
diff=qflx_infl(c)-infl_tmp
Expand Down
4 changes: 2 additions & 2 deletions components/clm/src/biogeophys/SoilWaterMovementMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -355,7 +355,7 @@ subroutine soilwater_zengdecker2009(bounds, num_hydrologyc, filter_hydrologyc, &
qflx_infl => waterflux_vars%qflx_infl_col , & ! Input: [real(r8) (:) ] infiltration (mm H2O /s)
qflx_tran_veg_col => waterflux_vars%qflx_tran_veg_col , & ! Input: [real(r8) (:) ] vegetation transpiration (mm H2O/s) (+ = to atm)
qflx_tran_veg_pft => waterflux_vars%qflx_tran_veg_patch , & ! Input: [real(r8) (:) ] vegetation transpiration (mm H2O/s) (+ = to atm)
qflx_rootsoi => waterflux_vars%qflx_rootsoi_col , & ! Output: [real(r8) (:,:) ] vegetation/soil water exchange (m H2O/s) (+ = to atm)
qflx_rootsoi => waterflux_vars%qflx_rootsoi_col , & ! Output: [real(r8) (:,:) ] vegetation/soil water exchange (mm H2O/s) (+ = to atm)

t_soisno => temperature_vars%t_soisno_col & ! Input: [real(r8) (:,:) ] soil temperature (Kelvin)
)
Expand Down Expand Up @@ -786,7 +786,7 @@ subroutine soilwater_zengdecker2009(bounds, num_hydrologyc, filter_hydrologyc, &
do j = 1, nlevsoi
do fc = 1, num_hydrologyc
c = filter_hydrologyc(fc)
qflx_rootsoi(c,j) = qflx_tran_veg_col(c) * rootr_col(c,j) * 1.e-3_r8 ![m H2O/s]
qflx_rootsoi(c,j) = qflx_tran_veg_col(c) * rootr_col(c,j)
enddo
enddo

Expand Down