[gradsusr] Creating a Binary File in Grads
Valeria Benson-Lira
vbensonl at asu.edu
Sat Jan 24 16:29:55 EST 2015
Hi Kishore,
Thank you for your suggestions. The problem is that I need to have the
script working for grads.
For instance, my precipitation control files open a .dat file, they are not
a netcdf.
In this case, I want to create a binary file (.dat) which has all the
precipitation values recorded in the 4 different files.
No calculations are needed.
Any suggestions?
Thanks,
Valeria
On Fri, Jan 23, 2015 at 11:11 PM, Kishore Ragi <kishoreragi at gmail.com>
wrote:
> After getting this, if you still want only binary files, please let me
> know. there are couple of options to make nc file to .dat file, even with
> grads.
>
> On Sat, Jan 24, 2015 at 11:37 AM, Kishore Ragi <kishoreragi at gmail.com>
> wrote:
>
>>
>>
>> *Hi Valeria Benson-Lira,*
>> If you don't have problem working with nc files, I will give you other
>> option to achieve your task.
>>
>> First, you make your ctl file into netcdf files as below
>>
>> cdo -f nc import_binary ENCB1_precip_apr2008.ctl precip.april.nc
>> cdo -f nc import_binary ENCB1_precip_may2008.ct precip.may.nc
>> cdo -f nc import_binary ENCB1_precip_jun2008.ctl precip.june.nc
>> cdo -f nc import_binary ENCB1_precip_july2008.ctl precip.july.nc
>>
>> now merge them all in a single file like below
>>
>> cdo mergetime precip.april.nc precip.may.nc precip.june.nc precip.july.nc
>> precip.AMJJ.nc
>>
>> Though, I did not help with grads, hope this helps.
>>
>> Cheers,
>>
>> Kishore
>>
>> On Sat, Jan 24, 2015 at 11:11 AM, Valeria Benson-Lira <vbensonl at asu.edu>
>> wrote:
>>
>>> Dear Grads users,
>>>
>>> Can someone please help me with my script? I am trying to write a binary
>>> file that contains all the values recorded in four different control files.
>>> My new binary file will have the hourly precipitation records from April
>>> through July. So far, the script works for the first month but not for the
>>> following months.
>>>
>>> This is the script:
>>>
>>>
>>> 'reinit'
>>>
>>> 'open ENCB1_precip_apr2008.ctl'
>>>
>>> 'open ENCB1_precip_may2008.ct'
>>>
>>> 'open ENCB1_precip_jun2008.ctl'
>>>
>>> 'open ENCB1_precip_july2008.ctl'
>>>
>>>
>>> 'set gxout fwrite'
>>>
>>> 'set fwrite ENCB1_AMJJ_precip.dat '
>>>
>>>
>>> *****Write April data First.
>>>
>>> *****Recall that April has 720 hours (30 days)
>>>
>>> *index*_t = 1
>>>
>>> while(*index*_t<=720)
>>>
>>> 'set t '*index*_t
>>>
>>> say *index*_t
>>>
>>> 'd precip.1 '
>>>
>>> *index*_t = *index*_t + 1
>>>
>>> endwhile
>>>
>>>
>>> *****Write May data Next.
>>>
>>> *****Recall that May has 744 hours (31 days)
>>>
>>> *index*_t = 1
>>>
>>> while(*index*_t<=744)
>>>
>>> 'set t '*index*_t
>>>
>>> say *index*_t
>>>
>>> 'd precip.2 '
>>>
>>> *index*_t = *index*_t + 1
>>>
>>> endwhile
>>>
>>>
>>> *****Write June data Next.
>>>
>>> *****Recall that June has 720 hours (30 days)
>>>
>>> *index*_t = 1
>>>
>>> while(*index*_t<=720)
>>>
>>> 'set t '*index*_t
>>>
>>> say *index*_t
>>>
>>> 'd precip.3 '
>>>
>>> *index*_t = *index*_t + 1
>>>
>>> endwhile
>>>
>>>
>>> *****Write July data Next.
>>>
>>> *****Recall that July has 744 hours (31 days)
>>>
>>> *index*_t = 1
>>>
>>> while(*index*_t<=744)
>>>
>>> 'set t '*index*_t
>>>
>>> say *index*_t
>>>
>>> 'd precip.4 '
>>>
>>> *index*_t = *index*_t + 1
>>>
>>> endwhile
>>>
>>>
>>> *************************
>>>
>>> 'disable fwrite'
>>>
>>> *************************
>>>
>>> *****Always reinitialize when finished,
>>>
>>> *****to avoid mistakingly writing extraneous data
>>>
>>> 'reinit'
>>>
>>>
>>> Do you have a suggestion of how I can fix it?
>>>
>>> Thanks,
>>>
>>> *Valeria Benson-Lira*
>>> *Graduate Teaching Assistant, Geography M.A.*
>>> *e-mail: vbensonl at asu.edu <vbensonl at asu.edu>*
>>>
>>>
>>>
>>> _______________________________________________
>>> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150124/15f2a9c3/attachment-0001.html
More information about the gradsusr
mailing list