[gradsusr] Writing a text output from Grads with header information
Arlindo da Silva
dasilva at alum.mit.edu
Thu Mar 3 20:44:05 EST 2011
On Tue, Mar 1, 2011 at 7:09 PM, Muhammad Rahiz <muhammad.rahiz at ouce.ox.ac.uk
> wrote:
> Dear Felix,
>
> The attached script may do what you want. What it does is to draw the
> values of lon-lat from a text file, and for each lon-lat value, it extracts
> value of variable at each time step and appends to an output txt file.
>
> The output file would contain the following (without headers);
>
> timestep lat lon var
> x.1.x 30.0 70 2.3
> ...
> x.120.x 30.0 70 2.1
> x.1.x 35.0 69.1 1.2
> ...
> x.120.x 35.0 69.1 1.3
>
> If you want to extract the timesteps, grep it;
>
> grep "x.1.x" output.txt > timestep.1.txt
>> grep "x.120.x" output.txt > timestep.120.txt
>>
>
> or in a shell script;
>
> #!/bin/bash
>
> for i in `seq 1 120`; do # where 120 = no. of timesteps
> grep "x.'$i'.x" output.txt > timestep.$i.txt
> done
>
> I've yet to test the script with multiple files according to your file
> specifications as I did not manage to download the files from the link you
> gave. Nonethless, this wouldn't be a problem as you can make the same
> "while-endwhile" loop to read in the files as I did for coordinates.txt. Or
> easier, concantenate the files.
>
> Do tell us how it goes...
>
>
> Muhammad Rahiz
> Researcher & DPhil Candidate (Climate Systems & Policy)
> School of Geography & the Environment
> University of Oxford
>
> On Tue, 1 Mar 2011, Felix Mutua wrote:
>
> Hi Muhammad,
>> I'm using TRMM 3B42 Rainfall data sample on this link.. The data is Netcdf
>> with each 3hrs as a different file..I.e 3*8*365 files(=timesteps)in a
>> year...
>>
>> I was thinking of doing a shell script to open and run the grads script
>> for
>> all the time steps producing the text output for
>> some modelling application.. I'm running Opengrads in Windows and also I
>> have access to a Unix server-based one..
>>
>> Thanks and Regards
>>
>> On Tue, Mar 1, 2011 at 6:32 PM, Muhammad Rahiz
>> <muhammad.rahiz at ouce.ox.ac.uk> wrote:
>> Hi Felix,
>>
>> At each lon & lat, you want to extract rainfall at every time
>> step i.e. 3 hours? Is that what you want to do? If yes, then
>> 1. how many time steps are you dealing with?
>> 2. what is the format of your file i.e. ctl or netcdf? If ctl,
>> then I need two samples of the filename.
>> 3. are you using GrADS in windows or linux?
>>
>> Off the top of my head, what you need is a double loop. If you
>> could give me the details above, then I could modify a similar
>> script that I have.
>>
>> --
>> Muhammad Rahiz
>> Researcher & DPhil Candidate (Climate Systems & Policy)
>> School of Geography & the Environment
>> University of Oxford
>>
>> On Tue, 1 Mar 2011, Felix Mutua wrote:
>>
>> Thanks Muhammad,
>> I've looked at it before...But my challenge is that I have
>> to do this for
>> many files (TRMM 3 hrly precipitation data) for one
>> year..Is there a way I
>> can automate or just get the formatted output in Grads?
>>
>> Thanks
>>
>> Regards
>>
>> On Tue, Mar 1, 2011 at 5:28 PM, Muhammad Rahiz
>> <muhammad.rahiz at ouce.ox.ac.uk> wrote:
>> You need the fprintf.gs script, availableathttp://
>> cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_v
>>
>> ar
>> iable_data_to_a_text_file.
>>
>>
>> Muhammad Rahiz
>> Researcher & DPhil Candidate (Climate Systems & Policy)
>> School of Geography & the Environment
>> University of Oxford
>>
>>
>> On Tue, 1 Mar 2011, Felix Mutua wrote:
>>
>> Hi All, Could someone kindly assist..I'm trying to export
>> a variable into an
>> text output in this format ( with header)
>>
>> Lat Lon Pre
>> ...........................................................
>> 1.5 34.5 4.56
>> 2.0 34.0 3.46
>> .. .... ......
>> .. .... ......
>>
>> For a specified period...
>>
>> Any help greatly appreciated
>>
>>
>> Regards
>>
>> Felix
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>>
>>
>>
>>
>> --
>> _________________________________________________________
>> Felix Mutua
>> School of Engineering, The University of Tokyo
>> River Engineering and Environment Laboratory
>> 7-3-1 Hongo Bunkyo-ku, Tokyo 113-8656, Japan
>> Phone +81-80-3553-0339,E-mail: Felix.mutua at gmail.com,
>> f-mutua at hydra.t.u-tokyo.ac.jp
>> Skype: Felix.mutua
>>
>>
>>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
If you are a python programmer, take a look at this pygrads recipe:
http://cookbooks.opengrads.org/index.php?title=Recipe-003:_Using_PyGrADS_to_save_variable_data_to_a_text_file
I think it does exactly what you want.
Arlindo
--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110303/0efd8274/attachment.html
More information about the gradsusr
mailing list