NCL绘制色斑图-多要素叠加

需求：下载GFS数据，绘制不同高度层（500hPa~950hPa）风温湿度；

GFS数据下载地址：https://www.cnblogs.com/tiandi/p/11045602.html

下载下来样例数据格式如：gfs.t06z.pgrb2.0p25.f000.grb2

代码如下：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"begin
;---Open file and read datafilename = "/cygdrive/d/learn/tonghang/gfs.t06z.pgrb2.0p25.f000.grb2"a    = addfile(filename,"r")lon = a->lon_0({100:125})lat = a->lat_0({15:40})lv=a->lv_ISBL0HGT= a->HGT_P0_L100_GLL0(:,{15:40},{100:125}) u = a->UGRD_P0_L100_GLL0(:,{15:40},{100:125})v = a->VGRD_P0_L100_GLL0(:,{15:40},{100:125})rh = a->RH_P0_L100_GLL0(:,{15:40},{100:125});---Start the graphics
do i = 0,11wks = gsn_open_wks("png","/cygdrive/d/learn/tonghang/975-500wind-height-rh/wind&height&rh"+lv(i)/100)gsn_define_colormap(wks,"MPL_Blues")
;---Set some resourcesres                   = Trueres@gsnDraw          = Falseres@gsnFrame         = Falseres@gsnMaximize       = True     ; maximize plot in frameres@cnFillOn          = True     ; turn on color fillres@gsnAddCyclic      = Falseres@mpOutlineOn           = Trueres@mpDataBaseVersion     = "MediumRes"res@mpDataSetName         = "Earth..4"  res@mpOutlineSpecifiers   = "China:states"res@mpOutlineBoundarySets = "AllBoundaries"res@mpFillOn =Falseres@mpMinLatF             =  15         ; min latres@mpMaxLatF             =  40         ; max latres@mpMinLonF             =  100        ; min lonres@mpMaxLonF             =  125        ; max lonres@mpGeophysicalLineThicknessF= 2.0 res@mpNationalLineThicknessF   = 2.0res@gsnLeftString              = lv(i)/100+"hPa  Wind&Height&RH"res@gsnRightString             = "2020/02/29/06(UTC)"res@gsnRightStringFontHeightF  = .02res@gsnLeftStringFontHeightF   = .02res@lbOrientation              = "vertical"res@lbTitleOn                  = True                  ; turn on titleres@lbTitleString              = "%"res@lbTitleFontHeightF         = .015                 ; make title smallerres@pmTickMarkDisplayMode      = "Always"res@tmXBTickSpacingF           = 40res@tiXAxisString              = "NCEP GFS 0.25~S~o~N~*0.25~S~o~N~       VALID TIME: 2020/02/29/06 +000hrs" res@tiXAxisFontHeightF         = .015res@cnLevelSelectionMode = "ExplicitLevels"          ;绘制等值线res@cnLevels             = (/80,90,100/)res@cnFillColors         = (/0,40,80,110,128/)plot0 = gsn_csm_contour_map(wks,rh(i,:,:),res)res1=Trueres1@gsnDraw              = False          res1@gsnFrame             = False res1@gsnAddCyclic      = False          res1@cnFillOn          = Falseres1@cnLinesOn         = Trueres1@cnLineThicknessF = 2.5res1@cnLineColor = "blue"res1@cnInfoLabelOn=Falseres1@gsnLeftString=""res1@gsnRightString=""plot1 = gsn_csm_contour(wks,HGT(i,:,:),res1)res2                         = True               ; plot mods desiredres2@gsnDraw                 = False          res2@gsnFrame                = False res2@gsnAddCyclic      = False;res2@vcRefMagnitudeF         = 10.                ; make vectors larger;res2@vcRefLengthF            = 0.050              ; ref vec lengthres2@vcGlyphStyle            = "WindBarb"         ; select wind barbs res2@vcMinDistanceF          = 0.025              ; thin out windbarbsres2@vcRefMagnitudeF         = 4.res2@vcRefLengthF            = 0.020res2@vcWindBarbScaleFactorF   = 2.5res2@vcWindBarbLineThicknessF  =2.res2@gsnLeftString=""res2@gsnRightString=""res2@vcRefAnnoString2On = Falseres2@lbLabelBarOn = Falseres2@vcRefAnnoOn = Falseplot2=gsn_csm_vector(wks,u(i,:,:),v(i,:,:),res2)overlay(plot0,plot2)overlay(plot0,plot1)draw(plot0)frame(wks)end do
end

