Zit je trouwens ook in Python te werken? Als je de NetCDF4 module gebruikt om de GFS OpenDAP data binnen te hengelen kan je daar ook een grid met latitudes vandaan halen.
zoiets:
lats = f.variables['lat'][:]
Dan zou je een oppervlak gewogen gemiddelde kunnen krijgen door de wortel van cosines van de latitude te gebruiken als gewicht.
Bijvoorbeeld:
mean_wind = np.average(wind, weights=np.sqrt(np.cos(np.deg2rad(lats))))
edit:
Moet de wortel vd cosines zijn ipv enkel de cosines...
Quote selectie