anchor.diff:  	ISET=NUESET+1
   	UESET(ISET).NAME=NAME
   	UESET(ISET).ENUM=KPELEMENT(1,N2)-KPELEMENT(2,N1)+1
  -    IF(KPELEMENT(1,N2)<KPELEMENT(2,N1)) THEN
  -        PAUSE
  -    ENDIF
  +	!if(ueset(ISET).enum==0) then
  +	!	pause
  +	!endif
  +    !IF(KPELEMENT(1,N2)<KPELEMENT(2,N1)) THEN
  +    !    PAUSE
  +    !ENDIF
   	ALLOCATE(UESET(ISET).ELEMENT(UESET(ISET).ENUM))
   	DO I=1,UESET(ISET).ENUM
   		UESET(ISET).ELEMENT(I)=KPELEMENT(2,N1)+I-1
  
positive.code:                 N2=SOILPROFILE(I).psoil(J).Z(2)
                  WRITE(CH1,'(I3)') I
                  WRITE(CH2,'(I3)') J
                  TITLE1="SP_"//TRIM(CH1)//"psoil_"//TRIM(CH2)
  				CALL Gen_USER_NSET_EXCA(N1,N2,TITLE1,SOILPROFILE(I).psoil(J).IUNSET)
  				CALL Gen_USER_ESET_EXCA(N1,N2,TITLE1,SOILPROFILE(I).psoil(J).IUESET)							 
  			else
  
negative.code:                 N2=SOILPROFILE(I).ASOIL(J).Z(2)
                  WRITE(CH1,'(I3)') I
                  WRITE(CH2,'(I3)') J
                  TITLE1="SP_"//TRIM(CH1)//"_ASOIL_"//TRIM(CH2)
  				CALL Gen_USER_NSET_EXCA(N1,N2,TITLE1,SOILPROFILE(I).ASOIL(J).IUNSET)
  				CALL Gen_USER_ESET_EXCA(N1,N2,TITLE1,SOILPROFILE(I).ASOIL(J).IUESET)							 
  			else
  