'pull ctsm_n04_clm4_5_16_r249 tags from svn'

HPSCTerrSys · Dec 28, 2017 · 487111d · 487111d
1 parent 35a76a9
commit 487111d
Show file tree

Hide file tree

Showing 14 changed files with 533 additions and 237 deletions.
diff --git a/src/biogeochem/CNDriverMod.F90 b/src/biogeochem/CNDriverMod.F90
@@ -36,7 +36,7 @@ module CNDriverMod
   use PhotosynthesisMod               , only : photosyns_type
   use ch4Mod                          , only : ch4_type
   use EnergyFluxType                  , only : energyflux_type
-  use SoilHydrologyType               , only : soilhydrology_type
+  use SurfRunoffSatMod                , only : surf_runoff_sat_type
   !
   ! !PUBLIC TYPES:
   implicit none
@@ -91,7 +91,7 @@ subroutine CNDriverNoLeaching(bounds,
        soilbiogeochem_nitrogenflux_inst, soilbiogeochem_nitrogenstate_inst,                &
        atm2lnd_inst, waterstate_inst, waterflux_inst,                                      &
        canopystate_inst, soilstate_inst, temperature_inst, crop_inst, ch4_inst,            &
-       dgvs_inst, photosyns_inst, soilhydrology_inst, energyflux_inst,                     &
+       dgvs_inst, photosyns_inst, surf_runoff_sat_inst, energyflux_inst,                   &
        nutrient_competition_method, cnfire_method)
     !
     ! !DESCRIPTION:
