anchor.diff:  	delta1=matproperty(soil.mat,8,istep)/180*PI()
   	if(abs(fi1)<1e-7) then
   		CDelta1=0.d0
  -	else
  -		CDelta1=asin(sin(delta1)/sin(fi1))
  +    else
  +        if(delta1>fi1) then
  +            stop "InputDataError. The wall friction angle > soil friction angle. sub=ep_Cal"
  +        else
  +            CDelta1=asin(sin(delta1)/sin(fi1))            
  +        endif
  +        
  +		
   	endif
   	
   	select case(soilprofile(isp).spm)
  
positive.code: 	
  	R1=0.0d0
  	do i=1,soilprofile(isoilprofile).nasoil
  		if(sf(soilprofile(isoilprofile).asoil(i).sf).factor(istep)==0) cycle
          !򻯿͸ǽ͸ˮǽˮѹȡ
  	    if(kpoint(ndimension,soilprofile(isoilprofile).asoil(i).z(2))<awl1) then!ˮλֲ
  
negative.code: 	enddo
  	
  	do i=soilprofile(isoilprofile).npsoil,1,-1
          soilprofile(isoilprofile).psoil(i).pw=0.d0
  		if(sf(soilprofile(isoilprofile).psoil(i).sf).factor(istep)==0) cycle
  		
          !򻯿͸ǽ͸ˮǽˮѹȡ
  