All,<div><br></div><div>   The most authoritative documentation on the hdivg function can be found here:</div><div><br></div><div><a href="http://opengrads.cvs.sourceforge.net/viewvc/opengrads/cola/grads2/src/gafunc.c?revision=1.1.1.9&amp;view=markup">http://opengrads.cvs.sourceforge.net/viewvc/opengrads/cola/grads2/src/gafunc.c?revision=1.1.1.9&amp;view=markup</a></div>
<div><br></div><div><a href="http://opengrads.cvs.sourceforge.net/viewvc/opengrads/cola/grads2/src/gafunc.c?revision=1.1.1.9&amp;view=markup"></a>That is, the source code :-). The hdivg() implementation can be found on line 2716. This file has the implementation of most grads functions. </div>
<div><br></div><div>BTW, the link above is the web interface to the opengrads CVS repository. We keep in this particular directory the original releases from COLA, before any opengrads patches are applied:</div><div><br></div>
<div><a href="http://opengrads.cvs.sourceforge.net/viewvc/opengrads/cola/grads2/src/">http://opengrads.cvs.sourceforge.net/viewvc/opengrads/cola/grads2/src/</a></div><div><br></div><div>There is a pulldown menu called &quot;Sticky Tag&quot; that you can use to select a particular grads version.</div>
<div><br></div><div>   Arlindo</div><div><br></div><div><br></div><div><br><br><div class="gmail_quote">On Mon, Nov 15, 2010 at 12:44 AM, Yaqiang Wang <span dir="ltr">&lt;<a href="mailto:yaqiang.wang@gmail.com">yaqiang.wang@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">The correct one should be:<br>
<br>
define du = cdiff(u,x)<br>
<div class="im">define dx = cdiff(lon,x)*3.1416/180<br>
</div>define dv = cdiff(v*cos(lat*3.1416/180),y)<br>
<div class="im">define dy = cdiff(lat,y)*3.1416/180<br>
display (du/dx+dv/dy)/(6.37e6*cos(lat*3.1416/180))<br>
<br>
Yaqiang<br>
<br>
</div><div class="im">Chinese Academy of Meteorological Sciences<br>
<a href="mailto:yaqiang.wang@gmail.com">yaqiang.wang@gmail.com</a><br>
<br>
</div><div><div></div><div class="h5">On Mon, Nov 15, 2010 at 12:47 PM, Yaqiang Wang &lt;<a href="mailto:yaqiang.wang@gmail.com">yaqiang.wang@gmail.com</a>&gt; wrote:<br>
&gt; Thanks for Jeff&#39;s reply! But it seems that in GrADS the hdivg function<br>
&gt; was calculated as following. Please mention the last sentence.<br>
&gt;<br>
&gt; define dv = cdiff(v,x)<br>
&gt; define dx = cdiff(lon,x)*3.1416/180<br>
&gt; define du = cdiff(u*cos(lat*3.1416/180),y)<br>
&gt; define dy = cdiff(lat,y)*3.1416/180<br>
&gt; display (du/dy+dv/dx)/(6.37e6*cos(lat*3.1416/180))<br>
&gt;<br>
&gt; Yaqiang<br>
&gt;<br>
&gt; Chinese Academy of Meteorological Sciences<br>
&gt; <a href="mailto:yaqiang.wang@gmail.com">yaqiang.wang@gmail.com</a><br>
&gt;<br>
&gt; On Mon, Nov 15, 2010 at 11:56 AM, Jeffrey Duda &lt;<a href="mailto:jdduda@iastate.edu">jdduda@iastate.edu</a>&gt; wrote:<br>
&gt;&gt; Yaqiang,<br>
&gt;&gt; You are exactly correct.<br>
&gt;&gt;<br>
&gt;&gt; Jeff Duda<br>
&gt;&gt;<br>
&gt;&gt; On Sun, Nov 14, 2010 at 9:22 PM, Yaqiang Wang &lt;<a href="mailto:yaqiang.wang@gmail.com">yaqiang.wang@gmail.com</a>&gt;<br>
&gt;&gt; wrote:<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; The GrADS documentation said that the cdiff function may be used to<br>
&gt;&gt;&gt; duplicate the calculation done by the hcurl function:<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; define dv = cdiff(v,x)<br>
&gt;&gt;&gt; define dx = cdiff(lon,x)*3.1416/180<br>
&gt;&gt;&gt; define du = cdiff(u*cos(lat*3.1416/180),y)<br>
&gt;&gt;&gt; define dy = cdiff(lat,y)*3.1416/180<br>
&gt;&gt;&gt; display (dv/dx-du/dy)/(6.37e6*cos(lat*3.1416/180))<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; But there is no description of hdivg function. I guess hdivg function<br>
&gt;&gt;&gt; could be duplicated using cdiff function as following. Please anybody<br>
&gt;&gt;&gt; tell me if it is correct or not. Thanks!<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; define dv = cdiff(v,x)<br>
&gt;&gt;&gt; define dx = cdiff(lon,x)*3.1416/180<br>
&gt;&gt;&gt; define du = cdiff(u*cos(lat*3.1416/180),y)<br>
&gt;&gt;&gt; define dy = cdiff(lat,y)*3.1416/180<br>
&gt;&gt;&gt; display (du/dx+dv/dy)/(6.37e6*cos(lat*3.1416/180))<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Yaqiang Wang<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Chinese Academy of Meteorological Sciences<br>
&gt;&gt;&gt; <a href="mailto:yaqiang.wang@gmail.com">yaqiang.wang@gmail.com</a><br>
&gt;&gt;&gt; _______________________________________________<br>
&gt;&gt;&gt; gradsusr mailing list<br>
&gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt;&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; --<br>
&gt;&gt; Jeff Duda<br>
&gt;&gt; Iowa State University<br>
&gt;&gt; Meteorology Graduate Student<br>
&gt;&gt; 3134 Agronomy Hall<br>
&gt;&gt; <a href="http://www.meteor.iastate.edu/~jdduda" target="_blank">www.meteor.iastate.edu/~jdduda</a><br>
&gt;&gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; gradsusr mailing list<br>
&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;<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><br clear="all"><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>
</div>