From 5ab19b667f5067ba935b87069d3b140c77a8ecab Mon Sep 17 00:00:00 2001
From: Sebastian Stammler <sebastian.stammler@gmail.com>
Date: Thu, 17 Oct 2024 15:22:46 +0200
Subject: [PATCH] Bugfix: Out-of-bounds error

---
 dustpy/std/dust.f90 | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/dustpy/std/dust.f90 b/dustpy/std/dust.f90
index d7fe60b..83ca935 100644
--- 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