anchor.diff:            nenl = nenl0
             nshl = nshl0
   
  -c          if(lcsyst.eq.itp_wedge_tri) lcsyst=nenbl ! may not be necessary
  +          if(lcsyst.eq.itp_wedge_tri) lcsyst=nenbl ! may not be necessary
             ngaussb = nintb(lcsyst)
   c
   c.... compute and assemble non-unit normal
  
positive.code:           nenl = nenl1
            nshl = nshl1
  
  c          if(lcsyst.eq.itp_wedge_tri) lcsyst=nenbl ! may not be necessary
            ngaussb = nintb(lcsyst)
  c
  c.... compute and assemble non-unit normal
  
negative.code:           ipord   = lcblkif(5, iblk)    ! polynomial order
            nenl0   = lcblkif(6, iblk)    ! number of vertices per element0
            nenl1   = lcblkif(7, iblk)    ! number of vertices per element1
            nenbl   = lcblkif(8, iblk)    ! number of vertices on boundary
            mater0  = lcblkif(9, iblk)
            mater1  = lcblkif(10,iblk)
            nshl0   = lcblkif(iblkif_nshl0,iblk)
  