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

Geeta Geeta geetag54 at yahoo.com
Tue Sep 2 05:10:32 MDT 2014


   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.nclwhich reads 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/for more 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140902/4127efc6/attachment-0001.html 


More information about the ncl-talk mailing list