750 DIFF=YS(I) -Y(I) 



760 QDIFF=Q(I+1) -Q < I ) 



770 <3<I)=Q<I> -DIFF/B2*<Q(I> /QDIFF) 



780 Q(I+1)=Q(I+1) -DIFF/B2*<6MI + 1) /QDIFF) 



790 Y(I)=YS(I) 



800 ENDIF 



810 C* 



820 C* Calc o-f Y starting -from element to the le-ft o-f minus 



830 C# point or boundary. Q is negative. 



840 C# 



850 30 DO 40 J=I-1,K,-1 



860 Y( J)=Y(J) -B2*(Q(J+1) -Q(J) ) 



870 IF<Y(J) .LT. YS(J) .AND. J.GE. YSBEG) THEN 



880 DIFF=YS(J) -Y(J) 



890 Q(J)=Q(J) +DIFF/B2 



900 Y(J)=YS(J) 



910 ENDIF 



920 40 CONTINUE 



930 1=1+1 



940 IFd.GE. YSEND+1) GOTO 100 



950 C* 



960 C# Calc o-f Y starting -from element tD the right of minus 



970 C* point or boundary. Q is positive. 



980 C* 



990 GOTO 10 



1000 100 CONTINUE 



1010 RETURN 



1020 END 



A5 



