%  Procedure for creating a convex array transducer
%
%  Calling:  Th = xdc_convex_focused_array (no_elements, width, height, kerf, Rconvex, Rfocus 
%                                           no_sub_x, no_sub_y, focus); 
%
%  Parameters:   no_elements -  Number of physical elements.
%                width       -  Width in x-direction of elements.
%                height      -  Width in y-direction of elements.
%                kerf        -  Width in x-direction between elements.
%                Rconvex     -  Convex radius.
%                Rfocus      -  Radius of elevation focus.
%                no_sub_x    -  Number of sub-divisions in x-direction of elements.
%                no_sub_y    -  Number of sub-divisions in y-direction of elements.
%                focus[]     -  Fixed focus for array (x,y,z). Vector with three elements.         
%
%  Return:    A handle Th as a pointer to this transducer aperture.
%
%  Version 1.01, June 26, 1998 by Joergen Arendt Jensen

function Th = xdc_convex_focused_array (no_elements, width, height, kerf, Rconvex, Rfocus, no_sub_x, no_sub_y, focus)

%  Check that all parameters are valid

  if (no_elements<1)
    error ('Field error: Illegal number of physical transducer elements')
    end

  if (width<=0) | (height<=0)
    error ('Field error: Width or height is negativ or zero')
    end

  if (kerf<0)
    error ('Field error: Kerf is negativ')
    end

  if (height > 2*Rfocus)
    error ('Field error: Illegal element height for the chosen elevation focus')
    end

  if ( (width*no_elements+kerf*(no_elements-1)) > (pi*Rconvex))
    error ('Field error: Illegal element width and kerf for the chosen convex radius')
    end

  if (Rconvex<0)
    error ('Field error: Convex radius is negative')
    end

  if (pi*Rconvex<=(kerf*(no_elements-1)+width*no_elements))
    error ('Field error: Width of elements is to large compared to Rconvex')
    end
    
  if (Rfocus<0)
    error ('Field error: Elevation focus radius is negative')
    end

  if (no_sub_x<1)
    error ('Field error: Number of mathematical elements must 1 or more in x-direction')
    end

  if (no_sub_y<2)
    error ('Field error: Number of mathematical elements in y-direction must 2 or more to model elevation focusing')
    end

  if (min(size(focus))~=1) | (max(size(focus))~=3)
    error ('Field error: Focus must be a vector with three elements')
    end

%  Call the C-part of the program to create aperture

  Th = Mat_field (1005,no_elements, width, height, kerf, Rconvex, Rfocus, no_sub_x, no_sub_y, focus);