爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 28423|回复: 49

[秀图] [求助]谁会画有地形这种图?

[复制链接]
发表于 2012-6-6 10:23:46 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
如题,我在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 为如下,能画大概如下图类似的地形覆盖的图吗?
dixing.jpg

我的数据的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
密码修改失败请联系微信:mofangbao
发表于 2012-6-6 11:24:19 | 显示全部楼层
其实我去做过这件事情,一种办法就是把2.5*2.5的数据插值为1*1,这样会存在一个问题,插值后地形的偏差较大。

另外,不插值,分别画图,然后覆盖,这样的缺点是等值线会出现断裂,特别是风场,会向地面吹。

最后发现,两种方法都不咋地
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-6-6 15:22:51 | 显示全部楼层

谢谢你的回复,目前我最重要的问题是WRF输出的数据里做叠加形成我附加的图那个效果,先对于这个数据不存在数据分辨率的问题!ctl 里有个变量为HGT,也就是26个层的地形数据,但是不知道怎么叠加形成附加的那个图的效果,知道怎么叠加吗?
密码修改失败请联系微信:mofangbao
发表于 2012-6-6 16:57:54 | 显示全部楼层
ARWpost是什么软件? 较高精度的地形数据可以用ETOPO
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-6-6 18:58:21 | 显示全部楼层
wlzhongouc 发表于 2012-6-6 16:57
ARWpost是什么软件? 较高精度的地形数据可以用ETOPO

转过WRF吗?处理模式输出数据生成ctl的软件。
密码修改失败请联系微信:mofangbao
发表于 2012-6-6 20:59:23 | 显示全部楼层
abd 发表于 2012-6-6 18:58
转过WRF吗?处理模式输出数据生成ctl的软件。

没有转过啊
密码修改失败请联系微信:mofangbao
发表于 2012-6-6 22:55:18 | 显示全部楼层
这种?我记得我当时grads画的时候用地形数据maskout画的
1.jpg
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-6-6 23:22:48 | 显示全部楼层
陌小默的 发表于 2012-6-6 22:55
这种?我记得我当时grads画的时候用地形数据maskout画的

是的是的!就这样,maskout 具体怎么用啊?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-6-6 23:26:20 | 显示全部楼层
本帖最后由 abd 于 2012-6-6 23:27 编辑
陌小默的 发表于 2012-6-6 22:55
这种?我记得我当时grads画的时候用地形数据maskout画的

我最多用:set z 1 26set lon 117.5
set lat 38 42
set gxout shaded
d height
set gxout shtream
d v;w*20
画出以下图:
地形图.png

请问用maskout 怎么画出你画的那种图啊???????


密码修改失败请联系微信:mofangbao
发表于 2012-6-6 23:27:17 | 显示全部楼层
abd 发表于 2012-6-6 10:26
我最多用:set z 1 26set lon 117.5
set lat 38 42
set gxout shaded

我的看下,好久没弄了,明天说,我要断电了,呜呜
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表