- 积分
- 757
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-12-30
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
输入的数据:把WRF的输出结果用CDO提取出累计降水并运算出某日降水后的nc文件
绘图方法:没有用NCL中专门对WRF后处理的函数,采用一般contour的绘图函数。
问题:将wrf结果中的XLONG,XLAT赋给要绘图的变量后,数据与底图的坐标不匹配。如果不对底图进行经纬度的设定,则把区域模拟的结果整个绘制在全球范围上。
根据模拟的区域对底图进行设定:
res@mpProjection = "LambertConformal"
res@mpLambertMeridianF = 115.
res@mpLimitMode = "LatLon"
res@mpMinLatF = 37.
res@mpMaxLatF = 42.
res@mpMinLonF = 114.
res@mpMaxLonF = 120.0
结果如图:
--------------------------------------------------------------------------------------------
但是,如果不对底图范围进行设置,就会出现下图的情况:
----------------------------------------------------------------------------------------------
下面是ncl——filedump对数据的描述:
dimensions:
x = 265
y = 255
XTIME = 1 // unlimited
variables:
float XLONG ( y, x )
standard_name : longitude
long_name : longitude
units : degree_east
_CoordinateAxisType : Lon
float XLAT ( y, x )
standard_name : latitude
long_name : latitude
units : degree_north
_CoordinateAxisType : Lat
double XTIME ( XTIME )
standard_name : time
units : minutes since 2012-07-10 00:00:00
calendar : standard
axis : T
float RAINNC ( XTIME, y, x )
units : mm
coordinates : XLONG XLAT
FieldType : 104
MemoryOrder : XY
description : ACCUMULATED TOTAL GRID SCALE PRECIPITATION
stagger :
-------------------------------------------------------------------------------------------------
完整代码如下:
begin
a = addfile("precip.nc","r")
precip = a->RAINNC
Lon = a->XLONG
Lat = a->XLAT
precip@lat2d = Lat
precip@lon2d = Lon
precip@units = "mm"
printVarSummary(precip)
wks = gsn_open_wks("png","precip1")
res = True
res@gsnMaximize = True
res@mpDataSetName = "Earth..4"
res@mpDataBaseVersion = "MediumRes"
res@mpOutlineBoundarySets = "National"
res@mpOutlineOn = True
res@mpOutlineBoundarySets = "NoBoundaries"
res@mpGeophysicalLineThicknessF= 4.
res@mpNationalLineThicknessF= 4.
res@mpOutlineSpecifiers = (/"China:states","taiwan"/)
res@mpFillOn = False
res@mpGeophysicalLineColor = "black"
res@mpNationalLineColor = "black"
res@mpUSStateLineColor = "black"
res@mpUSStateLineThicknessF = 3.
res@tfDoNDCOverlay = True ; Tell NCL you are doing a native plot
res@gsnAddCyclic = False
;res@mpProjection = "LambertConformal"
;res@mpLambertMeridianF = 115.
;res@mpLimitMode = "LatLon"
;res@mpMinLatF = 37.
;res@mpMaxLatF = 42.
;res@mpMinLonF = 114.
;res@mpMaxLonF = 120.0
res@cnFillOn = True
res@cnLinesOn = False
res@cnLevelSpacingF = 5. ; NCL chose 4.0
res@tfDoNDCOverlay = True
res@gsnAddCyclic = False
res@gsnLeftString = "Total Precipitation"
res@tiMainString = "WRF Simulation"
plot = gsn_csm_contour_map(wks,precip(0,:,:),res)
end
|
|