傲娇师兄开贴帮师妹调程序（暑假归来，继续）

vip545

哇哦，太感谢师兄的指导了

绯宜物语

师兄，可以帮我看一下这个程序吗？
program tem
integer*4,parameter:: x=144,y=73,z=16,t=12
integer*4 :: ix,iy,iz,it
real*4,parameter :: undef = -9.99E+33
real*4 var(x,y,z,t)
open(1,file='f:\t.grd',form='binary')
do it=1,t
do iz=1,z
do iy=1,y
do ix=1,x
read(1) var(ix,iy,iz,it)
enddo
enddo
enddo
enddo
open(2,file='f:\t0.grd',form='binary')
do it=1,t
do iy=1,y
do ix=1,x
iz=1
while ((var(ix,iy,it,iz) .GE. 0) .and. (var(ix,iy,it,iz+1) .LE. 0))
iz=iz+1
write(2) iz
enddo
enddo
enddo
close(1)
close(2)
end
我是要判断零度层临近两层的气压层，但是while那有问题。

虫儿飞

绯宜物语 发表于 2014-9-17 19:43
师兄，可以帮我看一下这个程序吗？
program tem
integer*4,parameter:: x=144,y=73,z=16,t=12

最好用 两个数相乘小于0作为判定条件

绯宜物语

虫儿飞 发表于 2014-9-17 20:34
最好用 两个数相乘小于0作为判定条件

嗯嗯，那这个地方除了这个怎么再修改呀？while var(ix,iy,it,iz)*var(ix,iy,it,iz+1) .LE. 0
提示E:\fortran\25\t0.f90(21): error FOR2934: lexical error: Unrecognized statement
E:\fortran\25\t0.f90(21): error FOR3852: syntax error

还有就是变量名和表示维度的量相同，怎么办？

kaidada

二师兄，师傅都被妖怪抓走了，你还有时间调戏妹子

金色牛神

哈哈……一起努力啦

atomlrhy

师兄给力啊

user_qxjy

这种“师兄精神”值得发扬广大，衷心赞一个！