[ncl-talk] How to get Monthly climatology of Level 2 TRMM swath data
Dennis Shea
shea at ucar.edu
Thu Apr 30 10:14:50 MDT 2020
*nmos = 12*
gbin = new ( (/*nmos,*nlat,mlon/), float )
gknt = new ( (/*nmos,*nlat,mlon/), integer)
Loop over each file and sum
nmo = ... ; determine the month [nmo=0,11] for the current file
* bin_sum*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/bin_sum.shtml>
(gbin(nmo,:,:),gknt(nmo,:,:),lon,lat,LON1D,LAT1D,ndtooned(x))
Later:
gknt = *where*(gknt.eq.0 , gknt at _FillValue, gknt)
gbin = gbin/gknt
=====================================
On Thu, Apr 30, 2020 at 5:06 AM Lyndz via ncl-talk <ncl-talk at ucar.edu>
wrote:
> Dear NCL-experts,
>
> I have a script for plotting the "DAILY" TRMM 2A25 swath data, which is
> attached to this email.
>
> I would like to get the MONTHLY CLIMATOLOGY of the files and plot them. I
> am not sure as to how I can modify the following script to generate these
> plots.
>
> I have TRMM 2A25 HDF files from 1998 to 2015. All located in the same
> directory.
>
>
> I'll appreciate any help on this matter.
>
> Sincerely,
> Lyndz
> _______________________________________________
> 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/20200430/84bcc361/attachment.html>
More information about the ncl-talk
mailing list