<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font face="Dax">Hi Jennifer, <br>
<br>
>> and my data file .dat weighs 32 Gb<br>
>A precise file size is required, down to the last byte. <br>
<br>
Now my data file .dat weighs 66 Gb<br>
68926752 - 66G 2013-09-22 11:11 test.dat<br>
<br>
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 ...<br>
<br>
My full GrADS descriptor file -><br>
<br>
dset ^test.dat<br>
options big_endian<br>
title MM5 data<br>
undef -9999.<br>
pdef 60 60 lcc 32.63 -112.54 1.00 1.00 30.00 60.00
-110.34 6000. 6000.<br>
xdef 149 linear -112.65 0.0270<br>
ydef 123 linear 32.63 0.0270<br>
zdef 25 levels<br>
0.99689<br>
0.99191<br>
0.98944<br>
0.98818<br>
0.97377<br>
0.94500<br>
0.91000<br>
0.87000<br>
0.82500<br>
0.77500<br>
0.72500<br>
0.67500<br>
0.62500<br>
0.57500<br>
0.52500<br>
0.47500<br>
0.42500<br>
0.37500<br>
0.32500<br>
0.27500<br>
0.22500<br>
0.17500<br>
0.12500<br>
0.07500<br>
0.02500<br>
tdef 1036 linear 00:00Z02JAN2012 60MN<br>
vars 33<br>
rc 0 99 accum conv pcn (cm)<br>
rn 0 99 accum non-c pcn (cm)<br>
ter 0 99 ter elevation (m)<br>
xlat 0 99 cross lat (degree)<br>
xlon 0 99 cross lon (degree)<br>
lu 0 99 land use<br>
t2m 0 99 2 m temperature (K)<br>
q2m 0 99 2m mix ratio (kg/kg)<br>
u10 0 99 10 m u wind (m/sec)<br>
v10 0 99 10 m v wind (m/sec)<br>
pslv 0 99 sea level prs (mb)<br>
u 25 99 u wind (m/s)<br>
v 25 99 v wind (m/s)<br>
w 25 99 vertical vel (m/s)<br>
pp 25 99 prs pert (Pa)<br>
t 25 99 temperature (C)<br>
q 25 99 mixing ratio (kg/kg)<br>
clw 25 99 cloud water (kg/kg)<br>
rnw 25 99 rain water (kg/kg)<br>
rtnd 25 99 rad tend (K/day)<br>
z 25 99 height AGL (m)<br>
h 25 99 geopot height (m)<br>
td 25 99 dewpoint temp (C)<br>
rh 25 99 rel humidity (%)<br>
th 25 99 potential temp (K)<br>
the 25 99 theta-e (K)<br>
prs 25 99 pressure (Pa)<br>
vor 25 99 vorticity (s-1)<br>
pv 25 99 potential vort (pvu)<br>
dbz 25 99 reflectivity (dbz)<br>
div 25 99 divergence (s-1)<br>
dir 25 99 dir (degrees)<br>
tadv 25 99 temp advection (K/s)<br>
endvars<br>
<br>
<br>
</font>
<div class="moz-cite-prefix">El 21/09/13 09:35, Leopoldo Alvarez
escribió:<br>
</div>
<blockquote cite="mid:523D5A35.3040202@gmail.com" type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<font face="Dax">Hi Jennifer, <br>
<br>
My binary data <i>grad.exe</i> file is 210 Kb. and data file
.dat weighs 32 Gb.<br>
<font color="#006600"><i>xsize*ysize*zsize*tsize*esize*4</i></font>
</font><font face="Dax">, Where I see these parameters? i don´t
know how to fix it, where ?, I need to recompile ?.<br>
</font><br>
<font face="Dax">best regards, <br>
<br>
Leopoldo ÁLVAREZ<br>
Tenerife - Canary Islands<br>
<br>
<br>
<br>
<br>
</font>
<div class="moz-cite-prefix">El 20/09/13 23:13, Leopoldo Alvarez
escribió:<br>
</div>
<blockquote cite="mid:523CC87B.5000601@gmail.com" type="cite">
<meta http-equiv="content-type" content="text/html;
charset=ISO-8859-1">
<font face="Dax">Hi jennifer, <br>
<br>
I am running a script (nuevo.gs) to extract info domain mm5
model data. </font><font face="Dax"><font face="Dax">I use </font><font
face="Dax"><font face="Dax">GrADS </font></font><font
face="Dax"><font face="Dax"> 2.0.2 </font>version.</font> I
use </font><font face="Dax">'options big_endian' but my data
file weighs 32 GB and<br>
shows multiple mutiple on screen- -><br>
<br>
<font color="#006600"><i>Low Level I/O Error: Seek error on
data file </i><i><br>
</i><i> Data file name = test.dat </i><i><br>
</i><i>650 rc=-1 fpos=-172300 pfi->fhdr=0</i></font></font><br>
<br>
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 <a moz-do-not-send="true" class="moz-txt-link-freetext"
href="http://gradsusr.org/pipermail/gradsusr/2008-April/006883.html">http://gradsusr.org/pipermail/gradsusr/2008-April/006883.html</a><br>
<br>
<font face="Dax">I'm on a machine that works -> uname -a<br>
<i>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></font><br>
<br>
<br>
nuevo.gs script<br>
***************<br>
***************<br>
<font color="#006600">'set lon -110.1829'<br>
'set lat 34.6518'<br>
'set lev 0.99199'<br>
<br>
#8737<br>
timeend = 4198<br>
i = 1<br>
while(i <= timeend)<br>
'set t 'i''<br>
'q time'<br>
timehr = substr(result, 8, 3)<br>
timedy = substr(result, 11, 2)<br>
timemo = substr(result, 13, 3)<br>
timeyr = substr(result, 16, 4)<br>
<br>
#velocidad viento<br>
'tg1=mag(u,v)'<br>
'd tg1'<br>
a.nr1 = subwrd(result,4)<br>
#direccion viento<br>
'd dir'<br>
a.nr11 = subwrd(result,10)<br>
# tv de Jeffrey iowa<br>
'tv=t+t*0.61*q'<br>
'rho=prs/(286.9*tv)'<br>
'd rho'<br>
a.nr = subwrd(result,4)<br>
'd t-273.15'<br>
a.nr2 = subwrd(result,10)<br>
'd rh'<br>
a.nr3 = subwrd(result,10)<br>
'd prs/100'<br>
a.nr4 = subwrd(result,10)<br>
<br>
val = timeyr%"-"%timemo%"-"%timedy"-"%timehr" "%a.nr1"
"%a.nr11" "%a.nr" "%a.nr2" "%a.nr3" "%a.nr4<br>
rec = write('drylake_80v.txt',val, append)<br>
i = i + 1<br>
endwhile</font><br>
*******************<br>
*******************<br>
<br>
best regards, <br>
<br>
Leopoldo ÁLVAREZ<br>
Tenerife - Canary Islands<br>
</blockquote>
<br>
</blockquote>
<br>
</body>
</html>