Groundwater irrigation

bld/CLMBuildNamelist.pm

@@ -2184,7 +2184,8 @@ sub setup_logic_crop {
         add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, "initial_seed_at_planting",
                     'use_crop'=>$nl->get_value('use_crop') );
      } else {
-        error_if_set( $nl, "Can NOT be set without crop on", "baset_mapping", "baset_latvary_slope", "baset_latvary_intercept" );
+	 error_if_set( $nl, "Can NOT be set without crop on", "baset_mapping", "baset_latvary_slope", "baset_latvary_intercept" );
+	 add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'crop_fsat_equals_zero' );
      }
   }
 }
@@ -2768,7 +2769,7 @@ sub setup_logic_irrigation_parameters {
      my $var;
      foreach $var ("irrig_min_lai", "irrig_start_time", "irrig_length",
                    "irrig_target_smp", "irrig_depth", "irrig_threshold_fraction",
-                   "limit_irrigation_if_rof_enabled") {
+                   "limit_irrigation_if_rof_enabled","use_groundwater_irrigation") {
         add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, $var);
      }
 

bld/namelist_files/namelist_defaults_clm4_5.xml

@@ -93,6 +93,9 @@ attributes from the config_cache.xml file (with keys converted to upper-case).
 <irrigate use_crop=".true." >.false.</irrigate>
 <irrigate use_crop=".false.">.true.</irrigate>
 
+<!-- Saturation excess runoff for crops  -->
+<crop_fsat_equals_zero>.false.</crop_fsat_equals_zero>
+
 <!-- MEGAN model -->
 <megan clm_accelerated_spinup="on" >0</megan>
 <megan clm_accelerated_spinup="off">1</megan>
@@ -273,6 +276,9 @@ attributes from the config_cache.xml file (with keys converted to upper-case).
 <limit_irrigation_if_rof_enabled phys="clm5_0">.false.</limit_irrigation_if_rof_enabled>
 <limit_irrigation_if_rof_enabled              >.false.</limit_irrigation_if_rof_enabled>
 
+<use_groundwater_irrigation>.false.</use_groundwater_irrigation>
+
+
 <!-- Snow veg treatment -->
 <snowveg_flag phys="clm4_5" >OFF</snowveg_flag>
 <snowveg_flag phys="clm5_0" >ON_RAD</snowveg_flag>

bld/namelist_files/namelist_definition_clm4_5.xml

@@ -437,11 +437,22 @@ Only applies if using an active runoff (ROF) model; otherwise, river storage-bas
 is turned off regardless of the setting of this namelist variable.
 </entry>
 
+<entry id="use_groundwater_irrigation" type="logical" category="clm_physics"
+       group="irrigation_inparm" valid_values="" >
+If TRUE, supply irrigation from groundwater (in addition to surface water).
+
+</entry>
+
 <entry id="irrigate" type="logical" category="clm_physics"
        group="clm_inparm"  >
 If TRUE, irrigation will be active.
 </entry> 
 
+<entry id="crop_fsat_equals_zero" type="logical" category="clm_physics"
+       group="clm_inparm"  >
+If TRUE, fsat will be set to zero for crop columns.
+</entry> 
+
 <entry id="maxpatch_glcmec" type="integer"  category="clm_physics"
        group="clm_inparm" valid_values="1,3,5,10,36" >
 Number of  multiple elevation classes over glacier points.

src/biogeophys/BalanceCheckMod.F90

