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