diff --git a/src/gray_core.f90 b/src/gray_core.f90
index 93009a7..5f3468b 100644
--- a/src/gray_core.f90
+++ b/src/gray_core.f90
@@ -553,7 +553,7 @@ contains
         write(msg, '(3x,a,g0.3," MW")') 'absoption: P=', pabs_beam
         call log_info(msg, mod='gray_core', proc='gray_main')
 
-        write(msg, '(3x,a,g0.3," MW")') 'current drive: I=', icd_beam*1.0e3_wp_
+        write(msg, '(3x,a,g0.3," kA")') 'current drive: I=', icd_beam*1.0e3_wp_
         call log_info(msg, mod='gray_core', proc='gray_main')
 
         if(ip < params%raytracing%ipass) then
@@ -1373,7 +1373,7 @@ contains
     dbvdc(3,1) = dbvcdc(3,1)
     dbvdc(1,2) = -bv(2)
     dbvdc(2,2) =  bv(1)
-    dbvdc(3,2) = dbvcdc(3,2)
+    dbvdc(3,2) = dbvcdc(3,2) ! = 0
     dbvdc(1,3) = dbvcdc(1,3)*csphi - dbvcdc(2,3)*snphi
     dbvdc(2,3) = dbvcdc(1,3)*snphi + dbvcdc(2,3)*csphi
     dbvdc(3,3) = dbvcdc(3,3)
@@ -1391,6 +1391,7 @@ contains
     ! B magnitude and derivatives
     btot  = norm2(bv)
 
+    ! dbtot(i) = d |B| / dxvcart(i)
     dbtot = matmul(bv, dbv)/btot
 
     yg = btot/Bres