24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

贺仪

铁杆木虫 (著名写手)

应助: 93 (初中生)
贵宾: 0.636
金币: 5937
散金: 1063
红花: 38
沙发: 2
帖子: 2499
在线: 588.7小时
虫号: 995125
注册: 2010-04-12
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] fortran编译 read出错完毕

错误提示
POSITION.F(31): error #6358: Constants and expressions are invalid in read-only I/O lists. [MSP]
   read(15,*) (msp(i),i=1,10)
前面msp有定义
character msp(10)
  read(15,*) (msp(i),i=1,10)
语法有问题吗？

一下是全部程序
   program getpos
ccccccccc  12 is SPPARKS OUTPUT  cccccccc
ccccccccc 15 is the molecure  ccccccccc
ccccccccc  13 is OUTPUT ccccccccccc
   integer mnum nsp
   real coord(2000,3)
   real species(10,3)
   integer nsp(10)
   real x y z
   character tile cellp cella cellb cellz msp(10)
   open(12,file='POSCAR',status='old')
   open(15,file='POSCAR',status='old')
   open(13,file='CONTCAR',status='new')
ccccccccccccccccccc sys param a b c cccccccccccccc
   read(12,*) tile
   write(13,*) tile
   read(12,*) cellp
   write(13,*) cellp
   read(12,*) cella ;
   write(13,*) cella
   read(12,*) cellb ;
   write(13,*) cellb
   read(12,*) cellc ;
   write(13,*) cellc
   read(15,*);
   read(15,*);
   read(15,*);
   read(15,*);
   read(15,*);
cccccccccccccccccccc atom char ccccccccccccccccccccccc
   read(15,*) (msp(i),i=1,10)
   write(13,*) (msp(i),i=1,10)
   read(12,*);
cccccccccccccccccccc atom num  cccccccccccccccccccccccc
   read(12,*) mnum
   read(15,*) (nsp(i),i=1,10)
   DO i=1, 10
      N_anum=mnum*nsp(i)
      write(13,N_anum)
   END DO
   read(12,*)
   read(15,*)
   write(13,*) 'R'
cccccccccccccccc  read molecue coord cccccccccccc
   DO i=1, mnum
      read(12,*) (coord(i,j),j=1,13)
   ENDDO
   close(12)
cc
   DO j=1, 10             !!! species cyc !!!!!!
         DO i=1, nsp(j)       !!! the j atom cyc !!!!!!
               read(15,*) x y z  !!!!!!!!!!!!!!!!!!!!!!!! molecure coord !!!!!!!!!!!!!
               DO k=1, mnum
                  nx=x+coord(k,1)
                  ny=y+coord(k,2)
                  nz=z+coord(k,3)
                  write(13,*) nx ,ny ,nz
               ENDDO
         ENDDO
   ENDDO
   close(15)
   close(13)
   END

[ Last edited by 贺仪 on 2011-8-23 at 14:53 ]

回复此楼

» 猜你喜欢

嘟啦啦嘟啦啦嘟

1楼 2011-08-23 09:57:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

★ ★
余泽成(金币+2): 辛苦了！ 2011-08-24 16:22:27

real x y z
改为
real x, y, z
或者
real :: x, y, z

你的文件之所以没有什么问题，是因为即使不声明，fortran 默认 x, y, z 是浮点数……

但如果加了，implicit none， fortran 则要求所有的变量必须声明……

赞一下(1人)

回复此楼

11楼2011-08-23 16:50:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 12 个回答

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

【答案】应助回帖

★ ★
贺仪(金币+5): 多谢指教！ 2011-08-23 14:53:23
余泽成(金币+2): 谢谢参与应助！ 2011-08-24 16:20:44

引用回帖:

1楼: Originally posted by 贺仪 at 2011-08-23 09:57:29:
错误提示
POSITION.F(31): error #6358: Constants and expressions are invalid in read-only I/O lists. [MSP]
read(15,*) (msp(i),i=1,10)
前面msp有定义
character msp(10)
read(15,*) (msp(i ...

character 类型数据不能用 read(15, *) 这种方式来读取吧。 * 只适合于数值，而且读出来的格式只能是浮点数……

如果你的 msp 是用来存单个字符（10个？），那么你可以用下面的代码试试（没有测试，不保证正确……）

CODE: