...
| Code Block | ||
|---|---|---|
| ||
afs = read('a.grib') # afs is a fieldset
bfs = read('b.grib') # bfs is a fieldset
a = values(afs) # a is a vector or a list of vectors
b = values(bfs) # b is a vector or a list of vectors
spd = sqrt(a*a + b*b) # spd is a vector or a list of vectors
spdfs = set_values(aafs, spd) # write the result back into a fieldset |
...