Added initialization of pores and resid

physics/sfc_drv.f

@@ -4,7 +4,9 @@
 !> This module contains the CCPP-compliant Noah land surface scheme driver.
       module lsm_noah
 
+      use machine,          only: kind_phys
       use set_soilveg_mod,  only: set_soilveg
+      use namelist_soilveg
 
       implicit none
 
@@ -20,11 +22,14 @@ module lsm_noah
 !! \htmlinclude lsm_noah_init.html
 !!
       subroutine lsm_noah_init(me, isot, ivegsrc, nlunit,
-     &                         errmsg, errflg)
+     &                         pores, resid, errmsg, errflg)
 
       implicit none
 
       integer,              intent(in)  :: me, isot, ivegsrc, nlunit
+
+      real (kind=kind_phys), dimension(:), intent(out) :: pores, resid
+
       character(len=*),     intent(out) :: errmsg
       integer,              intent(out) :: errflg
 
@@ -48,6 +53,9 @@ subroutine lsm_noah_init(me, isot, ivegsrc, nlunit,
       !--- initialize soil vegetation
       call set_soilveg(me, isot, ivegsrc, nlunit)
 
+      pores (:) = maxsmc (:)
+      resid (:) = drysmc (:)
+
       end subroutine lsm_noah_init
 
 
@@ -199,7 +207,7 @@ subroutine lsm_noah_run                                           &
      &       smcwlt2, smcref2, wet1, errmsg, errflg                     &
      &     )
 !
-      use machine , only : kind_phys
+      !use machine , only : kind_phys
       use funcphys, only : fpvs
 
       use surface_perturbation, only : ppfbet

physics/sfc_drv.meta

@@ -39,6 +39,22 @@
   type = integer
   intent = in
   optional = F
+[pores]
+  standard_name = maximum_soil_moisture_content_for_land_surface_model
+  long_name = maximum soil moisture for a given soil type for land surface model
+  units = m
+  dimensions = (30)
+  type = real
+  intent = out
+  kind = kind_phys
+[resid]
+  standard_name = minimum_soil_moisture_content_for_land_surface_model
+  long_name = minimum soil moisture for a given soil type for land surface model
+  units = m
+  dimensions = (30)
+  type = real
+  intent = out
+  kind = kind_phys
 [errmsg]
   standard_name = ccpp_error_message
   long_name = error message for error handling in CCPP

physics/sfc_drv_ruc.F90

@@ -32,7 +32,7 @@ subroutine lsm_ruc_init (me, master, isot, ivegsrc, nlunit,        &
                                tsfc_lnd, tsfc_wat,                       & ! in
                                tg3, smc, slc, stc,                       & ! in
                                zs, sh2o, smfrkeep, tslb, smois, wetness, & ! out
-                               tsice, errmsg, errflg)
+                               tsice, pores, resid, errmsg, errflg)
 
       implicit none
 !  ---  in
@@ -65,6 +65,8 @@ subroutine lsm_ruc_init (me, master, isot, ivegsrc, nlunit,        &
       real (kind=kind_phys), dimension(im,lsoil_ruc), intent(inout) :: tslb, smois
       real (kind=kind_phys), dimension(im,kice),      intent(out) :: tsice
 
+      real (kind=kind_phys), dimension(:), intent(out) :: pores, resid
+
       character(len=*),     intent(out) :: errmsg
       integer,              intent(out) :: errflg
 
@@ -113,6 +115,9 @@ subroutine lsm_ruc_init (me, master, isot, ivegsrc, nlunit,        &
       !--- initialize soil vegetation
       call set_soilveg_ruc(me, isot, ivegsrc, nlunit)
 
+      pores (:) = maxsmc (:)
+      resid (:) = drysmc (:)
+
       soiltyp(:) = 0
       vegtype(:) = 0
 
@@ -138,7 +143,7 @@ subroutine lsm_ruc_init (me, master, isot, ivegsrc, nlunit,        &
         endif
       enddo
 
-      call init_soil_depth_3 ( zs , dzs , lsoil_ruc )
+        call init_soil_depth_3 ( zs , dzs , lsoil_ruc )
 
         call rucinit   (flag_restart, im, lsoil_ruc, lsoil, nlev,   & ! in
                         me, master, lsm_ruc, lsm, slmsk,            & ! in

physics/sfc_drv_ruc.meta

@@ -263,6 +263,22 @@
   kind = kind_phys
   intent = out
   optional = F
+[pores]
+  standard_name = maximum_soil_moisture_content_for_land_surface_model
+  long_name = maximum soil moisture for a given soil type for land surface model
+  units = m
+  dimensions = (30)
+  type = real
+  intent = out
+  kind = kind_phys
+[resid]
+  standard_name = minimum_soil_moisture_content_for_land_surface_model
+  long_name = minimum soil moisture for a given soil type for land surface model
+  units = m
+  dimensions = (30)
+  type = real
+  intent = out
+  kind = kind_phys
 [errmsg]
   standard_name = ccpp_error_message
   long_name = error message for error handling in CCPP

physics/sfc_noahmp_drv.f

@@ -24,14 +24,19 @@ module noahmpdrv
 !! \section arg_table_noahmpdrv_init Argument Table
 !! \htmlinclude noahmpdrv_init.html
 !!
-      subroutine noahmpdrv_init(me, isot, ivegsrc, nlunit, errmsg,      &
-     &                          errflg)
+      subroutine noahmpdrv_init(me, isot, ivegsrc, nlunit, pores, resid,
+     &                          errmsg, errflg)
 
+        use machine,          only: kind_phys
         use set_soilveg_mod,  only: set_soilveg
+        use namelist_soilveg
 
         implicit none
 
         integer,              intent(in)  :: me, isot, ivegsrc, nlunit
+
+        real (kind=kind_phys), dimension(:), intent(out) :: pores, resid
+
         character(len=*),     intent(out) :: errmsg
         integer,              intent(out) :: errflg
 
@@ -54,6 +59,9 @@ subroutine noahmpdrv_init(me, isot, ivegsrc, nlunit, errmsg,      &
 
         !--- initialize soil vegetation
         call set_soilveg(me, isot, ivegsrc, nlunit)
+
+        pores (:) = maxsmc (:)
+        resid (:) = drysmc (:)
 
       end subroutine noahmpdrv_init
 

physics/sfc_noahmp_drv.meta

@@ -39,6 +39,22 @@
   type = integer
   intent = in
   optional = F
+[pores]
+  standard_name = maximum_soil_moisture_content_for_land_surface_model
+  long_name = maximum soil moisture for a given soil type for land surface model
+  units = m
+  dimensions = (30)
+  type = real
+  intent = out
+  kind = kind_phys
+[resid]
+  standard_name = minimum_soil_moisture_content_for_land_surface_model
+  long_name = minimum soil moisture for a given soil type for land surface model
+  units = m
+  dimensions = (30)
+  type = real
+  intent = out
+  kind = kind_phys
 [errmsg]
   standard_name = ccpp_error_message
   long_name = error message for error handling in CCPP