<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 Dj,<br>there are some tricks in building this procedure. I tested here with a GFS file and it worked (opengrads 2.0.a7.oga.3, but I think that the version is not a problem):<br><br>'open gfs.ctl'<br><br>* An area around your points. Note that, at least with the gfs model, the longitude have to be set in the positive way (0 to 360, or some range in between that):<br>'set lat 0 40'<br>'set lon 0 60'<br><br>'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(prmslmsl/100,x)/dx, cdiff(prmslmsl/100, y)/dy)'<br><br>*Below are three necessary steps:<br>'set x 1'<br>'set y 1'<br>'set gxout print'<br><br>* Here, longitude comes before latitude:<br>'d gr2stn(pg,44.8167,20.4667)'<br>* to print the result on the terminal or, <br>* alternatively in a gs script, you can capture the result:<br>resline=sublin(result,3)<br>pgvar=subwrd(resline,1)<br>say 'PGrad = 'pgvar<br><br>Hope this can solve your problem.<br><br>Regards,<br>Ricardo<br><br><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>Segunda-feira, 18 de Agosto de 2014 20:02:52<br><b>Assunto: </b>Re: [gradsusr] gr2stn doesn't work<br><br><div dir="ltr"><div><div><div><div>Hi Ricardo, <br><br></div>It didn't help. If I do that (or I set x and y to any other value), I get the following error message:<br>'Invalid station data dimension environment'. <br>
<br></div>I also tried to set station coordinates to a grid point (say lon = 20 and lat = 45), but it didn't work either. <br><br></div>By the way, my version of GrADS is v2.1.a1. Could that be an issue? <br><br></div>
Thanks,<br>Djordje<br></div>
<br>_______________________________________________<br>gradsusr mailing list<br>gradsusr@gradsusr.org<br>http://gradsusr.org/mailman/listinfo/gradsusr<br></blockquote><br></div></body></html>