From e9f5120db05a800cf492ee4b1c08ce28f130f2db Mon Sep 17 00:00:00 2001 From: Stefan Poll Date: Thu, 3 Sep 2026 17:12:30 +0200 Subject: [PATCH] eCLM-ParFlow: Couple urban landunits - bugfix: pfl_psi/pfl_h2osoi_liq were distributed over filter_soilc, but eCLM run over filter_hydrologyc, leading to soilwater_parflow wrote nan into h2osoi_liq for pervious road - distribute over filter_nolakec testing col%hydrologically_active, the predicate filterMod uses for filter_hydrologyc, originally leading to undefined TWS --- src/clm5/main/clm_driver.F90 | 15 ++++++++------- 1 file changed, 8 insertions(+), 7 deletions(-) diff --git a/src/clm5/main/clm_driver.F90 b/src/clm5/main/clm_driver.F90 index b7f36c920d..befbba7a45 100644 --- a/src/clm5/main/clm_driver.F90 +++ b/src/clm5/main/clm_driver.F90 @@ -1306,13 +1306,14 @@ subroutine clm_drv_init(bounds, & end do #ifdef COUP_OAS_PFL - ! TSMP/bldsva/intf_oas3/clm3_5/mct/receive_fld_2pfl.F90 - ! COUP_OAS_PFL - do f = 1, num_soilc - c = filter_soilc(f) - g = col%gridcell(c) - pfl_psi(c,:) = atm2lnd_inst%pfl_psi_grc(g,:) - pfl_h2osoi_liq(c,:) = atm2lnd_inst%pfl_h2osoi_liq_grc(g,:) + ! Cover every column that soilwater_parflow will later overwrite. + do f = 1, num_nolakec + c = filter_nolakec(f) + if (col%hydrologically_active(c)) then + g = col%gridcell(c) + pfl_psi(c,:) = atm2lnd_inst%pfl_psi_grc(g,:) + pfl_h2osoi_liq(c,:) = atm2lnd_inst%pfl_h2osoi_liq_grc(g,:) + end if end do #endif end associate