<div dir="ltr"><div><div>Hi Leopoldo,<br><br></div>One solution could be to split big MM5 output file into smaller files. You can do this with the following fortran program:<br><br></div><div><a href="ftp://ftp.ucar.edu/mesouser/MM5V3/Util/SplitMMFiles.f">ftp://ftp.ucar.edu/mesouser/MM5V3/Util/SplitMMFiles.f</a><br>
<br>Then you execute mm5_to_grads to each smaller files.<br><br>I hope this helps you.<br><br></div><div>Regards<br></div><div>Matias<br><br></div><div><br></div><br><br><div><br><br></div></div><div class="gmail_extra"><br>
<br><div class="gmail_quote">2013/9/23 Leopoldo Alvarez <span dir="ltr"><<a href="mailto:lap.bolsa@gmail.com" target="_blank">lap.bolsa@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div 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>El 21/09/13 09:35, Leopoldo Alvarez
escribió:<br>
</div><div><div class="h5">
<blockquote type="cite">
<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>El 20/09/13 23:13, Leopoldo Alvarez
escribió:<br>
</div>
<blockquote type="cite">
<font face="Dax">Hi jennifer, <br>
<br>
I am running a script (<a href="http://nuevo.gs" target="_blank">nuevo.gs</a>) 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 href="http://gradsusr.org/pipermail/gradsusr/2008-April/006883.html" target="_blank">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>
<a href="http://nuevo.gs" target="_blank">nuevo.gs</a> 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 href="http://a.nr" target="_blank">a.nr</a> = 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 href="http://a.nr" target="_blank">a.nr</a>" "%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>
</div></div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>