<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Jeff,<br>
<br>
My problem comes in using cdiff. Everything else should be fine (it is practically verbatim from the webpage). I have a suspicion it is because it uses x and y as arguments on two different lines but I don't really have much of a way to know without more Grads
 experience. Could this be causing the problem?<br>
<br>
Mason <br>
<div style="font-family: Times New Roman; color: rgb(0, 0, 0); font-size: 16px;">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF503599"><font color="#000000" face="Tahoma" size="2"><b>From:</b> gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of Jeff Duda [jeffduda319@gmail.com]<br>
<b>Sent:</b> Friday, November 11, 2011 3:17 PM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] Diagnostics in Ctl File?<br>
</font><br>
</div>
<div></div>
<div>Mason,<br>
I preferred to derive the formulas for spacings myself.&nbsp; Simple trigonometry reveals the following that you should use for dx and dy (i.e., meridional and zonal grid spacing):<br>
dx = Re*cos(phi)*d(theta)<br>
dy = Re*d(phi)<br>
<br>
where Re is the radius of the Earth (6371 to 6374 km), phi is latitude, theta is longitude, and the d operator means &quot;change in&quot;.&nbsp; To convert this to Grads scripting code, you use the cdiff command.&nbsp; Also note that all angle measures need to be in radians,
 so you'll have to convert them as you do the operation.&nbsp; Here are the literal Grads command lines I use for this:<br>
<br>
re = 6.371e3<br>
'define dx = 're'*cos(lat*(3.1419/180))*cdiff(lon*(3.14159/180))'<br>
'define dy = 're'*cdiff(lon*(3.14159/180))'<br>
<br>
Obviously you can assign the conversion factor pi/180 to a script variable to simplify the symbolism there:<br>
<br>
conv_factor = 3.14159/180<br>
e.g., 'define dy = 're'*cdiff(lon*'conv_factor')'<br>
<br>
If you try to display dx or dy, you should get a nice little image that contours or colors dx and dy (dx should decrease towards the poles and not change along a line of constant latitude, and dy should only be a function of latitude as well).&nbsp; Now use this
 to define temperature advection (assuming the temperature variable is tmpprs and wind is ugrdprs/vgrdprs):<br>
<br>
'define tadv = -ugrdprs*cdiff(tmpprs,x)/(2*dx) - vgrdprs*cdiff(tmpprs,y)/(2*dy)'<br>
<br>
The <a href="http://98cbar.gs" target="_blank">98cbar.gs</a> is indeed just a script that defines color levels and sets a color scheme for when you plot temperature advection.&nbsp; Changing the color tables is rather easy in Grads.&nbsp; You only need to define a color
 using 'set rgb # R G B', where # is anything between 16 and 255 (just a reference number) and R, G, and B are between 0 and 255.&nbsp; Then you set the contour intervals manually using 'set clevs , , , , , ...' where you insert literally the levels you want to
 use (separated by SPACES, not COMMAS).&nbsp; Then set the colors for those levels using 'set ccols , , , , , , ,&#43;1' again using the reference numbers for each color and separating by SPACES.&nbsp; The &quot;&#43;1&quot; is for the fact that you need one additional color for the number
 of levls you put.&nbsp; You can consult the <a href="http://www.iges.org/grads/gadoc/colorcontrol.html" target="_blank">
controlling colors</a> page for full info on what I just said.&nbsp; Also, I have a number of very simple color bar scripts that I've used in the past.&nbsp; I can provide them to all if anyone is interested.<br>
<br>
You can also create a simple color bar using the <a href="http://cbar.gs" target="_blank">
cbar.gs</a> script which is provided in the <a href="http://www.iges.org/grads/gadoc/library.html" target="_blank">
script library</a>.<br>
<br>
Let me know if you need any more help.<br>
<br>
Jeff<br>
<br>
<div class="gmail_quote">On Fri, Nov 11, 2011 at 2:50 PM, Rowell, Mason D. <span dir="ltr">
&lt;<a href="mailto:Mason.D.Rowell-1@ou.edu" target="_blank">Mason.D.Rowell-1@ou.edu</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>
<div style="direction: ltr; font-family: Tahoma; color: rgb(0, 0, 0); font-size: 10pt;">
Jeff and others,<br>
<br>
Strange. This is exactly what was suggested in the ctl file near the comments at the bottom (that which I pulled from nomads). The only difference is that they have also the inclusion of a
<a href="http://98cbar.gs" target="_blank">98cbar.gs</a> script right after the other expressions like that in this link you gave me, whose utility is a mystery since it is not necessary, at least suggested as so by the cdiff documenation given that it is not
 included. Perhaps it is just some coloring scheme. At any rate, the script fails at the first line (which attempts to define dtx), with each line coded exactly as it is shown (with the appropriate NARR temp as TMPprs and wind as UGRDprs and VGRDprs of course).
 I assume it will do so again at dty, and so forth. I originally didn't include define with these since I thought they weren't necessary but adding that in didn't solve the problem anyways. I suspect it is because cdiff calls x, and since I define my domain
 with lat long, it doesn't have a x variable to use? There definitely isn't a spacing issue this time...<br>
