module HydrologyDrainageMod !----------------------------------------------------------------------- ! !DESCRIPTION: ! Calculates soil/snow hydrology with drainage (subsurface runoff) ! use shr_kind_mod , only : r8 => shr_kind_r8 use shr_log_mod , only : errMsg => shr_log_errMsg use decompMod , only : bounds_type use clm_varctl , only : iulog, use_vichydro use clm_varcon , only : denh2o, denice, rpi, spval use atm2lndType , only : atm2lnd_type use glc2lndMod , only : glc2lnd_type use SoilHydrologyType , only : soilhydrology_type use SoilStateType , only : soilstate_type use TemperatureType , only : temperature_type use WaterFluxBulkType , only : waterfluxbulk_type use Wateratm2lndBulkType , only : wateratm2lndbulk_type use WaterStateBulkType , only : waterstatebulk_type use WaterDiagnosticBulkType , only : waterdiagnosticbulk_type use WaterBalanceType , only : waterbalance_type use GlacierSurfaceMassBalanceMod, only : glacier_smb_type use TotalWaterAndHeatMod, only : ComputeWaterMassNonLake use LandunitType , only : lun use ColumnType , only : col ! ! !PUBLIC TYPES: implicit none save ! ! !PUBLIC MEMBER FUNCTIONS: public :: HydrologyDrainage ! Calculates soil/snow hydrolog with drainage !----------------------------------------------------------------------- contains !----------------------------------------------------------------------- subroutine HydrologyDrainage(bounds, & num_hillslope, filter_hillslopec, & num_nolakec, filter_nolakec, & num_hydrologyc, filter_hydrologyc, & num_urbanc, filter_urbanc, & num_do_smb_c, filter_do_smb_c, & glc2lnd_inst, temperature_inst, & soilhydrology_inst, soilstate_inst, waterstatebulk_inst, & waterdiagnosticbulk_inst, waterbalancebulk_inst, waterfluxbulk_inst, & wateratm2lndbulk_inst, glacier_smb_inst) ! ! !DESCRIPTION: ! Calculates soil/snow hydrology with drainage (subsurface runoff) ! ! !USES: use landunit_varcon , only : istwet, istsoil, istice_mec, istcrop use column_varcon , only : icol_roof, icol_road_imperv, icol_road_perv, icol_sunwall, icol_shadewall use clm_varcon , only : denh2o, denice use clm_varctl , only : use_vichydro, use_hillslope use clm_varpar , only : nlevgrnd, nlevurb use clm_time_manager , only : get_step_size_real, get_nstep use SoilHydrologyMod , only : CLMVICMap, Drainage, PerchedLateralFlow, PerchedLateralFlowHillslope, LateralFlowPowerLaw, LateralFlowHillslope use SoilWaterMovementMod , only : use_aquifer_layer ! ! !ARGUMENTS: type(bounds_type) , intent(in) :: bounds integer , intent(in) :: num_nolakec ! number of column non-lake points in column filter integer , intent(in) :: filter_nolakec(:) ! column filter for non-lake points integer , intent(in) :: num_hydrologyc ! number of column soil points in column filter integer , intent(in) :: filter_hydrologyc(:) ! column filter for soil points integer , intent(in) :: num_urbanc ! number of column urban points in column filter integer , intent(in) :: filter_urbanc(:) ! column filter for urban points integer , intent(in) :: num_do_smb_c ! number of bareland columns in which SMB is calculated, in column filter integer , intent(in) :: filter_do_smb_c(:) ! column filter for bare land SMB columns integer , intent(in) :: num_hillslope ! number of soil hill cols. integer , intent(in) :: filter_hillslopec(:) ! column filter for designating all hill cols. type(glc2lnd_type) , intent(in) :: glc2lnd_inst type(temperature_type) , intent(in) :: temperature_inst type(soilhydrology_type) , intent(inout) :: soilhydrology_inst type(soilstate_type) , intent(inout) :: soilstate_inst type(waterstatebulk_type) , intent(inout) :: waterstatebulk_inst type(waterdiagnosticbulk_type) , intent(inout) :: waterdiagnosticbulk_inst type(waterbalance_type) , intent(inout) :: waterbalancebulk_inst type(waterfluxbulk_type) , intent(inout) :: waterfluxbulk_inst type(wateratm2lndbulk_type) , intent(inout) :: wateratm2lndbulk_inst type(glacier_smb_type) , intent(in) :: glacier_smb_inst ! ! !LOCAL VARIABLES: integer :: g,l,c,j,fc ! indices real(r8) :: dtime ! land model time step (sec) !----------------------------------------------------------------------- associate( & ! Input: layer thickness depth (m) dz => col%dz , & ! Input: column type ctype => col%itype , & ! Input: gridcell flux of flood water from RTM qflx_floodg => wateratm2lndbulk_inst%forc_flood_grc , & ! Input: rain rate [mm/s] forc_rain => wateratm2lndbulk_inst%forc_rain_downscaled_col , & ! Input: snow rate [mm/s] forc_snow => wateratm2lndbulk_inst%forc_snow_downscaled_col , & ! Input: water mass begining of the time step begwb => waterbalancebulk_inst%begwb_col , & ! Output:water mass end of the time step endwb => waterbalancebulk_inst%endwb_col , & ! Output:water mass end of the time step h2osoi_ice => waterstatebulk_inst%h2osoi_ice_col , & ! Output: ice lens (kg/m2) h2osoi_liq => waterstatebulk_inst%h2osoi_liq_col , & ! Output: liquid water (kg/m2) h2osoi_vol => waterstatebulk_inst%h2osoi_vol_col , & ! Output: volumetric soil water ! (0<=h2osoi_vol<=watsat) [m3/m3] qflx_evap_tot => waterfluxbulk_inst%qflx_evap_tot_col , & ! Input: qflx_evap_soi + qflx_evap_can + qflx_tran_veg qflx_snwcp_ice => waterfluxbulk_inst%qflx_snwcp_ice_col , & ! Input: excess solid h2o due to snow ! capping (outgoing) (mm H2O /s) [+] qflx_snwcp_discarded_ice => waterfluxbulk_inst%qflx_snwcp_discarded_ice_col, & ! excess solid h2o due to snow capping, ! which we simply discard in order to reset ! the snow pack (mm H2O /s) [+] qflx_snwcp_discarded_liq => waterfluxbulk_inst%qflx_snwcp_discarded_liq_col, & ! excess liquid h2o due to snow capping, ! which we simply discard in order to reset ! the snow pack (mm H2O /s) [+] qflx_drain_perched => waterfluxbulk_inst%qflx_drain_perched_col , & ! sub-surface runoff from perched zwt (mm H2O /s) qflx_rsub_sat => waterfluxbulk_inst%qflx_rsub_sat_col , & ! soil saturation excess [mm h2o/s] qflx_drain => waterfluxbulk_inst%qflx_drain_col , & ! sub-surface runoff (mm H2O /s) qflx_surf => waterfluxbulk_inst%qflx_surf_col , & ! surface runoff (mm H2O /s) qflx_infl => waterfluxbulk_inst%qflx_infl_col , & ! infiltration (mm H2O /s) qflx_qrgwl => waterfluxbulk_inst%qflx_qrgwl_col , & ! qflx_surf at glaciers, wetlands, lakes qflx_latflow_out => waterfluxbulk_inst%qflx_latflow_out_col , & ! lateral subsurface flow qflx_runoff => waterfluxbulk_inst%qflx_runoff_col , & ! total runoff (qflx_drain+qflx_surf+qflx_qrgwl) (mm H2O /s) qflx_runoff_u => waterfluxbulk_inst%qflx_runoff_u_col , & ! Urban total runoff (qflx_drain+qflx_surf) (mm H2O /s) qflx_runoff_r => waterfluxbulk_inst%qflx_runoff_r_col , & ! Rural total runoff (qflx_drain+qflx_surf+qflx_qrgwl) (mm H2O /s) qflx_ice_runoff_snwcp => waterfluxbulk_inst%qflx_ice_runoff_snwcp_col, & ! solid runoff from snow capping (mm H2O /s) qflx_sfc_irrig => waterfluxbulk_inst%qflx_sfc_irrig_col & ! surface irrigation flux (mm H2O /s) ) ! Determine time step and step size dtime = get_step_size_real() if (use_vichydro) then call CLMVICMap(bounds, num_hydrologyc, filter_hydrologyc, & soilhydrology_inst, waterstatebulk_inst) endif if (use_aquifer_layer()) then call Drainage(bounds, num_hydrologyc, filter_hydrologyc, & num_urbanc, filter_urbanc,& temperature_inst, soilhydrology_inst, soilstate_inst, & waterstatebulk_inst, waterfluxbulk_inst) else if(use_hillslope) then call PerchedLateralFlowHillslope(bounds, & num_hydrologyc, filter_hydrologyc, & soilhydrology_inst, soilstate_inst, & waterstatebulk_inst, waterfluxbulk_inst,& wateratm2lndbulk_inst) call LateralFlowHillslope(bounds, & num_hillslope, filter_hillslopec, & num_hydrologyc, filter_hydrologyc, & num_urbanc, filter_urbanc,& soilhydrology_inst, soilstate_inst, & waterstatebulk_inst, waterfluxbulk_inst, & wateratm2lndbulk_inst) else call PerchedLateralFlow(bounds, num_hydrologyc, filter_hydrologyc, & num_urbanc, filter_urbanc,& soilhydrology_inst, soilstate_inst, & waterstatebulk_inst, waterfluxbulk_inst) call LateralFlowPowerLaw(bounds, num_hydrologyc, filter_hydrologyc, & num_urbanc, filter_urbanc,& soilhydrology_inst, soilstate_inst, & waterstatebulk_inst, waterfluxbulk_inst) endif endif do j = 1, nlevgrnd do fc = 1, num_nolakec c = filter_nolakec(fc) if (ctype(c) /= icol_sunwall .and. ctype(c) /= icol_shadewall & .and. ctype(c) /= icol_roof) then h2osoi_vol(c,j) = h2osoi_liq(c,j)/(dz(c,j)*denh2o) + h2osoi_ice(c,j)/(dz(c,j)*denice) end if end do end do do j = 1, nlevurb do fc = 1, num_urbanc c = filter_urbanc(fc) if (col%itype(c) == icol_sunwall .or. col%itype(c) == icol_shadewall & .or. col%itype(c) == icol_roof) then h2osoi_vol(c,j) = h2osoi_liq(c,j)/(dz(c,j)*denh2o) + h2osoi_ice(c,j)/(dz(c,j)*denice) end if end do end do call ComputeWaterMassNonLake(bounds, num_nolakec, filter_nolakec, & waterstatebulk_inst, waterdiagnosticbulk_inst, & subtract_dynbal_baselines = .false., & water_mass = endwb(bounds%begc:bounds%endc)) ! Determine wetland and land ice hydrology (must be placed here ! since need snow updated from CombineSnowLayers) do fc = 1,num_nolakec c = filter_nolakec(fc) l = col%landunit(c) g = col%gridcell(c) if (lun%itype(l)==istwet .or. lun%itype(l)==istice_mec) then qflx_latflow_out(c) = 0._r8 qflx_drain(c) = 0._r8 qflx_drain_perched(c) = 0._r8 qflx_surf(c) = 0._r8 qflx_infl(c) = 0._r8 qflx_qrgwl(c) = forc_rain(c) + forc_snow(c) + qflx_floodg(g) - qflx_evap_tot(c) - qflx_snwcp_ice(c) - & qflx_snwcp_discarded_ice(c) - qflx_snwcp_discarded_liq(c) - & (endwb(c)-begwb(c))/dtime else if (lun%urbpoi(l) .and. ctype(c) /= icol_road_perv) then qflx_drain_perched(c) = 0._r8 qflx_rsub_sat(c) = spval qflx_infl(c) = 0._r8 end if qflx_ice_runoff_snwcp(c) = qflx_snwcp_ice(c) end do ! This call needs to be here so that it comes after the initial calculation of ! qflx_qrgwl and qflx_ice_runoff_snwcp, but before the ues of qflx_qrgwl in ! qflx_runoff. call glacier_smb_inst%AdjustRunoffTerms(bounds, num_do_smb_c, filter_do_smb_c, & waterfluxbulk_inst, & glc2lnd_inst, & qflx_qrgwl = qflx_qrgwl(bounds%begc:bounds%endc), & qflx_ice_runoff_snwcp = qflx_ice_runoff_snwcp(bounds%begc:bounds%endc)) do fc = 1,num_nolakec c = filter_nolakec(fc) l = col%landunit(c) qflx_runoff(c) = qflx_drain(c) + qflx_surf(c) + qflx_qrgwl(c) + qflx_drain_perched(c) if ((lun%itype(l)==istsoil .or. lun%itype(l)==istcrop) .and. col%active(c)) then qflx_runoff(c) = qflx_runoff(c) - qflx_sfc_irrig(c) end if if (lun%urbpoi(l)) then qflx_runoff_u(c) = qflx_runoff(c) else if (lun%itype(l)==istsoil .or. lun%itype(l)==istcrop) then qflx_runoff_r(c) = qflx_runoff(c) end if end do end associate end subroutine HydrologyDrainage end module HydrologyDrainageMod