About Fortrans program and a *.ctl in Linux systems

[M.Mahakur]

Try with:

 recl=4*33*17 instead of 33*17 in your open(2,....    statement.

Jim Bao wrote:

> Dear GrADS Users,
>      I have a data, and want to use GrADS to make a figure. I know
> the output has   problem.  But I don't know where is this problem.
> Please you help me check it. Thanks a lot!
>    Here, attached two files, one is Fortran program in Linux system by
> f77,
>    the other is *.ctl
>
>       Best regards
>
>         Zhenhao
>
>   File-1
>          program rw
> cccc
>         parameter (m=1, n1=33,n2=17)
>         real a(n1,n2,m)
>        open(2,file='/home/z3bao/dat1/dat2/read_test',
>      1 form='unformatted',
>      +      access='direct',recl=33*17)
>        open(1, file='/home/z3bao/dat1/dat2/olr_result',
>      1  form='formatted')
>        do ii=1,m
>          do it=1,n2
>           do ib=1,n1
>            read(1,*) a(ib,it,ii)
>            write(6,*) a(ib,it,ii)
>           end do
>          end do
>         end do
> cccc   from January 1950 to Sep. 1995 totally 549
>         ii=1
> cccc    do ii=1,m
> cccc        do it=1,n2
> ccccc         do ib=1,n1
>                write(2,rec=1) ((a(ib,it,ii),ib=1,n1),it=1,n2)
> cccccccc      enddo
> cccc        enddo
> cccc   enddo
>         stop
>         end
>
>   file-2
> DSET ^/home/z3dat/dat1/dat2/read_test
> TITLE Gridded Data
> UNDEF -999.0
> XDEF 33 linear 80 2.5
> YDEF 17 linear 20 2.5
> ZDEF  1 levels 700
> TDEF  1 linear 00z01apr85 6hr
> VARS  1
> hgt 1 99 U component
> endvars
>