Help with gridded data

[Yoshihiro Yamasaki]

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 <zmumba at YAHOO.COM> 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 <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!