用GRADS报错

南宫汗

首先通过ARWpost将WRFout文件转化为GRADS可以识别的CTL和DAT文件无报错，然后修改CTL文件夹中的路径，进行画图，ctl描述文件如下dset D:/test/20160224.dat

options byteswapped

undef 1.e30

title OUTPUT FROM WRF V3.6 MODEL

pdef  99  99 lcc  36.510 103.620   50.000   50.000 60.00000  30.00000  103.62000 27000.000  27000.000

xdef 321 linear   84.05481   0.12162163

ydef 226 linear   22.47371   0.12162163

zdef  59 linear 1 1  

tdef  61 linear 00Z24FEB2016     60MN      

VARS 105

XLAT           1 0  LATITUDE, SOUTH IS NEGATIVE(degree_north)

XLONG          1 0  LONGITUDE, WEST IS NEGATIVE(degree_east)

LU_INDEX       1 0  LAND USE CATEGORY (-)

VAR_SSO        1 0  variance of subgrid-scaleorography (m2)

LAP_HGT        1 0  Laplacian of orography (m)

U             59 0  x-wind component (m s-1)

V             59 0  y-wind component (m s-1)

W             59 0  z-wind component (m s-1)

PH            59 0  set

T             59 0  perturbation potentialtemperature (theta-t0) (K)

MU             1 0  perturbation dry air mass incolumn (Pa)

MUB            1 0  base state dry air mass incolumn (Pa)

NEST_POS       1 0  - (-)

P             59 0  perturbation pressure (Pa)

PB            59 0  BASE STATE PRESSURE (Pa)

P_HYD         59 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        59 0  Water vapor mixing ratio (kgkg-1)

QCLOUD        59 0  Cloud water mixing ratio (kgkg-1)

QRAIN         59 0  Rain water mixing ratio (kgkg-1)

QNWFA2D        1 0  Surface aerosol number concemission (kg-1 s-1)

SHDMAX         1 0  ANNUAL MAX VEG FRACTION (-)

SHDMIN         1 0  ANNUAL MIN VEG FRACTION (-)

SNOALB         1 0  ANNUAL MAX SNOW ALBEDO INFRACTION (-)

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 (Jm-2)

ACSNOM         1 0  ACCUMULATED MELTED SNOW (kgm-2)

SNOW          1  0  SNOWWATER EQUIVALENT (kg m-2)

SNOWH          1 0  PHYSICAL SNOW DEPTH (m)

CANWAT         1 0  CANOPY WATER (kg m-2)

SSTSK          1 0  SKIN SEA SURFACE TEMPERATURE(K)

COSZEN         1 0  COS of SOLAR ZENITH ANGLE(dimensionless)

LAI           1  0  LEAFAREA INDEX (m-2/m-2)

VAR            1 0  OROGRAPHIC VARIANCE (-)

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 onv-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 CUMULUSPRECIPITATION (mm)

RAINSH         1 0  ACCUMULATED SHALLOW CUMULUSPRECIPITATION (mm)

RAINNC         1 0  ACCUMULATED TOTAL GRID SCALEPRECIPITATION (mm)

SNOWNC         1 0  ACCUMULATED TOTAL GRID SCALESNOW AND ICE (mm)

GRAUPELNC      1 0  ACCUMULATED TOTAL GRID SCALEGRAUPEL (mm)

HAILNC         1 0  ACCUMULATED TOTAL GRID SCALEHAIL (mm)

REFL_10CM     59 0  Radar reflectivity (lamda = 10cm) (dBZ)

CLDFRA        59 0  CLOUD FRACTION (-)

SWDOWN         1 0  DOWNWARD SHORT WAVE FLUX ATGROUND SURFACE (W m-2)

GLW            1 0  DOWNWARD LONG WAVE FLUX ATGROUND SURFACE (W m-2)

SWNORM         1 0  NORMAL SHORT WAVE FLUX ATGROUND SURFACE (SLOPE-DEPENDENT) (W m-2)

OLR            1 0  TOA OUTGOING LONG WAVE (W m-2)

ALBEDO         1 0  ALBEDO (-)

CLAT           1 0  COMPUTATIONAL GRID LATITUDE,SOUTH IS NEGATIVE (degree_north)

ALBBCK         1 0  BACKGROUND ALBEDO (-)

EMISS          1 0  SURFACE EMISSIVITY (-)

NOAHRES        1 0  RESIDUAL OF THE NOAH SURFACEENERGY BUDGET (W m{-2})

TMN            1 0  SOIL TEMPERATURE AT LOWERBOUNDARY (K)

