作者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