pure function zones_2(n1,n2) result(A) integer, intent(in) :: n1,n2 complex(8) :: A(n1,n2) A = one end function zones_2