|
本帖最后由 LYCD 于 2023-11-14 10:02 编辑
[Fortran] 纯文本查看 复制代码 01 | DIMENSION GB 1 ( 2 , 2 ) , GB 2 ( 2 , 2 ) , JB ( 2 , 10 ) , GB ( 2 , 2 ) , |
02 | $ GE 1 ( 10 , 2 ) , GE 2 ( 10 , 2 ) , JE ( 10 , 10 ) , GE ( 2 , 2 ) , XB ( 2 ) |
03 | common / INV / GBx ( 70 , 70 ) , GBy ( 70 , 70 ) , GBai ( 70 , 70 ) |
04 | common / MATRIX / A ( 100 , 100 ) , B ( 100 , 100 ) , AA ( 100 , 100 ) |
05 | DOUBLE PRECISION GBX , GBY , GBAI , GB 1 , GB 2 , GB , GE 1 , GE 2 , GE , JE , A , B , AA |
07 | CHARACTER RI * 12 , R * 1 , ZZI * 10 , ZZO * 10 |
10 | WRITE ( * , '(1X,' 'Input the name of input file :' ')' ) |
12 | OPEN ( 1 , FILE = ZZI , STATUS = 'OLD' ) |
13 | WRITE ( * , '(1X,' 'Input the name of output file :' ')' ) |
15 | OPEN ( 2 , FILE = ZZO , STATUS = 'NEW' ) |
24 | c * * * * * * * * U 2 = 1 , U 8 = 2 , F 2 = 3 , F 8 = 4 * * * * * * * * * |
25 | c * * * * * * * K 1 & K 2 are the reserved variables * * * * * * * * |
35 | GB 2 ( 1 , 2 ) = -0.5 * L / ( ARES * E ) |
36 | GB 2 ( 2 , 1 ) = -0.5 * L / ( ARES * E ) |
43 | JB ( 1 , I ) = - RP * DX * 0.5 * X / E |
44 | JB ( 2 , I ) = - RP * DX * 0.5 * ( L - X ) / E |
47 | GE 2 ( I , 1 ) = 0.5 * X / ( ARES * E ) |
48 | GE 2 ( I , 2 ) = -0.5 * ABS ( L - X ) / ( ARES * E ) |
53 | JE ( I , J ) = -0.5 * ABS ( X 1 - X 2 ) * RP * DX / E |
|
|