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><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><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><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 that<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>> > > > 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>