<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: times new roman,new york,times,serif; font-size: 12pt; color: #000000'>Hi DjR,<br><br>if nobody has answared you yet...<br><br>After defining dx, dy and pg, try this:<br><br>'set x 1'<br>'set y 1'<br>'d gr2stn(pg,20.4667,44.8167)' <br><br>Hope this helps,<br>Ricardo<br><br><hr id="zwchr"><blockquote style="border-left:2px solid rgb(16, 16, 255);margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>De: </b>"Djordje Romanic" <djordje8@gmail.com><br><b>Para: </b>"GrADS Users Forum" <gradsusr@gradsusr.org><br><b>Enviadas: </b>Domingo, 17 de Agosto de 2014 17:35:08<br><b>Assunto: </b>[gradsusr] gr2stn doesn't work<br><br><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" target="_blank">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>
<br>_______________________________________________<br>gradsusr mailing list<br>gradsusr@gradsusr.org<br>http://gradsusr.org/mailman/listinfo/gradsusr<br></blockquote><br></div></body></html>