[請益]請救命~~我的fortran
Program CNT_dilute
implicit none
real:: cm(2394)
real::cor2(216)
real::a
integer::i,j,k,n
open(unit=10,file="D:\RECGLD2\cm.txt")
open(unit=70,file="D:\RECGLD2\cluster.txt")
read (10,*) cm
Write(70,'(1x,6F12.6)') cm
open(unit=100,file="D:\RECGLD2\cor2.txt")
open(unit=40,file="D:\RECGLD2\cnt.txt")
do i=1,n,6
do j=1,120,6
if((cm(i).GT.cor2(j)).and.(cm(i).LT.cor2(j+1))
&.and.(cm(i+1).GT.cor2(j+2)).and.(cm(i+1).LT.cor2(j+3))
&.and.(cm(i+2).GT.cor2(j+4)).and.(cm(i+2).LT.cor2(j+5)))then
write(40,'(1x,6F12.6)')cm
end if
end do
end do
stop
end
程式碼如上
我現在有一組座標的data
存成
x y z
cm(1) cm(2) cm(3) cm(4) cm(5) cm(6)
cm(7) cm(8) cm(9) cm(10) cm(11) cm(12)
cm(13) cm(14) cm(15) cm(16) cm(17) cm(18)
......
要判斷滿足
cor2(1)<cm(1)<cor2(2)
cor2(3)<cm(2)<cor2(4)
cor2(5)<cm(3)<cor2(6)
則寫入檔案
希望表示成---------------------------------------
(區間A)
滿足的座標
(區間B)
滿足的座標
寫再同一個檔案~~但是不同組空行
同組不空
不知道為什麼~~我寫的判斷式好像沒用
只會一直print一樣的data~~不會幫他們分組
有沒有人可以救救我...已經弄了一個多星期了
拜託~~~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.227.132.164
討論串 (同標題文章)
以下文章回應了本文:
完整討論串 (本文為第 1 之 3 篇):