Skip to content

Commit

Permalink
More prints
Browse files Browse the repository at this point in the history
  • Loading branch information
billsacks committed Apr 4, 2024
1 parent eff1138 commit 6c82971
Showing 1 changed file with 5 additions and 1 deletion.
6 changes: 5 additions & 1 deletion src/biogeophys/CanopyHydrologyMod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -549,6 +549,10 @@ subroutine BulkFlux_CanopyInterceptionAndThroughfall(bounds, num_nolakep, filter

! Canopy interception
qflx_intercepted_snow(p) = forc_snow(p) * fpisnow
if (p == 24) then
write(iulog,*) 'qflx_intercepted_snow, forc_snow, fpisnow, elai, esai = ', &
qflx_intercepted_snow(p), forc_snow(p), fpisnow, elai(p), esai(p)
end if
qflx_intercepted_liq(p) = qflx_liq_above_canopy(p) * fpiliq

else
Expand Down Expand Up @@ -704,7 +708,7 @@ subroutine UpdateState_AddInterceptionToCanopy(bounds, num_soilp, filter_soilp,

snocan(p) = max(0._r8, snocan(p) + dtime * qflx_intercepted_snow(p))
if (snocan(p) > 0 .and. snocan(p) < 1.e-290_r8) then
write(iulog,*) 'WJS: tiny snocan: AddInterceptionToCanopy: p, snocan = ', p, snocan(p)
write(iulog,*) 'WJS: tiny snocan: AddInterceptionToCanopy: p, snocan, qflx_intercepted_snow = ', p, snocan(p), qflx_intercepted_snow(p)
end if
liqcan(p) = max(0._r8, liqcan(p) + dtime * qflx_intercepted_liq(p))
end do
Expand Down

0 comments on commit 6c82971

Please sign in to comment.