爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 5995|回复: 4

[其他] NCL画图失败的问题

[复制链接]
发表于 2016-11-20 16:58:05 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 SG哭晕在厕所 于 2016-11-20 21:50 编辑

从官网上扒了一个例子自己试了试,但是一直显示fatal:(lat) is not a named dimension in variable (t).试着改了很多东西都没什么效果,被维度这个问题困扰挺久了,跪求大神帮忙解答
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin

;---Read in data
f     = addfile("wrfout_d01_2016-10-03_11:00:00","r")
t     = f->T2(0,:,:)                                    ; get t data

;---To plot multiple lines, you must put them into a mulidimensional array.
data      = new((/2,dimsizes(t&lat)/),float)

data(0,:) = u(0,:,{50})
data(1,:) = u(0,:,{51})

wks   = gsn_open_wks ("png","xy")                 ; send graphics to PNG file

;---Set plotting parameters

res                   = True                      ; plot mods desired
res@tiMainString      = "Two curve XY plot"       ; add title
;
; Similiar resources are xyLineThicknessF and xyLineColor,
; which will affect all lines in the array.
;
res@xyLineThicknesses = (/  1.0,   2.0/)          ; make second line thicker
res@xyLineColors      = (/"blue","red"/)          ; change line color

plot  = gsn_csm_xy (wks,t&lat,data,res) ; create plot

end
密码修改失败请联系微信:mofangbao
发表于 2016-11-20 22:26:19 | 显示全部楼层
wrfout的经纬度的变量名应该是XLONG和XLAT吧,你好好查查里面的变量。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-11-20 22:43:03 | 显示全部楼层
zly4814624 发表于 2016-11-20 22:26
wrfout的经纬度的变量名应该是XLONG和XLAT吧,你好好查查里面的变量。

我试过了,结果只是错误提示里的(lat)变成了(XLAT),还有south_north,degree_north什么的,都是一样的错误,搞不懂经纬变量名到底是哪个
密码修改失败请联系微信:mofangbao
发表于 2016-11-21 09:48:45 | 显示全部楼层
文件里面是不是直接就有维度的变量,你读一下呢
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-11-21 15:23:06 | 显示全部楼层
海盗船长 发表于 2016-11-21 09:48
文件里面是不是直接就有维度的变量,你读一下呢

这个物理量的信息是这样的
Variable: t
Type: float
Total Size: 18172 bytes
            4543 values
Number of Dimensions: 3
Dimensions and sizes:   [Time | 1] x [south_north | 59] x [west_east | 77]
Coordinates:
Number Of Attributes: 6
  FieldType :   104
  MemoryOrder : XY
  description : TEMP at 2 M
  units :       K
  stagger :
  coordinates : XLONG XLAT
到底应该怎么识别经度纬度我已经混乱了。。。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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