gelp for divergence

[Eduardo Agosta Scarel]

Thank you José and Heiner, I will try both. Then I will tell you!!
Eduardo.


2006/11/21, José Paulo Bonatti <bonatti at cptec.inpe.br>:
>
> Hi,
>
> May be this work:
>
> * k   is   the vertical index, km is the maximum
> k=1
> while (k <= km)
> k1=k-1
> k2=k+1
> if (k = 1)
>    k1=1
> endif
> if (k = km)
>    k2=km
> endif
> * Getting pressure in Pascal (if in zdef  the  pressure is given  in  mb
> = hPa)
> 'set z 'k1
> 'define p1=lev*100'
> 'set z 'k2
> 'define p2=lev*100'
> 'set z 'k
> * Vertical divergence in pressure coordinates considering omeg is the
> *               vertical pressure velocity in Pa/seg.
> *  vertdiv1, vertdiv2, ... will be the vertical divergence for levels 1,
> 2, ...
> 'define verdiv'k'=(omeg(z='k2')-omeg(z='k1'))/(p2-p1))'
> k=k+1
> endwhile
> *
> Bonatti
>
> Eduardo Agosta Scarel wrote:
>
> > Hi you all!!
> >
> > Does anybody how to estimate the vertical divergence?
> > I am trying to do so with the cdiff command but with no success.
> > Thanks for any help!!
> > Eduardo
>
>
> --
> *****************************************************************
>                     Jose Paulo Bonatti
>           Dr. em Meteorologia - Pesquisador Titular
>        Instituto Nacional de Pesquisas Espaciais - INPE
>    Centro de Previsao de Tempo e Estudos Climaticos - CPTEC
>          Divisao de Modelagem e Desenvolvimento - DMD
>                  (http://www.cptec.inpe.br/)
>              Rodovia Presidente Dutra, km 40 SPRJ
>           12630-000 Cachoeira Paulista - SP - Brasil
>       Tel: (+55) (12) 3186.8421  Fax: (+55) (12) 3101.2835
> *****************************************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20061121/97dc8a87/attachment.html