...
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 |
...