版塊導(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í)熱門版塊排行榜

返回列表

下一頁(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ù)此樓

» 猜你喜歡

深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿）已經(jīng)有10人回復(fù)
278求調(diào)劑已經(jīng)有5人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有15人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有11人回復(fù)
334求調(diào)劑已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有5人回復(fù)
293求調(diào)劑已經(jīng)有16人回復(fù)
302求調(diào)劑已經(jīng)有7人回復(fù)
312求調(diào)劑已經(jīng)有6人回復(fù)
281求調(diào)劑（0805）已經(jīng)有9人回復(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的小程序，編譯通不過(guò)，請(qǐng)大俠幫忙已經(jīng)有59人回復(fù)
【求助】fortran出錯(cuò)了已經(jīng)有3人回復(fù)
【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來(lái)吧已經(jīng)有105人回復(fù)
【求助】求Fortran90編譯器？已經(jīng)有11人回復(fù)
【求助】fortran出錯(cuò) 編譯通過(guò) 運(yùn)行溢出已經(jīng)有10人回復(fù)
【求助】FORTRAN編譯通過(guò)，但是運(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, *) 這種方式來(lái)讀取吧。 * 只適合于數(shù)值，而且讀出來(lái)的格式只能是浮點(diǎn)數(shù)……

如果你的 msp 是用來(lái)存單個(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

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

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

贊一下(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ò)了，呵呵，金幣送你了，常來(lái)程序語(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)通過(guò)編譯了
拿到自己想要的文件了就不管那么多了呵呵

贊一下

回復(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)通過(guò)編譯了
拿到自己想要的文件了就不管那么多了呵呵

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

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

這個(gè)就不清楚了…… 我不是專業(yè)編程的
real x y z 這樣的語(yǔ)句，你本來(lái)聲名了三個(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ā)表

[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	3/150	2026-03-18 00:13 by PY-Hzb
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	5/250	2026-03-17 22:01 by 睿08
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 344求調(diào)劑 +4	knight344 2026-03-16	4/200	2026-03-17 17:27 by ruiyingmiao
[考研] 本人考085602 化學(xué)工程專碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[基金申請(qǐng)] 有必要更換申報(bào)口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12	11啦啦啦 2026-03-11	35/1750	2026-03-13 21:25 by JourneyLucky
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長(zhǎng)
[考研] 283求調(diào)劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao