加载中…
个人资料
les-baleines
les-baleines
  • 博客等级:
  • 博客积分:0
  • 博客访问:4,683
  • 关注人气:2
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
相关博文
推荐博文
正文 字体大小:

NCL读取txt格式3列(经度,纬度,亮温)画图

(2018-04-05 16:02:44)
分类: NCL学习笔记
;***********************************
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"

begin
ascii_filename = "/nuist/p/workncl/yunyu.txt"
argu  = asciiread(ascii_filename,(/1324,3/),"float")
;1324行

obs = argu(:,2)
lat=argu(:,1)
lon=argu(:,0)
olon = fspan(120,146,256)
olat = fspan(10,30,153)

olon!0          = "lon"
olon@long_name  = "lon"
olon@units      = "degrees-east"
olon&lon        = olon
olat!0          = "lat"
olat@long_name  = "lat"
olat@units      = "degrees_north"
olat&lat        = olat
obs@_FillValue = 999999.000000
rscan = (/1/)   
dataobs = obj_anal_ic_Wrap(lon,lat,obs,olon,olat,rscan,False)
dataobs@units="K"

wks = gsn_open_wks ("png","obs")  
gsn_define_colormap( wks ,"Rainbow")         
res                         = True            
res@gsnAddCyclic            = False
res@tiMainString            ="obs" 
res@gsnMaximize             = True
;res@gsnDraw                 = False
;res@gsnFrame                = False

;>--------------------------------------------<</span>
          set for the map 
;>--------------------------------------------<</span>
res@mpMinLatF               = 10.                         
res@mpMaxLatF               = 30.
res@mpMinLonF               = 120.
res@mpMaxLonF               = 146.

res@mpFillOn                = True
res@mpOutlineOn             = True  ; Use outlines from shapefile
res@mpDataBaseVersion       = "MediumRes"
res@mpDataSetName           = "Earth..4"
res@mpOutlineSpecifiers     = (/"China:states","Taiwan"/)
;>--------------------------------------------<</span>
; set for the plot

res@cnFillOn                = True               
res@cnLinesOn               = False             
;res@cnLevelSpacingF         = 1.             
res@gsnSpreadColors         = True         
res@lbLabelAutoStride       = True
res@cnSmoothingOn           = True               
plot=gsn_csm_contour_map(wks,dataobs,res) 

end

0

阅读 评论 收藏 转载 喜欢 打印举报/Report
  • 评论加载中,请稍候...
发评论

    发评论

    以上网友发言只代表其个人观点,不代表新浪网的观点或立场。

      

    新浪BLOG意见反馈留言板 电话:4000520066 提示音后按1键(按当地市话标准计费) 欢迎批评指正

    新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 会员注册 | 产品答疑

    新浪公司 版权所有