@@ -175,7 +175,7 @@ subroutine CNDriverNoLeaching(bounds,
     type(ch4_type)                          , intent(in)    :: ch4_inst
     type(dgvs_type)                         , intent(inout) :: dgvs_inst
     type(photosyns_type)                    , intent(in)    :: photosyns_inst
-    type(soilhydrology_type)                , intent(in)    :: soilhydrology_inst
+    type(surf_runoff_sat_type)              , intent(in)    :: surf_runoff_sat_inst
     type(energyflux_type)                   , intent(in)    :: energyflux_inst
     class(nutrient_competition_method_type) , intent(inout) :: nutrient_competition_method
     class(cnfire_method_type)               , intent(inout) :: cnfire_method
@@ -737,7 +737,7 @@ subroutine CNDriverNoLeaching(bounds,
 
     call t_startf('CNFire')
     call cnfire_method%CNFireArea(bounds, num_soilc, filter_soilc, num_soilp, filter_soilp, &
-         atm2lnd_inst, energyflux_inst, soilhydrology_inst, waterstate_inst, &
+         atm2lnd_inst, energyflux_inst, surf_runoff_sat_inst, waterstate_inst, &
          cnveg_state_inst, cnveg_carbonstate_inst, &
          totlitc_col=soilbiogeochem_carbonstate_inst%totlitc_col(begc:endc), &
          decomp_cpools_vr_col=soilbiogeochem_carbonstate_inst%decomp_cpools_vr_col(begc:endc,1:nlevdecomp_full,1:ndecomp_pools), &

diff --git a/src/biogeochem/CNFireBaseMod.F90 b/src/biogeochem/CNFireBaseMod.F90
@@ -35,7 +35,7 @@ module CNFireBaseMod
   use CNVegNitrogenFluxType              , only : cnveg_nitrogenflux_type
   use SoilBiogeochemDecompCascadeConType , only : decomp_cascade_con
   use EnergyFluxType                     , only : energyflux_type
-  use SoilHydrologyType                  , only : soilhydrology_type  
+  use SurfRunoffSatMod                   , only : surf_runoff_sat_type
   use WaterstateType                     , only : waterstate_type
   use GridcellType                       , only : grc
   use ColumnType                         , only : col
@@ -263,7 +263,7 @@ end subroutine CNFireInterp
 
   !-----------------------------------------------------------------------
   subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_soilp, &
-       atm2lnd_inst, energyflux_inst, soilhydrology_inst, waterstate_inst, &
+       atm2lnd_inst, energyflux_inst, surf_runoff_sat_inst, waterstate_inst, &
        cnveg_state_inst, cnveg_carbonstate_inst, totlitc_col, decomp_cpools_vr_col, t_soi17cm_col)
     !
     ! !DESCRIPTION:
@@ -280,7 +280,7 @@ subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_
     integer                               , intent(in)    :: filter_soilp(:) ! filter for soil patches
     type(atm2lnd_type)                    , intent(in)    :: atm2lnd_inst
     type(energyflux_type)                 , intent(in)    :: energyflux_inst
-    type(soilhydrology_type)              , intent(in)    :: soilhydrology_inst
+    type(surf_runoff_sat_type)            , intent(in)    :: surf_runoff_sat_inst
     type(waterstate_type)                 , intent(in)    :: waterstate_inst
     type(cnveg_state_type)                , intent(inout) :: cnveg_state_inst
     type(cnveg_carbonstate_type)          , intent(inout) :: cnveg_carbonstate_inst

diff --git a/src/biogeochem/CNFireLi2014Mod.F90 b/src/biogeochem/CNFireLi2014Mod.F90
@@ -34,7 +34,7 @@ module CNFireLi2014Mod
   use CNVegNitrogenFluxType              , only : cnveg_nitrogenflux_type
   use SoilBiogeochemDecompCascadeConType , only : decomp_cascade_con
   use EnergyFluxType                     , only : energyflux_type
-  use SoilHydrologyType                  , only : soilhydrology_type  
+  use SurfRunoffSatMod                   , only : surf_runoff_sat_type
   use WaterstateType                     , only : waterstate_type
   use GridcellType                       , only : grc                
   use ColumnType                         , only : col                
@@ -83,7 +83,7 @@ end function constructor
 
   !-----------------------------------------------------------------------
   subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_soilp, &
-       atm2lnd_inst, energyflux_inst, soilhydrology_inst, waterstate_inst, &
+       atm2lnd_inst, energyflux_inst, surf_runoff_sat_inst, waterstate_inst, &
        cnveg_state_inst, cnveg_carbonstate_inst, totlitc_col, decomp_cpools_vr_col, t_soi17cm_col)
     !
     ! !DESCRIPTION:
@@ -106,7 +106,7 @@ subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_
     integer                               , intent(in)    :: filter_soilp(:) ! filter for soil patches
     type(atm2lnd_type)                    , intent(in)    :: atm2lnd_inst
     type(energyflux_type)                 , intent(in)    :: energyflux_inst
-    type(soilhydrology_type)              , intent(in)    :: soilhydrology_inst
+    type(surf_runoff_sat_type)            , intent(in)    :: surf_runoff_sat_inst
     type(waterstate_type)                 , intent(in)    :: waterstate_inst
     type(cnveg_state_type)                , intent(inout) :: cnveg_state_inst
     type(cnveg_carbonstate_type)          , intent(inout) :: cnveg_carbonstate_inst
@@ -164,7 +164,7 @@ subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_
          fd_pft             => pftcon%fd_pft                                   , & ! Input:
 
          btran2             => energyflux_inst%btran2_patch                    , & ! Input:  [real(r8) (:)     ]  root zone soil wetness                            
-         fsat               => soilhydrology_inst%fsat_col                     , & ! Input:  [real(r8) (:)     ]  fractional area with water table at surface       
+         fsat               => surf_runoff_sat_inst%fsat_col                   , & ! Input:  [real(r8) (:)     ]  fractional area with water table at surface       
          wf                 => waterstate_inst%wf_col                          , & ! Input:  [real(r8) (:)     ]  soil water as frac. of whc for top 0.05 m         
          wf2                => waterstate_inst%wf2_col                         , & ! Input:  [real(r8) (:)     ]  soil water as frac. of whc for top 0.17 m         
 

diff --git a/src/biogeochem/CNFireLi2016Mod.F90 b/src/biogeochem/CNFireLi2016Mod.F90
@@ -34,7 +34,7 @@ module CNFireLi2016Mod
   use CNVegNitrogenFluxType              , only : cnveg_nitrogenflux_type
   use SoilBiogeochemDecompCascadeConType , only : decomp_cascade_con
   use EnergyFluxType                     , only : energyflux_type
-  use SoilHydrologyType                  , only : soilhydrology_type  
+  use SurfRunoffSatMod                   , only : surf_runoff_sat_type
   use WaterstateType                     , only : waterstate_type
   use GridcellType                       , only : grc                
   use ColumnType                         , only : col                
@@ -84,7 +84,7 @@ end function constructor
 
   !-----------------------------------------------------------------------
   subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_soilp, &
-       atm2lnd_inst, energyflux_inst, soilhydrology_inst, waterstate_inst, &
+       atm2lnd_inst, energyflux_inst, surf_runoff_sat_inst, waterstate_inst, &
        cnveg_state_inst, cnveg_carbonstate_inst, totlitc_col, decomp_cpools_vr_col, t_soi17cm_col)
     !
     ! !DESCRIPTION:
