anchor.diff:          integer ifuncs(6), iarray(10)
   
           real*8 elDw(numel) ! element average of DES d variable
  +
  +        real*8, allocatable, dimension(:,:) :: HBrg
  +        real*8, allocatable, dimension(:) :: eBrg
  +        real*8, allocatable, dimension(:) :: yBrg
  +        real*8, allocatable, dimension(:) :: Rcos, Rsin
   c
   c  Here are the data structures for sparse matrix GMRES
   c
  
positive.code:      &                       acold,         x,
       &                       iBC,           BC,
       &                       res,           
       &                       BDiag,         a(mHBrg),      a(meBrg),
       &                       a(myBrg),      a(mRcos),      a(mRsin),
       &                       iper,          ilwork,
       &                       shp,           shgl,
       &                       shpb,          shglb,         solinc, 
  
negative.code:      &                       iBC,           BC,
       &                       colm,          rowp,          lhsk,
       &                       res,
       &                       BDiag,         hBrg,          eBrg,
       &                       yBrg,          Rcos,          Rsin,
       &                       iper,          ilwork,
       &                       shp,           shgl,
       &                       shpb,          shglb,         solinc,
  