版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

下一頁(yè)

賀儀

鐵桿木蟲 (著名寫手)

應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊(cè): 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] fortran編譯 read出錯(cuò) 完畢

錯(cuò)誤提示
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)
語(yǔ)法有問題嗎？

一下是全部程序
   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 ]

回復(fù)此樓

» 猜你喜歡

286求調(diào)劑已經(jīng)有9人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可已經(jīng)有10人回復(fù)
291求調(diào)劑已經(jīng)有5人回復(fù)
308求調(diào)劑已經(jīng)有7人回復(fù)
085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富已經(jīng)有4人回復(fù)
壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！已經(jīng)有3人回復(fù)
一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑已經(jīng)有4人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

急求可用的fortran編譯器已經(jīng)有13人回復(fù)
如何運(yùn)行Fortran編譯好的exe并自動(dòng)輸入文已經(jīng)有11人回復(fù)
請(qǐng)教一個(gè)fortran小程序編譯出錯(cuò)的問題，謝謝已經(jīng)有9人回復(fù)
寫了一個(gè)fortran90的小程序，編譯通不過，請(qǐng)大俠幫忙已經(jīng)有59人回復(fù)
【求助】fortran出錯(cuò)了已經(jīng)有3人回復(fù)
【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來吧已經(jīng)有105人回復(fù)
【求助】求Fortran90編譯器？已經(jīng)有11人回復(fù)
【求助】fortran出錯(cuò) 編譯通過運(yùn)行溢出已經(jīng)有10人回復(fù)
【求助】FORTRAN編譯通過，但是運(yùn)行不了已經(jīng)有19人回復(fù)
【求助】Fortran編譯目標(biāo)文件后的鏈接問題已經(jīng)有4人回復(fù)
求助學(xué)習(xí)linux下fortran編譯哪本書好！已經(jīng)有4人回復(fù)
【求助】win7 的用什么fortran 編譯器�。俊疽淹瓿伞� 已經(jīng)有14人回復(fù)

嘟啦啦嘟啦啦嘟

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲號(hào): 183750
注冊(cè): 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★ ★
賀儀(金幣+5): 多謝指教！ 2011-08-23 14:53:23
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-08-24 16:20:44

引用回帖:

1樓: Originally posted by 賀儀 at 2011-08-23 09:57:29:
錯(cuò)誤提示
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 類型數(shù)據(jù)不能用 read(15, *) 這種方式來讀取吧。 * 只適合于數(shù)值，而且讀出來的格式只能是浮點(diǎn)數(shù)……

如果你的 msp 是用來存單個(gè)字符（10個(gè)？），那么你可以用下面的代碼試試（沒有測(cè)試，不保證正確……）

CODE:

read(15, '(a1)') (msp(i), i=1,10)

贊一下(1人)

回復(fù)此樓

2樓2011-08-23 14:31:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲號(hào): 183750
注冊(cè): 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★ ★
賀儀(金幣+5): 多謝指教！ 2011-08-23 14:53:34
余澤成(金幣+2): 辛苦了！ 2011-08-24 16:20:55

看來我前面的說法沒有切入要點(diǎn)，你的問題在于
real x y z
character tile cellp cella cellb cellz msp(10)
上面兩行各個(gè)變量之間應(yīng)該用逗號(hào)來分隔……

當(dāng)然，字符的確不能用 read(15,*) 這種方式來讀……

贊一下(1人)

回復(fù)此樓

3樓2011-08-23 14:47:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

賀儀

鐵桿木蟲 (著名寫手)

應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊(cè): 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
余澤成(金幣+1): 謝謝參與應(yīng)助！ 2011-08-24 16:21:05

不是讀取的原因
在定義的時(shí)候用“：：”后就可以了

贊一下(1人)

回復(fù)此樓

嘟啦啦嘟啦啦嘟

4樓2011-08-23 14:55:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲號(hào): 183750
注冊(cè): 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

引用回帖:

4樓: Originally posted by 賀儀 at 2011-08-23 14:55:38:
不是讀取的原因
在定義的時(shí)候用“：：”后就可以了

不知道你這個(gè)是啥意思？

你給出的代碼中太多不規(guī)范的東西，呵呵，反正我用 gfortran 會(huì)提示很多錯(cuò)誤……

