- 积分
- 5402
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-4-16
- 最后登录
- 1970-1-1
![[NightWind陈] 粉丝数:425 微博数:2637 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
![](static/image/common//ico_lz.png)
楼主 |
发表于 2013-5-20 20:53:08
|
显示全部楼层
虫儿飞 发表于 2013-5-20 18:19
最好用风场和高度场结合起来比较好,个人体会
您好,可以帮我看下我的程序哪里有问题吗?我只用了风场的资料,03年6~8月份逐日资料。
program fugao
implicit none
integer,Parameter::x=37,y=25,mo=92
integer,Parameter::lon1=37,lon2=73,lat0=37,lat60=61
integer i,j,m,tt
real h(x,y,mo)
real gj(mo)
open(10,file='f:/fugao/hgt.2003.dat')
read(10) (((h(i,j,m),i=1,x),j=1,y),m=1,mo)
close(10)
!!!!!!!!!!!!!!!!!!!!!!!! 西脊点指数!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do m=1,mo
gj3(m)=0.0;tt=0
do i=lon1,lon2
do j=lat0,lat60
if(h(i,j,m)/10.0>=588 .and. h(i+1,j,m)/10.0>=588)then
tt=i
exit
end if
end do
if(tt/=0)then
exit
end if
end do
gj(m)=(tt-1)*2.5
if(gj(m)<=0)then
gj(m)=180
end if
open(11,file='f:/fugao/jidian.2003.dat')
write(11,*) gj3(m)
end do
stop
end
|
|