作者juangpeiyi (燒瓶中的小人 Homunculus)
看板Fortran
標題[請益] 這輸出格式的意義
時間Wed Dec 1 13:58:14 2010
擷取一部分放上來問
nit=nit+1
i1 = mod(nit,10)
i2 = mod(nit,100)/10
i3 = mod(nit,1000)/100
i4 = mod(nit,10000)/1000
i5 = mod(nit,100000)/10000
i6 = mod(nit,1000000)/100000
i7 = mod(nit,10000000)/1000000
i8 = mod(nit,100000000)/10000000
i9 = nit/100000000
ca1 = CHAR(i1+48)
ca2 = CHAR(i2+48)
ca3 = CHAR(i3+48)
ca4 = CHAR(i4+48)
ca5 = CHAR(i5+48)
ca6 = CHAR(i6+48)
ca7 = CHAR(i7+48)
ca8 = CHAR(i8+48)
ca9 = CHAR(i9+48)
!上面再補一些在完整一點
suffix= '.dat'
pt1 = 'shape'//ca9//ca8//ca7//ca6//ca5//ca4//ca3//ca2//ca1//suffix
pt2 = 'pc'//ca9//ca8//ca7//ca6//ca5//ca4//ca3//ca2//ca1//suffix
pt3 = 'Tc'//ca9//ca8//ca7//ca6//ca5//ca4//ca3//ca2//ca1//suffix
pt4 = 'Mc'//ca9//ca8//ca7//ca6//ca5//ca4//ca3//ca2//ca1//suffix
do i=1, 100, 1
t=s(i)
if (t<3) cycle
open(101,file=pt1,status='unknown')
write(101,*) 'VARIABLES= "r" , "Pc" , "z"'
write(101,*) 'ZONE I=',t+1,',J=1,f=point'
do l=1, t+1, 1
write(101,*) rn(i,l), Cn(i), zn(i,l)
write(*,*) rn(i,l), Cn(i), zn(i,l)
end do
end do
write(101,*) pn_0, Mn_1, C, c_1, c_2, A, Hi, gamma
close(101)
.
.
. (略)
F7 F5後
都不會跑出data
連把 pt1.dat 開(新增)出來都沒有
目前我比較疑惑的是
很多 雙斜線那幾行 我不是很確定那是在幹嘛
翻了書沒有翻到 (彭國倫f95那本)
不知道各位高手能指教一下嗎 ?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.117.58.26
1F：→ charlesdc:好像是相對還是絕對路徑的東東 不然就是為了方便而已 12/01 14:55
2F：→ gilocustom:雙除號是字串相連用的。 12/01 16:01
3F：→ charlesdc:我也不懂>< 可以再詳細說明一下他在這個範例中的用法嗎? 12/01 16:26
4F：推 Cypresslin:你的ca9~ca1應該是字串變數，//的功用就同二樓所說 12/01 19:07
5F：→ Cypresslin:若宣告character(len=5) :: StrA="Hello",StrB="World" 12/01 19:08
6F：→ Cypresslin:print*,StrA,StrB 跟 print*,StrA//StrB 效果是一樣的 12/01 19:09
7F：→ Cypresslin:但是這裡是要給定pt1變數一個檔案名稱，此時就不能使用 12/01 19:10
8F：→ Cypresslin:pt1=StrA,StrB的方式來讓pt1等於HelloWorld，要用// 12/01 19:11
9F：→ sjgau:// 的意思是 字串的串接，但是不好用。最好是 自己寫副程式 12/01 21:11
10F：→ juangpeiyi:為什麼這情況就必須用// 我還是不太了解 我再想想看 12/01 22:12
11F：→ juangpeiyi:謝謝喔^^" 12/01 22:12
12F：→ acyang:ca1紀錄個位數，ca2紀錄十位數，以此類推， 12/02 16:10
13F：→ acyang:最後檔名就會像是pc000003457.dat這樣，方便整理而已 12/02 16:12
14F：→ acyang:看這程式像是要用tecplot顯示等高線圖 12/02 16:12
15F：→ juangpeiyi:哇哇 acyang大好內行 我再研究看看 T^T 12/03 03:02
16F：→ juangpeiyi:但他不會出來dat擋 12/03 03:02
※ 編輯: juangpeiyi 來自: 140.117.59.134 (12/03 18:06)
17F：→ acyang:你應該檢查t=s(i)的值，if (t<3) cycle這行可能在搞鬼 12/04 16:16
18F：→ juangpeiyi:我檢查了!!!! 竟然t都是21這兩個值在跳 !!!! 12/04 16:52