- 积分
- 135
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-3-25
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 rolemodel 于 2017-9-23 17:14 编辑
下载的https://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=192&tid=61014&vid=18上的高度场nc资料,我用addfile看到HGT里面的单位是M,然后我画出来等高线上也是多少千米,请问这个用什么能转化成气压高度?这是代码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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
begin
hg = addfile("/cygdrive/d/yanjiusheng/hgt.2014.nc","r")
hgt = hg->hgt(947,5,:,:)
wks = gsn_open_wks("png","/cygdrive/c/Data/gaodu500")
res = True
res@tiMainString ="hgt500"
res@gsnMaximize = True
res@gsnDraw = False
res@gsnFrame = False
res@mpMinLatF = 17.
res@mpMaxLatF = 55.
res@mpMinLonF = 72.
res@mpMaxLonF = 136.
res@mpFillOn = True
res@mpOutlineOn = False ; Use outlines from shapefile
res@cnFillDrawOrder = "PreDraw"
res@mpDataBaseVersion = "MediumRes"
res@mpDataSetName = "Earth..4"
res@mpAreaMaskingOn = True
res@mpMaskAreaSpecifiers = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)
res@mpLandFillColor = "white"
res@mpInlandWaterFillColor = "white"
res@mpOceanFillColor = "white"
res@mpOutlineBoundarySets = "NoBoundaries"
res@cnFillOn = False
res@cnLinesOn = True
res@cnLineLabelsOn = True
res@cnLineDashSegLenF = 0.18
res@cnLineLabelInterval = 1
res@cnLineLabelPlacementMode = "constant"
res@cnLineLabelFontHeightF = 0.008
res@cnLevelSpacingF = 20.
res@gsnSpreadColors = True
res@lbLabelAutoStride = True
map=gsn_csm_contour_map(wks,hgt,res)
cnres = True
cnres@china = True ;draw china map or not
cnres@river = True ;draw changjiang&huanghe or not
cnres@province = True ;draw province boundary or not
cnres@nanhai = True ;draw nanhai or not
cnres@diqu = False ; draw diqujie or not
chinamap = add_china_map(wks,map,cnres)
draw(map)
frame(wks)
end
|
-
|