anchor.diff:  
   			
   			inode1=element(ieL1).node
  +			if(element(ieL1).ifreedof>0) then
  +				if(inode1(1)==freedof(element(ieL1).ifreedof).newnode) inode1(1)=freedof(element(ieL1).ifreedof).node
  +				if(inode1(2)==freedof(element(ieL1).ifreedof).newnode) inode1(2)=freedof(element(ieL1).ifreedof).node
  +			endif
   			iel2=pile(ipile).element_sp(aop,inode1)
   			element(iel2).isactive=1
   			
  
positive.code: 			
  			iel1=pile(IPILE).element(j)
  			nnode1=element(iel1).node
  			do k=1,2
  				!displacement
                  IF(ISEXCA2D==1) THEN
  
negative.code: 		integer::nseg=1,nnode=0,nel=0,system=0
          integer,allocatable::mat(:),kpoint(:)        
  		integer,allocatable::node(:),element(:)
  		integer,allocatable::element_sp(:,:) !ڵɵԪ
  		integer,allocatable::node2BCload(:)
          real(double),allocatable::Beamresult(:,:,:) !(node,13,istep),/dis,M,Q,Load/
  		INTEGER::NVA=14
  