<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<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 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>
</body>
</html>