Help with gridded data

[Ricardo Hallak]

Zilore,

sometimes it is only a problem of swap of words in the binary files.

Try use:

OPTIONS little_endian

or

OPTIONS big_endian

in the .ctl file.

The correct option will depend on the type of your machine.

Regards,

Ricardo

On Thu, 17 Mar 2005 09:18:59 -0800, zilore mumba wrote
> Dear Yoshihiro,
> I tried, as you advised;
> recl=4*l*m,
> write(10) cor, no loop,
> The resulting binary file looks much better than before, but 'd cor' gives 'display nan to nan interval inf.'
> Thanks all the same, i will keep trying.
> Zilore
>
>
> Yoshihiro Yamasaki <yamasaki at FIS.UA.PT> wrote:

Hi !
> Try the following :
>
> parameter ( L = , M = )
> dimension cor (L,M)
> open your file with ---> recl = 4 * L * M
> irec = 1
> write(10,rec=irec) cor
>
> no do loop is needed
>
> yyamazaki
>
> Em Thu, 17 Mar 2005 03:12:18 -0800
> zilore mumba escreveu:
> >
> >
> >
> > Dear Yoshihiro,
> >
> > Thanks very much for sparing your time to spot my silly
> >slip. Evidently the program was writing zeros. However it
> >is still not displaying as required.
> >
> >
> >
> >
> > Yoshihiro Yamasaki wrote:Dear
> >Zilone,
> >
> > Assuming that you have DIMENSION COR(55,57), you made a
> > mistake in your do loop index ( do j and not k )
> >
> > do J = 1 , m
> > irec=irec+1
> > write(10,rec=irec) (cor(i,j),i=1,l)
> > enddo
> >
> > good luck !
> >
> >
> > Em Wed, 16 Mar 2005 23:21:42 -0800
> > zilore mumba escreveu:
> >> I have an apparently very simple problem.
> >> I have correlation data cor(i,j) over an xy domain.
> >>Using a fortran program I have written the data into
> >>binary to display in Grads. Apparently Grads is able to
> >>open the data file.
> >> when I put recl=4*l, the display is blank (with only map
> >>of Africa of course)
> >> when I put recl=8*l, the display says constant field
> >>value=0.
> >>
> >> Here are snippets from my fortran program:
> >>
> >> open
> >>(10,file='zcor824.dat',form='unformatted',access='direct',
> >>recl=4*(l+1))
> >>
> >> irec=0
> >> do k=1,m
> >> irec=irec+1
> >> write(10,rec=irec) (cor(i,j),i=1,l)
> >> enddo
> >>
> >> and my ctl file
> >>
> >> dset zcor824.dat
> >> undef -1.00E+09
> >> title correlation
> >> xdef 55 linear -30 1.5
> >> ydef 57 linear -42 1.5
> >> zdef 1 levels 850
> >> tdef 999 linear 00Z29AUG97 6hr
> >> vars 1
> >> cor 1 99 linear correl coeff
> >> endvars
> >>
> >> Any help would be very welcome to resolve my problem
> >> Thanks all.
> >>
> >>
> >>
> >>
> >>
> >> ---------------------------------
> >> Do you Yahoo!?
> >> Yahoo! Small Business - Try our new resources site!
> >
> >
> >
> >
> >
> > ---------------------------------
> > Do you Yahoo!?
> > Yahoo! Small Business - Try our new resources site!
>
>
>
-----------------------------------------------------------------------
Do you Yahoo!?
> Yahoo! Small Business - Try our new resources site!

Hallak

Laboratório de Estudos de Tempestades Severas (LETS)
Departamento de Ciências Atmosféricas (DCA)
Instituto de Astronomia, Geofísica e Ciências Atmosféricas (IAG)
Universidade de São Paulo (USP)
Rua do Matão, 1226 - Cidade Universitária -São Paulo SP -05508-900
Fone.: +55 (0-XX-11) 3091-4735
Fax: +55 (0-XX-11) 3091-4714


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20050317/83107c10/attachment.html