@@ -108,7 +108,7 @@ subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_
     integer                               , intent(in)    :: filter_soilp(:) ! filter for soil patches
     type(atm2lnd_type)                    , intent(in)    :: atm2lnd_inst
     type(energyflux_type)                 , intent(in)    :: energyflux_inst
-    type(soilhydrology_type)              , intent(in)    :: soilhydrology_inst
+    type(surf_runoff_sat_type)            , intent(in)    :: surf_runoff_sat_inst
     type(waterstate_type)                 , intent(in)    :: waterstate_inst
     type(cnveg_state_type)                , intent(inout) :: cnveg_state_inst
     type(cnveg_carbonstate_type)          , intent(inout) :: cnveg_carbonstate_inst
@@ -171,7 +171,7 @@ subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_
          fd_pft             => pftcon%fd_pft                                   , & ! Input:
 
          btran2             => energyflux_inst%btran2_patch                    , & ! Input:  [real(r8) (:)     ]  root zone soil wetness                            
-         fsat               => soilhydrology_inst%fsat_col                     , & ! Input:  [real(r8) (:)     ]  fractional area with water table at surface       
+         fsat               => surf_runoff_sat_inst%fsat_col                   , & ! Input:  [real(r8) (:)     ]  fractional area with water table at surface       
          wf2                => waterstate_inst%wf2_col                         , & ! Input:  [real(r8) (:)     ]  soil water as frac. of whc for top 0.17 m         
 
          is_cwd             => decomp_cascade_con%is_cwd                       , & ! Input:  [logical  (:)     ]  TRUE => pool is a cwd pool                         

diff --git a/src/biogeochem/CNFireMethodMod.F90 b/src/biogeochem/CNFireMethodMod.F90
@@ -96,7 +96,7 @@ end subroutine CNFireInterp_interface
 
   !-----------------------------------------------------------------------
   subroutine CNFireArea_interface (this, bounds, num_soilc, filter_soilc, num_soilp, filter_soilp, &
-       atm2lnd_inst, energyflux_inst, soilhydrology_inst, waterstate_inst, &
+       atm2lnd_inst, energyflux_inst, surf_runoff_sat_inst, waterstate_inst, &
        cnveg_state_inst, cnveg_carbonstate_inst, totlitc_col, decomp_cpools_vr_col, t_soi17cm_col)
     !
     ! !DESCRIPTION:
