作者janewinnie (逃離)
看板Fortran
標題[問題] 實數轉字串 卻有錯誤
時間Fri Jun 26 23:59:35 2009
我想要在路徑下用迴圈開從-3.0, -2.8,...,0,...,2.8,3.0的.txt檔
我已經先建立好所需要用到的資料夾，包括"pk0.03","-3.0","-2.8".....
一開始都正常，但是到abil= -1.2以後，下一次的迴圈
就會變成abil=-0.9999996而不是-1.0，
到abil=0.0時就會出現錯誤，無法繼續跑下去
以下是我的程式碼
program main
implicit none
integer,parameter::ni=20
real::abil
character*7::r_name
character*2:: ni1
call int_to_cht(ni,ni1)
do abil=-3.0,3.0,0.2
call real_to_cht(abil,r_name,888)
open(250+abil,file="D:\pk0.03\"//TRIM(r_name)//"/SH_"&
//TRIM(r_name)//"_"//TRIM(ni1)//".txt",status='replace')
write(250+abil,"(g5.1,4(g10.1))") "迭代","Bias","RMSE","max_expo","not_use"
open(350+abil,file="D:\pk0.03\"//TRIM(r_name)//"/SH_1_"&
//TRIM(r_name)//"_"//TRIM(ni1)//".txt",status='replace')
write(350+abil,"(g5.1,3(g10.1))")"item","P(A)","P(S)","pk"
open(450+abil,file="D:\pk0.03\"//TRIM(r_name)//"/item_"&
//TRIM(r_name)//"_"//TRIM(ni1)//".txt",status='replace')
open(550+abil,file="D:\pk0.03\"//TRIM(r_name)//"/theta_hat_&
"//TRIM(r_name)//"_"//TRIM(ni1)//".txt",status='replace')
open(650+abil,file="D:\pk0.03\"//TRIM(r_name)//"/迭代次數 _&
"//TRIM(r_name)//"_"//TRIM(ni1)//".txt",status='replace')
end do !! abil
stop
end program
!-------------------------------------------
!實數轉字串
!-------------------------------------------
subroutine real_to_cht(abil,r_name,num)
implicit none
integer::num
real::abil
character*7::r_name
open(num,file="temp.txt")
write(num,"(f4.1)") abil
rewind(num)
read(num,"(a4)") r_name
close(num,status="delete")
return
end
!-------------------------------------------
!整數轉字串
!-------------------------------------------
subroutine int_to_cht(ni,ni1)
implicit none
integer::ni
character*7::ni1
open(987,file="temp.txt")
write(987,"(I2)") ni
rewind(987)
read(987,"(a2)") ni1
close(987,status="delete")
return
end
不知道是哪裡出了問題，
想請各位高手幫我看一下＠＠
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.117.239.136
1F：→ janewinnie:OPEN那邊因為版面關係,後面我用&換行,但fortran裡我都 06/27 00:06
2F：→ janewinnie:是打成一行 沒有分行 06/27 00:06
3F：→ snowanimal:把你的2.8*10吧,別用小數點 簡單的說就是用整數試試看 06/27 10:01
4F：→ janewinnie:解決了!!謝謝snowanimal :) 06/27 11:24