warning:MapSetTrans: map limits invalid - using maximal area

janweis

我的测试区域的namelist.wps设置如下，而后利用wps的geogrid.exe生成的geo_em_d01.nc，而后利用ncl查看HGT_M出现了一个错误，warning:MapSetTrans: map limits invalid - using maximal area。我把geo_em_d01.nc上传了，希望各位坛友帮助我，这个问题已经困扰了好久了，要是有坛友能抽空帮忙在WPS下跑一下，看看是不是我的WPS有问题，还是我生的geo_em_d01.nc有问题，还是我的NCL文件有问题，谢谢各位了，在此感激不尽！回帖就给金币啊！
[初步找到问题了]我在不更换其它参数，只修改投影时发现，我设置lambert投影时没有问题，设置mercator 投影时，就出现问题，应该是我的mercator投影设置有问题，但是不知道怎么解决，拜托各位了
namelist.wps
&share
wrf_core = 'ARW',
max_dom = 1,
io_form_geogrid = 2,
/

&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 13,
j_parent_start = 1, 18,
s_we = 1,
e_we = 82, 169,
s_sn = 1,
e_sn = 69, 91,
geog_data_res = '10m','2m',
dx = 30000,
dy = 30000,
map_proj = 'mercator',
ref_lat = 40,
ref_lon = 81,
truelat1 = 60,
geog_data_path = '/media/root/Media/WPS_GEOG/'
/
NCL文件
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"  ; Add special wrf functions

begin

a = addfile("../../WPS/geo_em.d01.nc","r")  


; type = "x11"
  type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"geo_mercator_5")       ; Create a plot workstation


  opts = True                                ; Set some Basic Plot options
  opts@MainTitle = "GEOGRID FIELDS"

  ter = a->HGT_M(0,:,:)                      ; Read the variable to memory

  res = opts                                 ; Use basic options for this field
  res@cnFillOn = True                        ; Create a color fill plot

  contour = wrf_contour(a,wks,ter,res)

  pltres = True                              ; Set plot options
  mpres = True                               ; Set map options
  plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)  ; Plot the data over a map background

end

longlivehj

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin

a = addfile("geo_em.d01.nc","r")  

  type = "pdf"
  wks = gsn_open_wks(type,"geo_mercator_5")

  ter = a->HGT_M(0,:,:)
  lat = a->XLAT_M(0, :, :)
  lon = a->XLONG_M(0, :, :)

  ter@lat2d = lat
  ter@lon2d = lon

  res = True
  res@cnFillOn = True
  res@cnLinesOn = False

  res@gsnAddCyclic = False
  res@mpLimitMode = "Corners"
  res@mpLeftCornerLatF = min(lat)
  res@mpLeftCornerLonF = min(lon)
  res@mpRightCornerLatF = max(lat)
  res@mpRightCornerLonF = max(lon)

  plt = gsn_csm_contour_map(wks, ter, res)

end

又是那隻貓

这应该是画图的问题

janweis

又是那隻貓 发表于 2015-3-30 12:36
这应该是画图的问题

版主，是我的那个NCL画图文件有错误对吗？还望您给予指导，如何修改，感激不尽啊版主
我刚刚是不是把那个回帖奖励没调好，不知道金币给您了没有

才仔

我就是来看看

janweis

才仔 发表于 2015-3-30 13:48
我就是来看看

版主，赚点金币不容易啊，还望帮忙解决问题啊

youyou

不懂楼主所说的东西 ，就是为了10块钱而来，不过还是希望楼主能早点解决疑问！！

evita

处于学习阶段中，纯顶楼（赚金币）而来

海中淡水鱼

意思是你要画的数据超过了你的底图设置的最大范围，两个可能：第一，你的底图设置的小了，比如你画全球的图但是经度只有359等等；第二，你画图的时候去的数据大了，例如你的100-200范围的图d的时候用了100-201.（没用过NCL，单纯从错误角度分析的，希望有所帮助。

子龙

好事，好事呀。