@@ -107,7 +107,7 @@ subroutine CNFireArea_interface (this, bounds, num_soilc, filter_soilc, num_soil
     use decompMod                          , only : bounds_type
     use atm2lndType                        , only : atm2lnd_type
     use EnergyFluxType                     , only : energyflux_type
-    use SoilHydrologyType                  , only : soilhydrology_type
+    use SurfRunoffSatMod                   , only : surf_runoff_sat_type
     use WaterstateType                     , only : waterstate_type
     use CNVegStateType                     , only : cnveg_state_type
     use CNVegCarbonStateType               , only : cnveg_carbonstate_type
@@ -122,7 +122,7 @@ subroutine CNFireArea_interface (this, bounds, num_soilc, filter_soilc, num_soil
     integer                               , intent(in)    :: filter_soilp(:) !  filter for soil patches
     type(atm2lnd_type)                    , intent(in)    :: atm2lnd_inst
     type(energyflux_type)                 , intent(in)    :: energyflux_inst
-    type(soilhydrology_type)              , intent(in)    :: soilhydrology_inst
+    type(surf_runoff_sat_type)            , intent(in)    :: surf_runoff_sat_inst
     type(waterstate_type)                 , intent(in)    :: waterstate_inst
     type(cnveg_state_type)                , intent(inout) :: cnveg_state_inst
     type(cnveg_carbonstate_type)          , intent(inout) :: cnveg_carbonstate_inst

diff --git a/src/biogeochem/CNFireNoFireMod.F90 b/src/biogeochem/CNFireNoFireMod.F90
@@ -16,7 +16,7 @@ module CNFireNoFireMod
   use CNVegNitrogenStateType             , only : cnveg_nitrogenstate_type
   use CNVegNitrogenFluxType              , only : cnveg_nitrogenflux_type
   use EnergyFluxType                     , only : energyflux_type
-  use SoilHydrologyType                  , only : soilhydrology_type  
+  use SurfRunoffSatMod                   , only : surf_runoff_sat_type
   use WaterstateType                     , only : waterstate_type
   use CNFireMethodMod                    , only : cnfire_method_type
   use CNFireBaseMod                      , only : cnfire_base_type
@@ -58,7 +58,7 @@ end function constructor
 
   !-----------------------------------------------------------------------
   subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_soilp, &
-       atm2lnd_inst, energyflux_inst, soilhydrology_inst, waterstate_inst, &
+       atm2lnd_inst, energyflux_inst, surf_runoff_sat_inst, waterstate_inst, &
        cnveg_state_inst, cnveg_carbonstate_inst, totlitc_col, decomp_cpools_vr_col, t_soi17cm_col)
     !
     ! !DESCRIPTION:
@@ -76,7 +76,7 @@ subroutine CNFireArea (this, bounds, num_soilc, filter_soilc, num_soilp, filter_
     integer                               , intent(in)    :: filter_soilp(:) ! filter for soil patches
     type(atm2lnd_type)                    , intent(in)    :: atm2lnd_inst
     type(energyflux_type)                 , intent(in)    :: energyflux_inst
-    type(soilhydrology_type)              , intent(in)    :: soilhydrology_inst
+    type(surf_runoff_sat_type)            , intent(in)    :: surf_runoff_sat_inst
     type(waterstate_type)                 , intent(in)    :: waterstate_inst
     type(cnveg_state_type)                , intent(inout) :: cnveg_state_inst
     type(cnveg_carbonstate_type)          , intent(inout) :: cnveg_carbonstate_inst

diff --git a/src/biogeochem/CNVegetationFacade.F90 b/src/biogeochem/CNVegetationFacade.F90
@@ -72,7 +72,7 @@ module CNVegetationFacade
   use CNDVType                        , only : dgvs_type
   use CNDVDriverMod                   , only : CNDVDriver, CNDVHIST
   use EnergyFluxType                  , only : energyflux_type
-  use SoilHydrologyType               , only : soilhydrology_type
+  use SurfRunoffSatMod                , only : surf_runoff_sat_type
   use FrictionVelocityMod             , only : frictionvel_type
   use SoilBiogeochemStateType         , only : soilBiogeochem_state_type
   use SoilBiogeochemCarbonStateType   , only : soilbiogeochem_carbonstate_type
@@ -783,7 +783,7 @@ subroutine EcosystemDynamicsPreDrainage(this, bounds, &
        soilbiogeochem_nitrogenflux_inst, soilbiogeochem_nitrogenstate_inst,     &
        atm2lnd_inst, waterstate_inst, waterflux_inst,                           &
        canopystate_inst, soilstate_inst, temperature_inst, crop_inst, ch4_inst, &
-       photosyns_inst, soilhydrology_inst, energyflux_inst,          &
+       photosyns_inst, surf_runoff_sat_inst, energyflux_inst,          &
        nutrient_competition_method, fireemis_inst)
     !
     ! !DESCRIPTION:
@@ -822,7 +822,7 @@ subroutine EcosystemDynamicsPreDrainage(this, bounds, &
     type(crop_type)                         , intent(inout) :: crop_inst
     type(ch4_type)                          , intent(in)    :: ch4_inst
     type(photosyns_type)                    , intent(in)    :: photosyns_inst
-    type(soilhydrology_type)                , intent(in)    :: soilhydrology_inst
+    type(surf_runoff_sat_type)              , intent(in)    :: surf_runoff_sat_inst
     type(energyflux_type)                   , intent(in)    :: energyflux_inst
     class(nutrient_competition_method_type) , intent(inout) :: nutrient_competition_method
     type(fireemis_type)                     , intent(inout) :: fireemis_inst
@@ -852,7 +852,7 @@ subroutine EcosystemDynamicsPreDrainage(this, bounds, &
          soilbiogeochem_nitrogenflux_inst, soilbiogeochem_nitrogenstate_inst,     &
          atm2lnd_inst, waterstate_inst, waterflux_inst,                           &
          canopystate_inst, soilstate_inst, temperature_inst, crop_inst, ch4_inst, &
-         this%dgvs_inst, photosyns_inst, soilhydrology_inst, energyflux_inst,          &
+         this%dgvs_inst, photosyns_inst, surf_runoff_sat_inst, energyflux_inst,          &
          nutrient_competition_method, this%cnfire_method)
 
     ! fire carbon emissions 

diff --git a/src/biogeophys/HydrologyNoDrainageMod.F90 b/src/biogeophys/HydrologyNoDrainageMod.F90
@@ -14,8 +14,9 @@ Module HydrologyNoDrainageMod
   use AerosolMod        , only : aerosol_type
   use EnergyFluxType    , only : energyflux_type
   use TemperatureType   , only : temperature_type
-  use SoilHydrologyType , only : soilhydrology_type  
+  use SoilHydrologyType , only : soilhydrology_type
   use SoilStateType     , only : soilstate_type
+  use SurfRunoffSatMod  , only : surf_runoff_sat_type
   use WaterfluxType     , only : waterflux_type
   use WaterstateType    , only : waterstate_type
   use CanopyStateType   , only : canopystate_type
@@ -42,7 +43,7 @@ subroutine HydrologyNoDrainage(bounds, &
        clm_fates, &
        atm2lnd_inst, soilstate_inst, energyflux_inst, temperature_inst, &
        waterflux_inst, waterstate_inst, &
-       soilhydrology_inst, aerosol_inst, &
+       soilhydrology_inst, surf_runoff_sat_inst, aerosol_inst, &
        canopystate_inst, soil_water_retention_curve)
     !
     ! !DESCRIPTION:
@@ -59,8 +60,9 @@ subroutine HydrologyNoDrainage(bounds, &
     use clm_time_manager     , only : get_step_size, get_nstep
     use SnowHydrologyMod     , only : SnowCompaction, CombineSnowLayers, DivideSnowLayers, SnowCapping
     use SnowHydrologyMod     , only : SnowWater, BuildSnowFilter 
-    use SoilHydrologyMod     , only : CLMVICMap, SetFracIce, SurfaceRunoff
+    use SoilHydrologyMod     , only : CLMVICMap, SetFracIce
     use SoilHydrologyMod     , only : SetQflxTopSoil, Infiltration, TotalSurfaceRunoff
+    use SoilHydrologyMod     , only : UpdateUrbanPonding
     use SoilHydrologyMod     , only : WaterTable, PerchedWaterTable
     use SoilHydrologyMod     , only : ThetaBasedWaterTable, RenewCondensation
     use SoilWaterMovementMod , only : SoilWater 
@@ -90,6 +92,7 @@ subroutine HydrologyNoDrainage(bounds, &
     type(waterstate_type)    , intent(inout) :: waterstate_inst
     type(aerosol_type)       , intent(inout) :: aerosol_inst
     type(soilhydrology_type) , intent(inout) :: soilhydrology_inst
+    type(surf_runoff_sat_type), intent(inout) :: surf_runoff_sat_inst
     type(canopystate_type)   , intent(inout) :: canopystate_inst
     class(soil_water_retention_curve_type), intent(in) :: soil_water_retention_curve
     !
@@ -177,17 +180,22 @@ subroutine HydrologyNoDrainage(bounds, &
       call SetFracIce(bounds, num_hydrologyc, filter_hydrologyc, &
            soilhydrology_inst, soilstate_inst, waterstate_inst)
 
-      call SurfaceRunoff(bounds, num_hydrologyc, filter_hydrologyc, num_urbanc, filter_urbanc, &
-           soilhydrology_inst, soilstate_inst, waterflux_inst, waterstate_inst)
+      call surf_runoff_sat_inst%SaturatedSurfaceRunoff(&
+           bounds, num_hydrologyc, filter_hydrologyc, col, &
+           soilhydrology_inst, soilstate_inst, waterflux_inst)
 
       call SetQflxTopSoil(bounds, num_hydrologyc, filter_hydrologyc, waterflux_inst)
 
       call Infiltration(bounds, num_hydrologyc, filter_hydrologyc, num_urbanc, filter_urbanc,&
-           energyflux_inst, soilhydrology_inst, soilstate_inst, temperature_inst, &
-           waterflux_inst, waterstate_inst)
+           energyflux_inst, soilhydrology_inst, soilstate_inst, surf_runoff_sat_inst, &
+           temperature_inst, waterflux_inst, waterstate_inst)
 
       call TotalSurfaceRunoff(bounds, num_hydrologyc, filter_hydrologyc, &
-           waterflux_inst)
+           num_urbanc, filter_urbanc, &
+           waterflux_inst, soilhydrology_inst, surf_runoff_sat_inst, waterstate_inst)
+
+      call UpdateUrbanPonding(bounds, num_urbanc, filter_urbanc, &
+           waterstate_inst, soilhydrology_inst, waterflux_inst)
 
       call Compute_EffecRootFrac_And_VertTranSink(bounds, num_hydrologyc, &
            filter_hydrologyc, soilstate_inst, canopystate_inst, waterflux_inst, energyflux_inst)