anchor.diff:            relres = dabs(res1-res)/res0
           endif
           if( hecMESH%my_rank==0 ) then
  -          write(*,'(a,i3,a,2e15.7)') ' - Resiual(',iter,') =',res,relres
  +          write(*,'(a,i3,a,2e15.7)') ' - Residual(',iter,') =',res,relres
             write(ISTA,'(''iter='',I5,''res/res0='',2E15.7)')iter,res,relres
           endif
   
  
positive.code:           relres = dabs(res1-res)/res0
          endif
          if( hecMESH%my_rank==0 ) then
            write(*,'(a,i3,a,2e15.7)') ' - Resiual(',iter,') =',res,relres
            write(ISTA,'(''iter='',I5,''res/res0='',2E15.7)')iter,res,relres
          endif
  
  
negative.code:         call hecmw_innerProduct_R(hecMESH,ndof,hecMAT%X,hecMAT%X,resX)
          resX = sqrt(resX)/n_node_global
          if( hecMESH%my_rank==0 ) then
            write(*,'(a,i3,a,e15.7)') ' - ResidualX    (',iter,') =',resX
            write(*,'(a,i3,a,e15.7)') ' - ResidualX+LAG(',iter,') =',sqrt(x_residual)/n_node_global
            write(*,'(a,i3,a,e15.7)') ' - ResidualQ    (',iter,') =',sqrt(q_residual)/n_node_global
          endif
  !   ----- update the small displacement and the displacement for 1step 
          do i=1,hecMESH%n_node*ndof
  