一个画地形高度的脚本

willnot

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/wrf/WRFUserARW.ncl"
begin
;
; The WRF ARW input file.  
; This needs to have a ".nc" appended, so just do it.
  a = addfile("geo_em.d01.nc","r")
; We generate plots, but what kind do we prefer?
type = "ps"
  wks = gsn_open_wks(type,"ter")

; Set some basic resources
  res = True
  res@MainTitle = "REAL-TIME WRF"
  res@Footer = False
  pltres = True
  ter_res = True
  opts_ter = ter_res
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  FirstTime = True
  FirstTimeMap = True
  mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for the file
  nd = dimsizes(mdims)
  xlat = wrf_user_getvar(a, "XLAT",0)
  xlon = wrf_user_getvar(a, "XLONG",0)
  ter = wrf_user_getvar(a, "HGT_M",0)
    dims = dimsizes(ter)
;---------------------------------------------------------------


          mpres = True
          pltres = True
          mpres@mpProvincialLineColor="black"         ; province line color
          mpres@mpProvincialLineThicknessF =3         ; province line thickness
                mpres@mpDataBaseVersion="MediumRes"
          mpres@mpDataSetName="./Earth..4"
          mpres@mpAreaMaskingOn = True
          mpres@mpMaskAreaSpecifiers = (/"China"/)
          mpres@mpOutlineSpecifiers = (/"China","China:Provinces"/)
          pltres@FramePlot = False
          optsM = res
          optsM@NoHeaderFooter = True
          optsM@cnFillOn = True
          optsM@lbTitleOn = False
          mpres@mpNationalLineThicknessF= 4         
          mpres@mpGeophysicalLineThicknessF= 4
                mpres@mpNationalLineColor="black"
                 mpres@mpGeophysicalLineColor="black"

      ; Plotting options for terrain               
               gsn_define_colormap(wks,"MPL_terrain")
              opts_ter@cnFillOn             = True
               opts_ter@cnLevelSelectionMode = "ExplicitLevels"       ; set explicit contour levels
         contour_ter  = wrf_contour(a,wks,ter,opts_ter)
    plot = wrf_map_overlays(a,wks,(/contour_ter/),pltres,mpres)
          frame(wks)

          pltres@FramePlot = True

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

end

煮茶的水

感谢分享代码

放逐之刃

谢谢楼主的分享

雨萌萌

楼主有成图可以看下吗？感谢！

342901wzx

geo_em.d01.nc地形数据，在哪可以下载？

大物小子

感谢分享代码

ly00

342901wzx 发表于 2020-5-24 08:55
geo_em.d01.nc地形数据，在哪可以下载？

你好，请问有找到地形图数据么

阿羊出没请注意

ly00 发表于 2021-10-10 14:47
你好，请问有找到地形图数据么

这个是WRF自己跑出来的地形数据