GrADS第6章变量和函数探究.ppt

第6章 变量和函数 ;6.1 变量名 ;例如： ga-open d:\model.ctl open d:\model.le.ctl 上例中，model.ctl文件的序列号为1，model.le.ctl文件的序列号为2。序列号1为缺省情况，可以省略。 ; 括号内的dimexpr：是对当前维数环境进行的局域维数设置表达式，该表达式仅对该变量的维数环境进行局域修正，因而仅对固定的维有效。 ;设置形式可以采用以下两种方式： ;例如： z.3（lev=500） tv.1（time-12hr） rh 缺省的当前文件中的变量rh q.2（t-1,lev=850）;dset ^model.le.dat options little_endian cray_32bit_ieee UNDEF -2.56E33 TITLE 5 Days of Sample Model Output XDEF 72 LINEAR 0.0 5.0 YDEF 46 LINEAR -90.0 4.0 ZDEF 7 LEVELS 1000 850 700 500 300 200 100 TDEF 5 LINEAR 02JAN1987 1DY vars 8 ps 0 99 Surface Pressure u 7 99 U Winds v 7 99 V Winds z 7 99 Geopotential Heights t 7 99 Temperature q 5 99 Specific Humidity ts 0 99 Surface Temperature p 0 99 Precipitation ENDVARS ;dset ^model.grb title Sample Model Data for lats4d Tutorial undef 1e+20 dtype grib index ^model.gmp xdef 72 linear 0.000000 5.000000 ydef 46 linear -90.000000 4.000000 zdef 7 levels 1000 850 700 500 300 200 100 tdef 5 linear 0Z1jan1987 1dy vars 8 ps 0 1, 1, 0, 0 Surface pressure [hPa] ua 7 33,100 Eastward wind [m/s] va 7 34,100 Northward wind [m/s] zg 7 7,100 Geopotential height [m] ta 7 11,100 Air Temperature [K] hus 7 51,100 Specific humidity [kg/kg] ts 0 11,105, 2 Surface (2m) air temperature [K] pr 0 59, 1, 0, 0 Total precipitation rate [kg/(m^2*s)] endvars ;;6.2 表达式;2. 运算规则;例如：;6.3 函数;2. 常用函数;例如： ga-open model.ctl set lev 500 d ave(z,t=1,t=5) 显示500hPa等压面上从第1时次到第5时次变量z的平均量;set parea 1 8 1 6 set grid off set grads off set xlint 20 set xlopts 1 2 0.18 set ylopts 1 2 0.18 set lon 60 180 set lat 0 60 set lev 200 set t 1 define a=ave(v,t=7,t=660,12) d a print c;例如： 求纬向平均： 如果ave(expr,lon=0,lon=360)在端点求 了两次平均；