XLAND          1 0  LAND MASK (1 FOR LAND, 2 FORWATER) (-)

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 THESURFACE (kg m-2 s-1)

LH             1 0  LATENT HEAT FLUX AT THE SURFACE(W m-2)

ACHFX          1 0  ACCUMULATED UPWARD HEAT FLUX ATTHE SURFACE (J m-2)

ACLHF          1 0  ACCUMULATED UPWARD LATENT HEATFLUX AT THE SURFACE (J m-2)

SNOWC          1 0  FLAG INDICATING SNOW COVERAGE(1 FOR SNOW COVER) (-)

SR             1 0  fraction of frozenprecipitation (-)

AFWA_CLOUD     1 0  AFWA Diagnostic: Cloud ceiling(m)

AFWA_CAPE      1 0  AFWA Diagnostic: ConvectiveAvail Pot Energy (J kg-1)

AFWA_ZLFC      1 0  AFWA Diagnostic: Level of FreeConvection (m)

AFWA_PLFC      1 0  AFWA Diagnostic: Pressure ofLFC (Pa)

MIDRH_MIN      1 0  Min Mid-level relative humidity(%)

AFWA_HAIL      1 0  AFWA Diagnostic: Hail Diameter(Weibull) (mm)

AFWA_LLWS      1 0  AFWA Diagnostic: 0-2000 ft windshear (m s-1)

AFWA_TORNA     1 0  AFWA Diagnostic: Tornado windspeed (Weibull) (m s-1)

LANDMASK       1 0  LAND MASK (1 FOR LAND, 0 FORWATER) (-)

LAKEMASK       1 0  LAKE MASK (1 FOR LAKE, 0 FORNON-LAKE) (-)

SST            1 0  SEA SURFACE TEMPERATURE (K)

SST_INPUT      1 0  SEA SURFACE TEMPERATURE FROMWRFLOWINPUT FILE (K)

pressure      59 0  Model pressure (hPa)

height        59 0  Model height (km)

tk            59 0  Temperature (K)

tc            59 0  Temperature (C)

ENDVARS

@ global String comment TITLE =  OUTPUT FROM WRF V3.6 MODEL

@ global String comment START_DATE =2016-02-24_00:00:00

@ global String commentSIMULATION_START_DATE = 2016-02-24_00:00:00

@ global String commentWEST-EAST_GRID_DIMENSION =   100

@ global String commentSOUTH-NORTH_GRID_DIMENSION =   100

@ global String commentBOTTOM-TOP_GRID_DIMENSION =    60

@ global String comment DX =     27000.00

@ global String comment DY =     27000.00

@ global String comment STOCH_FORCE_OPT=     0

@ global String comment GRIDTYPE = C

@ global String comment DIFF_OPT =     1

@ global String comment KM_OPT =     4

@ global String comment DAMP_OPT =     0

@ global String comment DAMPCOEF =         0.20

@ 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 SF_LAKE_PHYSICS=     0

@ global String commentSURFACE_INPUT_SOURCE =     1

@ global String comment SST_UPDATE =     0

@ global String comment GRID_FDDA =     0

@ global String comment GFDDA_INTERVAL_M=     0

@ global String comment GFDDA_END_H =     0

@ global String comment GRID_SFDDA =     0

@ global String comment SGFDDA_INTERVAL_M=     0

@ global String comment SGFDDA_END_H =     0

@ global String comment HYPSOMETRIC_OPT=     2

@ global String comment SF_URBAN_PHYSICS=     0

@ global String comment SHCU_PHYSICS =     0

@ global String comment MFSHCONV =     0

@ global String comment FEEDBACK =     1

@ global String comment SMOOTH_OPTION=     0

@ global String comment SWRAD_SCAT =         1.00

@ global String comment W_DAMPING =     0

@ global String comment DT =        90.00

@ global String comment RADT =        15.00

@ global String comment BLDT =         0.00

@ global String comment CUDT =         5.00

@ global String comment AER_OPT =     0

@ global String comment SWINT_OPT =     0

@ global String comment AER_TYPE =     1

@ global String comment AER_AOD550_OPT=     1

@ global String comment AER_ANGEXP_OPT=     1

@ global String comment AER_SSA_OPT =     1

@ global String comment AER_ASY_OPT =     1

@ global String comment AER_AOD550_VAL=         0.12

@ global String comment AER_ANGEXP_VAL=         1.30

@ global String comment AER_SSA_VAL =         0.00

@ global String comment AER_ASY_VAL =         0.00

@ global String comment MOIST_ADV_OPT=     1

@ global String comment SCALAR_ADV_OPT=     1

@ global String comment TKE_ADV_OPT =     1

