anchor.diff: -
  -
    subroutine readin(itype)
   	use SOLVERLIB
       use ds_hyjump
  
positive.code: 	!INITIALIZATION 
      
      IF(ISEXCA2D/=0) THEN
          CALL GenElement_EXCA2()
          solver_control.bfgm=continuum
          !CHECKDATA
          DO I=1,NSOILPROFILE
  		    do j=1,soilprofile(i).nasoil
  			    IF(KPOINT(NDIMENSION,soilprofile(i).asoil(j).z(1))<KPOINT(NDIMENSION,soilprofile(i).asoil(j).z(2))) THEN
                      PRINT *, "Z1 IS SMALLER THAN Z2.PLEASE CKECK. SOILPROFILE=, ACTIVE SIDE SOILLAYER=",I,J
                      STOP
                  ENDIF
  		    enddo
  		    do j=1,soilprofile(i).npsoil
                 IF(KPOINT(NDIMENSION,soilprofile(i).Psoil(j).z(1))<KPOINT(NDIMENSION,soilprofile(i).Psoil(j).z(2))) THEN
                      PRINT *, "Z1 IS SMALLER THAN Z2.PLEASE CHECK. SOILPROFILE=, PASSIVE SIDE SOILLAYER=",I,J
                      STOP
                  ENDIF
              enddo 
              
          ENDDO
          msg=INSERTMENUQQ (5, 0, $MENUENABLED, 'GraphSetting'c,NUL)
  		msg=INSERTMENUQQ (5, 1, $MENUENABLED, 'Gray'c, SetLineColor)
  		msg=INSERTMENUQQ (5, 2, $MENUENABLED, 'NoMarker'c, Marker)
  
negative.code: 	CHARACTER(1024)      ext
  	type(qwinfo) winfo
  	integer(4)::length,msg
      real(DPN)::AR1(100)
  	type(bc_tydef),allocatable::bf1(:),bf2(:)
      EXTERNAL::SetLineColor,Marker,LineStyle,SETBGCOLOR
      
  