<div dir="ltr"><div>Hi<br><br>z would refer to the lev (level), just as x would refer to lon and y refer to lat. If your descriptor is 'zdef levels 1000 900 800', that means set z 1 would mean set lev 1000. That is unless you have 1000 levels in your data. <br>
</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jul 18, 2013 at 10:20 AM, John D. Wilson <span dir="ltr"><<a href="mailto:jaydee.uu@ualberta.ca" target="_blank">jaydee.uu@ualberta.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Ron,<br>
<br>
I'd been using "set z 1 45" (for instance). It never occurred to me I<br>
could specify in mb. However I've just now tried<br>
<br>
set x 38<br>
set y 38<br>
set z 1000 900<br>
d ugrdhlev<br>
<br>
and got<br>
<br>
Data Request Error: Invalid grid coordinates<br>
Varying dimension 2 decreases: 1000 to 900<br>
Error ocurred getting variable 'ugrdhlev'<br>
Error ocurred at column 1<br>
DISPLAY error: Invalid expression<br>
Expression = ugrdhlev<br>
<br>
This is even though q ctlinfo returns<br>
<br>
ugrdhlev 58 0,2,2 ** (58 57 56 55 54.. 5 4 3 2 1) U-Component Of Wind<br>
[m/s]<br>
<br>
<br>
I also tried specifiying an interval for z in Pa -- similar complaint. John.<br>
<div class="HOEnZb"><div class="h5"><br>
On 17/07/13 18:02, Goodson,Ron [Edm] wrote:<br>
> Hi John<br>
><br>
> If you original levels were in pressure or height or sigma .. then I'm surprised it doesn't plot that as the vertical coord (if that is what you meant). Pretend the vertical coord is in mb in your original data. Did you do something like "set z 1 10" or did you do "set z 1000 500". The latter should cause the vertical coord to be mb.<br>
><br>
><br>
> If, however, you actually want all those different coord systems, then you will have to interpolate U to the various levels you want .. using something like the grads-script "<a href="http://zinterp.gs" target="_blank">zinterp.gs</a>" (in the script library)<br>
><br>
><br>
> Ron Goodson<br>
> Prairie and Northern Meteorological Service of Canada Science Section Unité des Sciences, Région des Prairies et du Nord, Service Météorologique du Canada<br>
> Environment Canada<br>
> Room 200 ; 4999 - 98 Avenue<br>
> Edmonton Alberta T6B 2X3 Environnement Canada<br>
> Pièce 200 ; 4999 - 98 Avenue<br>
> Edmonton Alberta T6B 2X3<br>
> Telephone | Téléphone 780 951 8791<br>
> Email | Courriel <a href="mailto:ron.goodson@ec.gc.ca">ron.goodson@ec.gc.ca</a><br>
><br>
><br>
><br>
> -----Original Message-----<br>
> From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On Behalf Of John D. Wilson<br>
> Sent: July 17, 2013 4:42 PM<br>
> To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> Subject: [gradsusr] plotting wind speed (or velocity) versus height rather than level number<br>
><br>
> Hi,<br>
><br>
> I'd like to plot a wind profile at fixed (X,Y,t). It's straightforward to plot (for instance) the u-component versus level number (``d<br>
> ugrdhlev'') but I'd like instead to display u versus height, or u versus pressure, or u versus sigma. I expect this is really easy but I haven't yet stumbled on the trick. Can anyone help? Thanks, John _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
<br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
J.D. Wilson (Prof)<br>
Rm. 3-277, CCIS Building<br>
Department of Earth & Atmospheric Sci.<br>
University of Alberta, Edmonton, AB T6G2E3<br>
<a href="tel:780%20492%200353" value="+17804920353">780 492 0353</a><br>
</font></span><div class="HOEnZb"><div class="h5"><br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</div></div></blockquote></div><br></div>