NCL代码分享-读取FNL资料画天气图

liuxiaoyue · 发表于 2017-11-17 22:51:15

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

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

x

本帖最后由 liuxiaoyue 于 2017-11-17 22:57 编辑

这是效果图~（槽线是用Photoshop手动添加的）
效果图.png

附上程序

500hPa-share.ncl (14.86 KB, 下载次数: 666)

2017-11-17 22:49 上传

点击文件名下载附件

这个是坐标轴刻度设定的一个函数，代码的一部分

undef("add_map_tickmarks2")
function add_map_tickmarks2(wks,plot,res)
local res2, bres, vpx, vpy, vpw, vph, xndc, yndc, npts, n, j, nlat, \
nlon, delta, bot_lon, top_lon, lft_lat, rgt_lat, xblabels, xbvalues, \
xtlabels, xtvalues, yllabels, ylvalues, yrlabels, yrvalues, xfix, \
xlat, xlon, yfix, annoid, anno_str
begin
;---Make a copy of the original resource list.
res2 = res
;---Retrieve edges of plot in NDC space.
getvalues plot
"vpXF" : vpx
"vpYF" : vpy
"vpWidthF" : vpw
"vpHeightF" : vph
end getvalues
;---Turn off tickmarks associated with map. We want to add our own.
setvalues plot
"pmTickMarkDisplayMode" : "Never"
end setvalues
;---Initialize resources for tickmark plot. User shouldn't change these.
bres = True
bres@vpXF = vpx
bres@vpYF = vpy
bres@vpWidthF = vpw
bres@vpHeightF = vph
bres@trXMinF = vpx
bres@trXMaxF = vpx + vpw
bres@trYMinF = vpy - vph
bres@trYMaxF = vpy
bres@tmEqualizeXYSizes = True
;---This resource the user can change in main code if desired.
bres@gsnTickMarksPointOutward = get_res_value(res2,"gsnTickMarksPointOutward",True)
;
; NDC Points to scan on X and Y axes. These arrays will be used to
; find the closest NDC pair that gets us close to the location where
; we want a tickmark.
;
npts = 100000 ; Increase to get closer match for tickmarks
xndc = fspan(vpx,vpx+vpw,npts)
yndc = fspan(vpy-vph,vpy,npts)
n = dimsizes(yndc)
xfix = new(n,float)
yfix = new(n,float)
xlon = new(n,float)
xlat = new(n,float)
delta = 0.001
;---Left axis tickmarks
if(isatt(res2,"tmYLValues")) then
lft_lat = get_res_value(res2,"tmYLValues",-1)
nlat = dimsizes(lft_lat)
ylvalues = new(nlat,float)
yllabels = new(nlat,string)
xfix = vpx + 0.0001 ; Just a smidge into the plot to make sure we don't
; get missing values returned.
;
; Loop across each left latitude value that we want a tickmark for,
; and try to find the closest X,Y NDC coordinate pair along this axis.
;
NhlNDCToData(plot,xfix,yndc,xlon,xlat)
do j=0,dimsizes(lft_lat)-1
NhlNDCToData(plot,xfix,yndc,xlon,xlat)
ii = minind(fabs(xlat-lft_lat(j)))
if(.not.any(ismissing(ii)).and.fabs(xlat(ii)-lft_lat(j)).le.delta)
yllabels(j) = fabs(lft_lat(j)) + ""
ylvalues(j) = yndc(ii(0))
if(lft_lat(j).lt.0) then
yllabels(j) = yllabels(j) + "~S~o~N~S"
end if
if(lft_lat(j).gt.0) then
yllabels(j) = yllabels(j) + "~S~o~N~N"
end if
end if
delete(ii)
end do
bres@tmYLMode = "Explicit"
bres@tmYLValues = ylvalues
bres@tmYLLabels = get_res_value(res2,"tmYLLabels",yllabels)
else
bres@tmYLOn = False
bres@tmYLLabelsOn = False
end if
;---Right axis tickmarks
if(isatt(res2,"tmYRValues")) then
rgt_lat = get_res_value(res2,"tmYRValues",-1)
nlat = dimsizes(rgt_lat)
yrvalues = new(nlat,float)
yrlabels = new(nlat,string)
xfix = vpx + vpw - 0.0001 ; Just a smidge into the plot to make sure we don't
; get missing values returned.
;
; Loop across each right latitude value that we want a tickmark for,
; and try to find the closest X,Y NDC coordinate pair along this axis.
;
do j=0,dimsizes(rgt_lat)-1
NhlNDCToData(plot,xfix,yndc,xlon,xlat)
ii = minind(fabs(xlat-rgt_lat(j)))
if(.not.any(ismissing(ii)).and.fabs(xlat(ii)-rgt_lat(j)).le.delta)
yrlabels(j) = fabs(rgt_lat(j)) + ""
yrvalues(j) = yndc(ii(0))
if(rgt_lat(j).lt.0) then
yrlabels(j) = yrlabels(j) + "~S~o~N~S"
end if
if(rgt_lat(j).gt.0) then
yrlabels(j) = yrlabels(j) + "~S~o~N~N"
end if
end if
delete(ii)
end do
bres@tmYROn = True
bres@tmYRLabelsOn = True
bres@tmYUseLeft = False
bres@tmYRMode = "Explicit"
bres@tmYRValues = yrvalues
bres@tmYRLabels = get_res_value(res2,"tmYRLabels",yrlabels)
else
bres@tmYUseLeft = False
bres@tmYROn = False
bres@tmYRLabelsOn = False
end if
;---Top axis tickmarks
if(isatt(res2,"tmXTValues")) then
top_lon = get_res_value(res2,"tmXTValues",-1)
nlon = dimsizes(top_lon)
xtvalues = new(nlon,float)
xtlabels = new(nlon,string)
yfix = vpy - 0.0001 ; Just a smidge into the plot to make sure we don't
; get missing values returned.
;
; Loop across each top longitude value that we want a tickmark for,
; and try to find the closest X,Y NDC coordinate pair along this axis.
;
do j=0,dimsizes(top_lon)-1
NhlNDCToData(plot,xndc,yfix,xlon,xlat)
ii = minind(fabs(xlon-top_lon(j)))
if(.not.any(ismissing(ii)).and.fabs(xlon(ii)-top_lon(j)).le.delta)
xtlabels(j) = fabs(top_lon(j)) + ""
xtvalues(j) = xndc(ii(0))
if(top_lon(j).lt.0) then
xtlabels(j) = xtlabels(j) + "~S~o~N~W"
end if
if(top_lon(j).gt.0) then
xtlabels(j) = xtlabels(j) + "~S~o~N~E"
end if
end if
delete(ii)
end do
bres@tmXTOn = True
bres@tmXTLabelsOn = True
bres@tmXUseBottom = False
bres@tmXTMode = "Explicit"
bres@tmXTValues = xtvalues
bres@tmXTLabels = get_res_value(res2,"tmXTLabels",xtlabels)
else
bres@tmXUseBottom = False
bres@tmXTOn = False
bres@tmXTLabelsOn = False
end if
;---Bottom axis tickmarks
if(isatt(res2,"tmXBValues")) then
bot_lon = get_res_value(res2,"tmXBValues",-1)
nlon = dimsizes(bot_lon)
xbvalues = new(nlon,float)
xblabels = new(nlon,string)
yfix = vpy-vph + 0.0001 ; Just a smidge into the plot to make sure
; we don't get missing values returned.
;
; Loop across each bottom longitude value that we want a tickmark for,
; and try to find the closest X,Y NDC coordinate pair along this axis.
;
do j=0,dimsizes(bot_lon)-1
NhlNDCToData(plot,xndc,yfix,xlon,xlat)
ii = minind(fabs(xlon-bot_lon(j)))
if(.not.any(ismissing(ii)).and.fabs(xlon(ii)-bot_lon(j)).le.delta)
xblabels(j) = fabs(bot_lon(j)) + ""
xbvalues(j) = xndc(ii(0))
if(bot_lon(j).lt.0) then
xblabels(j) = xblabels(j) + "~S~o~N~W"
end if
if(bot_lon(j).gt.0) then
xblabels(j) = xblabels(j) + "~S~o~N~E"
end if
end if
delete(ii)
end do
bres@tmXBMode = "Explicit"
bres@tmXBValues = xbvalues
bres@tmXBLabels = get_res_value(res2,"tmXBLabels",xblabels)
else
bres@tmXBOn = False
bres@tmXBLabelsOn = False
end if
;
; Now that we are done figuring out where to put tickmarks, and
; what labels to use, get any "tm" resources that might have been
; set by the user, and create a blank plot with thes new tickmarks.
;
;---Get rest of user resources that were set with "tm".
bres = get_res_eq(res2,"tm")
bres = True ; Above call will set bres to True if no "tm" resources, so
; make sure it is True still.
bres@gsnDraw = False
bres@gsnFrame = False
;
; Create blank plot with new tickmarks (don't use gsn_csm_blank_plot,
; because it wants to scale the size of your X and Y axes.)
;
blank = gsn_blank_plot(wks,bres)
;
; Attach new tickmarks to original plot. This will allow resizing
; if desired. The default is to attach one plot to the center of
; the other one. These two plots are already the same size.
;
annoid = gsn_add_annotation(plot,blank,False)
;
; Be sure to return the annotation id, otherwise the
; tickmarks will disappear.
;
anno_str = unique_string("annoid")
plot@$anno_str$ = annoid
return(plot)
end

