Bugfix: Out-of-bounds error

stammler · Oct 17, 2024 · 5ab19b6 · 5ab19b6
1 parent 51a3bf1
commit 5ab19b6
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/dustpy/std/dust.f90 b/dustpy/std/dust.f90
@@ -905,9 +905,11 @@ subroutine jacobian_hydrodynamic_generator(area, D, r, ri, SigmaGas, v, A, B, C,
   A(:, :)    = 0.d0
   B(:, :)    = 0.d0
   C(:, :)    = 0.d0
+  Vinv(:)    = 0.d0
+  w(:)       = 1.d0
 
   ! Grid cell volumes and distances
-  do ir=1, Nr
+  do ir=1, Nr-1
       Vinv(ir) = twopi / area(ir)
       w(ir) = r(ir+1) - r(ir)
   end do