help

[Omar Müller]

Aravequia,

If you run the code and then in Grads:
open bin.ctl
set gxout print
d val
You can confirm that the output is:
0 1 1 2 2 3

Omar.



On Wednesday 08 November 2006 13:32, José Antonio Aravéquia wrote:
> Omar,
>
> I noted that you didn't define nv parameter in your program.
>
> I did a test with the attached ctl and fortran code and it works well.
>
> I hope it helps.
>
> Aravequia
>
> 2006/11/8, Omar Müller <ovmuller at fich.unl.edu.ar>:
> > The problem is the 0 value, but I can't resolve them, when I have the
> > data:
> > 1
> > 2
> > 3
> > 1
> > 2
> > 3
> > I hope 1 1 2 2 3 3 in GRADS
> > but I obtain with options sequential
> > 0 1 1 2 2 3 3
> > and without options sequential
> > 3.36312e-44 0 1 1 2 2
> > Omar.
> >
> > > Omar,
> > >
> > > My best guess is that you forgot to include "options sequential" in
> > > your descriptor file.  That could account for the leading 0 in your
> > > output.
> > >
> > > If the remaining points are not in the order you expect, remember that
> > > grads loops through longitude first.
> > >
> > >     Diane
> > >
> > > Omar Müller wrote:
> > > > I prove the next code for create binary data files
> > > >
> > > >
> > > > program crea_dat
> > > >
> > > >         real dum(2,3)
> > > >         parameter (ni=1,nj=1,nk=1,nt=1)
> > > >         OPEN(8,FILE='bin.dat',status='replace',FORM='UNFORMATTED')
> > > >         OPEN(7,FILE='data.dat',STATUS='OLD')
> > > >         do t=1,nt !temp
> > > >           do v=1,nv !var
> > > >             do n=1,nk !level
> > > >                 do i=1,2 !long
> > > >                         do j=1,3 !lat
> > > >                                 read(7,*,END=150)d
> > > > 150                             dum(i,j)=d
> > > >                         end do
> > > >                 end do
> > > >                 write(8)dum
> > > >             end do
> > > >           end do
> > > >         end do
> > > >         CLOSE(8)
> > > >         CLOSE(7)
> > > > end
> > > >
> > > > where data.dat is:
> > > > 1
> > > > 2
> > > > 3
> > > > 1
> > > > 2
> > > > 3
> > > >
> > > > but I obtain in GRADS:
> > > > 0 1 1 2 2 3 when I use d command in gxout print.
> > > > Where is the error?
> > > > Plese help me, I'm from Argentina.
> > >
> > > --
> > > Diane Stokes
> > > Environmental Modeling Center
> > > National Weather Service/NOAA