<br>
Mason <br>
<div style="font-family: Times New Roman; color: rgb(0, 0, 0); font-size: 16px;">
<hr>
<div style="direction: ltr;"><font color="#000000" face="Tahoma" size="2"><b>From:</b>
<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] on behalf of Jeff Duda [<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>]<br>
<b>Sent:</b> Thursday, November 10, 2011 11:11 AM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] Diagnostics in Ctl File?<br>
</font><br>
</div>
<div></div>
<div>Mason,<br>
There is a good example of how to code temperature advection in the page for the cdiff command:
<a href="http://www.iges.org/grads/gadoc/gradfunccdiff.html" target="_blank">http://www.iges.org/grads/gadoc/gradfunccdiff.html</a><br>
The reason for all of the crazy manipulations is because the data points on your grid are not uniformly spaced on the globe.&nbsp; Thus corrections have to be made to get the correct distance between individual grid points.<br>
<br>
Jeff<br>
<br>
<div class="gmail_quote">On Thu, Nov 10, 2011 at 11:00 AM, Wesley Ebisuzaki <span dir="ltr">
&lt;<a href="mailto:Wesley.Ebisuzaki@noaa.gov" target="_blank">Wesley.Ebisuzaki@noaa.gov</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Mason,<br>
<br>
I think that you've discovered the &quot;private&quot; commands used by nomads<br>
plot software for defining new variables<br>
and custom plots. &nbsp; The nomads plot package is a wrapper for GrADS. &nbsp;It<br>
has the ability to run grads .gs files that<br>
are embedded within the ctl files.<br>
<br>
<br>
<a href="http://www.cpc.ncep.noaa.gov/products/wesley/ncep_data/plot_doc.htmlhttp://www.cpc.ncep.noaa.gov/products/wesley/ncep_data/plot_doc.html" target="_blank">http://www.cpc.ncep.noaa.gov/products/wesley/ncep_data/plot_doc.htmlhttp://www.cpc.ncep.noaa.gov/products/wesley/ncep_data/plot_doc.html</a><br>
<br>
Ask NCDC for a copy of the .gs files.<br>
<div><br>
&gt; Has anyone done this before with grads?<br>
</div>
How do you expect the developers to debug the scripts?<br>
<font color="#888888"><br>
&nbsp; &nbsp; Wesley<br>
</font>
<div>
<div></div>
<div><br>
<br>
<br>
Rowell, Mason D. wrote:<br>
&gt; All,<br>
&gt;<br>
&gt; I noticed ample comments toward the end of the ctl files available on the nomads servers from the NARR page, with a description of how to implement these in grads scripting. I tried to take this to plot temperature advection but it didn't work, b/c the last
 line of the comments calls for some function never provided. Has anyone done this before with grads?<br>
&gt;<br>
&gt; Mason<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
<font color="#888888"><br>
<br clear="all">
<br>
-- <br>
Jeff Duda<br>
Grad student - PhD, Meteorology<br>
University of Oklahoma School of Meteorology - Center for Analysis and Prediction of Storms<br>
<br>
</font></div>
</div>
</div>
</div>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
Jeff Duda<br>
Grad student - PhD, Meteorology<br>
University of Oklahoma School of Meteorology - Center for Analysis and Prediction of Storms<br>
<br>
</div>
</div>
</div>
</body>
</html>