Vanwege de soms onhandige dataset van GFS moet je alles per coordinaat (voor GFS iedere 0.5 graden) apart berekenen, daar gaat helaas veel tijd in zitten. Het plotten zelf kost nog geen seconde 
Je bedoelt dus dat je alles in een loop gooit? Dat is een bijzonder inefficiente manier van berekenen. De kracht van Python ligt juist in het feit dat je operaties op gehele arrays kunt uitvoeren. Probeer dus eerst de ruwe variabelen van GFS in een numpy array te zetten en vervolgens daar berekeningen op los te laten.
Quote selectie