Ulrike,<br><br>Also I see you may need to estimate the vertical gradient in a sistematic way and using a density wighted expression. I give you below the following script that you should adjust for your case, it may help (words and in Spanish, if needed, I can translate them later).<br>
That`s all.<br><br>Eduardo.<br><br><br><br><br>* k&nbsp;&nbsp; es el indice vertical, km es el maximo<br>km=12<br>k=1<br>while (k &lt;= km)<br>&nbsp;k1=k-1<br>&nbsp;k2=k+1<br>&nbsp;if (k = 1)<br>&nbsp;&nbsp; k1=1<br>&nbsp;endif<br>&nbsp;if (k = km)<br>&nbsp;&nbsp; k2=km<br>&nbsp;endif<br>
* los niveles p1 y p2 estan en hPa,<br>&nbsp;&#39;set z &#39;k1<br>&nbsp;&#39;define p1=lev&#39;<br>&nbsp;&#39;set z &#39;k2<br>&nbsp;&#39;define p2=lev&#39;<br>&nbsp;&#39;set z &#39;k<br><br>* Divergencia vertical en coordenadas de presión considerando omeg <br>
*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; la variable vertical.<br>*&nbsp; vertdiv1, vertdiv2, ... sera la divergencia vertical para los niveles 1,<br>2, ...<br>* en este caso la variable vertical está pesada por la densidad del aire<br>* de la manera siguiente<br>
* 1/(densidad)*Div/vertical(densidad*omeg)<br><br>* densidad del aire en funcion de alturas de nivel para amtosfera standard<br>&#39;roo=1.212&#39;<br>&#39;cte=0.9999&#39;<br>&#39;roz&#39;k1&#39;=roo*exp(-cte*p1)&#39;<br>
&#39;roz&#39;k2&#39;=roo*exp(-cte*p2)&#39;<br><br>&#39;define terma=(roz&#39;k2&#39;-roz&#39;k1&#39;)&#39;<br>&#39;define termb=(Euz(z=&#39;k2&#39;)+Euz(z=&#39;k1&#39;))/2&#39;<br>&#39;define termc=termb/terma&#39;<br><br>
&#39;p0=1000&#39;<br>&#39;C=p0*exp(log(p2/p0)-log(p1/p0)&#39;<br><br>&#39;define omeg(z=&#39;k2&#39;)=Euz(z=&#39;k2&#39;)&#39;<br>&#39;define omeg(z=&#39;k1&#39;)=Euz(z=&#39;k1&#39;)&#39;<br>&nbsp;<br>&nbsp;<br>&nbsp;&#39;define verdiv&#39;k&#39;=(termc+((omeg(z=&#39;k2&#39;)-omeg(z=&#39;k1&#39;))/(p2-p1))))*C&#39;<br>
&nbsp;k=k+1<br>endwhile<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>