[問題] 陣列中的變數

看板Fortran作者band33 (時鐘)時間13年前 (2011/12/12 15:28)推噓1(1推 0噓 5→)

留言6則, 3人參與討論串1/1

我想要將檔案中的座標讀入二維的陣列中程式碼如下: program RotStanOrien implicit none type :: Orientation real :: x real :: y real :: z end type Orientation type(Orientation) :: XYZ character(len=80) :: Infilename, Outfilename integer :: NAtoms, MaxAP3, AtmChg, error, i real :: InOri(NAtoms, 3) integer :: A(MaxAP3, 3) integer :: TrVec(3) NAtoms = 0 write(*,*) 'Please input the filename:' read(*,*) Infilename open(unit=99, file=Infilename, status="old", iostat=error ) if ( error /= 0 ) then write(*,*) 'Open file fail' stop end if do while( .true. ) read(99,*, end=500) AtmChg, XYZ%x, XYZ%y, XYZ%z NAtoms = NAtoms + 1 InOri(NAtoms, 3) = XYZ%z InOri(NAtoms, 2) = XYZ%y InOri(NAtoms, 1) = XYZ%x write(*,'F6.3,F6.3,F6.3') XYZ%x, XYZ%y, XYZ%z enddo 500 continue write(*,*) '---------------------------------' do i = 1, NAtoms write(*,'F6.3,F6.3,F6.3') InOri(i,1), InOri(i,2), InOri(i,3) enddo 這是輸出的結果: 0.050 0.060 0.040 0.010-1.000 0.020 虛線上半部是正確的座標，這是由第一個迴圈內 1.000 0.015 0.016 的XYZ%x、XYZ%y、XYZ%z三個變數所write出來的 --------------------------------- 0.050 0.050 0.050 下半部的座標則是錯的，這裡是第二個迴圈所write 0.010 0.010 0.010 出來的，每一列的三個值都是XYZ%x的值，我要的 1.000 1.000 1.000 是像上半部的值請問為什麼將陣列中的值設為XYZ的變數後就會跑掉？謝謝 -- -- ╭◢██◣═══╮ ╭══════╮ theanswer3 ╭═════╮ ║██◣◥ ◣ ╰══════╯◢█◣ ◢▌ ╰═════════╯ ω╴ ║ ║◥██◣ ◥█◣ ◢◢◣ ◢█◣ ◥◣◥ ██◣ ◢ ◢ ◢◢◣ ◢ ◢▌︿◤￣ο ║ ║◣◥theanswer3 █◤ █◢█ ◣◥◣theanswer3◢█ █◤ █▌█ ▄ ║ ║◥██◤ ◥◆ █ ◥◤◥ ◥█◤ ◤ ◤ ◥◤◥ █ ◥██ ╱ ﹨ ║ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.115.43.188