爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 21803|回复: 13

[其他] NCL坐标问题(解决)

[复制链接]
发表于 2014-2-20 09:51:57 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 BDRUSH 于 2014-2-22 00:24 编辑

大家看看我这个图,我有2个问题,第一个问题是里面的网格能去掉么,第二个问题是横纵坐标为啥经纬度不显示出来,需要设置什么吗? QQ截图20140220145354.jpg
下面是脚本。
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/csm/shea_util.ncl"

begin

  a1 = addfile("/cygdrive/D/Study/NCL/home/data/uwnd.1961.nc","r")
  a2 = addfile("/cygdrive/D/Study/NCL/home/data/vwnd.1961.nc","r")  

  u = short2flt(a1->uwnd(989,5,:,:))   
  v = short2flt(a2->vwnd(989,5,:,:))
  lon=a1->lon
  lat=a1->lat
wks=gsn_open_wks("pdf","1961290-1961300")

res=True

res@vcMinFracLengthF=0.25     
res@vcRefMagnitudeF=18.0     
res@vcRefLengthF=0.02  
res@vcMonoLineArrowColor=False
res@vcMinDistanceF             =0.015

res@vfXCStartV=lon(0)
res@vfXCEndV=lon(dimsizes(lon)-1)
res@vfYCStartV=lat(0)
res@vfYCEndV=lat(dimsizes(lat)-1)

res@mpLimitMode="LatLon"
res@mpCenterLonF=-180
res@mpMinLonF=-240
res@mpMaxLonF=-120
res@mpMinLatF=0
res@mpMaxLatF=60
res@mpPerimOn=True

res@tiMainString=":F26:Wind Feild---1961"
res@tiXAxisString=":F25:longitude"
res@tiYAxisString=":F25:latitude"
res@lbLabelFont="times-roman"

res@pmLabelBarDisplayMode="Always"
res@lbPerimOn=False
res@pmLabelBarSide="Bottom"
res@lbOrientation="Horizontal"

res@gsnScale=True

res@vpXF=0.2
res@vpYF=0.8
res@vpWidthF=0.6
res@vpHeightF=0.4

  cmap=(/(/1.00,1.00,1.00/),(/0.00,0.00,0.00/),(/.560,.500,.700/),(/.300,.300,.700/),(/.100,.100,.700/),(/.000,.100,.700/),(/.000,0.300,.700/),(/.000,0.500,.500/),(/.000,.700,.100/),(/.060,.680,.000/),(/.550,0.550,.000/),(/.570,0.420,.000/),(/.700,0.285,.000/),(/.700,0.180,.000/),(/.870,0.050,.000/),(/1.00,0.000,.000/),(/.700,0.700,.700/)/)
  gsn_define_colormap(wks,cmap)

  map=gsn_vector_map(wks,u,v,res)

end

1961290-1961300.pdf

248.86 KB, 下载次数: 4, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao
发表于 2014-2-20 11:01:43 | 显示全部楼层
楼主直接把图贴出来吧
密码修改失败请联系微信:mofangbao
发表于 2014-2-20 12:58:21 | 显示全部楼层
下载图还得花金币
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-2-20 14:56:50 | 显示全部楼层
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-2-20 14:57:28 | 显示全部楼层
尽头的尽头 发表于 2014-2-20 11:01
楼主直接把图贴出来吧

好了,谢谢提醒,之前忘了。
密码修改失败请联系微信:mofangbao
发表于 2014-2-20 15:35:30 | 显示全部楼层
网格线一般默认是关闭的,你再把脚本贴出来,让大家帮你看看
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-2-20 21:52:42 | 显示全部楼层
自己顶,丰衣足食,希望大家捧个人场
密码修改失败请联系微信:mofangbao
发表于 2014-2-20 22:01:01 | 显示全部楼层
{:5_231:}
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-2-20 22:14:06 | 显示全部楼层

你还真的就只捧捧人场啊,太客气了
密码修改失败请联系微信:mofangbao
发表于 2014-2-21 13:58:11 | 显示全部楼层
map=gsn_vector_map(wks,u,v,res)这句改为
map=gsn_csm_vector_map(wks,u,v,res)
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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