[gradsusr] gradsusr Digest, Vol 13, Issue 49

Zilore Mumba zmumba at gmail.com
Fri Mar 25 02:45:57 EDT 2011


Definition of lat/lon in a file.
With Excuses on people's time with a very basic problem I have, I
would like to request once more if someone can spare some of their
time to look at my definition of lat/lon in the ncl scripts attached.
1. I have binary data (not GRIB, but written into binary by fortran)
converted through Grads (Lats4d) into NetCdf. Both the original binary
file and the NetCdf versions work in Grads but not in ncl. I assumed
it did not work because the file was binary (not GRIB).
2. I took model output data in GRIB and made a ctl with grib2ctl and
converted to NetCdf in Grads. It does not work in ncl.
3. I downloaded 850hPa temperature data (in NetCdf format) from the
NCEP reanalysis page. It does not work in ncl, but works in Grads.

In all the three cases above the error is the same, as below
**************
 Copyright (C) 1995-2010 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 5.2.1
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
fatal:Either file (f1) isn't defined or variable (Air) is not a
variable in the file

fatal:Execute: Error occurred at or near line 20 in file native_2.ncl
***********************************
Here is how I introduced the file
!!!!!!!!!
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;*****************************************************
begin
 f1 = addfile("AirTempFeb11.nc","r")
 Air=f1->Air

 lat = f1->lat
 lon = f1->lon
 lat_len = dimsizes(lat)
 lon_len = dimsizes(lon)
; time = f1-> time
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

I would appreciate assistance, where am I making the mistake?
Appreciation
Zilore


