Binary to Ascii conversion

[Ben-Jei Tsuang]

Dear Yitbarek Kifle:

I have a utility "b2a", which convert GrADS binary data to ascii data.


Please http:/140.120.3.2/~bjtsuang/grads and search for b2a for details.
Once you download the file from http://140.120.30.2/~bjtsuang/grads/b2a.

Type "make" to create your own b2a utility.


Ben

-------------------

b2a: convert data from binary (REAL*4) (float) to ascii.
(http://140.120.30.2/~bjtsuang/grads/b2a)

b2a ver. 0.02

convert data from binary (REAL*4) (float) to ascii.

copyright by Chia Ying Tu - Revised from a2b.cpp by B. J. Tsuang, Jul., 2002

Usage: b2a infile <n>

Note : 1) n: number of data per row (default= 10).

     2) default output to 'standard output'

-------------------------------------



----- Original Message -----
From: "Diane Stokes" <Diane.Stokes at NOAA.GOV>
To: <GRADSUSR at LIST.CINECA.IT>
Sent: Thursday, March 08, 2007 6:08 AM
Subject: Re: Binary to Ascii conversion


> Kristian,
>
> Did you resolve this issue yet?
>
> If not, my best guess is that you don't have a row centered on 48N.
>
> There are few ways to check.  One would be to run the same script but,
> instead of 'd uwnd', do 'd lat' and check the latitudes that are printed.
>
> 48N likely sits somewhere between your rows 5 & 6, but closer to row 5.
>
>     Diane
>
>
> Kristian Horvath wrote:
>> Hello,
>>
>> I also benefited from the discussion and tried the following script,
>> which resulted in a peculiar result:
>> 'reinit'
>> 'sdfopen /doma/horvath/uwnd.sig995.1961.nc'
>> 'set t 1'
>> 'set lon 15 25'
>> 'set lat 38 48'
>> 'set gxout print'
>> 'set prnopts %10.3e 5 1'
>> 'd uwnd'
>> rc=write(ascii.txt,result)
>> rc=close(ascii.txt)
>> 'quit'
>> return
>>
>> The script works and the result is:
>>
>> Printing Grid -- 30 Values -- Undef = 32766
>>  3.100e+00  1.900e+00  1.400e+00  1.300e+00 -1.200e+00
>>  1.400e+00  1.100e+00  2.000e-01 -6.000e-01 -2.100e+00
>>  2.000e-01 -2.000e-01 -7.000e-01 -5.000e-01 -8.000e-01
>> -1.900e+00 -2.100e+00 -2.500e+00 -1.800e+00 -8.000e-01
>>  1.600e+00  1.400e+00 -1.000e-01 -8.000e-01  1.987e-06
>>  5.700e+00  3.900e+00  1.700e+00  6.000e-01  7.000e-01
>>
>> Now if I modify the script:
>> 'set lon 15 25'
>> 'set lat 38'    ; result is
>> Printing Grid -- 5 Values -- Undef = 32766
>>  3.100e+00  1.900e+00  1.400e+00  1.300e+00 -1.200e+00  ; which is ok.
>>
>> and the second modification:
>> 'set lon 15 25'
>> 'set lat 48'
>> Printing Grid -- 5 Values -- Undef = 32766
>>  1.600e+00  1.400e+00 -1.000e-01 -8.000e-01  1.987e-06  ; which is
>> not in accordance with the "area" example
>>
>> What is the reason for the above discrepancy?
>>
>> Thanks and regards,
>> Kristian
>>
>> At 23:05 22.2.2007, you wrote:
>>> If you don't need the data in binary (which is what you have in file
>>> yiki.dat), you can skip the "gxout fwrite" lines and use "gxout print"
>>> instead.  Eg, in a grads script:
>>>
>>> 'set gxout print'
>>> 'set t 613 697'
>>> 'display tloop(mean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))'
>>> rc=write(yiki.txt,result)
>>> say 'return code from write is: 'rc
>>> rc=close(yiki.txt)
>>> say 'return code from close is: 'rc
>>>
>>> In the above, the output from the display command gets stored in
>>> variable "result" and you write the contents of "result" to ascii file
>>> yiki.txt.
>>>
>>> See the doc for 'set prnopts' at
>>> http://grads.iges.org/grads/gadoc/gradcomdsetprnopts.html
>>> if you want to change the format of the text output.
>>>
>>>    Diane
>>>
>>>
>>> Yitbarek Kifle wrote:
>>>> Dear all,
>>>>
>>>> As a new and inexperienced user of Grads,i have come with a problem of
>>>> converting the Grads time series binary output to Ascii format.I have
>>>> read
>>>> almost all such problems on the users group list but that doesn't
>>>> seem to
>>>> help me. Here is the discription of what i did,
>>>>
>>>> 'set gxout fwrite'
>>>> 'set fwrite yiki.dat'
>>>> set t 613 697
>>>> d tloop(mean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))
>>>> disable fwrite
>>>>
>>>> i have the file Yiki.dat.
>>>>
>>>> Would somebody please tell me what to do next in step by step
>>>> basis,as I do
>>>> not have a programming skill?
>>>>
>>>> thanks in advance.
>>>
>>>
>>> ________ Information from NOD32 ________
>>> This message was checked by NOD32 Antivirus System for Linux Mail
>>> Server.
>>> http://www.nod32.com
>>>
>>> __________ NOD32 2072 (20070220) Information __________
>>>
>>> This message was checked by NOD32 antivirus system.
>>> http://www.eset.com
>>>
>>
>> _________________________________
>> Kristian Horvath                                email:
>> horvath at cirus.dhz.hr
>> Graduate Research Associate                     tel:+385-1-4565752
>> Meteorological and Hydrological Service         www:
>> http://www.dhmz.htnet.hr/index.php
>> Gric 3, 10000 Zagreb, Croatia                   personal www:
>> http://radar.dhz.hr/~horvath
>>