登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
篇幅可能不是很长，但是个人觉得还是蛮有用的（对于我这种初学者来说）
一、变量名 1、完全的变量名形式为： abbrev.file#(dimexpr,dimexpr,...) abbrev是数据描述文件中给出的变量名缩写。 file# 为包含此变量的已打开的文件序号，缺省为1，用"set dfile" 命令可定义当前的缺省文件。 dimexpr 为对当前维数环境进行的局域维数设置表达式，该表达式仅对该变量的维数环境进行局域修正，因而仅对固定的维有效。 2、绝对维数和相对维数： 绝对维数表达式为： X|Y|Z|T|LON|LAT|LEV|TIME =value 相对维数表达式为： X|Y|Z|T|LON|LAT|LEV|TIME +/-offset 例如： z.3(lev=500) 表示文件3中高度为500hPa等压面上的变量Z tv.1(time-12hr) 相对于当前时刻之前12小时时刻的1号文件中的变量tv rh 缺省的当前文件中的变量rh q.2(t-1,lev=850) 2号文件中相对于当前时刻前一时刻，高度为850hPa面上的变量q z(t+0) 其用法后面将介绍 注意：GrADS内部保留一些"预先定义" 的变量名，lat，lon ，lev，time，result,rc。 二、表达式 1、表达式的构成 GrADS中的表达式由运算符、运算域和括号组成。其中括号用于控制运算的次序； 2、运算符 + 加，- 减，* 乘，/ 除； 3、运算域 变量，函数和常数； 注意： 1） 运算是对相同网格点上的不同变量进行； 2） 运算时只要有一个变量在某格点的值为缺测，则该网格点的运算结果也为缺测值； 3） 被0除的结果也设为缺测。 4） 不能对不同维数尺度（即维数变动的范围不同）的格点进行运算操作，这个问题在同时打开多个描述文件时可能会碰到。 5） 如果某一格点的变动维数多于其他格点，则具有较少变动维数的格点其维数环境将被“拓展”，以便于运算操作。 例如： z - z(t-1) (变量Z的时间变化) t(lev=500)-t(lev=850)(500到850hPa等压面的变量t的变化) ave(z,t=1,t=5) (变量Z在第一到第五时刻的平均) z -ave(z,lon=0,lon=360,-b) (变量Z的纬向偏差)
内容提要一、新变量的定义二、气候变量的定义 一、定义一个新的变量1、语法： define varname = expr 特点：1）新变量可以用于随后的表达式中（它也可以用于随后的define和display命令中）。2）新变量不是储存在硬盘上而是在内存中，所以，要尽量避免定义维数太大的变量。 2、注意维数环境的设置1） 定义的变量可以有 0－4个变化的维数。2） 当4个维数同时变化时，define命令是GrADS中所有命令中 唯一可行的命令。3） 定义一个4维的变量，在define命令执行后，记住改变维数环境，使之 变化的维数小于4。4） 当Z和\或T维数变化时，define命令自动逐个固定Z和T的维数来计算表达式。你如何安排你的表达式，将影响计算的过程。 例如，假设定义一个4维的变量， set lon -180 0 set lat 0 90 set lev 1000 100 set t 1 10 define temp = rh 在define命令执行后，记住改变维数环境，使之变化的维数小于4。set t 5 set lev 500 d temp 3、固定维的理解如果你定义的一个变量具有某些固定的维数，随后使用这个变量，那么 固定的维数将被处理为 与该维数无关的变量。 例如：set lon -180 0 set lat 0 90 set lev 500 set t 10 define zave = ave(z,t=1,t=30) 这个变量的固定维数为Z和T，在任何的Z和T的设置下，它都可以被使用。如：set t 1 set lev 200 d zave 即使维数固定在200hPa，你仍然能够得到你定义它时的500hPa上的时间平均值。 更完整的请看附件
Grads讲义.doc
(47 KB, 下载次数: 92)
| 
窥视卡
雷达卡
发表于 2013-9-12 11:30:48
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
