关于计算 webster and yang的印度夏季风指数

LlovHY

才入门的新手，求教！最近一直被webster and yang定义印度夏季风指数计算结果所困扰，选取的是0——20N，40E——110E范围的850hPa与200hPa的纬向风差：U850-U200，可是利用ERA40和JRA55计算后得出的季风指数与最终的正确结果相关性总是很差。觉得问题应该出在了程序里可是又看不出错在哪里，所以求各位做过这方面的大神指点一下。此程序用的是1958年到2001年的ERA40的风场nc资料。


begin
   ymStrt  = 195806                           
   ymLast  = 200109                          
   yrStrt  = ymStrt/100
   yrLast  = ymLast/100
  latas   = 0
  latan   = 20
  lonaw   = 40
  lonae   = 110
f = addfile("wind.nc", "r")
  lev=f->level
  print(f)
  number=(/0,1/)   ;200,850hPa
  TIME   = f->time
  YYYY   = ut_calendar(TIME,-1)/100            
  iYYYY  = ind(YYYY.ge.yrStrt .and. YYYY.le.yrLast)
  uaa  =short2flt(f->u(iYYYY,number,{latas:latan},{lonaw:lonae}))
  printVarSummary(uaa)
uajjay=month_to_season (uaa, "JJA")
uajja1=uajjay(:,0,:,:)
uajja2=uajjay(:,1,:,:)
lata    = f->latitude({0:20})  
rad    = 4.0*atan(1.0)/180.0      
wgtya   = cos(lata*rad)
wgtx   =1.0   
uaave1 =wgt_areaave(uajja1,wgtya,wgtx,0)
uaave2 =wgt_areaave(uajja2,wgtya,wgtx,0)
index=uaave2-uaave1

sd_index= dim_standardize_n_Wrap(index,1,0)
print(sd_index)
asciiwrite("2002.ascii",sd_index)
asciiwrite("2002yuan.ascii",index)
fbindirwrite("2002.dat",sd_index)
fbindirwrite("2002yuan.dat",index)

lucinda筱橙

{:eb302:}{:eb302:}

lleoiu

你的数据从1958.6开始。下面算的是JJA，
uajjay=month_to_season (uaa, "JJA")
你仔细查一下数据，看看uajjay是不是JJA的平均。
JJA换成JFM试试

Mye011

lleoiu 发表于 2016-8-16 10:52
你的数据从1958.6开始。下面算的是JJA，
uajjay=month_to_season (uaa, "JJA")
你仔细查一下数据，看看u ...

我也有这样的问题，假如我的数据从198603-201602，我用了DJF的month_to_season函数，我需要2015年冬季的，那么输出aDJF(29,:,:），对吗

Mye011

lleoiu 发表于 2016-8-16 10:52
你的数据从1958.6开始。下面算的是JJA，
uajjay=month_to_season (uaa, "JJA")
你仔细查一下数据，看看u ...

从6月开始的数据都要减去6个月（JJA变成JFM），7月开始的数据都要-7来算，这个意思吗

yzzmc

学习了，谢谢

subtropical

楼主，JRA55资料可以分享下吗？谢谢了。

DABEIGE

求助一下楼主，我没看懂rad    = 4.0*atan(1.0)/180.0  是什么意思，我也去看了官网的例子。第一个例子里的rad    = 4.0*atan(1.0)/180.0，dlon   = abs(lon(2)-lon(1))*rr
  dx     = dlon*cos(lat*rad)
;                                    
  dy     = new ( jlat, typeof(dx))
我都没太看懂

帅蛋fly

楼主你好，想请问“最终的正确结果”是从哪里来的？