爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 17734|回复: 3

[脚本编辑] 求助:对高度场去除趋势后求WPSH强度指数

[复制链接]

新浪微博达人勋

发表于 2021-3-13 20:30:10 | 显示全部楼层 |阅读模式

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

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

x
在全球变暖的大背景下,西太副高增强,国家气候中心给的西太副高强度指数突变明显。我想对高度场去除趋势之后再自己算出夏季强度指数。希望有大神指点,毕业论文好难。。。。回帖或者加qq:952021247   附件里是1981-2020年全球500hPa月平均高度场数据。

1981_2020hgt.mon.mean.grd

19.38 MB, 下载次数: 2, 下载积分: 金钱 -5

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

新浪微博达人勋

发表于 2022-4-13 22:32:52 | 显示全部楼层
我写过用王斌老师(2013)定义计算的西太副高指数。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-13 22:35:00 | 显示全部楼层
clc,clear,clear off

ncdisp('F:\Data\hgt.mon.mean.nc')
hgt=ncread('F:\Data\hgt.mon.mean.nc','hgt');
  time=ncread('F:\Data\hgt.mon.mean.nc','time');
  lat=ncread('F:\Data\hgt.mon.mean.nc','lat');
  lon=ncread('F:\Data\hgt.mon.mean.nc','lon');
format = 'mm dd, yyyy HH:MM:SS.FFF AM';
time = double(time);
dstr = datestr((datenum('1800-01-01') + time./24),format);
TM = datevec(dstr);
hgt=squeeze(hgt(:,:,3,385:884));%数据从1980年1月开始截取;截取至2020年12月 ;
hgt(:,:,501:504)=nan;
%%
  hgt_year=reshape(hgt,[size(hgt,1),size(hgt,2),12,floor(size(hgt,3)/12)]);
   hgt_mean=nanmean(hgt_year,4);
   hgt_ano=hgt_year-hgt_mean;
%  截取数据
  lon_need1=find(lon>=115&lon<=150);
   lat_need1=find(lat>=15&lat<=25);
index=hgt_ano(lon_need1,lat_need1,6:8,:);
%%
%求出指数
WPSH_Osindex=squeeze(nanmean(nanmean(nanmean(index,1),2),3));
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-6-3 10:52:52 | 显示全部楼层
请问一下楼主在计算强度指数时用的纬向格距数值dx和经向格距数值dy是多大
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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