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 ]

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

急求可用的fortran编译器已经有13人回复
如何运行Fortran编译好的exe并自动输入文已经有11人回复
请教一个fortran小程序编译出错的问题，谢谢已经有9人回复
写了一个fortran90的小程序，编译通不过，请大侠帮忙已经有59人回复
【求助】fortran出错了已经有3人回复
【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧已经有105人回复
【求助】求Fortran90编译器？已经有11人回复
【求助】fortran出错编译通过运行溢出已经有10人回复
【求助】FORTRAN编译通过，但是运行不了已经有19人回复
【求助】Fortran编译目标文件后的链接问题已经有4人回复
求助学习linux下fortran编译哪本书好！已经有4人回复
【求助】win7 的用什么fortran 编译器啊？【已完成】已经有14人回复

嘟啦啦嘟啦啦嘟

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
专业: 污染生态化学

【答案】应助回帖

★ ★
贺仪(金币+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: