|
|
发表于 2003-4-11 13:07:40
|
显示全部楼层
SIMPLE全部源程序!!!!!FORTRAN语言!!!
[这个贴子最后由yuxin在 2003/04/11 01:10pm 编辑]
做如下修改:
***** ..\P1.FOR
*-----------------------------------------------------------------------
$INCLUDE:'SIMPLE.INC'
DIMENSION TH(22),THU(22),THDIF(22),THCV(22),THCVS(22)
***** P1.FOR
*-----------------------------------------------------------------------
INCLUDE 'SIMPLE.INC'
DIMENSION TH(22),THU(22),THDIF(22),THCV(22),THCVS(22)
*****
***** ..\SIMPLE.FOR
*-----------------------------------------------------------------------
$INCLUDE:'SIMPLE.INC'
*-----------------------------------------------------------------------
***** SIMPLE.FOR
*-----------------------------------------------------------------------
INCLUDE 'SIMPLE.INC'
*-----------------------------------------------------------------------
*****
***** ..\SIMPLE.FOR
*-----------------------------------------------------------------------
$INCLUDE:'SIMPLE.INC'
COMMON/CNTL/LSTOP
***** SIMPLE.FOR
*-----------------------------------------------------------------------
INCLUDE 'SIMPLE.INC'
COMMON/CNTL/LSTOP
*****
***** ..\SIMPLE.FOR
*-----------------------------------------------------------------------
$INCLUDE:'SIMPLE.INC'
*-----------------------------------------------------------------------
***** SIMPLE.FOR
*-----------------------------------------------------------------------
INCLUDE 'SIMPLE.INC'
*-----------------------------------------------------------------------
*****
后,用
Compaq Visual Fortran Optimizing Compiler Version 6.5
编译连接,运行并输入out10.dat
得到屏幕输出:
ITER T(4,4) T(5,3)
0 0.000000 0.000000
1 1.999978 1.720364
2 2.000000 1.720001
3 2.000000 1.720000
4 2.000000 1.720000
5 2.000000 1.720000
6 2.000000 1.720000
7 2.000000 1.720000
8 2.000000 1.720000
9 2.000000 1.720000
10 2.000000 1.720000
Stop - Program terminated.
文件out10.dat中
COMPUTATION IN CARTESIAN COORDINATES
****************************************
I = 2 3 4 5 6 7
XU = 0.00E+00 2.00E-01 4.00E-01 6.00E-01 8.00E-01 1.00E+00
J = 2 3 4 5 6 7
YV = 0.00E+00 4.00E-01 8.00E-01 1.20E+00 1.60E+00 2.00E+00
I = 1 2 3 4 5 6 7
X = 0.00E+00 1.00E-01 3.00E-01 5.00E-01 7.00E-01 9.00E-01 1.00E+00
J = 1 2 3 4 5 6 7
Y = 0.00E+00 2.00E-01 6.00E-01 1.00E+00 1.40E+00 1.80E+00 2.00E+00
************************** TEMP **************************
I = 1 2 3 4 5 6 7
J
7 2.00E+00 2.30E+00 2.90E+00 3.50E+00 4.10E+00 4.70E+00 5.00E+00
6 1.80E+00 2.08E+00 2.64E+00 3.20E+00 3.76E+00 4.32E+00 4.60E+00
5 1.40E+00 1.64E+00 2.12E+00 2.60E+00 3.08E+00 3.56E+00 3.80E+00
4 1.00E+00 1.20E+00 1.60E+00 2.00E+00 2.40E+00 2.80E+00 3.00E+00
3 6.00E-01 7.60E-01 1.08E+00 1.40E+00 1.72E+00 2.04E+00 2.20E+00
2 2.00E-01 3.20E-01 5.60E-01 8.00E-01 1.04E+00 1.28E+00 1.40E+00
1 0.00E+00 1.00E-01 3.00E-01 5.00E-01 7.00E-01 9.00E-01 1.00E+00
————————————————————————————————————————————
另外P5.FOR和P6.FOR有很多同名子目录(不同问题、参数),
使用时肯定只能摘出一个来用,P6.FOR中语句标号14和301没有,应去掉?:
***** ..\P6.for
2 CONTINUE
DO 14 I=1,4
LSOLVE(I)=.TRUE.
***** P6.for
2 CONTINUE
DO I=1,4
LSOLVE(I)=.TRUE.
*****
***** ..\P6.for
FL=0.
DO 301 I=2,L2
FL=FL+RHO(I,M1)*V(I,M2)*XCV(I)
***** P6.for
FL=0.
DO I=2,L2
FL=FL+RHO(I,M1)*V(I,M2)*XCV(I)
***** |
|