求助level程序计算一氧化碳报错

lqhenwunai

请教大家，我在用level16计算一氧化碳，输入文件如下
6 12 8 16 0 1
'co ground'
0.001 0.6 40 1.d-06
30 -1 0 89605.163975417
4 1 0 0 1
1  1 0
1.111  219.979422
1.112  199.041542
1.113  179.222983
1.114  160.479849
1.115  142.790194
1.116  126.197912
1.117  110.659108
1.118  96.173783
1.119  82.719988
1.12  70.319671
1.121  58.972833
1.122  48.613631
1.123  39.285959
1.124  30.989818
1.125  23.659365
1.126  17.360443
1.127  12.02721
1.128  7.681612
1.129  4.32365
1.13  1.931377
1.131  0.482844
1.132  0.0
1.133  0.460897
1.134  1.887482
1.135  4.213913
1.136  7.484085
1.137  11.67605
1.138  16.789809
1.139  22.803414
1.14  29.716865
4 1 2 3 1 1 2 -2
0 0 1 0 2 0 3 0
3 0 0.0
1.d0  -2.d-1  3.d-2  -4.d-3


报错：
At line 713 of file level16.f (unit = 5, file = 'stdin')
Fortran runtime error: End of file

Error termination. Backtrace:
#0  0x7fb84cf957ba
#1  0x7fb84cf96365
#2  0x7fb84cf96b2d
#3  0x7fb84d10d133
#4  0x7fb84d1060d4
#5  0x7fb84d107809
#6  0x560556c27337
#7  0x560556c2b040
#8  0x7fb84cc42082
#9  0x560556c00d39
#10  0xffffffffffffffff


得到的fort9 fort10文件只有表头，没有其它内容。

输出文件只输出以下内容：
Potential-1 for  C( 12)- O( 16)
================================
State has  OMEGA= 0   and energy asymptote:   Y(lim)= 89605.16398(cm-1)
Perform  4-point piecewise polynomial interpolation over   30 input points
Interpolation actually performed over modified input array:   Y(I) * r(I)**2
To make input points Y(i) consistent with  Y(lim),  add  Y(shift)=      0.0000
Scale input points:  (distance)* 1.000000000D+00  &  (energy)* 1.000000000D+00
             to get required internal units  [Angstroms & cm-1 for potentials]
      r(i)         Y(i)        r(i)         Y(i)        r(i)         Y(i)
   ----------------------   ----------------------   ----------------------
   1.11100000    219.9794   1.12100000     58.9728   1.13100000      0.4828
   1.11200000    199.0415   1.12200000     48.6136   1.13200000      0.0000
   1.11300000    179.2230   1.12300000     39.2860   1.13300000      0.4609
   1.11400000    160.4798   1.12400000     30.9898   1.13400000      1.8875
   1.11500000    142.7902   1.12500000     23.6594   1.13500000      4.2139
   1.11600000    126.1979   1.12600000     17.3604   1.13600000      7.4841
   1.11700000    110.6591   1.12700000     12.0272   1.13700000     11.6761
   1.11800000     96.1738   1.12800000      7.6816   1.13800000     16.7898
   1.11900000     82.7200   1.12900000      4.3236   1.13900000     22.8034
   1.12000000     70.3197   1.13000000      1.9314   1.14000000     29.7169
----------------------------------------------------------------------------
Extrapolate to   X .le. 1.1120  with
   Y=     -171.682  +1.264950D+29 * exp(-5.494114D+01*X)
Function for  X .GE.  1.1390   generated as
   Y=  89605.1640 - ( 8.191848D+09) * r** 13.038243  * exp{-(  11.519252*r)}
----------------------------------------------------------------------------

Calculate properties of the single potential described above
Potential-1 uses inner boundary condition of  zero value  at  RMIN

Eigenvalue convergence criterion is   EPS= 1.0D-06(cm-1)
Airy function at 3-rd turning point is quasibound outer boundary condition

Since state-1 has (projected) electronic angular momentum  OMEGA= 0
         eigenvalue calculations use centrifugal potential  [J*(J+1) - 0]/r**2


Solve for the   4 vibration-rotation levels of Potential-1:
   (v,J) =  (  0,  0)  (  1,  0)  (  2,  0)  (  3,  0)
   and automatically increment  J  in steps of  1 to a maximum value of   1

Matrix element arguments are powers of the distance  r (in Angstroms)
Coefficients of expansion for radial matrix element/expectation value argument:
       1.000000D+00 -2.000000D-01  3.000000D-02 -4.000000D-03


应该是我输入文件哪里写的不对，没有理解手册的意思。但不知道哪里写错了。求大家指教。谢谢。