| 查看: 1019 | 回复: 1 | ||
zhaoxj09295838铁杆木虫 (正式写手)
|
[求助]
求大神看看下面Fortran的例子，都是数组越界问题！在这里小弟先谢过了！！！ 已有1人参与
|
|
例题3-2 C WATER HAMMER FINITE DIMENSION V(100),H(100) READ(5,1) DX,DT,C,HO,VO,TOL,IMAX 1 FORMAT(6F7.0,I4) IMAX1=IMAX-1 V(I)=VO 2 H(I)=HO T=0 N=0 100 T=T+DT N=N+1 IF(T.GT.TOL) GO TO 3 SR=(TOL-T)/TOL GO TO 4 3 SR=0 4 CONTINUE V(IMAX)=SR*SQRT(2.*9.81*HO) DO 5 I=2,IMAX1 5 V(I)=V(I)-9.81*DT*(H(I)-H(I-1))/DX DO 6 I=2,IMAX1 6 H(I)=H(I)-DT*C**2/9.81*(V(I+1))-V(I)/DX WRITE (6,11) T 11 FORMAT(F10.3) WRITE(6,12) (V(I),I=1,IMAX) WRITE(6,12) (H(I),I=1,IMAX) 12 FORMAT (12F10.4) IF(N.LT.200) GO TO 100 STOP END 注：DX=600 DT=0.1 HO=5 VO=9.9 TOL=2 IMAX=12 例题4-3 C BACKWATER CURVES VARIABLE STEP METHOD DIMENSION H(50),V(50),DX(50) READ(5,1) CMAN,Q,SO,IMAX 1 FORMAT (3F7.0,I4) READ(5,2)H(I),H(IMAX) 2 FORMAT(2F7.0) IMAX1=IMAX-1 DO 3 I=1,IMAX H(I)=H(1)+(H(IMAX)-H(1))*(I-1)/IMAX1 3 V(I)=Q/H(I)/5 DO 4 I=1,IMAX1 HH=(H(I)+H(I+1))/2 A=5*HH P=2.*HH+5 R=A/P SF=Q**2/CMAN*2/R**1.333/A**2 4 DX(I)=(H(I+1)+V(I+1)**2/2./9.81-H(I)-V(I)**2/2./9.81)/(SO-SF) SUM=0 DO 5 I=1,IMAX1 5 SUM=SUM+DX(I) WRITE (6,6) SUM 6 FORMAT (F10.2) WRITE (6,7) (DX(I),I=1,IMAX1) 7 FORMAT (10F10.2) STOP END 注:CMAN=50——曼宁系数; Q=55.4——流量； SO=-0.001——底坡 DX(I)--流段长度 H（I）--水深 V(I)--流速 IMAX=30 例题4-4 C BCFS C METHOD DIMENSION H(100),V(100),B(100),HS(100),HST(100) READ(5,1) IMAX,DX,CMAN,Q,SO,H(1) 1 FORMAT(I4,5F7.0) READ(5,2)(B(I),I=1,IMAX) 2 FORMAT(10F7.0) V(1)=Q/B(1)/H(1) HS(1)=V(1)**2/2./9.81+H(1) DO 3 I=2,IMAX H(I)=10 DIFN=1 DH=1 20 DIF=DIFN NTEST=0 H(I)=H(I)-DH 15 V(I)=Q/B(I)/H(I) HH=(H(I-1)+H(I))/2 BB=(B(I-1)+B(I))/2 VV=Q/BB/HH R=HH**BB/(BB+2.*HH) SF=VV**2/CMAN**2/R**1.33 HS(I)=HS(I-1)+(SO-SF)*DX HST(I)=H(I)+V(I)**2/2./9.81 DIFN=HST(I)-HS(I) IF (NTEST.GT.0) GO TO 16 IF (DIF*DIFN.GT.0)GO TO 20 H(I)=H(I)+DH NTEST=1 GO TO 15 16 DH=DH/10 IF (DH.GT.0.001) GO TO 20 3 CONTINUE WRITE(6,4)(I,H(I),V(I),I=1,IMAX) 4 FORMAT (I6,2F10.3) STOP END 注：IMAX=30 DX=-1000 CMAN=30 Q=1000 SO=0.001 H(I)=5.6 |
回复此楼» 猜你喜欢
博士申请
已经有6人回复
-
评审有感
已经有11人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有10人回复
-
国自然上会要求
已经有9人回复
-
今年审到国自然15份，谈谈感受
已经有16人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fortran中的数组越界问题，
已经有3人回复
2楼2014-06-19 16:05:02
20