anchor.diff:          time   = 0
           yold   = y
           acold  = ac
  +
  +        allocate(HBrg(Kspace+1,Kspace))
  +        allocate(eBrg(Kspace+1))
  +        allocate(yBrg(Kspace))
  +        allocate(Rcos(Kspace))
  +        allocate(Rsin(Kspace))
  +
           if (mod(impl(1),100)/10 .eq. 1) then
   c
   c     generate the sparse data fill vectors
  
positive.code:      &                       iBC,           BC,
       &                       colm,          rowp,          lhsk,
       &                       res,
       &                       BDiag,         a(mHBrg),      a(meBrg),
       &                       a(myBrg),      a(mRcos),      a(mRsin),
       &                       iper,          ilwork,
       &                       shp,           shgl,
       &                       shpb,          shglb,         solinc,
  
negative.code:      &                       acold,         x,
       &                       iBC,           BC,
       &                       res,           
       &                       BDiag,         HBrg,          eBrg,
       &                       yBrg,          Rcos,          Rsin,
       &                       iper,          ilwork,
       &                       shp,           shgl,
       &                       shpb,          shglb,         solinc, 
  