...
| Code Block | ||
|---|---|---|
| ||
# Derive cross sections of temperature forecast and analysis
# in two netcdf variables, tfc_xs and tan_xs....
(...)
# derive the list of netcdf variables
var_list = variables(tfc_xs)
# create output netcdf
diff_xs = tfc_xs
# loop over variables and create fc-an difference cross-section
for i = 1 to count(var_list) do
setcurrent(tan_xs, i)
setcurrent(tfc_xs, i)
setcurrent(diff_xs, i)
diff_xs = tfc_xs - tan_xs
end for |
Extracting NetCDF values
If you need to have access to the data values of a netcdf current variable, you can use function values() :
val_list = values(netcdf)
which returns a vector with all the values for the current variable.You can then use the Vector functions to manipulate the data. This technique could even be used to create a new Geopoints variable with the netCDF data, or to insert the values into a GRIB field.
An alternative method for accessing individual values is to the use the function value() :
val = value(netcdf, n)
which returns the nth value from the current netcdf variable.