看板 Fortran 關於我們 聯絡資訊
程式碼如下,如果沒加pgplot的那些指令就很正常,而且可以跑 但一加下去後他就說我第五行宣告錯誤 讀的檔裡面的確都是浮點數沒錯.... program test integer :: i,pgopen istat=PGOPEN('raman.ps/vcps') if(istat<=0)stop' ERR opening for PS file!' real :: shift1(1423),count1(1423),shift2(305),count2(305) open(1,file='ethanol.txt',status='old') do i=1,1423 read(1,'(f7.2,1x,f7.3)') shift1(i),count1(i) end do do i=1,305 read(1,'(f7.3,1x,f7.3)') shift2(i),count2(i) end do call pgsci(1) call pgslw(4) call pgsch(2.0) call pgscf(3) call pgenv(500.323,3001.44,0,1000,0,1) call pglab('Raman shift(Reciprocal cm)','count','Raman spectrum of Ni-GEM') call pgslw(2) call pgsci(2) do i=1,1423,1 call pgline(1423,shift1(i),count1(i)) end do do i=1,305 call pgline(305,shift2(i),count2(i)) end do call pgend end 以上,感謝大家!! -- ◤ ◥ 塞魯表示:C字褲我14年前就在穿了! ﹨════│ ═════/ 2002十八號離子燙超直髮 ◥ ▂▂▃▅▆ 2003悟空兩件式混搭 ˙ | 2004達爾靴 ██▆▅▅▆ / 2005特南克斯短版外套 φJeans1020 ▅▃▂ ◣ ╱_▂▃ ████ 2009賽魯C字褲 鳥山明 引領流行20年 ◥◣ ▄▂▁ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.187.84
janhon:if 要放在宣告後面 12/03 22:32
meteorology:應該是你的ISTAT沒宣告~ 12/04 12:51
blc:一樓對了,宣告結束之後才能做其他事。 12/04 21:16
meteorology:ISTAT沒宣告也可以喔?? 學到了!! 12/06 01:06
blc:因為沒有寫implicit none,i開頭變數會預設成整數。 12/06 09:09
blc:好不好在其次,不會有錯誤訊息就是了。 12/06 09:10
meteorology:阿對齁 自從習慣加implicit none 之後 都忘了這個 12/10 23:49