anchor.diff:      ! nconb => node list body mesh
       ! ncondb => normal list body
       real(8),allocatable :: XYZE(:,:,:),DXYZE(:,:,:),TXYZE(:,:,:)
  -    ! xyze =>  used in convsb, new full mesh for node (3,8,elem_id)
  -    ! dxyze => used in convsb, new full mesh for normal (3,8,elem,id)
  +    ! xyze =>  before combine, new full mesh for node (3,8,elem_id)
  +    ! dxyze => before combine, new full mesh for normal (3,8,elem,id)
       real(8),allocatable :: XYZTP(:,:),DXYZTP(:,:)
       ! xyztp => convsb, combined node mesh, xyztp(3,node_id)
       ! dxyztp = > combined normal mesh, dxyztp(3,nrml_id)
  
positive.code: 200     CONTINUE
  
         NNF=M     
       write(11,  *) ' NNF=',NNF
         END IF    
  !
  ! -----------------------------------------
  
negative.code: 
         NNODE=M
       NNB=NNODE-NNF
      ! write(11,*)  ' NNB=',NNB,' NNF=',NNF,' NNODE=',NNODE
       !write(11,*) 
  !
  ! ===========================================================
  !
  