anchor.diff:          INTEGER::SHAPE=3
           INTEGER::V(4)=0
           INTEGER::HKEY=-1
  -        CHARACTER(64)::CKEY=""
  -        INTEGER::ISTRISURFACE=-1
  -        REAL(8)::NOR(3)=0.D0
  +        CHARACTER(64)::CKEY="" 
  +        INTEGER::ISTRISURFACE=0
  +        REAL(8)::UNORMAL(3)=0.D0
           INTEGER::NEL=0
           INTEGER,ALLOCATABLE::ELEMENT(:)        
       ENDTYPE
  
positive.code:     USE DS_Gmsh2Solver
      USE hashtbl
      IMPLICIT NONE
      INTEGER::I,J,K,IGP1,NFACE1,IEL1,NTRI1=0,N1=0,NQUA1=0,K1=0
      INTEGER,ALLOCATABLE::NODE1(:),TRI1(:,:)
  	REAL(8),ALLOCATABLE::XYZ1(:,:),FACENOR1(:,:)
      TYPE(DICT_DATA), ALLOCATABLE:: val1(:)
  
negative.code:             PHYSICALGROUP(IGP1).TRINODEL2G=NODE1(1:PHYSICALGROUP(IGP1).NTRINODEL2G)
              ALLOCATE(PHYSICALGROUP(IGP1).TRISURFACE(NTRI1),PHYSICALGROUP(IGP1).QUASURFACE(NQUA1))
              DO CONCURRENT (J=1:NFACE)
                  IF(ABS(FACE(J).ISTRISURFACE)==IGP1) THEN                    
                      IF(FACE(J).SHAPE==3) THEN
  						PHYSICALGROUP(IGP1).NTRISURFACE=PHYSICALGROUP(IGP1).NTRISURFACE+1
  						PHYSICALGROUP(IGP1).TRISURFACE(PHYSICALGROUP(IGP1).NTRISURFACE)=J
  