GRADS画图时出现问题

冰凌

我的用ARWpost转化后的ctl文件如下：
dset d:\a\2010010112wxl.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.3 MODEL
pdef  399 279 lcc  30.000  160.000  200.000  140.000  60.00000  30.00000  -98.00000  30000.000  30000.000
xdef 1379 linear   59.71485   0.13513514
ydef  791 linear  -21.70319   0.13513514
zdef   20 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
  50.00000
tdef    1 linear 06Z01JAN2010     360MN      
VARS   70
LU_INDEX       1  0  LAND USE CATEGORY (-)
U             20  0  x-wind component (m s-1)
V             20  0  y-wind component (m s-1)
W             20  0  z-wind component (m s-1)
PH            20  0  perturbation geopotential (m2 s-2)
PHB           20  0  base-state geopotential (m2 s-2)
T             20  0  perturbation potential temperature (theta-t0) (K)
T_INIT        20  0  initial potential temperature (K)
MU             1  0  perturbation dry air mass in column (Pa)
MUB            1  0  base state dry air mass in column (Pa)
P             20  0  perturbation pressure (Pa)
PB            20  0  BASE STATE PRESSURE (Pa)
P_HYD         20  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        20  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        20  0  Cloud water mixing ratio (kg kg-1)
QRAIN         20  0  Rain water mixing ratio (kg kg-1)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
TOPOSLPX       1  0  ELEVATION X SLOPE (-)
TOPOSLPY       1  0  ELEVATION Y SLOPE (-)
SHDMAX         1  0  ANNUAL MAX VEG FRACTION (-)
SHDMIN         1  0  ANNUAL MIN VEG FRACTION (-)
SNOALB         1  0  ANNUAL MAX SNOW ALBEDO IN FRACTION (-)
LANDUSEF      24  0  LANDUSE FRACTION BY CATEGORY (-)
SOILCTOP      16  0  SOIL CAT FRACTION (TOP) (-)
SOILCBOT      16  0  SOIL CAT FRACTION (BOTTOM) (-)
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 (-)
IVGTYP         1  0  DOMINANT VEGETATION CATEGORY (-)
ISLTYP         1  0  DOMINANT SOIL CATEGORY (-)
VEGFRA         1  0  VEGETATION FRACTION (-)
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)
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)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
CLAT           1  0  COMPUTATIONAL GRID LATITUDE, SOUTH IS NEGATIVE (degree_north)
CLONG          1  0  COMPUTATIONAL GRID LONGITUDE, WEST IS NEGATIVE (degree_east)
ALBBCK         1  0  BACKGROUND ALBEDO (-)
TMN            1  0  SOIL TEMPERATURE AT LOWER BOUNDARY (K)
XLAND          1  0  LAND MASK (1 FOR LAND, 2 FOR WATER) (-)
SNOWC          1  0  FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER) (-)
SR             1  0  fraction of frozen precipitation (-)
height        20  0  Model height (km)
tc            20  0  Temperature (C)
td            20  0  Dewpoint Temperature (C)
rh            20  0  Relative Humidity (%)
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)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.3 MODEL
@ global String comment SIMULATION_START_DATE = 2010-01-01_06:00:00
@ global String comment GRIDTYPE = C
@ global String comment MMINLU = USGS
@ global String comment WEST-EAST_GRID_DIMENSION =   400
@ global String comment SOUTH-NORTH_GRID_DIMENSION =   280
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    28
@ global String comment MAP_PROJ =     1
@ global String comment DX =     30000.00
@ global String comment DY =     30000.00
@ global String comment CEN_LAT =        30.00
@ global String comment CEN_LON =       160.00
@ global String comment TRUELAT1 =        30.00
@ global String comment TRUELAT2 =        60.00
@ global String comment MOAD_CEN_LAT =        30.00
@ global String comment STAND_LON =       -98.00
@ global String comment DIFF_OPT =     1
@ global String comment KM_OPT =     4
@ global String comment DAMP_OPT =     0
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =     3
@ global String comment RA_LW_PHYSICS =     1
@ global String comment RA_SW_PHYSICS =     1
@ global String comment SF_SFCLAY_PHYSICS =     1
@ global String comment SF_SURFACE_PHYSICS =     2
@ global String comment BL_PBL_PHYSICS =     1
@ global String comment CU_PHYSICS =     1
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     1
@ global String comment GRID_FDDA =     1
@ global String comment GRID_ID =     1
@ global String comment PARENT_ID =     0
@ global String comment I_PARENT_START =     1
@ global String comment J_PARENT_START =     1
@ global String comment PARENT_GRID_RATIO =     1
@ global String comment DT =       120.00
@ global String comment ISWATER =    16
@ global String comment ISICE =    24
@ global String comment ISURBAN =     1
@ global String comment ISOILWATER =    14




