750 DQ=(YS(I> -YCOLD(I) ) / <2#B) 



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



770 G(I)=Q<I)*DQ/QDIFF 



780 Q<I+1)=Q<I+1)*DQ/QDIFF 



790 ENDIF 



800 Y<I)=B*(GHI) -GKI + 1) ) +YCOLDCI) 



810 C* 



820 C# Calc o-f Y starting from element to the left o-f minus 



830 C# point or boundary. Q is negative. 



840 C* 



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



860 Q(J)=E<J)*GKJ + 1) +F< J) 



870 YC = 2*B*(GKJ) -Q ( J+l) ) +YC0LD(J) 



880 IF(YC.LT. YS(J) .AND. J.6E.YSBEG) THEN 



890 DQ=(YS(J) -YCOLD(J) ) / (2*B) 



900 Q(J)=Q(J+1)+DQ 



910 ENDIF 



920 Y ( J ) =B* <GKJ)-GKJ + 1>) + YC0LD ( J ) 



930 40 CONTINUE 



940 1=1+1 



950 IFd.GE.YSEND+1) GOTO 100 



960 C* 



970 C* Cp.Ic o+ Y starting from element to the right o-f minus 



980 C* point or boundary. Q is positive. 



990 C* 



1000 GOTO 10 



1010 100 CONTINUE 



1020 DO 110 I=YSBEG-1, 1, -1 



1030 Q(I)=E(I)*Q(I+1) +F(I) 



1040 Y(I)=B*(Q(I)-Q(I+1) ) +YC0LD(I) 



1050 110 CONTINUE 



1060 IF(YSEND.NE.N) THEN 



1070 DO 120 I=YSEND+1,N 



1080 GKI+1 )=P(I+1)*Q(I) +R(I+1) 



1090 Y(I)=B*< Q(I) -GKI+1) ) +YC0LD(I) 



1100 120 CONTINUE 



1110 ENDIF 



1120 RETURN 



1130 END 



A10 



