From a344c13aa4ff6e9d51ed23fb7de5088f28510e4b Mon Sep 17 00:00:00 2001
From: Austin Harris <jaharris87@gmail.com>
Date: Mon, 1 Jul 2024 16:45:14 -0400
Subject: [PATCH] Update t9rhofind calls to use new dummy argument bounds

---
 source/net.f90                | 3 ++-
 source/xnet_evolve.f90        | 2 +-
 source/xnet_integrate.f90     | 6 ++++--
 source/xnet_integrate_bdf.f90 | 3 ++-
 source/xnet_integrate_be.f90  | 3 ++-
 5 files changed, 11 insertions(+), 6 deletions(-)

diff --git a/source/net.f90 b/source/net.f90
index 59bd2b25..23b196f2 100644
--- a/source/net.f90
+++ b/source/net.f90
@@ -207,7 +207,8 @@ Program net
     Call read_thermo_file(thermo_file,thermo_desc,ierr)
 
     ! Determine thermodynamic conditions at tstart
-    Call t9rhofind(0,tstart,nstart,t9start,rhostart)
+    Call t9rhofind(0,tstart(zb_lo:zb_hi),nstart(zb_lo:zb_hi),t9start(zb_lo:zb_hi), &
+      & rhostart(zb_lo:zb_hi))
 
     ! Determine initial abundances
     Call load_initial_abundances(inab_file,abund_desc,ierr)
diff --git a/source/xnet_evolve.f90 b/source/xnet_evolve.f90
index aa837a4d..189d595b 100644
--- a/source/xnet_evolve.f90
+++ b/source/xnet_evolve.f90
@@ -65,7 +65,7 @@ Subroutine full_net(kstep)
     idiag0 = idiag
 
     ! Initialize trial time step abundances and conditions
-    Call t9rhofind(0,t,nt,t9,rho)
+    Call t9rhofind(0,t(zb_lo:zb_hi),nt(zb_lo:zb_hi),t9(zb_lo:zb_hi),rho(zb_lo:zb_hi))
     Do izb = zb_lo, zb_hi
       tdel_old(izb) = tdel(izb)
       tdel_next(izb) = tdel(izb)
diff --git a/source/xnet_integrate.f90 b/source/xnet_integrate.f90
index f93bd72e..bdcc94d2 100644
--- a/source/xnet_integrate.f90
+++ b/source/xnet_integrate.f90
@@ -169,7 +169,8 @@ Subroutine timestep(kstep,mask_in)
           tt(izb) = min(t(izb) + tdel(izb), tstop(izb))
         EndIf
       EndDo
-      Call t9rhofind(kstep,tt,ntt,t9t,rhot,mask_in = mask_profile)
+      Call t9rhofind(kstep,tt(zb_lo:zb_hi),ntt(zb_lo:zb_hi),t9t(zb_lo:zb_hi),rhot(zb_lo:zb_hi), &
+        & mask_in = mask_profile)
       Do izb = zb_lo, zb_hi
         If ( mask_profile(izb) .and. ntt(izb)-1 > nt(izb) ) Then
           Do j = nt(izb), ntt(izb)-1
@@ -199,7 +200,8 @@ Subroutine timestep(kstep,mask_in)
             tt(izb) = min(t(izb) + tdel(izb), tstop(izb))
           EndIf
         EndDo
-        Call t9rhofind(kstep,tt,ntt,t9t,rhot,mask_in = mask_profile)
+        Call t9rhofind(kstep,tt(zb_lo:zb_hi),ntt(zb_lo:zb_hi),t9t(zb_lo:zb_hi),rhot(zb_lo:zb_hi), &
+          & mask_in = mask_profile)
         Do izb = zb_lo, zb_hi
           If ( mask_profile(izb) .and. t9(izb) > 0.0 ) Then
             dtherm(izb) = 10.0*abs(t9t(izb)-t9(izb))/t9(izb) + abs(rhot(izb)-rho(izb))/rho(izb)
diff --git a/source/xnet_integrate_bdf.f90 b/source/xnet_integrate_bdf.f90
index b473020e..bf6dd650 100644
--- a/source/xnet_integrate_bdf.f90
+++ b/source/xnet_integrate_bdf.f90
@@ -680,7 +680,8 @@ Subroutine bdf_step(kstep)
         If ( iheat > 0 ) t9t(izb) = zt0(neq,0,izb)
       EndIf
     EndDo
-    If ( iheat == 0 ) Call t9rhofind(kstep,tt,ntt,t9t,rhot,mask_in = retry_ts(zb_lo:zb_hi))
+    If ( iheat == 0 ) Call t9rhofind(kstep,tt(zb_lo:zb_hi),ntt(zb_lo:zb_hi),t9t(zb_lo:zb_hi), &
+      rhot(zb_lo:zb_hi),mask_in = retry_ts(zb_lo:zb_hi))
 
     ! Initialize masks
     Do izb = zb_lo, zb_hi
diff --git a/source/xnet_integrate_be.f90 b/source/xnet_integrate_be.f90
index 8eb0faa7..6780bc37 100644
--- a/source/xnet_integrate_be.f90
+++ b/source/xnet_integrate_be.f90
@@ -98,7 +98,8 @@ Subroutine solve_be(kstep,its)
 
       ! Reset temperature and density for failed integrations
       lzstep = ( inr == 0 )
-      Call t9rhofind(kstep,tt,ntt,t9t,rhot,mask_in = lzstep)
+      Call t9rhofind(kstep,tt(zb_lo:zb_hi),ntt(zb_lo:zb_hi), &
+        & t9t(zb_lo:zb_hi),rhot(zb_lo:zb_hi),mask_in = lzstep)
       If ( iheat > 0 ) Then
         Do izb = zb_lo, zb_hi
           If ( inr(izb) == 0 ) Then