[ncl-talk] Request to replicate example script "gpm_1.ncl" for NetCDF data

[Dennis Shea]

[0] You will have to deal with HDF5 sooner or later.
    I think all new NASA based projects will be using this format.

[1] I do not think NCL people will be going to GIOVANNI to download files.
               *Our time is valuable too!!*
    If anything, this would be your task.

[2] You must use NCL 6.3.0 to use the GPM HDF5 data.
    This NCL version can read the GPM HDF5 file with groups.
    It is not possible to make it work with 6.2.*

[3] The change is trivial for an NCL user.  The HDF5 read section is
identical to what you would do for netCDF.

    fn = "3B-MO.MS.MRG.3IMERG.20140701-S000000-E235959.07.V03D.HDF5"
    f  = addfile(fn. "r")
;----------------------------------------------------------------------
; Access variable in group "/Grid"
; Note 1: Group access is via => syntax
; Note 2: The dimension order is (lon,lat) rather than the far more common
(lat,lon)
; Note 3: Currently, 6.3.0 does not support subsetting of group variables
upon input.
;         p = grp->$varName$({lonL:lonR},{latS:latN})   currently *not*
supported
;----------------------------------------------------------------------
    grp = f=>/Grid
    p   = grp->$varName$
    printVarSummary(p)        ; p(lon,lat)

    ---
    latS =   20.0
    latN =   50.0
    lonL = -110.0             ; lon  -180 to 180
    lonR =  -50.0
                              ; extract & reorder
    pregion = p({lat|latS:latN},{lon|lonL:lonR})   ; make regional and
(lat,lon)


On Tue, Apr 7, 2015 at 2:10 PM, <mmkamal at uwaterloo.ca> wrote:

> Hi,
>
> I was wondering whether it is possible to replicate the example script
> (https://www.ncl.ucar.edu/Applications/Scripts/gpm_1.ncl) for GPM
> NetCDF data (GIOVANNI offer to generate NetCDF GPM data set). It would
> be a great help for people like me, who are more familiar with
> manipulating NetCDF data using NCL.
>
>
> Two more subsidiary request:
>
> 1) Is it possible to create the example script that can be executed
> using at least NCL 2.6.2 for a wider audience.
>
> 2) Please write the script with a subset of data set instead of global
> map (as most of us use this product for RCM validation)
>
> Here is the link to generate NetCDF GPM data set.
>
>
>
> http://giovanni.gsfc.nasa.gov/giovanni/#service=ACCUMULATE&starttime=2014-07-01T04:00:00Z&endtime=2014-07-31T23:59:59Z&bbox=-97.1191,33.9653,-62.9473,52.106&data=GPM_3IMERGM_03_precipitation&variableFacets=dataFieldMeasurement%3APrecipitation%3BdataProductPlatformInstrument%3AGPM%3B
>
>
>
>
> Thanks
> Kamal
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> 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/20150407/7e130375/attachment.html