gelp for divergence

[José Paulo Bonatti]

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
*****************************************************************