奉献个给大家一个计算副高指数的程序，自己编的，请尊重他人劳动成果为谢，(*^__^*)！

123349818

支持楼主啊。正需要这方面的资料呢。呵呵。

Aaron1987

膜拜一下！学习学习

★翼★

很给力啊，以前要是找到了就好了

wuyejincao

太坏了，回帖可见。。。

阿呗

如果是闰年的话，一年就是366天，是不是还是需要做些改动的

丹玛斯.懋

原创的才是最经典的

阿呗

我心飞翔(ab) 发表于 2012-8-2 22:53

                               
登录/注册后可看大图

楼主，你好！我用该程序计算的出来东东，貌似不正确呀！如我用ncep 2.5*2.5 的资料 算出的 1980年的7月逐日 ...

请问你现在改好了么，我貌似算出来结果也不对

ayzqs

谢谢分享，呼呼

quwj

在论坛中看到lotina 分享的计算副高逐日指数的程序，个人觉得是否计算脊线位置指数和北界位置指数的部分做些改动才准确，像下面这样？多谢！

也就是在把格点数换算成纬度时，格点index减一在求出平均值之后进行。

请lotina 或者看过这个程序的朋友们给点建议，多谢啊！

!!!!!!!!!!!!!!!!!!!!!!!! 4 脊线位置指数 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do t=1,yr
do m=1,mo
gx4(m,t)=0.0;tt=0
do i=lon1,lon2
if(n588(i,m,t)/=0 .and. s588(i,m,t)/=0)then
tt=tt+1
gx4(m,t)=gx4(m,t)+(n588(i,m,t)+s588(i,m,t))/2.0
end if
end do
if(gx4(m,t)/=0)then
gx4(m,t)=(gx4(m,t)/float(tt)-1)*2.5-90.0
end if
end do
end do
print*, 'ok gx4'
!!!!!!!!!!!!!!!!!!!!!!!! 5 北界位置指数 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do t=1,yr
do m=1,mo
gb5(m,t)=0.0;tt=0
do i=lon1,lon2
if(n588(i,m,t)/=0.0)then
tt=tt+1
gb5(m,t)=gb5(m,t)+n588(i,m,t)
end if
end do
if(gb5(m,t)/=0)then
gb5(m,t)=(gb5(m,t)/float(tt)-1)*2.5-90.0
end if
end do
end do
print*, 'ok gb5'

carollee

谢谢楼主分享