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

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 64|回复: 2

[作图] wrfout画图问题

[复制链接] |关注本帖

新浪微博达人勋

1现金
想用ncl画wrfout的降水量图,得一个3*3的区域平均的降水量的值,然后逐小时输出一个值来做一个折线图
但就是在读取这个3*3的区域的时候没法读取
下面是我的程序
begin
  a  = addfile("../wrfout_d02_2015-08-18_00:00:00","r")
  rainc   =  wrf_user_getvar(a,"RAINC",-1)
  rainc@lat2d = wrf_user_getvar(a,"XLAT",0)   ; latitude/longitude
  rainc@lon2d = wrf_user_getvar(a,"XLONG",0)  ; required for plotting  
  rainnc  = wrf_user_getvar(a,"RAINNC",-1)
  rainnc@lat2d = wrf_user_getvar(a,"XLAT",0)   ; latitude/longitude
  rainnc@lon2d = wrf_user_getvar(a,"XLONG",0)  ; required for plotting      

  j = 0
  X = new(24, float)
  do i = 23,47
   raincj  = rainc(i+1,{28:31},{93:96})-rainc(i,{28:31},{93:96})
   rainncj  = rainnc(i+1,{28:31},{93:96}) -rainnc(i,{28:31},{93:96})
   rainTot = raincj + rainncj
   n = avg(rainTot)
   X(j) = n
   j=j+1
  end do
  res = True
  wks = gsn_open_wks("png","zxt")
  plot = gsn_csm_xy(wks,(/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24/),X,res)
end


这是报的错
fatal:Dimension (west_east) of (rainc) does not have an associated coordinate variable

密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-10-12 18:09:26 | 显示全部楼层 |取消关注该作者的回复
本帖最后由 男紫汗 于 2018-10-12 19:37 编辑

。。。。。。。。。。
密码修改失败请联系qq:937062711
回复

使用道具 举报

新浪微博达人勋

发表于 2018-10-12 22:23:10 | 显示全部楼层 |取消关注该作者的回复
男紫汗 发表于 2018-10-12 18:09
。。。。。。。。。。

raincj  = rainc(i+1,{28:31},{93:96})-rainc(i,{28:31},{93:96})
把大括号去掉
raincj  = rainc(i+1,28:31,93:96)-rainc(i,28:31,93:96)
密码修改失败请联系qq:937062711
回复

使用道具 举报

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

本版积分规则

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

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

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