anchor.diff:      expansile(num_elems+1:num_elems_new)=.false.
   
       call enter_exit(sub_name,2)
  -    
  +
     end subroutine reallocate_node_elem_arrays
  -
   
   !!!###############################################################
   
  
positive.code:     elems_at_node(1:num_nodes,0:3)=enodes_temp(1:num_nodes,0:3)
      deallocate(enodes_temp)
      elems_at_node(num_nodes+1:num_nodes_new,0:3)=0
      
      allocate(exp_temp(num_elems))
      exp_temp = expansile
      deallocate(expansile)
  
negative.code:     elem_cnct(-1:1,0:2,0:num_elems)=enodes_temp2(-1:1,0:2,0:num_elems)
      deallocate(enodes_temp2)
      elem_cnct(-1:1,0:2,num_elems+1:num_elems_new) = 0
  
      allocate(enodes_temp(num_ord,num_elems))
      enodes_temp=elem_ordrs
      deallocate(elem_ordrs)
  