<div dir="ltr"><div><div>Hello, <br><br></div>I am calculating the pressure gradients (PGs) from the NCEP/NCAR reanalysis I MSLP data. The resolution of the data is 2.5 x 2.5 degrees. I want to interpolate PGs to a particular coordinate that is not on the reanalysis grid. I was thinking to use gr2stn function for that but it does not work. This is my code: <br>

<br></div><div>&#39;sdfopen <a href="http://1948_1950.nc">1948_1950.nc</a>&#39;<br></div><div>&#39;define dx = 6.371e3 * cos(lat * 3.14159/180) * cdiff(lon * 3.14159/180, x)&#39;<br>&#39;define dy = 6.371e3 * cdiff(lat * 3.14159/180, y)&#39;<br>

&#39;define pg = mag(cdiff(slp,x)/dx, cdiff(slp, y)/dy)&#39;<br><br></div>&#39;d gr2stn(pg,20.4667,44.8167)&#39;  *This coordinate is just an example<br><br><br><div>This is the error message:<br><br>Error from GR2STN: 2nd argument is not a station expression<br>

Operation Error:  Error from gr2stn function<br>  Error ocurred at column 1<br>DISPLAY error:  Invalid expression <br>  Expression = gr2stn(pg,20.4667,44.8167)<br><br></div><div>Please could you help me with this?<br><br>

Thanks,<br>DjR<br></div></div>