[ncl-talk] ncl-talk Digest, Vol 130, Issue 1

H. Joe Lee hyoklee at hdfgroup.org
Tue Sep 2 08:50:53 MDT 2014


Hi, Geeta!

 I am sorry but I misunderstood your problem.

 I think you can match dimensions utilizing "clusterNumber" variable
for "ncluster" dimension and "freezingHeightIndex" variable for
"findex" dimension. They contain index values that can be used for
subsetting "cluster" variable.


--
Save the Earth. Save Earth data in HDF-EOS. Save Big data in HDF.


On Tue, Sep 2, 2014 at 8:47 AM, H. Joe Lee <hyoklee at hdfgroup.org> wrote:
> Hi, Geeta!
>
> 2A12 product doesn't contain "graupel" variable.
>
> Please download a correct product that has grapple variable like 2B31 product:
>
> http://disc.sci.gsfc.nasa.gov/precipitation/documentation/TRMM_README/TRMM_2B31_readme.shtml
>
>
>
> By the way, a complete 2A12 version 7 NCL example is available here:
>
> http://hdfeos.org/zoo/index_openGESDISC_Examples.php#TRMM
>
> --
> Save the Earth. Save Earth data in HDF-EOS. Save Big data in HDF.
>
>
> On Tue, Sep 2, 2014 at 6:27 AM, Geeta Geeta <geetag54 at yahoo.com> wrote:
>> I want to plot the GRAUPEL WATER CONTENT at a LAT/LON with height. So Y axis
>> will be the HEIGHT and on X axis it will Graupel
>> But the document of 2A12
>> (http://pps.gsfc.nasa.gov/Documents/filespec.TRMM.V7.pdf)%20says that The
>> information ablout this variable is CONTAINED in the two variables (Pg 119)
>> 7.5 2A12 - TMI Profiling 119
>>
>>
>> S = species(1-6)
>>
>> 1 = cloud water content
>>
>> 2 = rain water content
>>
>> 3 = cloud ice content
>>
>> 4 = snow water content
>>
>> 5 = graupel water content
>>
>> 6 = latent heat
>>
>> F = freezingHeightIndex
>>
>> C = clusterNumber
>>
>> In a Fortran program,
>>
>> Profile Value = clusterScale * cluster(S,F,L,C)
>>
>> In a C program,
>>
>> Profile Value = clusterScale * cluster[C-1][L-1][F-1][S-1]
>>
>> The two variables CLUSTER SCALE
>>  float clusterScale ( nscan, npixel, nspecies )
>>          hdf_name :     clusterScale
>>
>> AND CLUSTER
>>   float cluster ( ncluster, nlayer, nfindex, nspecies )
>>          hdf_name :     cluster
>>
>>
>> are of different dimensions.
>> I donot know how to proceed further with this.
>>
>> Geeta.
>>
>>
>> On Tuesday, 2 September 2014 4:45 PM, Alessandra Giannini
>> <alesall at iri.columbia.edu> wrote:
>>
>>
>>
>> you need to continue to edit your script to adapt to your specific need…
>> What is the variable that  you want to plot? If it's not called "latentHeat"
>> what is it called?
>> Then substitute just like you did for lat and lon…
>>
>>
>>
>>
>> On Sep 2, 2014, at 7:10 AM, Geeta Geeta <geetag54 at yahoo.com> wrote:
>>
>>>    Hi
>>> I have correctec the lines 34 and 35. of my file.
>>>
>>> 1 ;Reading HDF4-SDS files
>>>      2 ;*********** Load Libraries ************************************
>>>      3 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>      4 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>      5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>      6 ;**************************************************************
>>>      7 begin
>>>      8 ;***************************************************************
>>>      9 ; User Input
>>>      10 ;***************************************************************
>>>      11                                              ; INPUT
>>>      12    diri    = "./"                            ; input directory
>>>      13    fili    = "2A12.20100702.71931.7.HDF"      ;      file
>>>      14    var    = "latentHeat"                    ; desired variable
>>>      15
>>>      16                                              ; OUTPUT
>>>      17    PLOT    = True                            ; create a plot [?]
>>>      18    pltDir  = "./"                            ; directory for plot
>>> output
>>>      19  ;pltName = "2A12.100702.71931.6"          ; plot name
>>>      20    pltName = "hdf4sds"                      ; plot name
>>>      21    pltType = "ps"                            ; ps, pdf, x11, png
>>> [5.2.0]
>>>      22 ;***************************************************************
>>>      23 ; End User Input
>>>      24 ;***************************************************************
>>>      25 ; Read hdf
>>>      26 ;***************************************************************
>>>      27
>>>      28    f      = addfile (diri+fili, "r")        ; reading inputfile
>>>      29  ;print(f)
>>>      30
>>>      31    lat  = f->Latitude(:,:)                ; (scan, pixel)
>>>      32    lon  = f->Longitude(:,:)
>>>      33
>>>      34    lat at units = "degrees_north"            ; add units
>>>      35    lon at units = "degrees_east"
>>>      36
>>>          37  printMinMax(lat, True)
>>>      38  printMinMax(lon, True)
>>>      39
>>>      40    xs    = f->$var$                      ; ( scan, pixel, layer )
>>>      41    xs at _FillValue = inttoshort( -9999 )    ; visually look at
>>> file/data
>>>      42    xs at long_name  = xs at hdf_name
>>>      43  ;printVarSummary(xs)
>>>      44  ;printMinMax(xs, True)
>>>
>>>
>>> Now when I run this
>>> I get following error
>>> bash-3.2$ ncl hdf4sds_6.ncl
>>>  Copyright (C) 1995-2007 - All Rights Reserved
>>>  University Corporation for Atmospheric Research
>>>  NCAR Command Language Version 5.0.0
>>>  The use of this software is governed by a License Agreement.
>>>  See http://www.ncl.ucar.edu/for more details.
>>
>>> (0)
>>> (0)    min=-39.1053  max=39.0529
>>> (0)
>>> (0)    min=-179.999  max=180
>>> fatal:Either file (f) isn't defined or variable (latentHeat) is not a
>>> variable in the file
>>> fatal:Execute: Error occurred at or near line 40 in file hdf4sds_6.ncl
>>>
>>> bash-3.2$
>>>
>>> I am beginner to NCL. Started a week before as I have to plot the
>>> Hydrometeor profile. Kindly guide further.
>>>
>>>
>>> Geeta.
>>>
>>>
>>> On Tuesday, 2 September 2014 4:28 PM, Geeta Geeta <geetag54 at yahoo.com>
>>> wrote:
>>>
>>>
>>> Thanks.
>>> I corrected my mistake and this is how the NCL file looks like now.
>>> ;Reading HDF4-SDS files
>>> ;*********** Load Libraries ************************************
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>> ;**************************************************************
>>> begin
>>> ;***************************************************************
>>> ; User Input
>>> ;***************************************************************
>>>                                              ; INPUT
>>>    diri    = "./"                            ; input directory
>>>    fili    = "2A12.20100702.71931.7.HDF"      ;      file
>>>    var    = "latentHeat"                    ; desired variable
>>>
>>>                                              ; OUTPUT
>>>    PLOT    = True                            ; create a plot [?]
>>>    pltDir  = "./"                            ; directory for plot output
>>>  ;pltName = "2A12.100702.71931.6"          ; plot name
>>>    pltName = "hdf4sds"                      ; plot name
>>>    pltType = "ps"                            ; ps, pdf, x11, png [5.2.0]
>>> ;***************************************************************
>>> ; End User Input
>>> ;***************************************************************
>>> ; Read hdf
>>> ;***************************************************************
>>>
>>> "hdf4sds_6.ncl" 188 lines, 6725 charac
>>>
>>> I run NCL but got following message.
>>> bash-3.2$ ncl hdf4sds_6.ncl
>>>  Copyright (C) 1995-2007 - All Rights Reserved
>>>  University Corporation for Atmospheric Research
>>>  NCAR Command Language Version 5.0.0
>>>  The use of this software is governed by a License Agreement.
>>>  See http://www.ncl.ucar.edu/for more details.
>>
>>> fatal:Either file (f) isn't defined or variable (geolocation) is not a
>>> variable in the file
>>> fatal:Execute: Error occurred at or near line 31 in file hdf4sds_6.ncl
>>>
>>> bash-3.2$
>>>
>>>
>>> I am attaching the ncl file for reference.
>>>
>>> This is the header information that I extracted using
>>> bash-3.2$ ncl_filedump 2A12.20100702.71931.7.HDF >header2010
>>>
>>> I am attaching the header2010 as well. there is NO geolocation Variable in
>>> that. So what I should DO.
>>>
>>> Geeta.
>>>
>>>
>>> On Tuesday, 2 September 2014 3:57 PM, Alessandra Giannini
>>> <alesall at iri.columbia.edu> wrote:
>>>
>>>
>>>
>>> Hi Geeta,
>>>
>>> NCL is telling you what may be wrong…
>>>
>>> Either your filename is slightly off:
>>>
>>> in your script you call a file by the name "2A12.20100702.71931.6.HDF",
>>> but when you list, you have a file with a 7 at the end, instead of a 6:
>>> "2A12.20100702.71931.7.HDF"
>>>
>>> Or there is no variable called "geolocation" in the file you are reading
>>> from.
>>>
>>> What are latitude and longitude called in the file?
>>>
>>> best, alessandra
>>>
>>>
>>>
>>>
>>>
>>> --
>>> Alessandra Giannini
>>> IRI for Climate and Society - The Earth Institute at Columbia University
>>> P.O. Box 1000, Palisades NY 10964-8000
>>> phone/fax: +1 845 680-4473/4864 - email: alesall at iri.columbia.edu
>>>
>>> On Sep 2, 2014, at 6:21 AM, Geeta Geeta <geetag54 at yahoo.com> wrote:
>>>
>>> > Hi .
>>> > I am trying to plot the data in HDF format using NCL.
>>> >
>>> > I happen to see the ready script available
>>> > https://www.ncl.ucar.edu/Applications/Scripts/hdf4sds_6.nclwhichreads HDF
>>> > data.
>>
>>> >
>>> > On my system I downloaded this script and tried to run this but when I
>>> > run this, following message comes.
>>> > bash-3.2$ ncl hdf4sds_6.ncl
>>> >  Copyright (C) 1995-2007 - All Rights Reserved
>>> >  University Corporation for Atmospheric Research
>>> >  NCAR Command Language Version 5.0.0
>>> >  The use of this software is governed by a License Agreement.
>>> >  See http://www.ncl.ucar.edu/formore details.
>>
>>> > fatal:_NclCreateFile: Requested file does not exist as
>>> > (./2A12.20100702.71931.6.HDF) or as (./2A12.20100702.71931.6), can't add
>>> > file
>>> > fatal:Either file (f) isn't defined or variable (geolocation) is not a
>>> > variable in the file
>>> > fatal:Execute: Error occurred at or near line 31 in file hdf4sds_6.ncl
>>> >
>>> > bash-3.2$ ls -l *HDF*
>>> > -rw-r--r--    1 model    staff      55092053 Sep 01 15:54
>>> > 2A12.20100702.71931.7.HDF
>>> > -rw-r--r--    1 model    staff      55092079 Aug 12 11:13
>>> > 2A12.20130523.88390.7.HDF
>>> > -rw-r--r--    1 model    staff      55092085 Aug 12 15:27
>>> > 2A12.20140602.94247.7.HDF
>>> > bash-3.2$ ls -l *hdf*ncl*
>>> > -rw-r--r--    1 model    staff          6725 Sep 01 17:01 hdf4sds_6.ncl
>>> > bash-3.2$
>>> >
>>> > The files exist in my directory.
>>> > Can somebody tell me what's mistake???
>>> > Geeta.
>>> >
>>> >
>>> >
>>> >
>>> > _______________________________________________
>>> > ncl-talk mailing list
>>> > List instructions, subscriber options, unsubscribe:
>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>>
>>>
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>


More information about the ncl-talk mailing list