- 积分
- 24
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-5-3
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
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'
;
|
|