@@ -55,6 +55,11 @@ subroutine BeginWaterBalance(bounds, &
     !
     ! !DESCRIPTION:
     ! Initialize column-level water balance at beginning of time step
+
+    !
+    ! !USES:
+    use SoilWaterMovementMod , only : use_aquifer_layer
+
     !
     ! !ARGUMENTS:
     type(bounds_type)         , intent(in)    :: bounds     
@@ -78,14 +83,16 @@ subroutine BeginWaterBalance(bounds, &
          begwb        =>    waterbalancebulk_inst%begwb_col        & ! Output: [real(r8) (:)   ]  water mass begining of the time step    
          )
 
-   do fc = 1, num_nolakec
-       c = filter_nolakec(fc)
-       if (col%hydrologically_active(c)) then
-          if(zwt(c) <= zi(c,nlevsoi)) then
-             wa(c) = aquifer_water_baseline
-          end if
-       end if
-    end do
+      if(use_aquifer_layer()) then 
+         do fc = 1, num_nolakec
+            c = filter_nolakec(fc)
+            if (col%hydrologically_active(c)) then
+               if(zwt(c) <= zi(c,nlevsoi)) then
+                  wa(c) = aquifer_water_baseline
+               end if
+            end if
+         end do
+      endif
 
     call ComputeWaterMassNonLake(bounds, num_nolakec, filter_nolakec, &
          waterstatebulk_inst, waterdiagnosticbulk_inst, begwb(bounds%begc:bounds%endc))
@@ -199,7 +206,11 @@ subroutine BalanceCheck( bounds, &
           snow_sources            =>    waterfluxbulk_inst%snow_sources_col         , & ! Output: [real(r8) (:)   ]  snow sources (mm H2O /s)  
           snow_sinks              =>    waterfluxbulk_inst%snow_sinks_col           , & ! Output: [real(r8) (:)   ]  snow sinks (mm H2O /s)    
 
-          qflx_irrig              =>    waterfluxbulk_inst%qflx_irrig_col          , & ! Input:  [real(r8) (:)   ]  irrigation flux (mm H2O /s)             
+          qflx_sfc_irrig          =>    waterfluxbulk_inst%qflx_sfc_irrig_col       , & ! Input:  [real(r8) (:)   ]  irrigation flux (mm H2O /s)             
+          qflx_gw_uncon_irrig     =>    waterfluxbulk_inst%qflx_gw_uncon_irrig_col  , & ! Input:  [real(r8) (:)   ]  groundwater irrigation flux (mm H2O /s)
+          qflx_gw_con_irrig       =>    waterfluxbulk_inst%qflx_gw_con_irrig_col    , & ! Input:  [real(r8) (:)   ]  groundwater irrigation flux (mm H2O /s)
+          qflx_irrig_drip         =>    waterfluxbulk_inst%qflx_irrig_drip_col      , & ! Input:  [real(r8) (:)   ]  drip irrigation flux (mm H2O /s)
+          qflx_irrig_sprinkler    =>    waterfluxbulk_inst%qflx_irrig_sprinkler_col , & ! Input:  [real(r8) (:)   ]  sprinkler irrigation flux (mm H2O /s)
 
           qflx_glcice_dyn_water_flux => waterfluxbulk_inst%qflx_glcice_dyn_water_flux_col, & ! Input: [real(r8) (:)]  water flux needed for balance check due to glc_dyn_runoff_routing (mm H2O/s) (positive means addition of water to the system)
 
@@ -272,7 +283,10 @@ subroutine BalanceCheck( bounds, &
                   - (forc_rain_col(c)        &
                   + forc_snow_col(c)         &
                   + qflx_floodc(c)           &
-                  + qflx_irrig(c)            &
+                  + qflx_irrig_drip(c)       &
+                  + qflx_irrig_sprinkler(c)  &
+                  - qflx_gw_uncon_irrig(c)   &
+                  - qflx_gw_con_irrig(c)     &
                   + qflx_glcice_dyn_water_flux(c) &
                   - qflx_evap_tot(c)         &
                   - qflx_surf(c)             &
@@ -321,7 +335,11 @@ subroutine BalanceCheck( bounds, &
              write(iulog,*)'endwb                 = ',endwb(indexc)
              write(iulog,*)'begwb                 = ',begwb(indexc)
              write(iulog,*)'qflx_evap_tot         = ',qflx_evap_tot(indexc)*dtime
-             write(iulog,*)'qflx_irrig            = ',qflx_irrig(indexc)*dtime
+             write(iulog,*)'qflx_irrig_drip       = ',qflx_irrig_drip(indexc)*dtime
+             write(iulog,*)'qflx_irrig_sprinkler  = ',qflx_irrig_sprinkler(indexc)*dtime
+             write(iulog,*)'qflx_sfc_irrig        = ',qflx_sfc_irrig(indexc)*dtime
+             write(iulog,*)'qflx_gw_uncon_irrig   = ',qflx_gw_uncon_irrig(indexc)*dtime
+             write(iulog,*)'qflx_gw_con_irrig     = ',qflx_gw_con_irrig(indexc)*dtime
              write(iulog,*)'qflx_surf             = ',qflx_surf(indexc)*dtime
              write(iulog,*)'qflx_qrgwl            = ',qflx_qrgwl(indexc)*dtime
              write(iulog,*)'qflx_drain            = ',qflx_drain(indexc)*dtime
@@ -351,7 +369,11 @@ subroutine BalanceCheck( bounds, &
              write(iulog,*)'begwb                 = ',begwb(indexc)
 
              write(iulog,*)'qflx_evap_tot         = ',qflx_evap_tot(indexc)*dtime
-             write(iulog,*)'qflx_irrig            = ',qflx_irrig(indexc)*dtime
+             write(iulog,*)'qflx_irrig_drip       = ',qflx_irrig_drip(indexc)*dtime
+             write(iulog,*)'qflx_irrig_sprinkler  = ',qflx_irrig_sprinkler(indexc)*dtime
+             write(iulog,*)'qflx_sfc_irrig        = ',qflx_sfc_irrig(indexc)*dtime
+             write(iulog,*)'qflx_gw_uncon_irrig   = ',qflx_gw_uncon_irrig(indexc)*dtime
+             write(iulog,*)'qflx_gw_con_irrig     = ',qflx_gw_con_irrig(indexc)*dtime
              write(iulog,*)'qflx_surf             = ',qflx_surf(indexc)*dtime
              write(iulog,*)'qflx_qrgwl            = ',qflx_qrgwl(indexc)*dtime
              write(iulog,*)'qflx_drain            = ',qflx_drain(indexc)*dtime

src/biogeophys/CanopyHydrologyMod.F90

@@ -21,10 +21,10 @@ module CanopyHydrologyMod
   use AerosolMod      , only : aerosol_type
   use CanopyStateType , only : canopystate_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 WaterFluxBulkType       , only : waterfluxbulk_type
+  use Wateratm2lndBulkType    , only : wateratm2lndbulk_type
+  use WaterStateBulkType      , only : waterstatebulk_type
+  use WaterDiagnosticBulkType , only : waterdiagnosticbulk_type
   use ColumnType      , only : col                
   use PatchType       , only : patch                
   !
@@ -176,10 +176,10 @@ subroutine CanopyHydrology(bounds, &
      type(canopystate_type) , intent(in)    :: canopystate_inst
      type(temperature_type) , intent(inout) :: temperature_inst
      type(aerosol_type)     , intent(inout) :: aerosol_inst
-     type(waterstatebulk_type)  , intent(inout) :: waterstatebulk_inst
-     type(waterdiagnosticbulk_type)  , intent(inout) :: waterdiagnosticbulk_inst
-     type(waterfluxbulk_type)   , intent(inout) :: waterfluxbulk_inst
-     type(wateratm2lndbulk_type)   , intent(inout) :: wateratm2lndbulk_inst
+     type(waterstatebulk_type)      , intent(inout) :: waterstatebulk_inst
+     type(waterdiagnosticbulk_type) , intent(inout) :: waterdiagnosticbulk_inst
+     type(waterfluxbulk_type)       , intent(inout) :: waterfluxbulk_inst
+     type(wateratm2lndbulk_type)    , intent(inout) :: wateratm2lndbulk_inst
      !
      ! !LOCAL VARIABLES:
      integer  :: f                                            ! filter index
@@ -209,6 +209,7 @@ subroutine CanopyHydrology(bounds, &
      real(r8) :: qflx_through_snow(bounds%begp:bounds%endp)   ! direct snow throughfall [mm/s]
      real(r8) :: qflx_prec_grnd_snow(bounds%begp:bounds%endp) ! snow precipitation incident on ground [mm/s]
      real(r8) :: qflx_prec_grnd_rain(bounds%begp:bounds%endp) ! rain precipitation incident on ground [mm/s]
+     real(r8) :: qflx_liq_above_canopy(bounds%begp:bounds%endp) ! liquid water input above canopy (rain plus irrigation) [mm/s]
      real(r8) :: z_avg                                        ! grid cell average snow depth
      real(r8) :: rho_avg                                      ! avg density of snow column
      real(r8) :: temp_snow_depth,temp_intsnow                 ! temporary variables
@@ -267,13 +268,14 @@ subroutine CanopyHydrology(bounds, &
           qflx_prec_intr       => waterfluxbulk_inst%qflx_prec_intr_patch     , & ! Output: [real(r8) (:)   ]  interception of precipitation [mm/s]    
           qflx_prec_grnd       => waterfluxbulk_inst%qflx_prec_grnd_patch     , & ! Output: [real(r8) (:)   ]  water onto ground including canopy runoff [kg/(m2 s)]
           qflx_rain_grnd       => waterfluxbulk_inst%qflx_rain_grnd_patch     , & ! Output: [real(r8) (:)   ]  rain on ground after interception (mm H2O/s) [+]
-          qflx_irrig           => waterfluxbulk_inst%qflx_irrig_patch         , & ! Input: [real(r8) (:)    ]  irrigation amount (mm/s)           
+          qflx_irrig_drip      => waterfluxbulk_inst%qflx_irrig_drip_patch      , & ! Input: [real(r8) (:)    ]  drip irrigation amount (mm/s)           
+          qflx_irrig_sprinkler => waterfluxbulk_inst%qflx_irrig_sprinkler_patch , & ! Input: [real(r8) (:)    ]  sprinkler irrigation amount (mm/s)
+
           qflx_snowindunload   => waterfluxbulk_inst%qflx_snowindunload_patch , & ! Output: [real(r8) (:)   ]  canopy snow unloading from wind [mm/s]    
           qflx_snotempunload   => waterfluxbulk_inst%qflx_snotempunload_patch   & ! Output: [real(r8) (:)   ]  canopy snow unloading from temp. [mm/s]    
           )
 
        ! Compute time step
-
        dtime = get_step_size()
 
        ! Set status of snowveg_flag
@@ -287,7 +289,7 @@ subroutine CanopyHydrology(bounds, &
           g = patch%gridcell(p)
           l = patch%landunit(p)
           c = patch%column(p)
-
+          
           ! Canopy interception and precipitation onto ground surface
           ! Add precipitation to leaf water
 
@@ -308,12 +310,14 @@ subroutine CanopyHydrology(bounds, &
 
 
              if (col%itype(c) /= icol_sunwall .and. col%itype(c) /= icol_shadewall) then
+                if (frac_veg_nosno(p) == 1 .and. (forc_rain(c) + forc_snow(c) + qflx_irrig_sprinkler(p)) > 0._r8) then
 
-                if (frac_veg_nosno(p) == 1 .and. (forc_rain(c) + forc_snow(c)) > 0._r8) then
-
+                   ! total liquid water inputs above canopy
+                   qflx_liq_above_canopy(p) = forc_rain(c)+ qflx_irrig_sprinkler(p)
+
                    ! determine fraction of input precipitation that is snow and rain
-                   fracsnow(p) = forc_snow(c)/(forc_snow(c) + forc_rain(c))
-                   fracrain(p) = forc_rain(c)/(forc_snow(c) + forc_rain(c))
+                   fracsnow(p) = forc_snow(c)/(forc_snow(c) + qflx_liq_above_canopy(p))
+                   fracrain(p) = forc_rain(c)/(forc_snow(c) + qflx_liq_above_canopy(p))
 
                    ! The leaf water capacities for solid and liquid are different,
                    ! generally double for snow, but these are of somewhat less
@@ -339,21 +343,20 @@ subroutine CanopyHydrology(bounds, &
                       ! Direct throughfall
                       qflx_through_snow(p) = forc_snow(c) * (1._r8-fpi)
                    end if
-
                    ! Direct throughfall
-                   qflx_through_rain(p) = forc_rain(c) * (1._r8-fpi)
+                   qflx_through_rain(p) = qflx_liq_above_canopy(p) * (1._r8-fpi)
 
                    if (snowveg_on .or. snowveg_onrad) then
                       ! Intercepted precipitation [mm/s]
-                      qflx_prec_intr(p) = forc_snow(c)*fpisnow + forc_rain(c)*fpi
+                      qflx_prec_intr(p) = forc_snow(c)*fpisnow + (forc_rain(c) + qflx_irrig_sprinkler(p))*fpi
                       ! storage of intercepted snowfall, rain, and dew
                       snocan(p) = max(0._r8, snocan(p) + dtime*forc_snow(c)*fpisnow)    
-                      liqcan(p) = max(0._r8, liqcan(p) + dtime*forc_rain(c)*fpi)
+                      liqcan(p) = max(0._r8, liqcan(p) + dtime*qflx_liq_above_canopy(p)*fpi)
                    else
                       ! Intercepted precipitation [mm/s]
-                      qflx_prec_intr(p) = (forc_snow(c) + forc_rain(c)) * fpi
+                      qflx_prec_intr(p) = (forc_snow(c) + qflx_liq_above_canopy(p)) * fpi
                    end if
-
+                   
                    ! Water storage of intercepted precipitation and dew
                    h2ocan(p) = max(0._r8, h2ocan(p) + dtime*qflx_prec_intr(p))
 
@@ -419,7 +422,7 @@ subroutine CanopyHydrology(bounds, &
           if (col%itype(c) /= icol_sunwall .and. col%itype(c) /= icol_shadewall) then
              if (frac_veg_nosno(p) == 0) then
                 qflx_prec_grnd_snow(p) = forc_snow(c)
-                qflx_prec_grnd_rain(p) = forc_rain(c)
+                qflx_prec_grnd_rain(p) = forc_rain(c) + qflx_irrig_sprinkler(p)
              else
                 if (snowveg_on .or. snowveg_onrad) then
                    qflx_snowindunload(p)=0._r8 
@@ -452,8 +455,8 @@ subroutine CanopyHydrology(bounds, &
 
           ! Add irrigation water directly onto ground (bypassing canopy interception)
           ! Note that it's still possible that (some of) this irrigation water will runoff (as runoff is computed later)
-          qflx_prec_grnd_rain(p) = qflx_prec_grnd_rain(p) + qflx_irrig(p)
-
+          qflx_prec_grnd_rain(p) = qflx_prec_grnd_rain(p) + qflx_irrig_drip(p)
+             
           qflx_prec_grnd(p) = qflx_prec_grnd_snow(p) + qflx_prec_grnd_rain(p)
 
           qflx_snow_grnd_patch(p) = qflx_prec_grnd_snow(p)           ! ice onto ground (mm/s)