急！！！！求帮忙关于降水率折线图的一些问题，谢谢大家

Maziy

想用这个脚本来绘制出的降水率随时间变化的折线图，但是run的时候总是会提示出错，以下是脚本
'reinit'
'open h:\zcd3\c50.ctl'
*'e:\13colors.gs'
'set mpdset hires'
'set display color white'
'set grid off'
'set grads off'
'set x 0 121'
'set y 0 121'
'set z 1'
frame = 1
endtime=48
runscript = 1
dist =2
while(runscript)
'set t 'dist
'set x 60'
'set y 60'
'set z 1'
'define aa1=aave(RAINC(t='dist')-RAINC(t='dist'-1),x=0,x=121,y=0,y=121)*1000'
'define aa2=aave(RAINNC(t='dist')-RAINNC(t='dist'-1),x=0,x=121,y=0,y=121)*1000'
'd aa1+aa2'
if(ans)
*'printim 'frame'.gif gif'
*'printim e:\rainrate'frame'.gif gif white'
frame=frame+1
endif
pull dummy
if(dis_t=endtime)
runscript=0
endif
dist = dist+1
endwhile
;
请问是什么问题，他总是提示define这个语句define的那个变量有问题

超台十八

每次都这样，资料来源都没提供，你以为人家这么容易能看出来呀。。。

Maziy

超台十八 发表于 2015-2-12 13:15
每次都这样，资料来源都没提供，你以为人家这么容易能看出来呀。。。

不好意思，忘记附上ctl了

Maziy

这是脚本对应的ctl
dset h:\zcd3\c75.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.4.1 MODEL
xdef   121 linear 0 0.0001
ydef   121 linear 0 0.0001
zdef   79 levels  
   0.25000
   0.50000
   0.75000
   1.00000
   1.25000
   1.50000
   1.75000
   2.00000
   2.25000
   2.50000
   2.75000
   3.00000
   3.25000
   3.50000
   3.75000
   4.00000
   4.25000
   4.50000
   4.75000
   5.00000
   5.25000
   5.50000
   5.75000
   6.00000
   6.25000
   6.50000
   6.75000
   7.00000
   7.25000
   7.50000
   7.75000
   8.00000
   8.25000
   8.50000
   8.75000
   9.00000
   9.25000
   9.50000
   9.75000
  10.00000
  10.25000
  10.50000
  10.75000
  11.00000
  11.25000
  11.50000
  11.75000
  12.00000
  12.25000
  12.50000
  12.75000
  13.00000
  13.25000
  13.50000
  13.75000
  14.00000
  14.25000
  14.50000
  14.75000
  15.00000
  15.25000
  15.50000
  15.75000
  16.10000
  16.50000
  16.75000
  17.00000
  17.25000
  17.50000
  17.75000
  18.00000
  18.25000
  18.50000
  18.75000
  19.00000
  19.25000
  19.50000
  19.75000
  20.00000
tdef   49 linear 00Z01JAN2000       5MN      
VARS  103
LU_INDEX       1  0  LAND USE CATEGORY (-)
VAR_SSO        1  0  variance of subgrid-scale orography (m2)
LAP_HGT        1  0  Laplacian of orography (m)
U             79  0  x-wind component (m s-1)
V             79  0  y-wind component (m s-1)
W             79  0  z-wind component (m s-1)
PH            79  0  perturbation geopotential (m2 s-2)
PHB           79  0  base-state geopotential (m2 s-2)
T             79  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             79  0  perturbation pressure (Pa)
PB            79  0  BASE STATE PRESSURE (Pa)
P_HYD         79  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        79  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        79  0  Cloud water mixing ratio (kg kg-1)
QRAIN         79  0  Rain water mixing ratio (kg kg-1)
QICE          79  0  Ice mixing ratio (kg kg-1)
QSNOW         79  0  Snow mixing ratio (kg kg-1)
QGRAUP        79  0  Graupel mixing ratio (kg kg-1)
QNICE         79  0  Ice Number concentration (  kg-1)
QNRAIN        79  0  Rain Number concentration (  kg(-1))
SHDMAX         1  0  ANNUAL MAX VEG FRACTION (-)
SHDMIN         1  0  ANNUAL MIN VEG FRACTION (-)
SNOALB         1  0  ANNUAL MAX SNOW ALBEDO IN FRACTION (-)
TSLB           5  0  SOIL TEMPERATURE (K)
SMOIS          5  0  SOIL MOISTURE (m3 m-3)
SH2O           5  0  SOIL LIQUID WATER (m3 m-3)
SMCREL         5  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)
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)
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)

lqouc

我觉得不是define的问题，前面还有维度的错误哈。你的ctl指出空间格点121*121
gs里面你写的是0-121，自己想一下0-121是121个点么。

Maziy

lqouc 发表于 2015-2-12 15:20
我觉得不是define的问题，前面还有维度的错误哈。你的ctl指出空间格点121*121
gs里面你写的是0-121，自己 ...

我修改了一下，都改成1 121，但是还是出来一样的错误提示

lqouc

Maziy 发表于 2015-2-12 15:49
我修改了一下，都改成1 121，但是还是出来一样的错误提示

RAINC(t='dist')-RAINC(t='dist'-1),
把-1放在引号里，或者aave之前先单独define下这两个变量，多尝试下啦

南南偶卖糕Oo

请问是否有解决呢？》