<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>'sdfopen <a href="http://1948_1950.nc">1948_1950.nc</a>'<br></div><div>'define dx = 6.371e3 * cos(lat * 3.14159/180) * cdiff(lon * 3.14159/180, x)'<br>'define dy = 6.371e3 * cdiff(lat * 3.14159/180, y)'<br>
'define pg = mag(cdiff(slp,x)/dx, cdiff(slp, y)/dy)'<br><br></div>'d gr2stn(pg,20.4667,44.8167)' *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>