计算物理2.docVIP

数 组 定义格式 DIMENSION A1(N1), A2(N1, N2), ASD(N1, N2, N3)，K(N,M) 定义了三个数组。 说明： 1 定义语句，放在执行语句之前； 2 数组类型遵循I-N规则。可强行说明。 例 DIMENSION A1(4,5), I(10), M(5,5), SS(100) REAL I, M 不能写成REAL I(10), M(5,5) INTEGER A1 LOGICAL LA 或直接写成： REAL I(10), M(100) INTEGER A1(4,5), SS(100) 3 同一程序块中，数组名和变量名不能相重，也不能用同一数组名代表不同的数组。例如： DIMENSION X(5), X(5,5) X=3. 4 标准fortran 中N1 ,N2, N3 为正整数，定义数组的大小。 很多版本允许如下格式定义： DIMENSION X(2 : 5), XX(4 : 5, 5 : 6) 定义了X数组有X(2 )，X(3 )，X(4 )，X(5 )四个元素。 XX有XX(4 ,5 )， XX(4，6) ，XX( 5, 5)，XX(5, 6) 四个元素。 还允许下标为负整数，例如： DIMENSION X(-2 :2) 定义了X数组有X(-2 )，X(-1 )，X(0 )，X(1 ) ，X(2 )五个元素。 特别注意：限制的放宽往往会有一定的隐患： 例如：DIMESION A( 2 , 2 ) DO I= - 2,2 DO J= - 2,2 A( J, I )=I ENDDO ENDDO WRITE(*,*) A STOP END 发现结果的值不对，所以特别注意数组越界的问题。 数组的输入和输出 多维数组元素在计算机内存中按列排列， DIMENSION A( 3 , 3 ) 顺序为： A( 1 , 1 ) A( 2 , 1 ) A( 3 , 1 ) A( 1 , 2 ) A( 2 , 2 ) A( 3 , 2 ) A( 1 , 3 ) A( 2 , 3 ) A( 3 , 3 ) 输入输出方式 1 A( I, J ) = e 表达式 2 输入 READ( N , *) A 从文件或键盘读，特别注意顺序，按列。 输出 WRITE( N, *) A 写到文件或显示器。 3 输入 DO 100 I=1,3 DO 100 J=1,3 100 READ( N , * ) A( I, J ) 该循环式按行顺序读入。 输出 DO 100 I=1,3 DO 100 J=1,3 100 WRITE( N , * ) A( I, J ) 该循环式按行顺序输出。 4 隐含循环 =========================== WRITE( N , * ) (A( I ), I=1,3) 顺序等效于 DO 100 I=1 , 3 100 WRITE(N,*) A(I) =========================== WRITE( N , * ) (A( I, J ), J=1,3) DO 100 I=1 , 3 100 WRITE(N,*) A( I, J ) 3=========================== WRITE( N , * ) ( (A ( I, J ) , J = 1,3), I = 1, 3) 顺序等效于 DO 100 I= 1, 3 DO 100 J= 1, 3 100 WRITE( N , * ) A ( I, J ) =========================== WRITE( N , * ) ( (A（i） ，b（i）, I = 1, 3) 顺序等效于 DO 100 I= 1, 3 DO 100 J= 1, 3 100 WRIT