Program Main
  Implicit None
  Integer :: i,j,n=3
  Integer, Dimension(3,3) :: a=RESHAPE((/1,2,3,4,5,6,7,8,9/),(/3,3/))

  Write(*,*) ' Old  a = '
  Do i=1,n ; Write(*,*) a(i,:) ; Enddo
  Forall (i=1:n, j=1:n, j>i)  a(i,j)=a(j,i) 
  Write(*,*) ' New  a = '
  Do i=1,n ; Write(*,*) a(i,:) ; Enddo

  Forall (i=1:n, j=1:n, a(i,j)>=5) 
      a(i,j)=0
  End Forall
  Write(*,*) ' New  a = '
  Do i=1,n ; Write(*,*) a(i,:) ; Enddo

End Program