贊一下

回復(fù)此樓

5樓2011-08-23 15:00:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

賀儀

鐵桿木蟲 (著名寫手)

應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊(cè): 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

余澤成: 剛才發(fā)錯(cuò)了，呵呵，金幣送你了，常來程序語(yǔ)言版！ 2011-08-24 16:22:03

引用回帖:

5樓: Originally posted by snoopyzhao at 2011-08-23 15:00:35:
不知道你這個(gè)是啥意思？

你給出的代碼中太多不規(guī)范的東西，呵呵，反正我用 gfortran 會(huì)提示很多錯(cuò)誤……

是的錯(cuò)誤很多這是最開始的現(xiàn)在調(diào)好了，程序如下：
   program getpos
ccccccccc  12 is SPPARKS OUTPUT  cccccccc
ccccccccc 15 is the molecure  ccccccccc
ccccccccc  13 is OUTPUT ccccccccccc
   integer mnum Nanum
   real::coord(2000,3)
   integer::nsp(4)
   real x y z
   character tile cellp
   character::msp(4)
   real::cella(3)
   real::cellb(3)
   real::cellc(3)
   open(12,file='POSCARS',status='old')
   open(15,file='POSCARV',status='old')
   open(13,file='CONTCAR',status='new')
ccccccccccccccccccc sys param a b c cccccccccccccc
   read(12,*) tile
   write(13,"(2XA10)"

tile
   read(12,*) cellp
   write(13,*) cellp
   read(12,*) (cella(i),i=1,3)
   write(13,"(3(2XF10.7))"

(cella(i),i=1,3)
read(12,*) (cellb(i),i=1,3)
write(13,"(3(2XF10.7))"

(cellb(i),i=1,3)
read(12,*) (cellc(i),i=1,3)
write(13,"(3(2XF10.7))"

(cellc(i),i=1,3)
   read(15,*)
   read(15,*)
   read(15,*)
   read(15,*)
   read(15,*)
cccccccccccccccccccc atom char ccccccccccccccccccccccc
   read(15,*) (msp(i),i=1,4)
   write(13,"(4(2XA2))"

(msp(i),i=1,4)
   read(12,*)
cccccccccccccccccccc atom num  cccccccccccccccccccccccc
   read(12,*) mnum
   read(15,*) (nsp(i),i=1,4)
   DO i=1, 4
         Nanum=mnum*nsp(i)
      write(13,10) Nanum
   END DO
   read(12,*)
   read(15,*)
   write(13,"(/,A1)"

'D'
   DO i=1, mnum
      read(12,*) (coord(i,j),j=1,3)
   END DO
   close(12)
cc
   DO j=1, 4             !!! species cyc !!!!!!
         DO i=1, nsp(j)       !!! the j species atom cyc !!!!!!
               read(15,*) x, y, z  !!!!!!!!!!!!!!!!!!!!!!!! molecure coord !!!!!!!!!!!!!
               DO k=1, mnum
            write(13,100) x+coord(k,1), y+coord(k,2), z+coord(k,3)
               END DO
         END DO
   END DO
   close(15)
   close(13)
10 FORMAT(2XI5.0,\)
100 FORMAT(3(2XF10.6))
   END

贊一下(1人)

回復(fù)此樓

嘟啦啦嘟啦啦嘟

6樓2011-08-23 15:06:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲號(hào): 183750
注冊(cè): 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★
賀儀(金幣+5): 多謝指教！ 2011-08-23 15:29:55
余澤成(金幣+1): 謝謝參與應(yīng)助！ 2011-08-24 16:22:51

至少
integer mnum Nanum
character tile cellp
還是不正確的吧……
另外，\ 似乎也有問題……

贊一下(1人)

回復(fù)此樓

7樓2011-08-23 15:10:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

賀儀

鐵桿木蟲 (著名寫手)

應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊(cè): 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

7樓: Originally posted by snoopyzhao at 2011-08-23 15:10:46:
至少
integer mnum Nanum
character tile cellp
還是不正確的吧……
另外，\ 似乎也有問題……

我用ifort -O2 已經(jīng)通過編譯了
拿到自己想要的文件了就不管那么多了呵呵

贊一下

回復(fù)此樓

嘟啦啦嘟啦啦嘟

8樓2011-08-23 15:14:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲號(hào): 183750
注冊(cè): 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★
賀儀(金幣+5): 多謝指教！ 2011-08-23 15:30:10
余澤成(金幣+1): 謝謝應(yīng)助到底，呵呵！ 2011-08-24 16:23:08

引用回帖:

8樓: Originally posted by 賀儀 at 2011-08-23 15:14:08:
我用ifort -O2 已經(jīng)通過編譯了
拿到自己想要的文件了就不管那么多了呵呵

你運(yùn)行一下看看得到的結(jié)果是不是你想要的，呵呵……

在 Fortran 中，像 real x y z 這樣的語(yǔ)句，你本來聲名了三個(gè)變量，但編譯之后是一個(gè)變量 xyz。你或許會(huì)說，為什么沒有報(bào)錯(cuò)，你試著在 program 那句之后加一句 implicit none 看看……

贊一下(1人)

回復(fù)此樓

9樓2011-08-23 15:16:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

賀儀

鐵桿木蟲 (著名寫手)

應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊(cè): 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

9樓: Originally posted by snoopyzhao at 2011-08-23 15:16:59:
你運(yùn)行一下看看得到的結(jié)果是不是你想要的，呵呵……

在 Fortran 中，像 real x y z 這樣的語(yǔ)句，你本來聲名了三個(gè)變量，但編譯之后是一個(gè)變量 xyz。你或許會(huì)說，為什么沒有報(bào)錯(cuò)，你試著在 program 那句之后加 ...

這個(gè)就不清楚了…… 我不是專業(yè)編程的
real x y z 這樣的語(yǔ)句，你本來聲名了三個(gè)變量，但編譯之后是一個(gè)變量 xyz 應(yīng)該怎么聲明呢？

我要處理的是坐標(biāo)文件，得到的文件打開是沒有問題的，我看了結(jié)構(gòu)也沒什么問題

贊一下

回復(fù)此樓

嘟啦啦嘟啦啦嘟

10樓2011-08-23 15:26:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主賀儀的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 086502化學(xué)工程342求調(diào)劑 +6	阿姨復(fù)古不過 2026-03-27	6/300	2026-03-28 07:06 by wangy0907
[考研] 085600 286分材料求調(diào)劑 +5	麻辣魷魚 2026-03-27	6/300	2026-03-28 04:51 by wxiongid
[考研] 085701環(huán)境工程求調(diào)劑 +9	多久上課 2026-03-27	9/450	2026-03-28 03:58 by fmesaito
[考研] 求調(diào)劑 +8	張zz111 2026-03-27	9/450	2026-03-28 03:41 by fmesaito
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 一志愿211院校 344分東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5	丶風(fēng)雪夜歸人丶 2026-03-26	8/400	2026-03-27 19:22 by 丶風(fēng)雪夜歸人丶
[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調(diào)劑 +3	Micky11223 2026-03-25	3/150	2026-03-27 14:52 by 醉在風(fēng)里
[考研] 329求調(diào)劑 +4	星野? 2026-03-26	4/200	2026-03-27 12:00 by 不吃魚的貓
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 336材料求調(diào)劑 +7	陳瀅瑩 2026-03-26	9/450	2026-03-27 00:20 by wxiongid
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid
[考研] 求調(diào)劑 +5	蘆lty 2026-03-25	6/300	2026-03-26 20:49 by 不吃魚的貓
[考研] 求調(diào)劑 +8	Auroracx 2026-03-22	8/400	2026-03-26 19:55 by 不吃魚的貓
[考研] 286求調(diào)劑 +13	Faune 2026-03-21	13/650	2026-03-26 19:52 by peike
[考研] 279求調(diào)劑 +6	紅衣隱官 2026-03-21	6/300	2026-03-26 18:32 by 不吃魚的貓
[考研] 309求調(diào)劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 311求調(diào)劑 +3	冬十三 2026-03-24	3/150	2026-03-24 21:31 by peike
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原

24小時(shí)熱門版塊排行榜

[求助] fortran編譯 read出錯(cuò) 完畢

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: