help

[Omar Müller]

Aravequia
Yes, I can do this correctly using DIRECT='ACCESS', RECL=6*4
Maybe the problem is the compiler type or the machine, I don't know, but I can
work now with DIRECT ACCESS files!!.
Thanks
Omar

On Thursday 09 November 2006 10:58, José Antonio Aravéquia wrote:
> Omar,
>
> The result in GrADS is the same as when I run the program.
>
> See my output is :
>
> $ crea_dat.exe
>   1.  1.  2.  2.  3.  3.
>
>
> $ gradsc -l
>
> Grid Analysis and Display System (GrADS) Version 1.8SL11
> Copyright (c) 1988-2001 by Brian Doty
> Center for Ocean-Land-Atmosphere Studies
> Institute for Global Environment and Society
> All Rights Reserved
>
> Config: v1.8SL11 32-bit little-endian readline printim
>
> Issue 'q config' command for more information.
>
> GX Package Initialization: Size = 11 8.5
> ga-> open bin.ctl
> Scanning description file:  bin.ctl
> Data file bin.dat is open as file 1
> LON set to 0 2.5
> LAT set to -50.2 -45.2
> LEV set to 1000 1000
> Time values set: 2000:1:1:0 2000:1:1:0
> ga-> set gxout print
> ga-> d val
> Printing Grid -- 6 Values -- Undef = 1e+20
> 1 1 2 2 3 3
> ga->
> ---------------------------------------------------------------------------
>--
>
> 2006/11/9, Omar Müller <ovmuller at fich.unl.edu.ar>:
> > 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