c volts.inc
      common /sigma2/sigma2(njp1,nk2p3)
      common/volt/psi(njp1,nk2p3),curnt(njp1,nk2p3),
     $		sigmap(njp1,nk2p3),sigmah(njp1,nk2p3),psii(njp1,nk2p1),
     $		vxint(nj,nk2),vyint(nj,nk2),vzint(nj,nk2),
     $		dflux(njp1,nk2p3),denergy(njp1,nk2p3)
      common/wts/eta1(njp1,nk2p1),chi1(njp1,nk2p1),
     $		etahf(nj,nk2),chihf(nj,nk2),
     $		psieta(nj,nk2),psichi(nj,nk2),
     $		vxeta(nj,nk2),vyeta(nj,nk2),vzeta(nj,nk2),
     $		vxchi(nj,nk2),vychi(nj,nk2),vzchi(nj,nk2)