ncl如何mask掉青藏高原地形呢

yang69can

画再分析资料，我想把青藏高原部分mask掉，我现在基本上可以去掉，但是经纬度加不进来，求指教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"

begin
    fnl = addfile("./fnl_20140522_12_00.grb2","r")
    pw = fnl->PWAT_P0_L200_GLL0
    hgt = fnl->HGT_P0_L1_GLL0
    z = fnl->HGT_P0_L100_GLL0({85000},:,:)

    qMask = (/mask(pw,hgt.gt.z,False)/)


    wks = gsn_open_wks("png","MASK")
    res = True
    res@gsnDraw =  True
    res@gsnFrame = True
    res@cnFillOn = True
    res@cnLinesOn = False
    res@mpOutlineOn = True
    res@mpLimitMode = "LatLon"
    res@mpMinLatF = 10
    res@mpMaxLatF = 40
    res@mpMinLonF = 90
    res@mpMaxLonF = 125
    res@cnLevelSelectionMode = "ExplicitLevels"
    res@cnLevels = (/30,35,40,45,50,55/)
    res@cnFillColors = (/"RoyalBlue","CornflowerBlue","Lightblue","LightCyan","Khaki","Gold","Orange"/)


    plot = gsn_csm_contour_map_ce(wks,qMask,res)

end
然后会提示错误：

(0)     check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.
(0)     A valid latitude coordinate array should have a 'units' attribute equal to one of the following values:
(0)         'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0)     check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn't contain one at all.
(0)     A valid longitude coordinate array should have a 'units' attribute equal to one of the following values:
(0)         'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'

stepdance

是这样的图吗？

yang69can

麦田_smile 发表于 2014-12-8 16:41
({85000},:,:)里面的85000是啥意思呢 谢谢

850hPa哈！不同的资料这个不一定一样，像NECP第一类再分析资料里高度的值直接是用hPa表示的。第二类再分析资料是用pa表示的！EC资料也是用hPa表示的。我这里是用的grib2所以是pa。
在ncl中表示高度有两种基本方式，第一种就是以上这种，大括号里边加入实际高度！第二种就是用排序表示。在再分析资料中850hPa是第21层，那么表示为（20，：，：）
希望给你讲清楚了

yang69can

longlivehj

qMask = (/mask(pw,hgt.gt.z,False)/)
(/.../)一加，坐标什么的就木有了！

yang69can

longlivehj 发表于 2014-8-10 17:49
qMask = (/mask(pw,hgt.gt.z,False)/)
(/.../)一加，坐标什么的就木有了！

请问具体怎么加呢？能举个例子不，谢谢

longlivehj

yang69can 发表于 2014-8-10 18:28
请问具体怎么加呢？能举个例子不，谢谢

意思是，让你把(/.../)去掉！

yang69can

longlivehj 发表于 2014-8-10 18:35
意思是，让你把(/.../)去掉！

去掉还是那个错误，没有变

yang69can

longlivehj

yang69can 发表于 2014-8-10 18:41
去掉还是那个错误，没有变

嘿嘿，在mask前加qMask=pw

yang69can

longlivehj 发表于 2014-8-10 20:45
嘿嘿，在mask前加qMask=pw

哇，您好厉害！！ 能告诉我原因吗？谢谢

yang69can

yang69can 发表于 2014-8-10 21:47
哇，您好厉害！！ 能告诉我原因吗？谢谢

还有在问下您，这个图边缘好像不是很平滑，是正常的吗？还是有什么解决方法