<html><head></head><body><div>it seems your request is outside file domain</div><div><br></div><div>open file:</div><div>X is varying   Lon = -58.489 to -54.793   X = 1 to 463</div><div>Y is varying   Lat = 20.948 to 24.548   Y = 1 to 451</div><div><br></div><div>request:</div><div><div>lon1 = -53.56</div>lon2 = -52.745<div>lat1 = 19.532</div><div>lat2 = 18.927</div><div><br></div></div><div><br></div><div>bye</div><div>Davide</div><div><br></div><div style="margin: 0px; padding: 0px; box-sizing: content-box; color: rgb(34, 34, 34);"><i style="box-sizing: content-box;"><font size="2" style="box-sizing: content-box;"><span id="mcntmcntext-gen4583" style="box-sizing: content-box; font-family: Calibri, sans-serif;">Dott. Davide Sacchetti</span></font></i></div><div style="margin: 0px; padding: 0px; box-sizing: content-box; color: rgb(34, 34, 34);"><i style="box-sizing: content-box;"><font size="2" style="box-sizing: content-box;"><span id="mcntmcntext-gen4583" style="box-sizing: content-box; font-family: Calibri, sans-serif;"></span></font></i></div><i style="box-sizing: content-box; color: rgb(34, 34, 34);"><font size="2" style="box-sizing: content-box;"><span id="mcntmcntext-gen4583" style="box-sizing: content-box; font-family: Calibri, sans-serif;">Dipartimento Stato dell’Ambiente e Tutela dai Rischi Naturali<br style="box-sizing: content-box;">U.O. Clima Meteo Idro<br style="box-sizing: content-box;">Viale Brigate Partigiane 2, 16129 Genova (I)<br style="box-sizing: content-box;">tel: +39 0106437535<br style="box-sizing: content-box;">mail: <a href="mailto:mauro.damonte@arpal.liguria.it" class="mailto" title="Send email to mauro.damonte@arpal.liguria.it" target="_blank" style="box-sizing: content-box;">davide.sacchetti@arpal.liguria.it</a>  web: <a href="http://www.meteoliguria.it/" title="http://www.meteoliguria.it" target="_blank" style="box-sizing: content-box;">https://www.arpal.liguria.it/homepage/meteo.html</a></span></font></i><br><br><br><div><strong>

Da:

</strong>

Sim Aberson <sim.aberson@noaa.gov>

<br>

<strong>

A:

</strong>

GrADS Users Forum <gradsusr@gradsusr.org>

<br>

<strong>

Inviato:

</strong>

24/03/2022 15:50

<br>

<strong>

Oggetto:

</strong>

[gradsusr] arbitrary cross section

<blockquote class="mori" style="margin:0 0 0 .8ex;border-left:1px solid #CCC;padding-left:1ex;">Attenzione: Questa e` la prima volta che ricevi un'email da questo mittente.

Attenzione: Assicurati che sia qualcuno di cui ti fidi.

<br>

<br>I am trying to plot a cross section of wind speed through the center of

<br>a hurricane.  I adapted a script I found many places online, but I only

<br>get errors.  The script is:

<br>

<br>'open TEDDY20L_202009181830_Postr_Mean.ctl'

<br>'set grads off'

<br>'set zlog on'

<br>'set x 1'

<br>'set y 1'

<br>'set lev 1000 100'

<br>lat1 = 19.532

<br>lon1 = -53.56

<br>lat2 = 18.927

<br>lon2 = -52.745

<br>lon = lon1

<br>'collect 1 free'

<br>while (lon <= lon2)

<br>  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)

<br>  say ' lat,lon = 'lat' 'lon

<br> 'collect 1

<br>gr2stn(mag(maskout(ugrdprs,hgtprs),maskout(vgrdprs,hgtprs)),'lon','lat')'

<br>  lon = lon + 0.1

<br>endwhile

<br>

<br>'set x 14 16'

<br>'set xaxis 'lon1' 'lon2

<br>'set clab on'

<br>'set gxout shaded'

<br>'d coll2gr(1,-u)'

<br>'set gxout contour'

<br>'set cint 5'

<br>'d coll2gr(1,-u)'

<br>

<br>When I run this, I get:

<br>

<br>

<br> lat,lon = 19.532 -53.56

<br>Data Request Warning:  Request is completely outside file limits

<br>Data Request Warning:  Request is completely outside file limits

<br>Data Request Warning:  Request is completely outside file limits

etc.

<br>

All the points are within the domain, but I get all undefined values.

<br>The file dimensions are

<br>

<br>ga-> q dims

<br>Default file number is: 1

<br>X is fixed     Lon = -58.489  X = 1

<br>Y is fixed     Lat = 20.948  Y = 1

<br>Z is varying   Lev = 1000 to 100   Z = 1 to 43

<br>T is fixed     Time = 18Z18SEP2020  T = 1

<br>E is fixed     Ens = 1  E = 1

<br>

<br>Obviously I'm doing something wrong, but I have no idea what.  I tried

<br>it without the maskout function.  I tried it on a single variable.  Same

results, so it's not those.

<br>

<br>You can find the data and control file at

<br> https://urlsand.esvalabs.com/?u=https%3A%2F%2Fwww.aoml.noaa.gov%2Fftp%2Fhrd%2Faberson%2FHEDASanalyses%2F&e=7312d17f&h=3fbe6db9&f=y&p=n .

<br>

<br>

<br>-- 

<br>Sim Aberson (he/him/his)

<br>NOAA/AOML/Hurricane Research Division   Miami, FL

<br> https://urlsand.esvalabs.com/?u=https%3A%2F%2Fwww.aoml.noaa.gov%2Fhrd%2Fpeople%2Fsimdaberson%2F&e=7312d17f&h=4a5afe79&f=y&p=n 

<br>

<br>AOML stands on the land of the Tequesta, Calusa, Miccosukee and Seminole

<br>peoples.  This acknowledgment is one way to educate ourselves about this

land, its history, and our relationships with the land and each other.

<br>_______________________________________________

<br>gradsusr mailing list

<br>gradsusr@gradsusr.org

<br> https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr&e=7312d17f&h=5ee7e1e3&f=y&p=n 

<br></blockquote></div></body></html>