help

José Antonio Aravéquia araveq at GMAIL.COM
Thu Nov 9 09:38:09 EST 2006


Omar,

I did that test with g77 compiler (Gnu) on a Cygwin terminal in a Pentium IV
machine.

g77 crea_dat.f -o crea_dat.exe

Aravéquia

2006/11/9, Omar Müller <ovmuller at fich.unl.edu.ar>:
>
> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20061109/ae588a0e/attachment.html 


More information about the gradsusr mailing list