program ch0800a
implicit none
integer, parameter:: n = 3
integer:: i,j
integer, dimension(1:n)     :: colsum = 0
integer, dimension (1:n,1:n):: a = 0
! work with “a” column by column for efficiency
   print *,'input elements of  array "a" size = ',n
   print*,'1 element per line column-wise'
   do j =1,n 
      do i=1,n
        read*, a(i,j)
      end do
   end do
  ! calculate column sums
   do j=1,n
      colsum(j) = 0
      do i=1,n
         colsum(j) = colsum(j) + a(i,j)
      end do
   end do
! print out a and column sums below
   do i =1,n
      print*,a(i,1:n)
   end do
   print*, colsum(1:n)
end program ch0800a