作者balloom (叶子的离开...)
看板Programming
标题[请益]请救命～～我的ｆｏｒｔｒａｎ
时间Sun Jun 15 20:19:00 2008
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