read shf options block; will need to circle back and test TS functionality later

Author: emorway-usgs

src/Model/GroundWaterEnergy/PbstBase.f90

@@ -112,7 +112,7 @@ subroutine read_options(this)
     character(len=MAXCHARLEN) :: fname
     logical :: isfound
     logical :: endOfBlock
-    logical(LGP) :: foundchildclassoption
+    logical(LGP) :: found
     integer(I4B) :: ierr
     ! -- formats
     character(len=*), parameter :: fmtts = &
@@ -153,11 +153,11 @@ subroutine read_options(this)
         case default
           !
           ! -- Check for child class options
-          call this%pbst_options(keyword, foundchildclassoption)
+          call this%pbst_options(keyword, found)
           !
           ! -- Defer to subtype to read the option;
           ! -- if the subtype can't handle it, report an error
-          if (.not. this%read_option(keyword)) then
+          if (.not. found) then
             write (errmsg, '(a,3(1x,a),a)') &
               'Unknown', trim(adjustl(this%packName)), "option '", &
               trim(keyword), "'."

src/Model/GroundWaterEnergy/gwe-shf.f90

@@ -108,16 +108,51 @@ subroutine shf_allocate_arrays(this)
   !! This routine overrides PbstBaseType%bnd_options
   !<
   subroutine shf_options(this, option, found)
-    
+    ! -- dummy
+    class(ShfType), intent(inout) :: this
+    character(len=*), intent(inout) :: option
+    logical, intent(inout) :: found
+    !
     found = .true.
-    foundgcclassoption = .false.
     select case (option)
-    case ('FLOW_PACKAGE_NAME')
-      call this%parser%GetStringCaps(this%flowpackagename)
-      write (this%iout, '(4x,a)') &
-        'THIS '//trim(adjustl(this%text))//' PACKAGE CORRESPONDS TO A GWF &
-        &PACKAGE WITH THE NAME '//trim(adjustl(this%flowpackagename))
-
+    case ('DENSITY_AIR')
+      this%rhoa = this%parser%GetDouble()
+      if (this%rhoa <= 0.0) then
+        write (errmsg, '(a)') 'Specified value for the density of &
+          &the atmosphere must be greater than 0.0.'
+        call store_error(errmsg)
+        call this%parser%StoreErrorUnit()
+      else
+        write (this%iout, '(4x,a,1pg15.6)') &
+          "The density of the atmosphere has been set to: ", this%rhoa
+      end if
+    case ('HEAT_CAPACITY_AIR')
+      this%cpa = this%parser%GetDouble()
+      if (this%cpa <= 0.0) then
+        write (errmsg, '(a)') 'Specified value for the heat capacity of &
+          &the atmosphere must be greater than 0.0.'
+        call store_error(errmsg)
+        call this%parser%StoreErrorUnit()
+      else
+        write (this%iout, '(4x,a,1pg15.6)') &
+          "The heat capacity of the atmosphere has been set to: ", this%cpa
+      end if
+    case ('DRAG_COEFFICIENT')
+      this%cd = this%parser%GetDouble()
+      if (this%cd <= 0.0) then
+        write (errmsg, '(a)') 'Specified value for the drag coefficient &
+          &must be greater than 0.0.'
+        call store_error(errmsg)
+        call this%parser%StoreErrorUnit()
+      else
+        write (this%iout, '(4x,a,1pg15.6)') &
+          "The heat capacity of the atmosphere has been set to: ", this%cpa
+      end if
+    case default
+      write (errmsg, '(a,a)') 'Unknown SHF option: ', trim(option)
+      call store_error(errmsg)
+      call this%parser%StoreErrorUnit()
+    end select
   end subroutine shf_options
 
   !> @brief Deallocate package memory