|
编译环境:IVF 2013,compliers XE 2013
用QUICK WIN作图
代码如下:
[Fortran] 纯文本查看 复制代码 02 | integer , PARAMETER :: IR = 400 |
03 | integer , parameter :: JR = 400 |
04 | INTEGER IS ( 0 : IR +1 , 0 : JR +1 ) , TMAX , ISN ( 1 : 8 ) , NSTATE , T , NR , IX , IY |
05 | WRITE ( * , * ) "PLEASE INPUT THE TIME STEP " |
14 | DISTANCE = SQRT ( 1.0 * ( I - IRC ) * * 2 +1.0 * ( J - JRC ) * * 2 ) |
15 | IF ( DISTANCE .LT. R ) IS ( I , J ) = 2 |
20 | OPEN ( 1 , FILE = "E:\LUKE.DAT" ) |
27 | ISN = ( / IS ( IX -1 , JY -1 ) , IS ( IX -1 , JY ) , IS ( IX -1 , JY +1 ) , IS ( IX , JY -1 ) , IS ( IX , JY +1 ) , IS ( IX +1 , JY -1 ) , IS ( IX +1 , JY ) , IS ( IX +1 , JY +1 ) / ) |
28 | E 0 = COUNT ( ISN .NE. IS ( IX , JY ) ) |
33 | E = COUNT ( ISN .NE. NSTATE ) |
35 | DE = E - E 0 + NSTATE - IS ( IX , JY ) +2.5 * RD -1.25 |
36 | IF ( DE .LT. 0.0 ) IS ( IX , JY ) = NSTATE |
37 | ISRE = SETCOLOR ( IS ( IX , JY ) ) |
41 | WRITE ( 1 , * ) T , COUNT ( IS .EQ. 2 ) |
最后出来这个东西:
急啊,这到底什么回事呢?我代码哪里写错了?
|
|