第11章程序设计347 READ(S, N !土条总数 READ(5,)RW,UWL!水容重,坡外水位 READ5,*)IN!线段总数 DO 3021=L IN READ(5,*)(C(D)J=1,3)!线段信息 302 CONTINUE !点总数 READ,*)(L,XN(D),YN①)=1,NN)!点坐标 DOI=LN READ(,)I,(D,Y ENDDO IF(PH EQ.O)GOTO DO I=LIPH READ5,*)(ICPH(,D),J=1,3)!线段信息 18 READ(S, )IN, ( L, PF(), PC), PF1(), PCI(), PDW(I), PDS) %LRU IFl, IN) CONTINUE 对土条信息重新处理,保证IC,1)的x坐标比IC(1,2)的小 DO 37=l.INI M2=C(I,2) IF(XN(MI). LTXN(M2)GO TO 37 IC(, 2FMI 37 CONTINUE F(PH EQ0)GO TO 672 C REARRANGE ICPH(I )IN THE SAME WAY AS IC(, J) DO 38 I=LIPH MI=ICPH(I, 1) M2=ICPH(, 2) IF(XN(MD).LTXN(M2))GO TO 38 PH(, IFICPH(, 2) H(,2)=MI 38 CONTINUE 672 CONTINUE GTENSION=O FIND THE LINES WHICH INTERCEPT THE CONSIDERED SLICE MX:IC(. YYY=YN(MX) DO 813=lINI DO813Il=1,2 Y=IC(, D IF(YNGY)GT.YYYYYY=YNY 813 CONTINUE YYY=YYY+100 X3(1=X2(1) Y3(1)=Y2(1) DO211K=1 QN=0. W(k=0 DHREO WL(KFYYY第 11 章 程序设计 347 READ(5,*) N ! 土条总数 READ(5,*) RW,UWL ! 水容重 坡外水位 READ(5,*) IN1 ! 线段总数 DO 302 I=1,IN1 READ(5,*) (IC(I,J),J=1,3) !线段信息 302 CONTINUE READ(5,*)NN ! 点总数 READ(5,*) (L,XN(I),YN(I),I=1,NN) ! 点坐标 DO I=1,N READ(5,*) I,X(I),Y(I) ENDDO READ(5,*)IPH IF(IPH.EQ.0)GOTO 18 DO I=1,IPH READ(5,*) (ICPH(I,J),J=1,3) !线段信息 ENDDO 18 READ(5,*)IN,(L,PF(I),PC(I),PF1(I),PC1(I),PDW(I),PDS(I), %LRU(I),I=1,IN) CONTINUE C 对土条信息重新处理 保证 IC(I,1)的 x 坐标比 IC(I,2)的小 DO 37 I=1,IN1 M1=IC(I,1) M2=IC(I,2) IF(XN(M1).LT.XN(M2))GO TO 37 IC(I,1)=IC(I,2) IC(I,2)=M1 37 CONTINUE IF(IPH.EQ.0)GO TO 672 C REARRANGE ICPH(I,J) IN THE SAME WAY AS IC(I,J) DO 38 I=1,IPH M1=ICPH(I,1) M2=ICPH(I,2) IF(XN(M1).LT.XN(M2))GO TO 38 ICPH(I,1)=ICPH(I,2) ICPH(I,2)=M1 38 CONTINUE 672 CONTINUE GTENSION=0 C FIND THE LINES WHICH INTERCEPT THE CONSIDERED SLICE MX=IC(1,1) YYY=YN(MX) DO 813 I=1,IN1 DO 813 II=1,2 IY=IC(I,II) IF(YN(IY).GT.YYY)YYY=YN(IY) 813 CONTINUE YYY=YYY+100. X3(1)=X2(1) Y3(1)=Y2(1) DO 211 K=1,N MF=0 QN=0. W(K)=0. DHR=0 WL(K)=YYY