[ncl-talk] Area average of a variable from WRF output

Bill Ladwig ladwig at ucar.edu
Tue Dec 27 10:41:49 MST 2016


Hi Jose,

The function you are using (wgt_areaave) is for rectilinear coordinates.
WRF uses curvilinear coordinates (latitude and longitude are
two-dimensional arrays), so you need to use the wgt_areaave2 function (see
http://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_areaave2.shtml).

Also, for the weights, there is no GWGT variable in a WRF file, so you'll
need to come up with your own weights.  Dennis has a suggestion from this
archive that might be a good starting point:

https://www.ncl.ucar.edu/Support/talk_archives/2011/2789.html

Hope this helps,

Bill


On Thu, Dec 15, 2016 at 4:50 PM, Jose Augusto Paixão Veiga <
veiga.uea at gmail.com> wrote:

> Dear all,
>
> I would appreciate if someone could share a NCL script example that
> calculates the area average of precipitation from a WRF domain (or a
> specific area average inside a domain).
>
> I've used an example from http://www.ncl.ucar.edu/
> Document/Functions/Built-in/wgt_areaave.shtml website. However, I've
> reached no success. After that, I "googled" for help in tis issue but with
> no exit again.
>
> The script I've used as an initial test is attached. The error is
> presented bellow:
>
> fatal:["Execute.c":5773]:variable (GWGT) is not in file (f)
>
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 20 in
> file media-area2.ncl
>
> Any help will be very appreciated.
>
> Thanks
>
> PS: My objective is to plot area averaged precipitation total versus time
> integration.
>
>
>
> ======================================================
> ​*Prof. José Augusto Paixão Veiga (CREA n.: 24161)*
>
> *Coordenador de Projetos Institucionais*
> Universidade do Estado do Amazonas
> ​(UEA)​
>
> Pró-Reitoria de Pesquisa e Pós-Graduação - PROPESP​
> Av. Djalma Batista, 3578, Manaus-AM Brasil
> (92)
> ​3646-7297
>
> Departamento de Meteorologia
> ​ -​
>
> Escola Superior de Tecnologia (EST)
> Av. Darcy Vargas, 1200, Manaus-AM Brasil
> (92) ​
> 3878
> ​-​
> 4317
> =====================================================
> Skype: veiga_j.a.p.
> *CV*: http://lattes.cnpq.br/4027612512091565
> *LinkedIn*: https://www.linkedin.com/in/josé-augusto-veiga-a9a28130
> *Site*: https://scientificmet.wordpress.com
> =====================================================
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161227/1f74f7e3/attachment.html 


More information about the ncl-talk mailing list