底图区域不对 经纬度不显示

追风的阳光

下面是我的脚本文件，请大家帮我指正错误，谢谢！
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/wrf/WRFUserARW.ncl"
begin
a=addfile("/home/zssapr/WRF3.5.1/20110710/WRFV3/test/em_real/wrfout_d01_2011-07-10_00:00:00","r")
r =a->RAINNC
e =a->RAINC

lat =a->XLAT
lon =a->XLONG
;printVarSummary(lat)
print(lat)
dimll =dimsizes(lat)
nlat =dimll(0)
mlon =dimll(1)
wks =gsn_open_wks("ps","rain_contour2d02")
gsn_define_colormap(wks,"precip4_11lev")


res  =True
res@tiMainString      ="rain"
res@cnFillOn            = True          ; turn on color
res@cnLinesOn           = True
res@gsnAddCyclic          = False       ; regional data: not cyclic
res@tfDoNDCOverlay        = True        ; set True for native (direct) mapping
res@mpProjection          = "Mercator"  ;"LambertConformal"

res@mpLimitMode           = "Corners"
res@mpLeftCornerLatF      = lat
res@mpLeftCornerLonF      = lon

; res@mpRightCornerLatF     = lat(:,nlat-1,mlon-1)
; res@mpRightCornerLonF     = lon(:,nlat-1,mlon-1)
res@mpOutlineDrawOrder    = "PostDraw"  ; draw continental outline last
res@mpFillDrawOrder       = "PreDraw"
res@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries
res@mpFillOn              = False       ; turn off map fill

rain =r(0,:,:)+r(1,:,:)+r(2,:,:)+r(3,:,:)+r(4,:,:)+r(5,:,:)+r(6,:,:)+r(7,:,:)+e(0,:,:)+e(1,:,:)+e(2,:,:)+e(3,:,:)+e(4,:,:)+e(5,:,:)+e(6,:,:)+e(7,:,:)plot = gsn_csm_contour_map(wks,rain,res)draw(plot)
end

freekiller

res@pmTickMarkDisplayMode = "Always"

river

ncl为什么标题里要加上grads呢

一水天下

还在想跟楼主请教grads出图怎么这么好看……

追风的阳光

river 发表于 2014-4-18 19:47
ncl为什么标题里要加上grads呢

不知道。。难道是在grads版发的。。。我记得我点的是ncl，请问这个脚本有错在何处呢？

追风的阳光

freekiller 发表于 2014-4-18 16:29
res@pmTickMarkDisplayMode = "Always"

我明天拿我电脑试试。。

river

追风的阳光 发表于 2014-4-19 15:13
不知道。。难道是在grads版发的。。。我记得我点的是ncl，请问这个脚本有错在何处呢？

我不用NCL所以不清楚哎······

longlivehj

建议使用wrf_map_resources从文件中读取绘图区域信息。具体而言，注释掉下面的行：
res@mpProjection          = "Mercator"  ;"LambertConformal"
res@mpLimitMode           = "Corners"
res@mpLeftCornerLatF      = lat
res@mpLeftCornerLonF      = lon
res@mpRightCornerLatF     = lat(:,nlat-1,mlon-1)
res@mpRightCornerLonF     = lon(:,nlat-1,mlon-1)
用res = wrf_map_resources(a,res)这句替换！

追风的阳光

longlivehj 发表于 2014-4-19 21:02
建议使用wrf_map_resources从文件中读取绘图区域信息。具体而言，注释掉下面的行：
res@mpProjection      ...

最后用了http://www.ncl.ucar.edu/Applications/Scripts/WRF_pcp_1.ncl的这个脚本，成功出图了。后来发现wrf输出的U、V量是c网格的，计算散度等变量时需要把他们插值到一点上去，请问具体用哪一个脚本比较好呢？
linint2_points还是其他？谢谢！

longlivehj

追风的阳光 发表于 2014-4-21 20:42
最后用了http://www.ncl.ucar.edu/Applications/Scripts/WRF_pcp_1.ncl的这个脚本，成功出图了。后来发现 ...

用wrf_user_unstagger来做插值