Operators are:
+
Addition -
Subtraction *
Multiplication /
Division =
Equal ( ==
is also acceptable)!=
Not Equal>
Greater than>=
Greater than or equal <
Less than <=
Less than or equal&
And ( &&
is also acceptable)|
Or ( ||
is also acceptable)Operands are:
variable specifications, functions, and constants
.Operations cannot be done between grids that have different scaling in their varying dimensions -- i.e., grids that have different rules for converting the varying dimensions from grid space to world coordinate space. This can only be encountered when you are attempting operations between grids from different files that have different scaling rules.
If one grid has more varying dimensions than the other, the grid with fewer varying dimensions is 'expanded' and the operation is performed.
Expression evaluation in GrADS is recursive, so that multiple expressions may be nested together.
(GrADS version 2.0.a7+) Variable specifications can include a dimension expression to set time as an offset from the variable's initial time.
(GrADS version 2.1.1.b0+) The result of a logical operation is boolean -- an answer to a yes/no question. If the expression is true the result is 1, if the expression is false the answer is -1 (instead of zero).
This is slightly different from the usual convention, but it is implemented this way in GrADS to make it easier to use logical operators with the maskout()
function. An if()
function has also been implemented to use logical operators in expressions of the form if-then-else. Note: the logical operators have not yet been implemented for station data.
slp/100
(Convert sea level pressure units from hPa to mb)z-z(t-1)
(Height change over one time step)z-z(offt=0)
(Height change since initial time)t(lev=500)-t(lev=850)
(Temp difference between 500
and 850mb)ave(z,t=1,t=5)
(Average of z over first 5 times in file) sum(prec(offt+0),t=1,t=4)
(Accumulated precipitation -- sum of 2nd through 5th time steps in file) z-ave(z,lon=0,lon=360,-b)
(Remove zonal mean)tloop(aave(p,global))
(Time series of globally averaged precip)tsfc>=0
(Surface temperatures greater than or equal to zero)maskout(tsfc,tsfc>=0)
(Set all negative surface temperatures to be undefined)(tsfc>30)|(tsfc<-30)
(Surface temperature extremes)(lat>=-5)&(lat<=5)&(lon>=190)&(lon<=240)
(Lat/Lon bounding box -- Nino3.4)maskout(tsfc,(lat>=-5)&(lat<=5)&(lon>=190)&(lon<=240))
(Surface temperature only defined within the Lat/Lon box)