已编写好求青藏高原上上空200hpa纬向风30m/s出现的最西经度的fortran程序

adeline

fortran程序，我去验证了一些，发现误差很大
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_1.grd',form='binary')
open(2,file='e:\jingdu\uwnd30_1.grd',form='binary')
open(3,file='e:\jingdu\uwnd30_1.txt')
do t=1,nt
    do j=1,ny
    do i=1,nx
    k=14*(j-1)+i
   read(1)uwnd(i,j,t)
   write(*,*)'ok'
  if (uwnd(i,j,t)>=30.0) then
   a(k,t)=(i-1)*2.5+72.5
   write(*,*)a(k,t)
  else
    a(k,t)=0
  endif
  enddo
enddo
enddo
! 100 format()
  !!!!!!!!!!!!!!!!!!!!!!!!!!!
do t=1,nt
   do i=1,97
    if(a(i,t)/=0) m(t)=a(i,t)
   enddo
  do j=i+1,97
   if (m(t)>=a(i+1,t).and.a(i+1,t)>0.) m(t)=a(i+1,t)
   enddo
   enddo
write(2)(m(i),i=1,nt)
write(3,100)(m(i),i=1,nt)
write(*,*)m
100 format(f7.2)
end
*生成grd文件的gs文件，提取青藏高原6月份资料
'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'
;

言深深

你已经发过两次类似的帖子，自己去主动编写程序，这很好。
如下几个问题，或许对你有益：
1、你说的纬向风，最西边，是不是就是你14个格点中当纬向风速高于30就标记，那么有两种方案：
a)从西边往东边找，找到大于30的，就结束这个循环进行下一个纬度寻找；
b)从东边往西边找，找到所有大于30的，标记最后一个。
似乎你是用a进行寻找的，但是标记的时候记录了所有的，不知道你能不能分清楚哪一个对应哪一个？
2、a(k,t)=(i-1)*2.5+72.5，将经纬度全都合到一起，有点儿眉毛胡子一把抓的感觉，你难道不希望标记一下，到底是哪一个纬度上面的哪一个经度是你需要的吗？
3、你最终输出的是什么，不是很明白。