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

[mmkamal at uwaterloo.ca]

Hi Dannis,


Thanks for the code and the information that "new NASA based projects  
will be using this format".


Thanks
Kamal




Quoting Dennis Shea <shea at ucar.edu>:

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