Merge pull request #439 from danieljprice/moddump-sink

(moddump_sink) correct displayed value of sink luminosity
danieljprice · Jun 29, 2023 · d9393a1 · d9393a1
2 parents b020e27 + d43597b
commit d9393a1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/utils/moddump_sink.f90 b/src/utils/moddump_sink.f90
@@ -77,7 +77,7 @@ subroutine modify_dump(npart,npartoftype,massoftype,xyzh,vxyzu)
     xyzmh_ptmass(1,isinkpart) = newx
     print*,'x-coordinate changed to ',xyzmh_ptmass(1,isinkpart)
 
-    Lnuc = xyzmh_ptmass(1,ilum)
+    Lnuc = xyzmh_ptmass(ilum,isinkpart)
     Lnuc_cgs = Lnuc * unit_energ / utime
     call prompt('Enter new sink heating luminosity in erg/s:',Lnuc_cgs,0.)
     xyzmh_ptmass(ilum,isinkpart) = Lnuc_cgs / unit_energ * utime