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

Mon Sep 23 16:19:53 EDT 2013

Hi Leopoldo,

One solution could be to split big MM5 output file into smaller files.  You
can do this with the following fortran program:

ftp://ftp.ucar.edu/mesouser/MM5V3/Util/SplitMMFiles.f

Then you execute mm5_to_grads to each smaller files.

I hope this helps you.

Regards
Matias

2013/9/23 Leopoldo Alvarez <lap.bolsa at gmail.com>

>  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
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130923/75fe94e2/attachment.html 