anchor.diff:            call fstr_rcap_send( fstrCPL )
           endif
   
  -        if( fstrPARAM%fg_couple_first==5 .and. revocap_flag==0 ) then
  -          do j = 1, hecMAT%NP * ndof
  -            hecMAT%B(j) = prevB(j)
  -          enddo
  -          go to 1000
  -        endif
  -
  -        if( fstrPARAM%fg_couple_first==6 ) then
  +        if( fstrPARAM%fg_couple_first==5 .or. &
  +            fstrPARAM%fg_couple_first==6 ) then
             call fstr_get_convergence( revocap_flag )
  -          if( revocap_flag == 0 ) then
  +          if( revocap_flag==0 ) then
               do j = 1, hecMAT%NP * ndof
                 hecMAT%B(j) = prevB(j)
               enddo
  +            if( fstrPARAM%fg_couple_first==5 ) then
  +              go to 1000
  +            else
  +              call fstr_rcap_get( fstrCPL )
  +              call dynamic_mat_ass_couple( hecMESH, hecMAT, fstrSOLID, fstrCPL )
  +              go to 2000
  +            endif
  +          else
  +            if( fstrPARAM%fg_couple_first==6 .and. &
  +                i /= fstrDYNAMIC%n_step ) then
  +              call fstr_rcap_get( fstrCPL )
  +              call dynamic_mat_ass_couple( hecMESH, hecMAT, fstrSOLID, fstrCPL )
  +            endif
             endif
           endif
   
  -        if( fstrPARAM%fg_couple_first==4 .or. &
  -            fstrPARAM%fg_couple_first==6 ) then
  +        if( fstrPARAM%fg_couple_first==4 ) then
             call fstr_rcap_get( fstrCPL )
             call dynamic_mat_ass_couple( hecMESH, hecMAT, fstrSOLID, fstrCPL )
  -        endif
  -
  -        if( fstrPARAM%fg_couple_first==6 .and. revocap_flag==0 ) then
  -          go to 2000
           endif
         endif
   !C *****************************************************
  
positive.code:             end if
      call hecMAT_finalize( MAT0 )
  
  
  !C-- end of implicit dynamic analysis
  !C
  
negative.code: 
  1000  continue
  
          if( fstrPARAM%fg_couple_first==5 .or. &
              fstrPARAM%fg_couple_first==6 ) then
            do j = 1, hecMAT%NP * ndof
              prevB(j) = hecMAT%B(j)
            enddo
          endif
  
          if( fstrPARAM%fg_couple_first==1 .or. &
  