复制代码

本人也是初学NCL，有什么不对的地方希望大家多多指正！

小其其格 · 发表于 2017-11-17 23:38:50

很棒，下下来学习学习！

zongqing123 · 发表于 2017-11-18 09:47:03

小其其格发表于 2017-11-17 23:38
很棒，下下来学习学习！

很棒，很想学习一下如何写函数

river · 发表于 2017-11-18 13:49:16

很漂亮，感谢分享

WZZ809525579 · 发表于 2017-11-18 14:06:22

挺不错，谢谢分享

rabin_xu · 发表于 2017-11-21 10:36:58

这图画的挺漂亮的，学习一下

成吉思汗__杜 · 发表于 2017-11-21 11:17:50

{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}

贫道敬孔 · 发表于 2017-11-24 17:58:21

多谢楼主分享

{:5_238:}

sunjyi · 发表于 2017-11-26 15:32:51

{:5_213:}{:5_213:}

liuxiaoyue · 发表于 2017-12-7 08:58:01

转角转好运发表于 2017-12-6 18:04
会写函数的牛逼楼主，求带，求飞.

呀这个函数不是我写的，官网上提供的，我只会用不会写呢

		自动登录	找回密码
密码			立即注册

[作图] NCL代码分享-读取FNL资料画天气图

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

点评

评分

浏览过的版块