ncl中怎么叠加shp文件（未解决）

小莹子

利用WRFDomain画出青藏高原的范围，在这里选择WRFDomain的原因是自己还不会设置格点数，所以一般就是通过可视化确定单层或者多层嵌套，然后根据生成的namelist.wps修改namelist.input，生成wrfout。
我想添加地形，然后在气象家园中下载的青藏高原的shp文件，然后查了说使用gsn_add_shapefile_polylines。
在ncl官网下载gsn_add_shapefile_polylines脚本，然后根据文献中给出的分析范围25-40度N，75-105度E。然后修改filename = "/home/syrus/data/WRFDomain/shp/Qinghai-Tibet_Plateau.shp"
但是只出现了范围的空图，没有实质内容。在摸索

就一直在纠结是自己代码是不是写错的问题，后来就在别人的指导意见下，重新下载了青藏高原的shp文件，显示成功，那可能就是第一次下载的文件有问题。
但是precip的图和shp图是独立的，没有叠加在一起，还没有解决

filename="/home/syrus/data/Tibet/tibetan/tibetan.shp"
  pres             = True
  pres@gsLineColor = "blue"
  pres@gsLineThicknessF =2
poly=gsn_add_shapefile_polylines(wks,plot,filename,pres)

draw(plot)
frame(wks);   Example script to produce plots for a WRF real-data run,


http://www.ncl.ucar.edu/Applications/wrf.shtml

sfz13851581045

顶一下，求解决！

jiangjinyue

请问能上传一下青藏高原的shp文件吗，谢谢

小莹子

jiangjinyue 发表于 2019-5-14 22:00
请问能上传一下青藏高原的shp文件吗，谢谢

已上传，这个也是我在家园的帖子里下载的，建议你可以多找找。

深度狙杀

楼主最后怎么解决的？

猎人

雪中送碳，非常感谢！