low level I/O error

[Roop]

Hi,

I am getting the following error beyond t=169:

Low Level I/O Error: Read  error on data file
Data File name = boxave_gpcp.r4
Error reading 13 bytes at location 24180
Data Request Error: Error for variable 'gpcp'
Error occurred at column 1
Display error: Invalid expression
Expression = gpcp

I am using GrADS 1.9.0-rc1 on a PC.

My control file looks like:

dset ^boxave_gpcp.r4
undef -99999.0
title GPCP data
xdef 13 linear 71.25 2.5
ydef 11 linear 6.25 2.5
tdef 264 linear Jan1960 1mo
zdef 01 levels 1
vars 1
gpcp 0 99 monsoon onset GPCP
endvars

I am trying to average "precipitation" over a certain area.

In case, you would like to know what my Fortran program looks like:

       Program Box Average

       parameter(nx=144)
       parameter(ny=72)
       parameter(nt=264)
!       parameter(nlon=13)
!       parameter(nlat=11)
       parameter(lon1=30)
       parameter(lon2=42)
       parameter(lat1=39)
       parameter(lat2=49)
       parameter(undef=-99999.0)

       integer ix
       integer iy
       integer it

       real gpcp(nx,ny,nt)
       real gpcp_ave(nt)
       real count(nt)

       open(13,file="anomaly_gpcp.r4",form="unformatted",
     & recl=nx*ny*nt*4)
       read(13,rec=1)gpcp

       do it=1,nt

             count(it)=0.0
             gpcp_ave(it)=0.0

       enddo

       do ix=lon1,lon2
         do iy=lat1,lat2
           do it=1,nt

             if(gpcp(ix,iy,it).ne.undef) then

             gpcp_ave(it)=gpcp_ave(it)+gpcp(ix,iy,it)
             count(it)=count(it)+1.0

             endif

           enddo
         enddo
       enddo

       do it=1,nt

           gpcp_ave(it)=gpcp_ave(it)/count(it)

       enddo

       open(24,file="boxave_gpcp.r4",form="unformatted",
     &  recl=nt*4)
       write(24,rec=1)gpcp_ave

       write(*,*)"Completed GPCP Box Average"

       stop
       end

Can anyone help me as to what is the problem?

Thanks,
Roop.