Help with gridded data

[zilore mumba]

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 <yamasaki at FIS.UA.PT> 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!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20050317/e7efed79/attachment.html