- 积分
- 647
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-11-21
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 janweis 于 2015-3-21 23:44 编辑
我的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
|
-
|