Update t9rhofind calls to use new dummy argument bounds

starkiller-astro · Jul 1, 2024 · a344c13 · a344c13
1 parent 96626e7
commit a344c13
Show file tree

Hide file tree

Showing 5 changed files with 11 additions and 6 deletions.
diff --git 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
@@ -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
@@ -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
@@ -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
@@ -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