@ global String comment DIFF_6TH_OPT =     0

@ global String comment DIFF_6TH_FACTOR=         0.12

@ global String comment OBS_NUDGE_OPT=     0

@ global String comment BUCKET_MM =        -1.00

@ global String comment BUCKET_J =        -1.00

@ global String comment PREC_ACC_DT =         0.00

@ global String comment SF_OCEAN_PHYSICS=     0

@ global String comment ISFTCFLX =     0

@ global String comment ISHALLOW =     0

@ global String comment ISFFLX =     1

@ global String comment ICLOUD =     1

@ global String comment ICLOUD_CU =     0

@ global String comment TRACER_PBLMIX=     1

@ global String comment SCALAR_PBLMIX=     0

@ global String comment GRAV_SETTLING=     0

@ global String comment DFI_OPT =     0

@ global String commentWEST-EAST_PATCH_START_UNSTAG =     1

@ global String commentWEST-EAST_PATCH_END_UNSTAG =    99

@ global String commentWEST-EAST_PATCH_START_STAG =     1

@ global String commentWEST-EAST_PATCH_END_STAG =   100

@ global String commentSOUTH-NORTH_PATCH_START_UNSTAG =     1

@ global String commentSOUTH-NORTH_PATCH_END_UNSTAG =    99

@ global String comment SOUTH-NORTH_PATCH_START_STAG=     1

@ global String commentSOUTH-NORTH_PATCH_END_STAG =   100

@ global String commentBOTTOM-TOP_PATCH_START_UNSTAG =     1

@ global String commentBOTTOM-TOP_PATCH_END_UNSTAG =    59

@ global String commentBOTTOM-TOP_PATCH_START_STAG =     1

@ global String commentBOTTOM-TOP_PATCH_END_STAG =    60

@ 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 CEN_LAT =        36.51

@ global String comment CEN_LON =       103.62

@ global String comment TRUELAT1 =        30.00

@ global String comment TRUELAT2 =        60.00

@ global String comment MOAD_CEN_LAT =        36.51

@ global String comment STAND_LON =       103.62

@ global String comment POLE_LAT =        90.00

@ global String comment POLE_LON =         0.00

@ global String comment GMT =         0.00

@ global String comment JULYR =  2016

@ global String comment JULDAY =    55

@ global String comment MAP_PROJ =     1

@ global String comment MAP_PROJ_CHAR =Lambert Conformal

@ global String comment MMINLU = USGS

@ global String comment NUM_LAND_CAT =    24

@ global String comment ISWATER =    16

@ global String comment ISLAKE =    -1

@ global String comment ISICE =    24

@ global String comment ISURBAN =     1

@ global String comment ISOILWATER =    14

再用GRADS画垂直方向上的图片时，在GRAS窗口下运行 d u 结果报错 invalid dimension environment looping dimension does not vary 求大神指点

xzllzx

U=>u，ctl中的变量需要写成这种形式

xzllzx

诶 报错是说你的维数没有设置变化画不了剖面，应该是gs有问题？我以往ctl都是写U=>u这样的，下意识回复了一个，不好意思！

南宫汗

xzllzx 发表于 2018-2-27 15:38
诶 报错是说你的维数没有设置变化画不了剖面，应该是gs有问题？我以往ctl都是写U=>u这样的，下意识回复了一 ...

谢谢 ，我是直接在GRADS中输入的命令报错的，然后通过编写GS脚本，垂直方向上可以出图了。顺便想请教一下垂直方向上的维数应该如何正确设置啊，其中Z变量0 59始终有点摸不着头脑啊，求大神不吝赐教啊

南宫汗

xzllzx 发表于 2018-2-27 15:38
诶 报错是说你的维数没有设置变化画不了剖面，应该是gs有问题？我以往ctl都是写U=>u这样的，下意识回复了一 ...

zdef  59 linear 1 1，这条语句是垂直方向上有59个分层的意思吗？但是我在WRF模式中设置的顶高为2000米啊

南宫汗

xzllzx 发表于 2018-2-27 15:35
U=>u，ctl中的变量需要写成这种形式

我以前用ARWpost处理WRFout文件时，生成的Ctl文件，只需要单独修改一下路径，GRADS就可以识别的，不晓得，这次是哪里出错了啊，但是编辑GS脚本却是可以出图，不晓得啥问题啊

xzllzx

南宫汗 发表于 2018-2-28 10:00
我以前用ARWpost处理WRFout文件时，生成的Ctl文件，只需要单独修改一下路径，GRADS就可以识别的，不晓得 ...

应该是59层的意思，但是哪里错了我也不懂...