<div>Thank you José and Heiner, I will try both. Then I will tell you!!</div>
<div>Eduardo.<br><br>&nbsp;</div>
<div><span class="gmail_quote">2006/11/21, José Paulo Bonatti &lt;<a href="mailto:bonatti@cptec.inpe.br">bonatti@cptec.inpe.br</a>&gt;:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi,<br><br>May be this work:<br><br>* k&nbsp;&nbsp; is&nbsp;&nbsp; the vertical index, km is the maximum<br>k=1<br>while (k &lt;= km)
<br>k1=k-1<br>k2=k+1<br>if (k = 1)<br>&nbsp;&nbsp; k1=1<br>endif<br>if (k = km)<br>&nbsp;&nbsp; k2=km<br>endif<br>* Getting pressure in Pascal (if in zdef&nbsp;&nbsp;the&nbsp;&nbsp;pressure is given&nbsp;&nbsp;in&nbsp;&nbsp;mb<br>= hPa)<br>'set z 'k1<br>'define p1=lev*100'<br>'set z 'k2
<br>'define p2=lev*100'<br>'set z 'k<br>* Vertical divergence in pressure coordinates considering omeg is the<br>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; vertical pressure velocity in Pa/seg.<br>*&nbsp;&nbsp;vertdiv1, vertdiv2, ... will be the vertical divergence for levels 1,
<br>2, ...<br>'define verdiv'k'=(omeg(z='k2')-omeg(z='k1'))/(p2-p1))'<br>k=k+1<br>endwhile<br>*<br>Bonatti<br><br>Eduardo Agosta Scarel wrote:<br><br>&gt; Hi you all!!<br>&gt;<br>&gt; Does anybody how to estimate the vertical divergence?
<br>&gt; I am trying to do so with the cdiff command but with no success.<br>&gt; Thanks for any help!!<br>&gt; Eduardo<br><br><br>--<br>*****************************************************************<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Jose Paulo Bonatti
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Dr. em Meteorologia - Pesquisador Titular<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Instituto Nacional de Pesquisas Espaciais - INPE<br>&nbsp;&nbsp; Centro de Previsao de Tempo e Estudos Climaticos - CPTEC<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Divisao de Modelagem e Desenvolvimento - DMD
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (<a href="http://www.cptec.inpe.br/">http://www.cptec.inpe.br/</a>)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Rodovia Presidente Dutra, km 40 SPRJ<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;12630-000 Cachoeira Paulista - SP - Brasil<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Tel: (+55) (12) 
3186.8421&nbsp;&nbsp;Fax: (+55) (12) 3101.2835<br>*****************************************************************<br></blockquote></div>