行家们能不能帮我看看我程序问题呢，需要求出青藏高原上空200hpa大于等于30m/s的最...

adeline

fortran程序如下，给的是格点数据，在grd文件中已将青藏高原的格点范围处理好了，附上处理数据，需要求出66个月200hpa上纬向风大于等于30m/s的格点的最西经度，然后画出时间与经度折线图，观察急流中心的变化，但是程序运行过程中无错误，结果运行不出来，如图所示
program main
implicit none
integer i,j,t,k,ny,nx,nt
parameter (nx=14,ny=7,nt=66)
real uwnd(nx,ny,nt),a(98,nt),b(nt),m(nt)
open(1,file='e:\jingdu\uwnd200.grd',form='binary')
open(2,file='e:\jingdu\uwnd30.dat',form='binary')
do t=1,nt
    do j=1,ny
    do i=1,nx
   read(1,100)uwnd(i,j,t)
  if (uwnd(i,j,t)>=30.0) then
  k=i*j
   a(k,t)=(i-1)*2.5+72.5
     k=k+1
  else
    a(k,t)=0
    k=k+1
  endif
  enddo
enddo
enddo
100 format()
  !!!!!!!!!!!!!!!!!!!!!!!!!!!
do t=1,nt
    m(t)=a(1,1)
  do i=1,98
   if (m(t)>=a(i+1,t)) m(t)=a(i+1,t)
   enddo
   enddo
write(2,*)(m(i),i=1,nt)
end
*grd的gs文件
'reinit'
'sdfopen e:\data\uwnd.mon.mean.nc'
'set gxout fwrite'
'set fwrite e:\data\uwnd200.grd'
it=6
while(it<=794)
'set lon 72.5 105'
'set lat 25 40'
'set z 10'
'set t 'it''
'd uwnd'
it=it+12
endwhile
'disable fwrite'
;

lqouc

请补图

adeline

运行的结果如图

adeline

lqouc 发表于 2015-5-5 12:40
请补图

不好意思哈，图在下面

lqouc

精简了一下，顺便给你加上了缺省的设置，试试吧
program main
implicit none
integer i,j,t,k,ny,nx,nt
parameter (nx=14,ny=7,nt=66)
real uwnd,a
open(1,file='e:\jingdu\uwnd200.grd',form='binary')
open(2,file='e:\jingdu\uwnd30.dat',form='binary')
do t=1,nt
    do j=1,ny
    a=-999.0
    do i=1,nx
    read(1)uwnd
    if (uwnd>=30.0) then
     a=(i-1)*2.5+72.5
     write(2)a
     goto 100
    endif
   enddo
write(2)a
100 continue
enddo
enddo
close(2)
close(1)
end

adeline

lqouc 发表于 2015-5-5 12:49
精简了一下，顺便给你加上了缺省的设置，试试吧
program main
implicit none

还是不行，出现了一样的问题

lqouc

adeline 发表于 2015-5-5 12:58
还是不行，出现了一样的问题

你直接用的我给改过的程序？我把你的format删掉了。不应该还出现一样的报错啊。

adeline

lqouc 发表于 2015-5-5 13:08
你直接用的我给改过的程序？我把你的format删掉了。不应该还出现一样的报错啊。

重试了下，可以运行，但是写出的dat如何识别呢？

lqouc

adeline 发表于 2015-5-5 13:29
重试了下，可以运行，但是写出的dat如何识别呢？

你不是要给出时间纬度剖面么，输出的直接就是啊，ctl部分你完善下就行
undef -999.0
xdef  1 linear 0 1
ydef  14 linear 这里写你自己的起点和步长
zdef  1 levels 1000
tdef  66 linear 这里写你自己的起点和步长