From a5085e6313f8cbe73bc81a7350709bca31e3eacf Mon Sep 17 00:00:00 2001
From: Eric Morway <emorway@usgs.gov>
Date: Thu, 10 Aug 2023 14:12:47 -0700
Subject: [PATCH] fixes the following comments:
 https://github.com/MODFLOW-USGS/modflow6/pull/1306#discussion_r1272274220
 https://github.com/MODFLOW-USGS/modflow6/pull/1306#discussion_r1272275972
 https://github.com/MODFLOW-USGS/modflow6/pull/1306#discussion_r1272318323

---
 src/Model/TransportModel/tsp1adv1.f90 | 15 ++-------------
 src/Model/TransportModel/tsp1ssm1.f90 |  2 +-
 2 files changed, 3 insertions(+), 14 deletions(-)

diff --git a/src/Model/TransportModel/tsp1adv1.f90 b/src/Model/TransportModel/tsp1adv1.f90
index 7e3b25bf1ed..ea431151e9f 100644
--- a/src/Model/TransportModel/tsp1adv1.f90
+++ b/src/Model/TransportModel/tsp1adv1.f90
@@ -18,9 +18,7 @@ module TspAdvModule
     integer(I4B), pointer :: iadvwt => null() !< advection scheme (0 up, 1 central, 2 tvd)
     integer(I4B), dimension(:), pointer, contiguous :: ibound => null() !< pointer to model ibound
     type(TspFmiType), pointer :: fmi => null() !< pointer to fmi object
-    real(DP), dimension(:), pointer, contiguous :: cpw => null() ! pointer to GWE heat capacity of water
-    real(DP), dimension(:), pointer, contiguous :: rhow => null() ! fixed density of water
-    real(DP), pointer :: eqnsclfac => null() !< governing equation scale factor; =1. for solute; =rhow*cpw for energy
+    real(DP), pointer :: eqnsclfac => null() !< governing equation scale factor; =1.0 for solute; =rhow*cpw for energy
 
   contains
 
@@ -114,14 +112,12 @@ end subroutine adv_df
   !!
   !!  Method to allocate and read static data for the ADV package.
   !<
-  subroutine adv_ar(this, dis, ibound, cpw, rhow)
+  subroutine adv_ar(this, dis, ibound)
     ! -- modules
     ! -- dummy
     class(TspAdvType) :: this
     class(DisBaseType), pointer, intent(in) :: dis
     integer(I4B), dimension(:), pointer, contiguous, intent(in) :: ibound
-    real(DP), dimension(:), pointer, contiguous, optional, intent(in) :: cpw
-    real(DP), dimension(:), pointer, contiguous, optional, intent(in) :: rhow
     ! -- local
     ! -- formats
 ! ------------------------------------------------------------------------------
@@ -130,10 +126,6 @@ subroutine adv_ar(this, dis, ibound, cpw, rhow)
     this%dis => dis
     this%ibound => ibound
     !
-    ! -- if part of a GWE simulation, need heat capacity(cpw) and density (rhow)
-    if (present(cpw)) this%cpw => cpw
-    if (present(rhow)) this%rhow => rhow
-    !
     ! -- Return
     return
   end subroutine adv_ar
@@ -165,7 +157,6 @@ subroutine adv_fc(this, nodes, matrix_sln, idxglo, cnew, rhs)
         if (this%dis%con%mask(ipos) == 0) cycle
         m = this%dis%con%ja(ipos)
         if (this%ibound(m) == 0) cycle
-!!        qnm = this%fmi%gwfflowja(ipos)
         qnm = this%fmi%gwfflowja(ipos) * this%eqnsclfac
         omega = this%adv_weight(this%iadvwt, ipos, n, m, qnm)
         call matrix_sln%add_value_pos(idxglo(ipos), qnm * (DONE - omega))
@@ -370,8 +361,6 @@ subroutine adv_da(this)
     !
     ! -- nullify pointers
     this%ibound => null()
-    nullify (this%cpw)
-    nullify (this%rhow)
     !
     ! -- Scalars
     call mem_deallocate(this%iadvwt)
diff --git a/src/Model/TransportModel/tsp1ssm1.f90 b/src/Model/TransportModel/tsp1ssm1.f90
index 4d3973b70ce..20c96c837d9 100644
--- a/src/Model/TransportModel/tsp1ssm1.f90
+++ b/src/Model/TransportModel/tsp1ssm1.f90
@@ -367,7 +367,7 @@ subroutine ssm_term(this, ipackage, ientry, rrate, rhsval, hcofval, &
     ! -- set requested values
     if (present(hcofval)) hcofval = hcoftmp
     if (present(rhsval)) rhsval = rhstmp
-    if (present(rrate)) rrate = (hcoftmp * ctmp - rhstmp)
+    if (present(rrate)) rrate = hcoftmp * ctmp - rhstmp
     if (present(cssm)) cssm = ctmp
     if (present(qssm)) qssm = qbnd
     !