<div dir="ltr">Hello,<div><br></div><div>I am using GFS model to generate wind forecasts for different locations. To do an interpolation on a station level, I am using the gr2stn function. However I came across a problem with wind direction there.</div><div><br></div><div>With winds turning around northerly directions, the computed wind direction sometimes displays completely irregular values, for example switching to southeasterly winds in one forecast horizon and then getting back to north. I think there might be an error in handling of ugrd/vgrd data in gr2stn() or I am doing something wrong. </div><div><br></div><div>Here is the way I extract wind direction:</div><div>&#39;d gr2stn(atan2(ugrd10m,vgrd10m),&lt;lon&gt;,&lt;lat&gt;)*57.2957795+180.0)&#39;</div><div><br></div><div><br></div><div>Can you give me a hint what I am doing wrong? Is there another way to interpolate wind direction to stations? Thank you for your help!</div><div><br></div><div><br></div><div>Best regards,</div><div><br></div><div>Jonas</div></div>