implicit none

include 'DimArrays.inc'
include 'comread.inc'
integer flin,flout,ntre,ios,chkgo,I,J,K,II,JJ,I1,I2,J1,J2,npf
integer npfi(npun),npfj(npun),numc,totq,MedQ
real Z,ZI,PEN,P,MAXQUO
character*30 NameFl
!elegxos timis dim,jdim
write(*,*) 'dim=',dim,'jdim=',jdim
pause
numc=0
totq=0

do I=1,idim
read(flin,*) (QUO(I,J),j=1,jdim)
do J=1,jdim
C if the elevation is less than 0, then the cell is
C outside the catchment, therefore QUO=-1
!elegxos timis upsometrou
write(*,*) 'To upsometro Q(i,j) einai ',QUO(I,J),'me I=',I,'kai J=',J
if (QUO(I,J).LT.0) then
QUO(I,J)=-1.
else
totq=totq+QUO(I,J)
numc=numc+1
end if
end do
end do
!elegxos timis totq,numc
write(*,*) 'totq=',totq,'numc=',numc
pause