grads提取wrf模式后处理的高度场，画图出错，请指点！

渭水听风

是这样的，我用arwpost处理完wrf数据，生成ctl和dat文件：
ctl文件：dset g:/0604/hs/0604_hs_pre_namelist.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.2.1 MODEL
pdef  258 195 lcc  38.800  100.500  129.500   98.000  38.80000  38.80000  100.80000   9000.000   9000.000
xdef  753 linear   85.19873   0.04054054
ydef  423 linear   29.86214   0.04054054
zdef   19 levels  
1000.00000
950.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
450.00000
400.00000
350.00000
300.00000
250.00000
200.00000
150.00000
100.00000
tdef   25 linear 00Z04JUN2012     180MN      
VARS  109
LU_INDEX       1  0  LAND USE CATEGORY (-)
U             19  0  x-wind component (m s-1)
V             19  0  y-wind component (m s-1)
W             19  0  z-wind component (m s-1)
PH            19  0  perturbation geopotential (m2 s-2)
PHB           19  0  base-state geopotential (m2 s-2)
T             19  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             19  0  perturbation pressure (Pa)
PB            19  0  BASE STATE PRESSURE (Pa)
P_HYD         19  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        19  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        19  0  Cloud water mixing ratio (kg kg-1)
QRAIN         19  0  Rain water mixing ratio (kg kg-1)
QICE          19  0  Ice mixing ratio (kg kg-1)
QSNOW         19  0  Snow mixing ratio (kg kg-1)
QGRAUP        19  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)
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)
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        19  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)
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)
pressure      19  0  Model pressure (hPa)
geopt         19  0  Geopotential (m2/s2)
height        19  0  Model height (km)
tk            19  0  Temperature (K)
tc            19  0  Temperature (C)
theta         19  0  Potential Temperature (K)
rh            19  0  Relative Humidity (%)
clflo          1  0  Low Cloud Fraction (%)
clfmi          1  0  Mid Cloud Fraction (%)
clfhi          1  0  High Cloud Fraction (%)
rh2            1  0  Relative Humidity at 2m (%)
wspd          19  0  Wind Speed (m s-1)
wdir          19  0  Wind Direction (Degrees)
ws10           1  0  Wind Speed at 10 M (m s-1)
wd10           1  0  Wind Direction at 10 M (Degrees)
u10m           1  0  Rotated wind component (m s-1)
v10m           1  0  Rotated wind component (m s-1)
slp            1  0  Sea Levelp Pressure (hPa)
max_dbz        1  0  Max Reflectivity (-)
cape          19  0  CAPE (J/kg)
cin           19  0  CIN (J/kg)
mcape          1  0  MCAPE (J/kg)
mcin           1  0  MCIN (J/kg)
lcl            1  0  LCL (meters AGL)
lfc            1  0  LFC (meters AGL)
ENDVARS
提取高度场的gs文件是：
'reinit'
'open g:/0604/hs/0604_hs_pre_namelist.ctl'
'set gxout fwrite'
'set fwrite g:/0604/hs/grd2sta/geo/500/geo500.dat'
'set t 9'
'set z 11'
'd (ph+phb)/98.1'
'disable fwrite'
'reinit'
在写一个ctl文件：
dset g:/0604/hs/grd2sta/geo/500/geo500.dat
*options  byteswapped
undef 1.e30
*pdef  258 195 lcc  38.800  100.500  129.500   98.000  38.80000  38.80000  100.80000   9000.000   9000.000
xdef  753 linear   85.19873   0.04054054
ydef  423 linear   29.86214   0.04054054
zdef   1 levels 500  
tdef   1 linear 00Z05JUN2012     180MN      
VARS  1
h  1  99
endvars
然后写个gs文件画图：
'reinit'
'open g:/0604/hs/grd2sta/geo/500/geo500.ctl'
'set gxout contour'
'set mpdset cnworld'
'set lev 500'
'set t 1'
'd h'
'printim  g:/0604/hs/grd2sta/geo/500/geo500-1.png white'
画出的图是：geo500-1.png
file:///C:/Users/asus/AppData/Roaming/Tencent/Users/380420721/QQ/WinTemp/RichOle/H44EPOQ4J0~V)QO%7BGT~S16N.jpg

实际应该是geo500.png
[img]file:///C:/Users/asus/AppData/Roaming/Tencent/Users/380420721/QQ/WinTemp/RichOle/XH[%7D_[G13_(JIX~3C3]~RNT.jpg[/img]
请教各位大侠，哪儿出问题了，帮忙看看!谢谢了！

渭水听风

自己先顶一下，高手帮帮忙指点一下！

abd

你对原始ctl文件直接打开后
'set t 9'
'set z 11'
'd (ph+phb)/98.1'
画出来看看不行吗？

abd

你先一步一步查错，你先画出那两个变量，再画出对应高度和事件的
'd (ph+phb)/98.1'

兰溪之水

顺便检查下缺省值~

渭水听风

abd 发表于 2012-8-16 22:41

                               
登录/注册后可看大图

你对原始ctl文件直接打开后
'set t 9'
'set z 11'

这样可以的，画出来的图是正确的。

Forcast

3楼的是一种解决办法。最主要的原因是楼主
pdef  258 195 lcc  38.800  100.500  129.500   98.000  38.80000  38.80000  100.80000   9000.000   9000.000
xdef  753 linear   85.19873   0.04054054
ydef  423 linear   29.86214   0.04054054
没有对它们的理解和进行试验，想当然的去转存数据造成的。
这个数据实际只有pdef  258 195 个点，但GRADS画图时可以内插到753*423个点上再画图。数据能否转存fwrite为753*423个点没有试过。

渭水听风

兰溪之水 发表于 2012-8-16 22:54

                               
登录/注册后可看大图

顺便检查下缺省值~

缺省值难道不是模式结果输出后默认的吗？

渭水听风

Forcast 发表于 2012-8-17 08:25

                               
登录/注册后可看大图

3楼的是一种解决办法。最主要的原因是楼主
pdef  258 195 lcc  38.800  100.500  129.500   98.000  38.80 ...

请教牛人，模式结果用arwpost转出来的ctl文件里面直接就是这样的，我看前面也有人用类似方法提取过降水资料！如果这样不行，那应该怎么弄，求赐教！谢谢

Forcast

渭水听风 发表于 2012-8-17 08:37

                               
登录/注册后可看大图

请教牛人，模式结果用arwpost转出来的ctl文件里面直接就是这样的，我看前面也有人用类似方法提取过降水资 ...

你可以把你后面给输出量的CTL里的PDF前*去掉试试，也就是你的CTL与原始的CTL一样，只是变量从多个变成1个而已