[gradsusr] Low Level I/O Error: Seek error on data file [GrADS] -

[Leopoldo Alvarez]

Hi Jennifer,

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

Now my data file .dat weighs 66 Gb
68926752 - 66G 2013-09-22 11:11 test.dat

I can see the data well, until the time of step 1064 (set t 1064, then 
'd t' it`s ok), but when i exceeded from 1064 time of step gives me that 
error ...

My full GrADS descriptor file ->

dset ^test.dat
options big_endian
title MM5 data
undef -9999.
pdef   60   60 lcc   32.63 -112.54    1.00    1.00   30.00   60.00 
-110.34   6000.   6000.
xdef  149 linear -112.65  0.0270
ydef  123 linear   32.63  0.0270
zdef  25 levels
  0.99689
  0.99191
  0.98944
  0.98818
  0.97377
  0.94500
  0.91000
  0.87000
  0.82500
  0.77500
  0.72500
  0.67500
  0.62500
  0.57500
  0.52500
  0.47500
  0.42500
  0.37500
  0.32500
  0.27500
  0.22500
  0.17500
  0.12500
  0.07500
  0.02500
tdef 1036 linear 00:00Z02JAN2012   60MN
vars 33
rc       0 99 accum conv pcn (cm)
rn       0 99 accum non-c pcn (cm)
ter      0 99 ter elevation (m)
xlat     0 99 cross lat (degree)
xlon     0 99 cross lon (degree)
lu       0 99 land use
t2m      0 99 2 m temperature (K)
q2m      0 99 2m mix ratio (kg/kg)
u10      0 99 10 m u wind (m/sec)
v10      0 99 10 m v wind (m/sec)
pslv     0 99 sea level prs (mb)
u         25 99 u wind (m/s)
v         25 99 v wind (m/s)
w         25 99 vertical vel (m/s)
pp        25 99 prs pert (Pa)
t         25 99 temperature (C)
q         25 99 mixing ratio (kg/kg)
clw       25 99 cloud water (kg/kg)
rnw       25 99 rain water (kg/kg)
rtnd      25 99 rad tend (K/day)
z         25 99 height AGL (m)
h         25 99 geopot height (m)
td        25 99 dewpoint temp (C)
rh        25 99 rel humidity (%)
th        25 99 potential temp (K)
the       25 99 theta-e (K)
prs       25 99 pressure (Pa)
vor       25 99 vorticity (s-1)
pv        25 99 potential vort (pvu)
dbz       25 99 reflectivity (dbz)
div       25 99 divergence (s-1)
dir       25 99 dir (degrees)
tadv      25 99 temp advection (K/s)
endvars


El 21/09/13 09:35, Leopoldo Alvarez escribió:
> Hi Jennifer,
>
> My binary data /grad.exe/ file is 210 Kb. and data file .dat weighs 32 Gb.
> /xsize*ysize*zsize*tsize*esize*4/ , Where I see these parameters? i 
> don´t know how to fix it, where ?, I need to recompile ?.
>
> best regards,
>
> Leopoldo ÁLVAREZ
> Tenerife - Canary Islands
>
>
>
>
> 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
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130923/73626ef4/attachment.html