diff --git a/bbb/bbb.v b/bbb/bbb.v
index a9290cef..217f131d 100755
--- a/bbb/bbb.v
+++ b/bbb/bbb.v
@@ -1726,7 +1726,6 @@ volpsor(0:nx+1,0:ny+1,1:nisp) _real [1/s]+maybeinput #current src into ions in c
 volmsor(0:nx+1,0:ny+1,1:nisp) _real [kg m/s**2] +maybeinput #up mom src in cell ix,iy
 voljcsor(0:nx+1,0:ny+1) _real [A]      +maybeinput #uniform core-region curr sor. in ix,iy
 volpsorg(0:nx+1,0:ny+1,1:ngsp) _real [1/s]+maybeinput #curr source for gas in cell ix,iy
-pwrsorg(0:nx+1,0:ny+1,1:ngsp) _real [1/s]+maybeinput #power source for gas in cell ix,iy; only for manufactured solution verification
 pondpot(0:nx+1,0:ny+1)  _real [V] /0./ +maybeinput #elec ponderomotive potential
 psgov_use(0:nx+1,0:ny+1,1:ngsp) _real [1/m**3 s]+maybeinput #user-specified gas source
 jcvsor                  real [A]  /0./ +maybeinput #total core-region current for voljcsor
diff --git a/bbb/oderhs.m b/bbb/oderhs.m
index db6de4b2..8c562c19 100755
--- a/bbb/oderhs.m
+++ b/bbb/oderhs.m
@@ -4583,7 +4583,7 @@ ccc             call volave(nx, ny, j2, j5, i2, i5, ixp1(0,0), ixm1(0,0),
      .                       ( -mi(1)*t1*t2*(psor(ix,iy,1)+tv)
      .                         +0.5*mi(1)*t1*t1*
      .                          (psor(ix,iy,1)+psorrg(ix,iy,1)+2*tv) )
-               reseg(ix,iy,1) = reseg(ix,iy,1) + pwrsorg(ix,iy,1)
+               reseg(ix,iy,1) = reseg(ix,iy,1)
      .                            - t0+0.5*mg(1) * ( (t1-t2)*(t1-t2)
      .                                              +temp3+temp4 )
      .                            * (psorrg(ix,iy,1)+tv)