阻塞高压识别，由于从来没有接触过，不理解~！求源程序

飞儿ＬＯＶＥ

基于Tibaldi and Molteni (1990)的阻高识别技术方法，建立阻高指数序列，该阻塞高压识别方法在业务上较成熟，在南北半球中高纬度具有普遍的适用性。阻塞高压客观确定是局地和瞬时的，其判别标准如下：对500hPa位势高度取5天滑动平均，然后对每个经度进行位势高度的梯度计算。向南方向的500hPa高度梯度GHGS和向北方向的500hPa高度梯度GHGN计算如下：
GHGS=(Z(φo)--Z(φs))/φo-φs
GHGN=(Z(φn)--Z(φo))/φn-φo
Z：位势高度（gpm）
φ：纬度
φn=80+δ，φ0=60+δ，φs=40+δ，δ=5,0,-5
对某时某经度任意一个δ值，如果条件满足：（i）GHGS>0;（ii）GHGN<-10即即确认该时（取5天中间的一天为该时刻）该经度受阻塞高压影响。
program main
implicit none
integer i,j,k,l,m,n,it,nx,ny,nz,nt,dy,dyy,yr,cir,yrr
parameter(nx=45,ny=25,nz=1,dy=365,yrr=15,nt=5479)
real hgt4(nx,ny,nz,dy,yrr),hgt5(nx,ny,nz,nt)
!位势高度五天滑动平均
open(10,file='E:\lwshuju\exdata\11.grd',form='binary',status='old')
open(11,file='E:\lwshuju\111.txt')
read(10) (((((hgt4(i,j,k,l,m),i=1,nx),j=1,ny),k=1,nz),l=1,dy),m=1,yrr)
do m=1,yrr
  do l=1,dy
   do k=1,nz
    do j=1,ny
     do i=1,nx
     hgt5(i,j,k,l+(m-1)*365)=hgt4(i,j,k,l,m)
     enddo
    enddo
   enddo
  enddo
enddo
do it=3,nt-2
  do k=1,nz
    do j=1,ny
     do i=1,nx
     hgt5(i,j,k,it)=(hgt5(i,j,k,it-2)+hgt5(i,j,k,it-1)+hgt5(i,j,k,it+1)+hgt5(i,j,k,it+2))/4
     enddo
    enddo
  enddo
enddo
do it=3,nt-2
  do k=1,nz
    do j=1,ny
     do i=1,nx
        write(11,*) hgt5(i,j,k,it)
     enddo
    enddo
  enddo
enddo
end

!!!!!!!!!!!!!!!!!!!!!!!!!!!阻塞高压识别
real ghgs(i,j,k,it),ghgn(i,j,k,it)
do it=1,nt
   do k=1,nz
    do j=11,15,2
     do i=1,nx
    ghgs(i,j,k,it)=(hgt5(i,j,k,it)-hgt5(i,j-8,k,it))/20
    ghgn(i,j,k,it)=(hgt5(i,j+8,k,it)-hgt5(i,j,k,it))/20
      enddo
     enddo
    enddo
enddo
下面就不会写了~

qxtlyf

Tibaldi and Molteni (1990)的阻高识别技术方法   可以介绍一下

承德气象

同求啊，我现在也在做这个，可是不会弄，不知道楼主现在弄好没啊，能不能发给我呢?

承德气象

承德气象 发表于 2012-5-2 09:56

                               
登录/注册后可看大图

同求啊，我现在也在做这个，可是不会弄，不知道楼主现在弄好没啊，能不能发给我呢?

我的邮箱是805171539@qq。com

惜梦逐水

楼主的原程序编好了吗，求共享，还有Tibaldi and Molteni (1990)的阻高识别技术的文章能共享一下嘛，569669529@qq.com

松下风

谢谢楼主，好东西，学习了

惜梦逐水

松下风 发表于 2012-12-14 21:48

                               
登录/注册后可看大图

谢谢楼主，好东西，学习了

求发送啊，谢谢啦569669529@qq.com

ayzqs

哇，不错的东西

Journing90

好东西，求共享，谢楼主啊！

ZOE

不知道在这个基础上修改 能识别冷涡么