请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9094|回复: 13

[经验总结] T&M方法的阻塞高压指数(美国气候预测中心采用方法)

[复制链接]

新浪微博达人勋

发表于 2018-11-9 17:46:36 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
program dd
parameter(nx=144,ny=73,nz=21,nt=60)
real::h(nx,ny,nt),a,b,c,d,mmx(nx,nt),gs1(nx,nt),gn1(nx,nt),gs2(nx,nt),gn2(nx,nt),gs3(nx,nt),gn3(nx,nt),index1(nx,nt),index2(nx,nt),index3(nx,nt)
integer i,j,k
external smax
!!!读取500hPa的数据
open(20,file='G:\coldwave\blockingindex/500hgt.grd',form='binary')
do k = 1,nt
  do j = 1,ny
read(20)(h(i,j,k),i=1,nx)
  enddo
enddo
print*,h
do t = 1,nt
  do x = 1,nx
!!!计算北指数和南指数,三个扰动值都计算了
d    =   0
a    = 68 + d
b    = 60 + d
c    = 52 + d
gs1(x,t) = (h(x,b,t)-h(x,c,t))/20   
gn1(x,t) = (h(x,a,t)-h(x,b,t))/20
d    =   -2
a    = 68 + d
b    = 60 + d
c    = 52 + d
gs2(x,t) = (h(x,b,t)-h(x,c,t))/20   
gn2(x,t) = (h(x,a,t)-h(x,b,t))/20  
d    =   2
a    = 68 + d
b    = 60 + d
c    = 52 + d
gs3(x,t) = (h(x,b,t)-h(x,c,t))/20
gn3(x,t) = (h(x,a,t)-h(x,b,t))/20

!!将符合不同时满足两个条件的三个指数都取为0,这样处理缺测值可以设置为0
if(gs1(x,t) > 0 .and. gn1(x,t) < -10)then
index1(x,t) = gs1(x,t)
else
index1(x,t)  = 0
endif
if(gs2(x,t) > 0 .and. gn2(x,t) < -10)then
index2(x,t)  = gs2(x,t)
else
index2(x,t)  = 0
endif
if(gs3(x,t) > 0 .and. gn3(x,t) < -10)then
index3(x,t)   = gs3(x,t)
else
index3(x,t)  = 0
endif                           
!!!!!=====================输出三个指数最大的一个
mmx(x,t)      = smax(index1(x,t),index2(x,t),index3(x,t))

   enddo     !对应经度循环
enddo        !对应时间循环
!print*,mmx
open(22,file='G:\coldwave\blockingindex/2018.HGT.blocking.txt',status='replace')
open(44,file='G:\coldwave\blockingindex/2018.HGT.blocking.grd',status='replace',form='binary')
do t=1,nt
do x=1,nx
write(22,*)mmx(x,t)
write(44)mmx(x,t)  
enddo
enddo
close(22)
close(44)

end program

!!三个数,返回其中最大值的外部子程序
function smax(a,b,c)
real a,b,c
if(a > b)then
   if(a>c)then
    e = a
   else
    e = c
    endif
else
    if(b>c)then
e = b
else
e = c
endif
endif
smax = e
end function smax

评分

参与人数 1金钱 +10 收起 理由
Jerman5 + 10 很给力!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-11-10 12:11:20 | 显示全部楼层
C:\Users\Administrator\Desktop\1panel_ghgs_obs_nh_nrm.gif
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-11-10 12:13:46 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-11-10 18:02:41 | 显示全部楼层
谢谢楼主分享!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-12-23 10:26:30 | 显示全部楼层
非常好用,感谢分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-5-23 21:16:48 | 显示全部楼层
楼主,文献里北中南三个纬度选择的是40,60,80,delta选的-5,0,5,而你选的并不一样,想问问你这样改可以的吗,是参考的哪篇文献呢?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-5-23 22:00:38 | 显示全部楼层
littleqi 发表于 2019-5-23 21:16
楼主,文献里北中南三个纬度选择的是40,60,80,delta选的-5,0,5,而你选的并不一样,想问问你这样改可以的 ...

经纬度转换成格点看,68对应80N
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-5-23 22:10:53 | 显示全部楼层
董伟是我 发表于 2019-5-23 22:00
经纬度转换成格点看,68对应80N

哦哦明白了,多谢!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-6-29 14:06:03 | 显示全部楼层
谢谢分享!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2022-7-19 20:35:17 | 显示全部楼层
楼主你好,我想问一下这个阻高指数做出来之后怎么计算阻高的频次和强度呢?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表