[求助]WRF模拟输出的数据如何画地形阴影图？

abd

如题，我在http://bbs.06climate.com/forum.php?mod=viewthread&tid=5363 贴中也问了差不多的问题，至今尚未得到解决。虽然http://bbs.06climate.com/forum.php?mod=viewthread&tid=2663 贴中luoziwuhui虽然提供了一个方法，但是一直忽略了后面一大堆人回帖提问的关于地形数据的分辨率的问题。即，他提供的地形分辨率太低，才2.5°×2.5°，垂直方向才12层，无法满足目前较高分辨率数据画图的需求。更不用说我用WRF模拟出来的4公里分辨率的数据了！
所以，我想问问论坛里的高手，我用ARWpost来生成的数据ctl 为如下，能画大概如下图类似的地形覆盖的图吗？
dset ^t20-3.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.3 MODEL
pdef  192 192 lcc  39.500  116.500   96.500   96.500  60.00000  30.00000  117.00000   4000.000   4000.000
xdef  559 linear  111.43380   0.01801802
ydef  418 linear   35.67949   0.01801802
zdef   26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
tdef   49 linear 00Z14JUL2011      30MN
VARS  100
LU_INDEX       1  0  LAND USE CATEGORY (-)
U             26  0  x-wind component (m s-1)
V             26  0  y-wind component (m s-1)
W             26  0  z-wind component (m s-1)
PH            26  0  perturbation geopotential (m2 s-2)
PHB           26  0  base-state geopotential (m2 s-2)
T             26  0  perturbation potential temperature (theta-t0) (K)
MU             1  0  perturbation dry air mass in column (Pa)
MUB            1  0  base state dry air mass in column (Pa)
NEST_POS       1  0  - (-)
P             26  0  perturbation pressure (Pa)
PB            26  0  BASE STATE PRESSURE (Pa)
P_HYD         26  0  hydrostatic pressure (Pa)
Q2             1  0  QV at 2 M (kg kg-1)
T2             1  0  TEMP at 2 M (K)
TH2            1  0  POT TEMP at 2 M (K)
PSFC           1  0  SFC PRESSURE (Pa)
U10            1  0  U at 10 M (m s-1)
V10            1  0  V at 10 M (m s-1)
QVAPOR        26  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        26  0  Cloud water mixing ratio (kg kg-1)
QRAIN         26  0  Rain water mixing ratio (kg kg-1)
QICE          26  0  Ice mixing ratio (kg kg-1)
QSNOW         26  0  Snow mixing ratio (kg kg-1)
QGRAUP        26  0  Graupel mixing ratio (kg kg-1)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
TSLB           4  0  SOIL TEMPERATURE (K)
SMOIS          4  0  SOIL MOISTURE (m3 m-3)
SH2O           4  0  SOIL LIQUID WATER (m3 m-3)
SMCREL         4  0  RELATIVE SOIL MOISTURE (-)
SEAICE         1  0  SEA ICE FLAG (-)
XICEM          1  0  SEA ICE FLAG (PREVIOUS STEP) (-)
SFROFF         1  0  SURFACE RUNOFF (mm)
UDROFF         1  0  UNDERGROUND RUNOFF (mm)
IVGTYP         1  0  DOMINANT VEGETATION CATEGORY (-)
ISLTYP         1  0  DOMINANT SOIL CATEGORY (-)
VEGFRA         1  0  VEGETATION FRACTION (-)
GRDFLX         1  0  GROUND HEAT FLUX (W m-2)
ACGRDFLX       1  0  ACCUMULATED GROUND HEAT FLUX (J m-2)
SNOW           1  0  SNOW WATER EQUIVALENT (kg m-2)
SNOWH          1  0  PHYSICAL SNOW DEPTH (m)
CANWAT         1  0  CANOPY WATER (kg m-2)
SST            1  0  SEA SURFACE TEMPERATURE (K)
SSTSK          1  0  SKIN SEA SURFACE TEMPERATURE (K)
LAI            1  0  Leaf area index (area/area)
H_DIABATIC    26  0  MICROPHYSICS LATENT HEATING (K s-1)
MAPFAC_M       1  0  Map scale factor on mass grid (-)
MAPFAC_MX      1  0  Map scale factor on mass grid, x direction (-)
MAPFAC_MY      1  0  Map scale factor on mass grid, y direction (-)
MF_VX_INV      1  0  Inverse map scale factor on v-grid, x direction (-)
F              1  0  Coriolis sine latitude term (s-1)
E              1  0  Coriolis cosine latitude term (s-1)
SINALPHA       1  0  Local sine of map rotation (-)
COSALPHA       1  0  Local cosine of map rotation (-)
HGT            1  0  Terrain Height (m)
TSK            1  0  SURFACE SKIN TEMPERATURE (K)
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINSH         1  0  ACCUMULATED SHALLOW CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
SNOWNC         1  0  ACCUMULATED TOTAL GRID SCALE SNOW AND ICE (mm)
GRAUPELNC      1  0  ACCUMULATED TOTAL GRID SCALE GRAUPEL (mm)
HAILNC         1  0  ACCUMULATED TOTAL GRID SCALE HAIL (mm)
CLDFRA        26  0  CLOUD FRACTION (-)
SWDOWN         1  0  DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE (W m-2)
GLW            1  0  DOWNWARD LONG WAVE FLUX AT GROUND SURFACE (W m-2)
SWNORM         1  0  NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT) (W m-2)
OLR            1  0  TOA OUTGOING LONG WAVE (W m-2)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
ALBEDO         1  0  ALBEDO (-)
ALBBCK         1  0  BACKGROUND ALBEDO (-)
EMISS          1  0  SURFACE EMISSIVITY (-)
NOAHRES        1  0  RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET (W m{-2})
TMN            1  0  SOIL TEMPERATURE AT LOWER BOUNDARY (K)
XLAND          1  0  LAND MASK (1 FOR LAND, 2 FOR WATER) (-)
UST            1  0  U* IN SIMILARITY THEORY (m s-1)
PBLH           1  0  PBL HEIGHT (m)
HFX            1  0  UPWARD HEAT FLUX AT THE SURFACE (W m-2)
QFX            1  0  UPWARD MOISTURE FLUX AT THE SURFACE (kg m-2 s-1)
LH             1  0  LATENT HEAT FLUX AT THE SURFACE (W m-2)
ACHFX          1  0  ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE (J m-2)
ACLHF          1  0  ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE (J m-2)
SNOWC          1  0  FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER) (-)
SR             1  0  fraction of frozen precipitation (-)
POTEVP         1  0  accumulated potential evaporation (W m-2)
SNOPCX         1  0  snow phase change heat flux (W m-2)
SOILTB         1  0  bottom soil temperature (K)
AVG_FUEL_F     1  0  fuel remaining averaged to atmospheric grid (1)
UAH            1  0  wind at fire_wind_height (m/s)
VAH            1  0  wind at fire_wind_height (m/s)
pressure      26  0  Model pressure (hPa)
height        26  0  Model height (km)
tk            26  0  Temperature (K)
tc            26  0  Temperature (C)
theta         26  0  Potential Temperature (K)
rh            26  0  Relative Humidity (%)
dbz           26  0  Reflectivity (-)
max_dbz        1  0  Max Reflectivity (-)
cape          26  0  CAPE (J/kg)
cin           26  0  CIN (J/kg)
ENDVARS

hzwjy

WRF输出结果有地形高度的，就是你这ctl里的 HGT 变量。

abd

hzwjy 发表于 2012-6-6 15:07

                               
登录/注册后可看大图

WRF输出结果有地形高度的，就是你这ctl里的 HGT 变量。

嗯，我也发现了，但是不知道怎么叠加形成我附加的那个图的效果，知道怎么叠加吗？

7iong

楼主问题解决没，同求

abd

7iong 发表于 2013-3-19 14:33

                               
登录/注册后可看大图

楼主问题解决没，同求

没呢，郁闷！

宁小才

求楼主出解决方法！谢谢

boyq

楼主问题解决没？同求啊

abd

boyq 发表于 2013-4-17 10:13

                               
登录/注册后可看大图

楼主问题解决没？同求啊

没呢，等大侠出来帮忙解决...

七夕路

求大神解救！

七夕路

abd 发表于 2013-4-17 15:29

                               
登录/注册后可看大图

没呢，等大侠出来帮忙解决...

我杂感觉我们做的东西差不多，我留下名字，王妮妮。你有解决方法了恳请告诉我啊。先谢了！