Omar,<br>
<br>
I did that test with g77 compiler (Gnu) on a Cygwin terminal in a Pentium IV machine.<br>
<br>
g77 crea_dat.f -o crea_dat.exe <br>
<br>
Aravéquia<br><br><div><span class="gmail_quote">2006/11/9, Omar Müller <<a href="mailto:ovmuller@fich.unl.edu.ar">ovmuller@fich.unl.edu.ar</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Aravequia<br>Yes, I can do this correctly using DIRECT='ACCESS', RECL=6*4<br>Maybe the problem is the compiler type or the machine, I don't know, but I can<br>work now with DIRECT ACCESS files!!.<br>Thanks<br>Omar<br><br>
On Thursday 09 November 2006 10:58, José Antonio Aravéquia wrote:<br>> Omar,<br>><br>> The result in GrADS is the same as when I run the program.<br>><br>> See my output is :<br>><br>> $ crea_dat.exe<br>
> 1. 1. 2. 2. 3. 3.<br>><br>><br>> $ gradsc -l<br>><br>> Grid Analysis and Display System (GrADS) Version 1.8SL11<br>> Copyright (c) 1988-2001 by Brian Doty<br>> Center for Ocean-Land-Atmosphere Studies
<br>> Institute for Global Environment and Society<br>> All Rights Reserved<br>><br>> Config: v1.8SL11 32-bit little-endian readline printim<br>><br>> Issue 'q config' command for more information.<br>>
<br>> GX Package Initialization: Size = 11 8.5<br>> ga-> open bin.ctl<br>> Scanning description file: bin.ctl<br>> Data file bin.dat is open as file 1<br>> LON set to 0 2.5<br>> LAT set to -50.2 -45.2
<br>> LEV set to 1000 1000<br>> Time values set: 2000:1:1:0 2000:1:1:0<br>> ga-> set gxout print<br>> ga-> d val<br>> Printing Grid -- 6 Values -- Undef = 1e+20<br>> 1 1 2 2 3 3<br>> ga-><br>
> ---------------------------------------------------------------------------<br>>--<br>><br>> 2006/11/9, Omar Müller <<a href="mailto:ovmuller@fich.unl.edu.ar">ovmuller@fich.unl.edu.ar</a>>:<br>> > Aravequia,
<br>> ><br>> > If you run the code and then in Grads:<br>> > open bin.ctl<br>> > set gxout print<br>> > d val<br>> > You can confirm that the output is:<br>> > 0 1 1 2 2 3<br>> >
<br>> > Omar.<br>> ><br>> > On Wednesday 08 November 2006 13:32, José Antonio Aravéquia wrote:<br>> > > Omar,<br>> > ><br>> > > I noted that you didn't define nv parameter in your program.
<br>> > ><br>> > > I did a test with the attached ctl and fortran code and it works well.<br>> > ><br>> > > I hope it helps.<br>> > ><br>> > > Aravequia<br>> > >
<br>> > > 2006/11/8, Omar Müller <<a href="mailto:ovmuller@fich.unl.edu.ar">ovmuller@fich.unl.edu.ar</a>>:<br>> > > > The problem is the 0 value, but I can't resolve them, when I have the<br>> > > > data:
<br>> > > > 1<br>> > > > 2<br>> > > > 3<br>> > > > 1<br>> > > > 2<br>> > > > 3<br>> > > > I hope 1 1 2 2 3 3 in GRADS<br>> > > > but I obtain with options sequential
<br>> > > > 0 1 1 2 2 3 3<br>> > > > and without options sequential<br>> > > > 3.36312e-44 0 1 1 2 2<br>> > > > Omar.<br>> > > ><br>> > > > > Omar,
<br>> > > > ><br>> > > > > My best guess is that you forgot to include "options sequential" in<br>> > > > > your descriptor file. That could account for the leading 0 in your
<br>> > > > > output.<br>> > > > ><br>> > > > > If the remaining points are not in the order you expect, remember<br>> ><br>> > that<br>> ><br>> > > > > grads loops through longitude first.
<br>> > > > ><br>> > > > > Diane<br>> > > > ><br>> > > > > Omar Müller wrote:<br>> > > > > > I prove the next code for create binary data files
<br>> > > > > ><br>> > > > > ><br>> > > > > > program crea_dat<br>> > > > > ><br>> > > > > > real dum(2,3)<br>> > > > > > parameter (ni=1,nj=1,nk=1,nt=1)
<br>> > > > > ><br>> > > > > > OPEN(8,FILE='bin.dat',status='replace',FORM='UNFORMATTED')<br>> > > > > > OPEN(7,FILE='data.dat',STATUS='OLD')<br>> > > > > > do t=1,nt !temp
<br>> > > > > > do v=1,nv !var<br>> > > > > > do n=1,nk !level<br>>
> > > >
>
do i=1,2 !long<br>> > > > >
>
do j=1,3 !lat<br>> > > > >
>
read(7,*,END=150)d<br>> > > > > >
150
dum(i,j)=d<br>> > > > >
>
end do<br>> > > > >
>
end do<br>> > > > >
>
write(8)dum<br>> > > > > > end do<br>> > > > > > end do<br>> > > > > > end do<br>> > > > > > CLOSE(8)<br>> > > > > > CLOSE(7)
<br>> > > > > > end<br>> > > > > ><br>> > > > > > where data.dat is:<br>> > > > > > 1<br>> > > > > > 2<br>> > > > > > 3
<br>> > > > > > 1<br>> > > > > > 2<br>> > > > > > 3<br>> > > > > ><br>> > > > > > but I obtain in GRADS:<br>> > > > > > 0 1 1 2 2 3 when I use d command in gxout print.
<br>> > > > > > Where is the error?<br>> > > > > > Plese help me, I'm from Argentina.<br>> > > > ><br>> > > > > --<br>> > > > > Diane Stokes
<br>> > > > > Environmental Modeling Center<br>> > > > > National Weather Service/NOAA<br></blockquote></div><br>