anchor.diff:      if( fstrPARAM%fg_couple == 1) then
         if( fstrPARAM%fg_couple_first==5 .or. &
             fstrPARAM%fg_couple_first==6 ) then
  -        allocate( prevB(0:hecMAT%NP)      ,STAT=ierror )
  +        allocate( prevB(hecMAT%NP*ndof)      ,STAT=ierror )
               if( ierror /= 0 ) then
                 write(idbg,*) 'stop due to allocation error <fstr_solve_LINEAR_DYNAMIC, prevB>'
                 write(idbg,*) '  rank = ', hecMESH%my_rank,'  ierror = ',ierror
  
positive.code:         endif
  
          if( fstrPARAM%fg_couple_first==6 ) then
            do j = 1, hecMAT%NP * 3
              prevB(j) = hecMAT%B(j)
            enddo
          endif
  
negative.code:         if( fstrPARAM%fg_couple_first==6 ) then
            call fstr_get_convergence( revocap_flag )
            if( revocap_flag == 0 ) then
              do j = 1, hecMAT%NP * ndof
                hecMAT%B(j) = prevB(j)
              enddo
            endif
  