[gradsusr] arbitrary cross section

[Davide Sacchetti]

it seems your request is outside file domain


open file:
X is varying   Lon = -58.489 to -54.793   X = 1 to 463
Y is varying   Lat = 20.948 to 24.548   Y = 1 to 451


request:

lon1 = -53.56lon2 = -52.745
lat1 = 19.532
lat2 = 18.927




bye
Davide


Dott. Davide Sacchetti
Dipartimento Stato dell’Ambiente e Tutela dai Rischi Naturali
U.O. Clima Meteo Idro
Viale Brigate Partigiane 2, 16129 Genova (I)
tel: +39 0106437535
mail: davide.sacchetti at arpal.liguria.it  web: https://www.arpal.liguria.it/homepage/meteo.html



 Da:   Sim Aberson <sim.aberson at noaa.gov> 
 A:   GrADS Users Forum <gradsusr at gradsusr.org> 
 Inviato:   24/03/2022 15:50 
 Oggetto:   [gradsusr] arbitrary cross section 

Attenzione: Questa e` la prima volta che ricevi un'email da questo mittente. 
Attenzione: Assicurati che sia qualcuno di cui ti fidi. 
 
I am trying to plot a cross section of wind speed through the center of 
a hurricane.  I adapted a script I found many places online, but I only 
get errors.  The script is: 
 
'open TEDDY20L_202009181830_Postr_Mean.ctl' 
'set grads off' 
'set zlog on' 
'set x 1' 
'set y 1' 
'set lev 1000 100' 
lat1 = 19.532 
lon1 = -53.56 
lat2 = 18.927 
lon2 = -52.745 
lon = lon1 
'collect 1 free' 
while (lon <= lon2) 
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1) 
  say ' lat,lon = 'lat' 'lon 
 'collect 1 
gr2stn(mag(maskout(ugrdprs,hgtprs),maskout(vgrdprs,hgtprs)),'lon','lat')' 
  lon = lon + 0.1 
endwhile 
 
'set x 14 16' 
'set xaxis 'lon1' 'lon2 
'set clab on' 
'set gxout shaded' 
'd coll2gr(1,-u)' 
'set gxout contour' 
'set cint 5' 
'd coll2gr(1,-u)' 
 
When I run this, I get: 
 
 
 lat,lon = 19.532 -53.56 
Data Request Warning:  Request is completely outside file limits 
Data Request Warning:  Request is completely outside file limits 
Data Request Warning:  Request is completely outside file limits 
etc. 
 
All the points are within the domain, but I get all undefined values. 
The file dimensions are 
 
ga-> q dims 
Default file number is: 1 
X is fixed     Lon = -58.489  X = 1 
Y is fixed     Lat = 20.948  Y = 1 
Z is varying   Lev = 1000 to 100   Z = 1 to 43 
T is fixed     Time = 18Z18SEP2020  T = 1 
E is fixed     Ens = 1  E = 1 
 
Obviously I'm doing something wrong, but I have no idea what.  I tried 
it without the maskout function.  I tried it on a single variable.  Same 
results, so it's not those. 
 
You can find the data and control file at 
 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 . 
 
 
--  
Sim Aberson (he/him/his) 
NOAA/AOML/Hurricane Research Division   Miami, FL 
 https://urlsand.esvalabs.com/?u=https%3A%2F%2Fwww.aoml.noaa.gov%2Fhrd%2Fpeople%2Fsimdaberson%2F&e=7312d17f&h=4a5afe79&f=y&p=n  
 
AOML stands on the land of the Tequesta, Calusa, Miccosukee and Seminole 
peoples.  This acknowledgment is one way to educate ourselves about this 
land, its history, and our relationships with the land and each other. 
_______________________________________________ 
gradsusr mailing list 
gradsusr at gradsusr.org 
 https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr&e=7312d17f&h=5ee7e1e3&f=y&p=n  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20220325/2f6192eb/attachment.html>