[gradsusr] gradsusr Digest, Vol 24, Issue 37
Zilore Mumba
zmumba at gmail.com
Mon Feb 27 10:03:44 EST 2012
Thank you very much to Jayakrishnan. I did install cdo and tried it on
a completely different file and it works perfect. Thanks at least I
have something new to use.
However it looks like the WRF output files are written in a different
Netcdf format. Below is the error I get with cdo.
Thank you very much for your time and for something new I have learnt
Zilore
[rma at localhost run]$ cdo -f grb copy wrfout_d01_2012-01-27_00:00:00 output.grb
Warning (cdfSetVar) : Inconsistent variable definition for XLAT_U!
Warning (cdfSetVar) : Inconsistent variable definition for XLONG_U!
Warning (cdfSetVar) : Inconsistent variable definition for XLAT_V!
Warning (cdfSetVar) : Inconsistent variable definition for XLONG_V!
Warning (cdfInqContents) : time variable not found!
Segmentation fault (core dumped)
On Mon, Feb 27, 2012 at 9:27 AM, <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. Clip data with a shapefile (Youssef Bennouna)
> 2. Re: Interpolation (Alexandre Santos)
> 3. Grads 2.0.1 Compilation Problem (Spencer Ogden)
> 4. Meteorological field in circular area
> (Rodrigo de Souza Barreto Mathias)
> 5. Re: Meteorological field in circular area (Arlindo da Silva)
> 6. netcdf to grib (Zilore Mumba)
> 7. Re: netcdf to grib (Jayakrishnan PR)
> 8. Reg: how to convert station data to gridded data
> (malleswararao maguluri)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 24 Feb 2012 09:35:20 -0000
> From: "Youssef Bennouna" <ybennouna at gmail.com>
> Subject: [gradsusr] Clip data with a shapefile
> To: <gradsusr at gradsusr.org>
> Message-ID: <001401ccf2d7$a147f600$e3d7e200$@com>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi all grads users,
>
> I would like to clip data (in grads format) by a polygon (shapefile). To be
> more precise, I want to have 0's (or no data) outside an area defined by a
> shapefile and original values inside the area, this is in the purpose of
> doing some statistics. Is there any easy way to do this ?
>
> Thank you for your help.
>
>
>
> Youssef BENNOUNA
>
> ETUDES et MESURES les 5 DOMAINES
>
> www.les5domaines.com
>
> 80, Rue Oued Ziz - AGDAL, RABAT
>
> Tel/Fax : +212 (0) 537 68 03 60
>
> GSM : +212 (0) 661 10 97 43
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120224/80bd6d85/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Fri, 24 Feb 2012 09:52:35 -0300
> From: Alexandre Santos <puffy78 at gmail.com>
> Subject: Re: [gradsusr] Interpolation
> To: gradsusr at gradsusr.org
> Message-ID:
> <CAMh9esWcOhw-Vb_kEXAQWxWnh61w9FkR+cJb5gJV5z_fGik9tQ at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> ello,
>>
>> Have a file at format text and like of the open in grads with your ctl.
>> How do this program in fortran for interpolation of the field precipitation?
>>
>> file:
>>
>> lat lon precipitation
>> -5.5 -35.6 10.00 mm
>> -5.4 -37.3 50.3 mm
>> -6.2 -34.3 67.0 mmm
>> .....
>> .....
>> ....
>>
>> thanks
>>
>> Alexandre
>>
>>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120224/ddab4003/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Fri, 24 Feb 2012 12:51:14 -0500
> From: Spencer Ogden <spencer at spencerogden.com>
> Subject: [gradsusr] Grads 2.0.1 Compilation Problem
> To: gradsusr at gradsusr.org
> Message-ID: <4F47CE12.3020009 at spencerogden.com>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> I'm trying to get into grads for the analysis of ocean wind patterns.
> I've worked through a few compilation problems, but have hit a road
> block and I wasn't able to turn up any similar problems on Google.
>
> I'm running on Gentoo and I have manually installed libpng 1.2.34 and
> g2clib 1.2.3. This has allowed configure to enable grib2 and netcdf. I'm
> running in to a problem with the compile though. When compiling gaio.c,
> NC_CHUNKED throws an error as an undeclared symbol. Seems like this
> could be a case of having the wrong version of a header, but I haven't
> been able to track down where I should start looking.
>
> Thanks for any advice,
>
> Spencer
>
> Here's the last bit of ./configure:
> +=========================================================================+
> GrADS 2.0.1 : Built Fri Feb 24 11:52:02 EST 2012 for i686-pc-linux-gnu
>
> +=========================================================================+
>
> +-------------------------------+
> | |
> | GrADS will be built with: |
> | |
> | - GUI disabled |
> | + readline enabled |
> | + geotiff enabled |
> | - shapefile disabled |
> | + printim enabled |
> | + grib2 enabled |
> | - hdf4 disabled |
> | + hdf5 enabled |
> | + netcdf-4 enabled |
> | - OPeNDAP disabled |
> | |
> +-------------------------------+
>
> `make` ouput:
> gcc -DHAVE_CONFIG_H -I. -I. -I. -I/usr/include -I/usr/include
> -I/usr/include -I/usr/include/hdf -g -O2 -rdynamic -c gaio.c
> gaio.c: In function 'garead':
> gaio.c:753:3: warning: format '%ld' expects type 'long int', but
> argument 4 has type 'size_t'
> gaio.c:753:3: warning: format '%ld' expects type 'long int', but
> argument 4 has type 'size_t'
> gaio.c: In function 'gancsetup':
> gaio.c:2446:17: error: 'NC_CHUNKED' undeclared (first use in this function)
> gaio.c:2446:17: note: each undeclared identifier is reported only once
> for each function it appears in
> gaio.c:2462:3: warning: format '%ld' expects type 'long int', but
> argument 4 has type 'size_t'
> gaio.c:2462:3: warning: format '%ld' expects type 'long int', but
> argument 4 has type 'size_t'
> make[2]: *** [gaio.o] Error 1
> make[2]: Leaving directory `/home/spencer/grads-2.0.1/src'
> make[1]: *** [all] Error 2
> make[1]: Leaving directory `/home/spencer/grads-2.0.1/src'
> make: *** [all-recursive] Error 1
>
>
>
>
>
> ------------------------------
>
> Message: 4
> Date: Fri, 24 Feb 2012 21:05:32 -0800 (PST)
> From: Rodrigo de Souza Barreto Mathias <rodrigowind at yahoo.com.br>
> Subject: [gradsusr] Meteorological field in circular area
> To: "gradsusr at gradsusr.org" <gradsusr at gradsusr.org>
> Message-ID:
> <1330146332.86104.YahooMailNeo at web161505.mail.bf1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Good morning!
> I wanna calculate a meteorological field in a circular area with radius R starting ?at a point (lat, lon). For example: average temperature in a circular area with 200 km of radius from a point 20N/020W. I don't know a specific function in Grads to do what I want. Could you help me with some function or a script using a mask?
>
> Thanks a lot!
> Rodrigo Mathias
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120224/0ff260c7/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Sat, 25 Feb 2012 10:27:30 -0500
> From: Arlindo da Silva <dasilva at alum.mit.edu>
> Subject: Re: [gradsusr] Meteorological field in circular area
> To: Rodrigo de Souza Barreto Mathias <rodrigowind at yahoo.com.br>, GrADS
> Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <CALz7xZeQeTSsbo6=y4qq9_94fjSvcwYGuLPZA=Vfxn9AOkoJ6Q at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> On Sat, Feb 25, 2012 at 12:05 AM, Rodrigo de Souza Barreto Mathias <
> rodrigowind at yahoo.com.br> wrote:
>
>> Good morning!
>> I wanna calculate a meteorological field in a circular area with radius R
>> starting at a point (lat, lon). For example: average temperature in a
>> circular area with 200 km of radius from a point 20N/020W. I don't know a
>> specific function in Grads to do what I want. Could you help me with some
>> function or a script using a mask?
>>
>>
> The attached script computes a circular mask with distance R from a
> reference (lon,lat). For your point (20W,20N) and R=200 do this:
>
> ga-> circ_mask -20 20 200
>
> This will define a variable called "mask" with "1' inside the circle, undef
> elsewhere. To get the average temperature inside the circle:
>
> ga-> d aave(mask*temperature,global)
>
> Good Luck,
>
> Arlindo
>
>
> --
> Arlindo da Silva
> dasilva at alum.mit.edu
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120225/28361791/attachment-0001.html
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: circ_mask.gs
> Type: application/octet-stream
> Size: 788 bytes
> Desc: not available
> Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120225/28361791/attachment-0001.obj
>
> ------------------------------
>
> Message: 6
> Date: Mon, 27 Feb 2012 07:55:41 +0200
> From: Zilore Mumba <zmumba at gmail.com>
> Subject: [gradsusr] netcdf to grib
> To: gradsusr at gradsusr.org
> Message-ID:
> <CAKS2qbL6k6LpevD_JHU52S5R9CtOs5XUKcpS7JUxWWMFNXXW0g at mail.gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> Excuses for posting what may be a very simple issue or an issue
> already tackled. I have searched and have not found the answer. I have
> WRf output in NetCdf format. WRF have a software ARWpost to convert
> the files for working in grads. This software is too heavy for my
> system.
> With sdfopen I get the error below.
> Is there a solution to this?
> Thanks for assistance
> Zilore
> ga-> sdfopen wrfout_d01_2011-11-19_00:00:00
> Scanning self-describing file: wrfout_d01_2011-11-19_00:00:00
> gadsdf: SDF file has no discernable X coordinate.
> ga->
>
>
> ------------------------------
>
> Message: 7
> Date: Mon, 27 Feb 2012 12:52:45 +0530
> From: Jayakrishnan PR <prjayakrishnan at gmail.com>
> Subject: Re: [gradsusr] netcdf to grib
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <CA+P5ZEFXR0Xj6jOeMo+ZR67AoXcGtvhpqOiTSQ9tinOuCdOzFQ at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear Zilore,
> Try this one. First install CDO in your system. then run the
> following commands.
>
>
> cdo -f grb copy filename.nc output.grb
> cdo gradsdes2 output.grb
>
> This will produce a grib file corresponding to your data and also the
> corresponding ctl file. Then you can open the ctl file and do processing.
> may be helpful. Best of Luck.
>
>
> --
> Sincerely
> ***********************************************
> Jayakrishnan.P.R
> CSIR-Senior Research Fellow
> Department of Atmospheric Sciences
> Cochin University of Science and Technology (CUSAT), Cochin-682 016
> Kerala, India.
> Mob: 09895417565
>
> ***********************************************
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120227/e4ec7051/attachment-0001.html
>
> ------------------------------
>
> Message: 8
> Date: Mon, 27 Feb 2012 12:57:25 +0530
> From: malleswararao maguluri <mallesh4science at gmail.com>
> Subject: [gradsusr] Reg: how to convert station data to gridded data
> To: gradsusr at gradsusr.org
> Message-ID:
> <CAM+Tnx1ovEhJ9jpjYBp3yf9nTp9VxyOU_k73AeDW9WP5oCL1CQ at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> dear all,
>
> i want to convert station data to gridded data. can u give any
> suggestions or programs related to that using fortran and grads.
>
> i done some sample cases with the help of grads cook books, but i am facing
> the following problems
> 1) how to convert sample data to binary file using fortran and how we will
> apply ocares (cressman analysis technique on that) using grads
> 2) how to install stnmap utility in windows and how to operate it.
>
>
> pls help me regarding this.
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120227/192b2040/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 24, Issue 37
> ****************************************
More information about the gradsusr
mailing list