爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5272|回复: 1

[作图] ncl画风玫瑰图

[复制链接]

新浪微博达人勋

发表于 2018-4-17 15:14:51 | 显示全部楼层 |阅读模式

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

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

x
用站点数据画风玫瑰图,老是报错,我把程序和报错信息贴出来,请大侠们帮忙看看是程序哪里要改?
程序:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
begin
f1="E:\cygwin\home\Sta288\2017wind.csv"
wrData = asciiread(f1,-1,"string")
  delim = ";"
  wspd = str_get_field(wrData,0,delim)
  wdir = str_get_field(wrData,1,delim)                             
  wspd@long_name = "Wind Speed"
  wspd@units     = "m/s"
  wdir@long_name = "Wind Direction"

; specify plot arguments
  numPetals =  8                    ; N, NE, E, SE, S, SW, W, NW
  circFr    = 10.
  spdBounds         = (/ 10., 20., 30., 100. /)
; generate color plot
  wks = gsn_open_wks("png","rose")  ; send graphics to PNG file
  res               = True
  res@tiMainString  = "Wind Rose: Color + Variable Thickness"
  colorBounds       = (/ "blue", "green", "yellow", "red" /)
  wrColor           = WindRoseColor (wks,wspd,wdir,numPetals,circFr,spdBounds,colorBounds,res)
  end

报错信息截图:
QQ截图20180417150838.png


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

新浪微博达人勋

发表于 2018-5-3 18:15:49 | 显示全部楼层
是不是 需要把wdir 和 wspd 赋值时 加上tofloat 也就是wspd = tofloat(str_get_field(wrData,0,delim))
  wdir =tofloat( str_get_field(wrData,1,delim) )  就是需要让wspd 和wdir类型变为浮点型  。我也不是很懂啦 你试试
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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