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

阿呗

如果是闰年的话，一年就是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

谢谢楼主分享

我心飞翔(ab)

carollee 发表于 2012-10-26 22:25

                               
登录/注册后可看大图

谢谢楼主分享

实际副高的情况很复杂（比如副高存双脊线，又比如副高断裂，等等），可能还要加诸多判断，才能反映真实的副高！

carollee

我心飞翔(ab) 发表于 2012-10-27 14:41

                               
登录/注册后可看大图

实际副高的情况很复杂（比如副高存双脊线，又比如副高断裂，等等），可能还要加诸多判断，才能反映真实的 ...

好的，谢谢

诗琳

思路清晰，适合菜鸟学习，好贴

一水天下

多谢楼主！