给出完整代码，完整输入文件。

我贴出的是输入部分的代码，但是输入只能输入X-DIRECTION GRID SPACING (M)
    0.0  100.0   10.0    4.0    1.0    1.0    1.0    1.0    1.0    1.0
    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0
    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    4.0   10.0
  100.0    0.0

这是需要输入的

X-DIRECTION GRID SPACING (M)
    0.0  100.0   10.0    4.0    1.0    1.0    1.0    1.0    1.0    1.0
    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0
    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    4.0   10.0
  100.0    0.0
Y-DIRECTION GRID SPACING (M)
    0.0  100.0   10.0    4.0    1.0    1.0    1.0    1.0    1.0    1.0
    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0
    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    4.0   10.0
  100.0    0.0
Z-DIRECTION GRID SPACING (M)
    0.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0    1.0
    2.0    2.0    2.0    4.0   10.0  100.0    0.0
REFERENCE POINT
     16     16    0.0    0.0
RESISTIVITY BLOCK (OHM*M)
   4
   1  15   1  15   1  15     100.0
  16  30   1  15   1  15     200.0
   1  15  16  30   1  15     300.0
  16  30  16  30   1  15     400.0
NUM. OF CURRENT NODE AND CURRENT VALUES
   2   1.0
CURRNT NODE POINT (X,Y,Z)
       0.0       0.0       0.0
       0.0       0.0       5.0
WRITE CONTROL PARAM. ( X-START,X-END,Y-START,Y-END,RHOa:1 or V:2 )
   4  28   4  28   1

C---------------------------------------
C  INPUT GRID SEPARATION in free format
C---------------------------------------
      READ(5,9)  DUMMY
    9 FORMAT(A50)
      READ(5,*) (X(I),I=0,IE)
      READ(5,9)  DUMMY
      READ(5,*) (Y(J),J=0,JE)
      READ(5,9)  DUMMY
      READ(5,*) (Z(K),K=0,KE)
      READ(5,9)  DUMMY
      READ(5,*) IXREF,JYREF,VXREF,VYREF
C--------------------------------------------------
C  CALC DISTANCE from the origin: DX(I),DY(J),DZ(K)

具体代码，具体分析

贴出代码。