可是用grads画图却是这样的




这是为什么呢？请各位指教


ARWpost应该是没问题的，GRADS也应该没问题，对于其他文件都是可以用的。可是对于这个文件却出现这样的问题。我的WRF的output是在大型机上转的，但是由于大型机上ARWpost装不起来就下载到自己的电脑上转换的，但是就出现这样的问题

冰凌

图没有显示出来  在附件中的文件中有

MeteoInfo

冰凌 发表于 2012-3-26 21:24

                               
登录/注册后可看大图

图没有显示出来  在附件中的文件中有

pdef  399 279 lcc  30.000  160.000  200.000  140.000  60.00000  30.00000  -98.00000  30000.000  30000.000
xdef 1379 linear   59.71485   0.13513514
ydef  791 linear  -21.70319   0.13513514

这个明显有问题，不太清楚你的模式的区域是怎样的，不过pdef中的参数肯定不对，（200, 140）是模式区域的中心点，从pdef里的参数来看该点的经/纬度为160E/30N（假设这个是正确的），那么标准经度应该也是160E，而不是-98。再者两个标准纬度的顺序应该是先南后北，也就是（30, 60）。在Lambert投影下格点间隔为30000可能也太大了些。

MeteoInfo

两个标准纬度的顺序在实际应用中并不重要，30000的间隔也是有可能的（你的模式区域好像跨度很大），不过标准经度应该是需要改的，所以你可以把pdef改为如下试试：

pdef  399 279 lcc  30.000  160.000  200.000  140.000  60.00000  30.00000 160.00000  30000.000  30000.000

冰凌

MeteoInfo 发表于 2012-3-26 22:33

                               
登录/注册后可看大图

两个标准纬度的顺序在实际应用中并不重要，30000的间隔也是有可能的（你的模式区域好像跨度很大），不过标准 ...

谢谢你的详细解答，我看看

冰凌

MeteoInfo 发表于 2012-3-26 22:19

                               
登录/注册后可看大图

pdef  399 279 lcc  30.000  160.000  200.000  140.000  60.00000  30.00000  -98.00000  30000.000  30 ...

确实是标准经度的问题  谢谢

Eudora99

楼主，您可不可以介绍一下pdef后每一个数据的意义，或者给一下这方面资料的链接，我是刚在学这个，劳驾您了

Eudora99

Eudora99 发表于 2012-9-5 21:21

                               
登录/注册后可看大图

楼主，您可不可以介绍一下pdef后每一个数据的意义，或者给一下这方面资料的链接，我是刚在学这个，劳驾您了

问题已经得到解决：

茉莉花

我用ARWpost处理wrf输出的数据，然后用grads显示，发现只有fields中list出的诊断变量（ fields = 'height,pressure,tk,tc'）可以正常显示，其他变量均不正常。这是什么原因？

杉杉来迟

MeteoInfo 发表于 2012-3-26 22:33
两个标准纬度的顺序在实际应用中并不重要，30000的间隔也是有可能的（你的模式区域好像跨度很大），不过标 ...

标准经度是什么意思呀 我在设wrf的namelist里也遇到这个了 不知道啥意思 是随便设吗