[gradsusr] GrADS executable NOT good [Leopoldo Alvarez] -

Jennifer Adams jma at cola.iges.org
Mon Sep 23 13:56:06 EDT 2013 

On Sep 23, 2013, at 1:46 PM, Leopoldo Alvarez wrote:

> Hi Jennifer, 
> 
> I am writing again because I have urgent need to solve this problem. Thanks!.
> 
> My binary data grad.exe file is 210 Kb
This is not important
> (mm5_to_grads package)
I am not familiar with the mm5_to_grads package and how it creates binary output … however …

> and my data file .dat weighs 32 Gb
A precise file size is required, down to the last byte. 

> AND grads file (GrADS 2.0.2) is 1.5 Mb.
The version number is important, but the size of the executable isn't. 

> xsize*ysize*zsize*tsize*esize*4 , Where I see these parameters? i don´t know how to fix it!!, I need to recompile ?.
These are the dimensions of your grid, taken from the XDEF, YDEF, ZDEF, TDEF, and EDEF (if it's there) statements in your GrADS descriptor file. The 4 is for the four bytes per floating point number. For each variable, calculate the grid size and add them up to see if the size of your binary file is correct. It is only slightly more difficult if some of your variables are z-varying and others aren't. 

> The answer and solution I will post it on the list of GrADS users. Sure!!.
If you need more help, you should include your full GrADS descriptor file. 
--Jennifer

> 
> best regards, 
> 
> Leopoldo ÁLVAREZ
> Tenerife - Canary Islands
> 
> 
> 
> 
> On Sep 20, 2013, at 6:13 PM, Leopoldo Alvarez wrote:
> 
> > Hi jennifer, 
> > 
> > I am running a script (nuevo.gs) to extract info domain mm5 model data. I use GrADS 2.0.2 version. I use 'options big_endian' but my data file weighs 32 GB and
> > shows multiple mutiple on screen- ->
> > 
> > Low Level I/O Error:  Seek error on data file 
> 
> 
> This error generally occurs when your binary data file is too small. Check your file size -- it should be xsize*ysize*zsize*tsize*esize*4 (assuming 4 bytes per data point and you're not writing out sequential data).
> --Jennifer
> 
> 
> El 20/09/13 23:13, Leopoldo Alvarez escribió:
>> Hi jennifer, 
>> 
>> I am running a script (nuevo.gs) to extract info domain mm5 model data. I use GrADS 2.0.2 version. I use 'options big_endian' but my data file weighs 32 GB and
>> shows multiple mutiple on screen- ->
>> 
>> Low Level I/O Error:  Seek error on data file 
>>   Data file name = test.dat 
>> 650 rc=-1 fpos=-172300 pfi->fhdr=0
>> 
>> I think the problem is when running the mm5_to_grads.csh and creates the .ctl and .dat and not script. I can see the data well, until the time of step 1062 (set t 1062, then 'd t' it`s ok), but when i exceeded from 1062 time of step gives me that error ... How i can run my script saving this problem???. One clue http://gradsusr.org/pipermail/gradsusr/2008-April/006883.html
>> 
>> I'm on a machine that works -> uname -a
>> Linux login1 2.6.16.53-0.16-ppc64 #1 SMP Tue Oct 2 16:57:49 UTC 2007 ppc64 ppc64 ppc64 GNU/Linux
>> 
>> 
>> nuevo.gs script
>> ***************
>> ***************
>> 'set lon -110.1829'
>> 'set lat 34.6518'
>> 'set lev 0.99199'
>> 
>> #8737
>> timeend = 4198
>> i = 1
>> while(i <= timeend)
>> 'set t 'i''
>> 'q time'
>> timehr = substr(result, 8, 3)
>> timedy = substr(result, 11, 2)
>> timemo = substr(result, 13, 3)
>> timeyr = substr(result, 16, 4)
>> 
>> #velocidad viento
>> 'tg1=mag(u,v)'
>> 'd tg1'
>> a.nr1 = subwrd(result,4)
>> #direccion viento
>> 'd dir'
>> a.nr11 = subwrd(result,10)
>> # tv de Jeffrey iowa
>> 'tv=t+t*0.61*q'
>> 'rho=prs/(286.9*tv)'
>> 'd rho'
>> a.nr = subwrd(result,4)
>> 'd t-273.15'
>> a.nr2 = subwrd(result,10)
>> 'd rh'
>> a.nr3 = subwrd(result,10)
>> 'd prs/100'
>> a.nr4 = subwrd(result,10)
>> 
>> val = timeyr%"-"%timemo%"-"%timedy"-"%timehr" "%a.nr1" "%a.nr11" "%a.nr" "%a.nr2" "%a.nr3" "%a.nr4
>> rec = write('drylake_80v.txt',val, append)
>> i = i + 1
>> endwhile
>> *******************
>> *******************
>> 
>> best regards, 
>> 
>> Leopoldo ÁLVAREZ
>> Tenerife - Canary Islands
> 
> 
> 
> 
> 
> 
> <logo.png>
> Leopoldo ÁLVAREZ
> CEO - Managing-Partner
> Harmattan Solutions SLU
> Tenerife - Canary Islands
> Mobile: (+34) 663 380 856
> <skp.png>Skype: leopoldo.alvarez El 20/09/13 23:15, Leopoldo Alvarez escribió:
>> I send to gradsusr at gradsusr.org from lap.bolsa at gmail.com email.
>> i don´t know how to fix it.
>> Leopoldo
>> 
>> <Mail Attachment.png>
>> Leopoldo ÁLVAREZ
>> CEO - Managing-Partner
>> Harmattan Solutions SLU
>> Tenerife - Canary Islands
>> Mobile: (+34) 663 380 856
>> <Mail Attachment.png>Skype: leopoldo.alvarez
>> El 20/09/13 22:59, Jennifer Adams escribió:
>>> It will be better for you and me if you would post your questions to the users forum instead of writing directly to me. --Jennifer
>>> 
>>> 
>>> On Sep 20, 2013, at 5:12 PM, Leopoldo Alvarez wrote:
>>> 
>>>> Hi jennifer, 
>>>> 
>>>> I am running a script (nuevo.gs) to extract info grads domain mm5 model data. I use 'options big_endian' but my data file weighs 32 GB and
>>>> shows multiple on screen- ->
>>>> 
>>>> Low Level I/O Error:  Seek error on data file 
>>>>   Data file name = test.dat 
>>>> 650 rc=-1 fpos=-172300 pfi->fhdr=0
>>>> 
>>>> Can i can run my script saving this problem???
>>>> 
>>>> 
>>>> nuevo.gs script
>>>> ***************
>>>> ***************
>>>> 'set lon -110.1829'
>>>> 'set lat 34.6518'
>>>> 'set lev 0.99199'
>>>> 
>>>> #8737
>>>> timeend = 4198
>>>> i = 1
>>>> while(i <= timeend)
>>>> 'set t 'i''
>>>> 'q time'
>>>> timehr = substr(result, 8, 3)
>>>> timedy = substr(result, 11, 2)
>>>> timemo = substr(result, 13, 3)
>>>> timeyr = substr(result, 16, 4)
>>>> 
>>>> #velocidad viento
>>>> 'tg1=mag(u,v)'
>>>> 'd tg1'
>>>> a.nr1 = subwrd(result,4)
>>>> #direccion viento
>>>> 'd dir'
>>>> a.nr11 = subwrd(result,10)
>>>> # tv de Jeffrey iowa
>>>> 'tv=t+t*0.61*q'
>>>> 'rho=prs/(286.9*tv)'
>>>> 'd rho'
>>>> a.nr = subwrd(result,4)
>>>> 'd t-273.15'
>>>> a.nr2 = subwrd(result,10)
>>>> 'd rh'
>>>> a.nr3 = subwrd(result,10)
>>>> 'd prs/100'
>>>> a.nr4 = subwrd(result,10)
>>>> 
>>>> val = timeyr%"-"%timemo%"-"%timedy"-"%timehr" "%a.nr1" "%a.nr11" "%a.nr" "%a.nr2" "%a.nr3" "%a.nr4
>>>> rec = write('drylake_80v.txt',val, append)
>>>> i = i + 1
>>>> endwhile
>>>> *******************
>>>> *******************
>>>> 
>>>> best regards, 
>>>> 
>>>> Leopoldo ÁLVAREZ
>>>> 
>>>> 
>>>> <logo.png>
>>>> Leopoldo ÁLVAREZ
>>>> CEO - Managing-Partner
>>>> Harmattan Solutions SLU
>>>> Tenerife - Canary Islands
>>>> Mobile: (+34) 663 380 856
>>>> <skp.png>Skype: leopoldo.alvarez El 18/09/13 21:29, Jennifer Adams escribió:
>>>>> Oh, your descriptor file already had options little_endian. Change that to 'options big_endian' and then remove 'options byteswapped'. 
>>>>> --Jennifer
>>>>> 
>>>>> On Sep 18, 2013, at 3:54 PM, Leopoldo Alvarez wrote:
>>>>> 
>>>>>> I have to change the header .ctl ?? How?? 
>>>>>> 
>>>>>> <fcbcafhb.png>
>>>>>> 
>>>>>> 
>>>>>> <logo.png>
>>>>>> Leopoldo ÁLVAREZ
>>>>>> CEO - Managing-Partner
>>>>>> Harmattan Solutions SLU
>>>>>> Tenerife - Canary Islands
>>>>>> Mobile: (+34) 663 380 856
>>>>>> <skp.png>Skype: leopoldo.alvarez
>>>>>> 
>>>>>> El 18/09/13 20:38, Jennifer Adams escribió:
>>>>>>> Looks like a big v. little endian issue. Try putting 'options byteswapped' in the descriptor file and see if that changes anything. --Jennifer
>>>>>>> 
>>>>>>> 
>>>>>>> On Sep 18, 2013, at 3:32 PM, Leopoldo Alvarez wrote:
>>>>>>> 
>>>>>>>> Dear Jennirfer Adams, 
>>>>>>>> 
>>>>>>>> I am using the GrADS graphics package for plotting MM5 model outputs...
>>>>>>>> 
>>>>>>>> I'm on a machine that works -> uname -a
>>>>>>>> Linux login1 2.6.16.53-0.16-ppc64 #1 SMP Tue Oct 2 16:57:49 UTC 2007 ppc64 ppc64 ppc64 GNU/Linux
>>>>>>>> 
>>>>>>>> I think the problem happens when I run mm5tograds.csh. The grads.ese executable compiled NOT good, I think.
>>>>>>>> If I paint the terrain domain or the temperature parameter or other... decimal values ​​shows me always ... example below.
>>>>>>>> 
>>>>>>>> These are the compilation options I use ->
>>>>>>>> <ghccdiaj.png>
>>>>>>>> 
>>>>>>>> <cfefadda.png>
>>>>>>>> 
>>>>>>>> 
>>>>>>>> <hgjhefia.png>
>>>>>>>> 
>>>>>>>> 
>>>>>>>> and if I use readv3 to view MMOUTS_DOMAINx output is OK ->
>>>>>>>> 
>>>>>>>> 
>>>>>>>> <ebjefjie.png>
>>>>>>>> 
>>>>>>>> 
>>>>>>>> Best regards, 
>>>>>>>> 
>>>>>>>> Leopoldo ÁLVAREZ
>>>>>>>> 
>>>>>>>> -- 
>>>>>>>> <logo.png>
>>>>>>>> Leopoldo ÁLVAREZ
>>>>>>>> CEO - Managing-Partner
>>>>>>>> Harmattan Solutions SLU
>>>>>>>> Tenerife - Canary Islands
>>>>>>>> Mobile: (+34) 663 380 856
>>>>>>>> <skp.png>Skype: leopoldo.alvarez
>>>>>>> 
>>>>>>> --
>>>>>>> Jennifer M. Adams
>>>>>>> Center for Ocean-Land-Atmosphere Studies (COLA)
>>>>>>> 111 Research Hall, Mail Stop 2B3
>>>>>>> George Mason University
>>>>>>> 4400 University Drive
>>>>>>> Fairfax, VA 22030 
>>>>>>> 
>>>>>>> 
>>>>>>> 
>>>>>>> 
>>>>>>> 
>>>>>> 
>>>>> 
>>>>> --
>>>>> Jennifer M. Adams
>>>>> Center for Ocean-Land-Atmosphere Studies (COLA)
>>>>> 111 Research Hall, Mail Stop 2B3
>>>>> George Mason University
>>>>> 4400 University Drive
>>>>> Fairfax, VA 22030 
>>>>> 
>>>>> 
>>>>> 
>>>>> 
>>>>> 
>>>> 
>>> 
>>> --
>>> Jennifer M. Adams
>>> Center for Ocean-Land-Atmosphere Studies (COLA)
>>> 111 Research Hall, Mail Stop 2B3
>>> George Mason University
>>> 4400 University Drive
>>> Fairfax, VA 22030 
>>> 
>>> 
>>> 
>>> 
>>> 
>> 
> 

--
Jennifer M. Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
111 Research Hall, Mail Stop 2B3
George Mason University
4400 University Drive
Fairfax, VA 22030 





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130923/9a2c9e02/attachment-0001.html

More information about the gradsusr mailing list