｛补充｝NCL处理数值时画出来的图范围不对，而且格网也不对

janweis

我的WRFout文件的基本参数如下,
TITLE :  OUTPUT FROM WRF V3.6 MODEL
      START_DATE : 2013-07-20_00:00:00
      SIMULATION_START_DATE : 2013-07-20_00:00:00
      WEST-EAST_GRID_DIMENSION : 103
      SOUTH-NORTH_GRID_DIMENSION : 55
      BOTTOM-TOP_GRID_DIMENSION : 30
      DX : 20000
      DY : 20000
；
WEST-EAST_PATCH_START_UNSTAG : 1
      WEST-EAST_PATCH_END_UNSTAG : 102
      WEST-EAST_PATCH_START_STAG : 1
      WEST-EAST_PATCH_END_STAG : 103
      SOUTH-NORTH_PATCH_START_UNSTAG : 1
      SOUTH-NORTH_PATCH_END_UNSTAG : 54
      SOUTH-NORTH_PATCH_START_STAG : 1
      SOUTH-NORTH_PATCH_END_STAG : 55
      BOTTOM-TOP_PATCH_START_UNSTAG : 1
      BOTTOM-TOP_PATCH_END_UNSTAG : 29
      BOTTOM-TOP_PATCH_START_STAG : 1
      BOTTOM-TOP_PATCH_END_STAG : 30
；
I_PARENT_START : 1
      J_PARENT_START : 1
      PARENT_GRID_RATIO : 1
      CEN_LAT : 40.99999
      CEN_LON : 87
      TRUELAT1 : 30
      TRUELAT2 : 1e+20
      MOAD_CEN_LAT : 40.99999
      STAND_LON : 1e+20
      POLE_LAT : 90
      POLE_LON :  0
      GMT :  0
      JULYR : 2013
      JULDAY : 201
      MAP_PROJ : 3
      MAP_PROJ_CHAR : Mercator

另外，我的ncl文件是根据官网上修改的，出的图挺奇怪的,
第一、他的范围不对，第二，他的网格线也不对
刚刚学习也不是太懂，望大神指导，

; load functions and procedures
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "./WRFUserARW.ncl"

begin
;
; The WRF ARW input file.
; This needs to have a ".nc" appended, so just do it.
; WRF ARW input file (NOTE, your wrfout file does not need the .nc,
; but NCL needs it so make sure to add it in the line below)
  a = addfile("./wrfout_d01_2013-07-20_00:00:00.nc","r")

; Output on screen (type="x11"). Output will be called "plt_Surface1"
; We generate plots, but what kind do we prefer?
  type = "x11"
; type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"plt_PW")
  gsn_define_colormap(wks,"gsdtol")

; Set some basic resources
  res = True
  res@MainTitle                   = "REAL-TIME WRF"     ; Give plot a main title
  res@Footer = False                                    ; Set Footers off
                 
  pltres = True                                         ; Plotting resources
  mpres = True                                          ; Map resources
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  pw = wrf_user_getvar(a,"pw",(/1,6,1/))            ; get only a couple of times
  times = wrf_user_getvar(a,"times",(/1,6,1/))      ; get the same time information
  ntimes = dimsizes(times)                     ; get number of times extracted
  
  opts = res                                            ; Add basic resources                          
  opts@cnFillOn = True                                  ; Shaded plot
  opts@lbLabelBarOn = True                              

  do it = 0,ntimes-1                      ; loop over all times we got from file
    opts@TimeLabel = times(it)                  ; Set Valid time to use on plots
    contour  = wrf_contour(a,wks,pw(it,:,:),opts)
    plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
  end do  
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
end
根本问题
我利用我的wrfout数据一步步往下弄，最后发现了为什么出现我上述的格网不对，是因为我在添加map background的时候出现了问题，即使用plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)将覆盖层放上的时候，就有问题了
warning:MapSetTrans: map limits invalid - using maximal area

longlivehj

janweis 发表于 2015-3-19 21:15
您说的是那个XLAT和XLON么？要怎样才能提取出来与他对应呢？感觉wrf_user_getvar做不到

用读取nc文件变量的方式啊，一个箭头就可以了。坐标本来就是和数据点一一对应的。

longlivehj

图呢？

janweis

file:///root/%E5%9B%BE%E7%89%87/2015-03-18%2017:35:17%20%E7%9A%84%E5%B1%8F%E5%B9%95%E6%88%AA%E5%9B%BE.png

janweis

longlivehj 发表于 2015-3-18 19:07
图呢？

不好意思啦版主，在图书馆，网有点卡，没传上来

janweis

longlivehj 发表于 2015-3-18 19:07
图呢？

版主，我在终端上使用ncl命令直接输出后，发现，我提取出来的数据的坐标信息，并不是按照我的namelist里面所设计的一样，而是我的we和sN值，就是从０开始的，而不是我想要的地理坐标值，这一点应该怎么调整啊

longlivehj

janweis 发表于 2015-3-19 20:53
版主，我在终端上使用ncl命令直接输出后，发现，我提取出来的数据的坐标信息，并不是按照我的namelist里 ...

是这样子的。有专门的变量存储坐标。

janweis

longlivehj 发表于 2015-3-19 20:59
是这样子的。有专门的变量存储坐标。

您说的是那个XLAT和XLON么？要怎样才能提取出来与他对应呢？感觉wrf_user_getvar做不到

janweis

longlivehj 发表于 2015-3-19 21:27
用读取nc文件变量的方式啊，一个箭头就可以了。坐标本来就是和数据点一一对应的。

版主您好，不好意思，又要打扰您了，按照您的指导，我利用我的wrfout数据一步步往下弄，最后发现了为什么出现我上述的格网不对，是因为我在添加map background的时候出现了问题，即使用plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)将覆盖层放上的时候，就有问题了，请问这个问题如何解决
我弄了一个简单的只有HGT数据的图，当不加那个时是这样的，没有网格，起始坐标是０，单位没有
加了之后就有了网格，并且偏移了，变成了单位是度的了

longlivehj

janweis 发表于 2015-3-21 17:53
版主您好，不好意思，又要打扰您了，按照您的指导，我利用我的wrfout数据一步步往下弄，最后发现了为什么 ...

没遇到过这种问题。从wrfout中读取经纬度，范围是正确的吗？wrfout有没有被修改过？