Ulrike,<br><br>The function cdiff works as a horizontal plane differental, not vertical, for qhich it needs equal space grid points. Also when you use it probably you should need to convert the x and y direction into spherical&nbsp; differential arcs&nbsp; as:<br>
<br>&#39;define dy = cdiff(lat*1000.0,y) * 3.1416/180&#39;<br>&#39;define dx = cdiff(lon*1000.0,x) * 3.1416/180&#39;<br>&#39;define dvx = cdiff(vinterp,x)&#39;<br>

 &nbsp;&#39;define duy = cdiff(uinterp,y)&#39;<br>&#39;define zeta = ((dvx / (cos(lat*3.1416/180) * dx)) -<br>(duay / dy))/(6.37e6)&#39;<br>&#39;define zeta = &nbsp;(dvx/dx-duy/dy)&#39;<br>
<br>for instance.<br><br><br><br><br><div class="gmail_quote">2008/6/23 Ulrike Wissmeier &lt;<a href="mailto:ulrike@meteo.physik.uni-muenchen.de">ulrike@meteo.physik.uni-muenchen.de</a>&gt;:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear all,<br>
in order to compute/plot the vertical gradient of vertical vorticity zeta, I<br>
couldn&#39;t use the cdiff-function, as there seem to be problems with<br>
cdiff(Expr,z). In the past (see archive) it was suggested to use "define<br>
dExpr=Expr(z+1)-Expr(z-1)" instead of cdiff, however, that does not work for<br>
my calculated vertical vorticity zeta either:<br>
<br>
grads-script:<br>
 &nbsp;&#39;set lon 0 60&#39;<br>
 &nbsp;&#39;set lat 0 60&#39;<br>
 &nbsp;&#39;set z 0 30&#39;&#39;<br>
<br>
 &nbsp;&#39;define dx = cdiff(lon*1000.0,x)&#39;<br>
 &nbsp;&#39;define dy = cdiff(lat*1000.0,y)&#39;<br>
 &nbsp;&#39;define dvx = cdiff(vinterp,x)&#39;<br>
 &nbsp;&#39;define duy = cdiff(uinterp,y)&#39;<br>
 &nbsp;&#39;define zeta = &nbsp;(dvx/dx-duy/dy)&#39;<br>
<br>
 &nbsp;&#39;set z 1 29&#39;<br>
<br>
 &nbsp;&#39;define dzetaz = zeta(z+1)-zeta(z-1)&#39;<br>
 &nbsp;&#39;define dzhz = zh(z+1)-zh(z-1)&#39;<br>
<br>
 &nbsp;&#39;set lon &#39;lon&#39;&#39;<br>
 &nbsp;&#39;display dzetaz/dzhz&#39;<br>
<br>
(where zh is the position of the vertical grid points z=1,2,3,...)<br>
<br>
This gives the following error message:<br>
 &nbsp;Syntax Error: &nbsp;Invalid Operand<br>
 &nbsp;&#39;dzetaz&#39; not a variable or function name<br>
 &nbsp;Error ocurred at column 1<br>
 &nbsp;DISPLAY error: &nbsp;Invalid expression<br>
 &nbsp;Expression = dzetaz/dzhz<br>
<br>
I would be very glad if anyone can help,<br>
thanks in advance,<br>
<font color="#888888">Ulrike<br>
</font></blockquote></div><br>