...
Note that the average() function does not return the physically correct average of the field but merely the average of all field values using the following formula:
| Mathdisplay |
|---|
average = \frac {1}{N} \sum_{i}^{N}f_{i} |
. To get the physically correct average value of a field, use the integrate() function (which employs a cosine latitude weighting)
...