效果如下：（就放三张图吧）

NCL绘制色斑图-多要素叠加相关推荐

使用NCL绘制安徽省的轮廓图
使用NCL绘制安徽省的轮廓图,其它省可以参照设置: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load &q ...
Python绘制极化图时坐标轴的旋转及两散点图的叠加绘制
Python绘制极化图时坐标轴的旋转及两散点图的叠加绘制我的sen_angle是n*2的数组,第一列是卫星方位角,第二列是卫星天顶角,sun_angle第一列是太阳方位角,第二列是太阳天顶角. im ...
python简单代码画曲线图教程-Python绘制折线图和散点图的详细方法介绍（代码示例）...
本篇文章给大家带来的内容是关于Python绘制折线图和散点图的详细方法介绍(代码示例),有一定的参考价值,有需要的朋友可以参考一下,希望对你有所帮助. 1.绘制折线图和散点图要用到matplotlib ...
python画折线图代码-Python绘制折线图和散点图的详细方法介绍（代码示例）
本篇文章给大家带来的内容是关于Python绘制折线图和散点图的详细方法介绍(代码示例),有一定的参考价值,有需要的朋友可以参考一下,希望对你有所帮助. 1.绘制折线图和散点图要用到matplotlib ...
3、绘制E-R图：数据库概要设计阶段
需求分析阶段过后,就进入了概要设计阶段. 在概要设计阶段,我们需要和项目团队中的其他成员及客户进行沟通,讨论数据库的设计是否满足客户的业务和数据处理需求. 数据库设计和机械.建筑等行业一样,机械行业需 ...
tableau绘制雷达图（4步法）
为什么要用tableau? 我称之为(静态.动态图的区别) 我们可以用python,excel进行雷达图的绘制,那么为什么要用tableau? 因为,tableau的图可以是动态的多个视图的整合,你可 ...
利用Python的folium包绘制城市道路图！专业地图！
写在前面很长一段时间内,我都在研究在线地图的开发者文档,百度地图和高德地图的开发者中心提供了丰富的在线地图服务,虽然有一定的权限限制,但不得不说,还是给我的科研工作提供了特别方便的工具,在博客前面我 ...
matlab怎么画两个自变量的图_关系图怎么画？一款实用的绘制关系图设计软件
关系图是指实体-联系图,是用来描述现实世界的概念模型.关系图应用范围很广.比如人物关系图.零件关系图等等.关系图主要由三部分构成:矩形框里写实体名.椭圆符号表示属性.菱形框中标明何种联系.同时还用线条 ...
R-forestplot包| HR结果绘制森林图
本文首发于"生信补给站"微信公众号,https://mp.weixin.qq.com/s/2W1W-8JKTM4S4nml3VF51w 更多关于R语言,ggplot2绘图,生信分析 ...
R统计绘图-corrplot绘制热图及颜色、字体等细节修改1
有师妹想要更改热图的颜色和字体,想着之前相关性绘图等推文只是使用corrplot默认的颜色绘图,为了帮师妹解惑,今天就写一篇,怎么设置热图颜色和字体等细节到推文.其实看一遍R语言实战|入门3:图形初阶 ...

NCL绘制色斑图-多要素叠加

NCL绘制色斑图-多要素叠加相关推荐

最新文章

热门文章