On Mon, Mar 21, 2011 at 6:49 PM,  <gradsusr-request at gradsusr.org> wrote:
> Send gradsusr mailing list submissions to
>        gradsusr at gradsusr.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
>        http://gradsusr.org/mailman/listinfo/gradsusr
> or, via email, send a message with subject or body 'help' to
>        gradsusr-request at gradsusr.org
>
> You can reach the person managing the list at
>        gradsusr-owner at gradsusr.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gradsusr digest..."
>
>
> Today's Topics:
>
>   1. WRF output (Dr.Gamal El Afandi)
>   2. Re: WRF output (hachelaf rabah)
>   3. Regards Recipe-020: Spaghetti and Box-and-whisker plots   and
>      sorting in Grads (Becker, Bernd)
>   4. Re: Regards Recipe-020: Spaghetti and Box-and-whisker     plots
>      and sorting in Grads (Jennifer Adams)
>   5. Re: curvilinear data (Jennifer Adams)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 18 Mar 2011 13:11:03 -0700 (PDT)
> From: "Dr.Gamal El Afandi" <gamalafandy at yahoo.com>
> Subject: [gradsusr] WRF output
> To: gradsusr at gradsusr.org
> Message-ID: <144313.75860.qm at web113403.mail.gq1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi All,
> Could you please let me know if anyone has some scripts to visualize wrf output.
> ?
> With my kind regards
> Gamal
> _____________________________________________________
> Dr.Gamal Salah El Afandi,Postdoctoral Fellow.
> Data to Insight Center?of Pervasive Technology Institute
> Indiana University Bloomington USA
> 2719 E 10th St.
> Innovation Center
> Room 130C.?
> Bloomington, IN ?47408,USA
> Phone:+1 812-856-1953
> Fax:?+1 812-856-1735
> ?
> http://www.indiana.edu.
> http://pervasivetechnologylabs.iu.edu/d2i/people
> ?
> Al Azhar?University
> Faculty of science?at Cairo?
> Department of Astronomy and Meteorology
> Nasr City,Postal Code 11884
> Cairo,Egypt
> http://www.uazhar.edu.eg/astmet/Gamal.htm
> http://www.itap.purdue.edu/pto/NATO_KEEO/partCoDir_en.html
> http://www.ucar.edu/governance/iap/iapreps.shtml
> Tel. 0020238701259 Home???Tel. 0020222618163 Work
> Tel. 0020103770593 Cell??????Fax 0020222629356 Work
> ________________________________________________________
>
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110318/06a13ef9/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Fri, 18 Mar 2011 22:34:48 +0100
> From: hachelaf rabah <r.hachelaf at gmail.com>
> Subject: Re: [gradsusr] WRF output
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
>        <AANLkTinsONf81SQ9=FdXLUaWu+WepBXXBieyrf6OAPO5 at mail.gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> Hello Dr.Gamal
>
> you can use ARWpost to generat .ctl and .data files
> or with WPPV3 to convert netcdf WRF output to grib format and plot it
> also with GRADS.
> So what parameters you want to visualize.
>
>
> 2011/3/18, Dr.Gamal El Afandi <gamalafandy at yahoo.com>:
>> Hi All,
>> Could you please let me know if anyone has some scripts to visualize wrf
>> output.
>>
>> With my kind regards
>> Gamal
>> _____________________________________________________
>> Dr.Gamal Salah El Afandi,Postdoctoral Fellow.
>> Data to Insight Center of Pervasive Technology Institute
>> Indiana University Bloomington USA
>> 2719 E 10th St.
>> Innovation Center
>> Room 130C.
>> Bloomington, IN  47408,USA
>> Phone:+1 812-856-1953
>> Fax: +1 812-856-1735
>>
>> http://www.indiana.edu.
>> http://pervasivetechnologylabs.iu.edu/d2i/people
>>
>> Al Azhar University
>> Faculty of science at Cairo
>> Department of Astronomy and Meteorology
>> Nasr City,Postal Code 11884
>> Cairo,Egypt
>> http://www.uazhar.edu.eg/astmet/Gamal.htm
>> http://www.itap.purdue.edu/pto/NATO_KEEO/partCoDir_en.html
>> http://www.ucar.edu/governance/iap/iapreps.shtml
>> Tel. 0020238701259 Home   Tel. 0020222618163 Work
>> Tel. 0020103770593 Cell      Fax 0020222629356 Work
>> ________________________________________________________
>>
>>
>>
>>
>
>
> --
> Cordialement,
> Best regards,
> Rabah Hachelaf
>
>     ______
>    (          )
>   (            )
>  (________)
>     /////////
>
>
> ------------------------------
>
> Message: 3
> Date: Mon, 21 Mar 2011 15:45:13 -0000
> From: "Becker, Bernd" <bernd.becker at metoffice.gov.uk>
> Subject: [gradsusr] Regards Recipe-020: Spaghetti and Box-and-whisker
>        plots   and sorting in Grads
> To: "GrADS Users Forum" <gradsusr at gradsusr.org>
> Message-ID:
>        <E209E2066317644FB35E5AA0E8C1FB7F29EF40 at EXXMAIL02.desktop.frd.metoffice.com>
>
> Content-Type: text/plain;       charset="us-ascii"
>
>
>
> Hello,
>
> what have we available in Grads that can sort by size?
>
> In recipe-020, box and whisker plot ranges are calculated by
> fraction of standard deviation.
>
> Assume, one were to plot quantiles after ranking:
>
> how would one achieve the ranking?
>
> of a stack of 2-d fields?
>
> This can be a compute intensive and time consuming affair.
> But it would certainly be nice to have.
>
> (in metview, macro "distribution" does the trick)
>
> many thanks for your suggestions.
>
> Bernd.
>
>
> --
> Bernd Becker   Climate Impacts Product Development
> Met Office Hadley Centre    FitzRoy Road  Exeter   Devon EX1 3PB  United
> Kingdom
> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050
> E-mail:bernd.becker at metoffice.gov.uk -
> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/
> bernd-becker
>
>
>
>
> ------------------------------
>
> Message: 4
> Date: Mon, 21 Mar 2011 12:15:19 -0400
> From: Jennifer Adams <jma at cola.iges.org>
> Subject: Re: [gradsusr] Regards Recipe-020: Spaghetti and
>        Box-and-whisker plots and sorting in Grads
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <17CAFB17-2803-4B28-A343-CE9613A06D44 at cola.iges.org>
> Content-Type: text/plain; charset="us-ascii"
>
> There's no way to do sorting in GrADS right now, but it's high on our
> list of features we'd like to implement. It will fall under the
> category of 'operations on define objects' and will require a fair bit
> of available memory, depending on how much data you want to sort.
> --Jennifer
>
>
> On Mar 21, 2011, at 11:45 AM, Becker, Bernd wrote:
>
>>
>>
>> Hello,
>>
>> what have we available in Grads that can sort by size?
>>
>> In recipe-020, box and whisker plot ranges are calculated by
>> fraction of standard deviation.
>>
>> Assume, one were to plot quantiles after ranking:
>>
>> how would one achieve the ranking?
>>
>> of a stack of 2-d fields?
>>
>> This can be a compute intensive and time consuming affair.
>> But it would certainly be nice to have.
>>
>> (in metview, macro "distribution" does the trick)
>>
>> many thanks for your suggestions.
>>
>> Bernd.
>>
>>
>> --
>> Bernd Becker   Climate Impacts Product Development
>> Met Office Hadley Centre    FitzRoy Road  Exeter   Devon EX1 3PB
>> United
>> Kingdom
>> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050
>> E-mail:bernd.becker at metoffice.gov.uk -
>> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/
>> bernd-becker
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> --
> Jennifer M. Adams
> IGES/COLA
> 4041 Powder Mill Road, Suite 302
> Calverton, MD 20705
> jma at cola.iges.org
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110321/6d2f37e9/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Mon, 21 Mar 2011 17:01:05 -0400
> From: Jennifer Adams <jma at cola.iges.org>
> Subject: Re: [gradsusr] curvilinear data
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <6D6ACA6F-1D68-48C8-9662-94783407FAE0 at cola.iges.org>
> Content-Type: text/plain; charset="us-ascii"
>
> It is possible to use PDEF to map curvilinear coordinates to a regular
> lat/lon grid, but I have never tried, and I know of nobody who has
> already done this with GrADS. It requires some knowledge of how to
> convert any lat/lon value into real-valued i,j coordinates in the
> curvilinear grid. It has always been easier to do this interpolation
> outside of GrADS.
> --Jennifer
>
>
> On Mar 18, 2011, at 2:00 PM, James T. Potemra wrote:
>
>> As far as I know, GrADS will only be able to handle curvilinear grids
>> that are curvilinear in one dimension, i.e., lat(i)/lon(i), not
>> lat(i,j)/lon(i,j). Of course you can always make the control file
>> whatever you want, but the display will not be correct.
>>
>> Jim
>>
>> Edward wrote:
>>>
>>> hi, gradsuser,
>>>
>>> Does anybody has an idea to deal with the curvilinear data in grads?
>>> Do I need list all long/lat in ctl file? or  interpolate outside
>>> grads
>>> before grads?
>>>
>>> Thanks for any help.
>>>
>>>
>>> Dawn
>>> ------------------------------------------------------------------------
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> --
> Jennifer M. Adams
> IGES/COLA
> 4041 Powder Mill Road, Suite 302
> Calverton, MD 20705
> jma at cola.iges.org
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110321/3bf686b6/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 13, Issue 49
> ****************************************
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: lat_time_1.ncl
Type: application/octet-stream
Size: 2684 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20110325/9e7fed50/attachment-0006.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: native_2.ncl
Type: application/octet-stream
Size: 2711 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20110325/9e7fed50/attachment-0007.obj 


More information about the gradsusr mailing list