第八章Grads定义新变量.doc

GRADS (Grid Analysis and Display System)讲义（8） 内容提要一、新变量的定义二、气候变量的定义一、定义一个新的变量1、语法： define varname = expr 特点：1）新变量可以用于随后的表达式中（它也可以用于随后的define和display命令中）。2）新变量不是储存在硬盘上而是在内存中，所以，要尽量避免定义维数太大的变量。2、注意维数环境的设置1）定义的变量可以有04个变化的维数。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上的时间平均值。4、变化维的理解当定义变量具有变化的维数环境时，然后再将这个维数固定起来，这时所得的变量将被固定在这个维数上：set lon -180 0 set lat 0 90 set lev 500 set t 10 define temp = z set lat 40 d temp 将会显示定义格点中400N上的资料。set lat -40 d temp 定义变量400S上的资料将被显示，当400S超过了原来变量定义的维数环境，这时，资料值将显示为缺测。6、局地维数环境的设置和理解1）如果局地环境是在定义变量时变化的维数环境内，则显示局地维数环境下的变量值。2）如果这个维数在定义变量时是固定的，则局地维数环境将失效。d temp(lat=50) d temp(t=15) 以上显示命令中，变量temp在T维数是固定的，因此t=15的局地维数环境将失效。注意：目前的define命令仅仅可用于格点资料。 二、定义气候变量1、语法： 如果你定义了格点变量，你可以告诉GrADS新变量是一个气候值，即：你想在更宽的时间维数环境下使用新变量。其命令是：modify varname varname：是定义变量名（define在这个命令前使用）。Seasonal： GrADS假定格点资料是月平均的气候值，Diurnal：格点资料被假定为平均的日变化值（这时资料的时间间隔应小于一天）。2、举例 以下是一个计算十年平均气候值的例子：set lat 0 40set lon 60 150 set z 1 set t 1 12 define have = ave(h,t+0,t=48,1yr)气候场（12 have），距平场（48 hano） 下面我们再来看一下如何计算一个变量的年际异常值：modify have seasonal set t 1 48define hano=h-haved hano 在这个例子中，所求的是10年120个时次的异常值，但zave变量仅在时次1到12有定义，通过modify命令设置它为气候值后，当计算时次13、25（第二、三年的一月份）时，自动使用时次1的zave值（一月份的气候值），其余月份类推。合起来：set lat 0 40set lon 60 150 set z 1 set t 1 12 define have = ave(h,t+0