

%FD_coefs = ReadData('hemisphere.out');

F = FD_coefs.F;
w = FD_coefs.w;
k = FD_coefs.k;
T = FD_coefs.T;
pts = FD_coefs.pts;
conn = FD_coefs.conn;
water_depth = FD_coefs.water_depth;
wave_amp = FD_coefs.wave_amp;
wave_angle = FD_coefs.wave_angle;
sym_planes = FD_coefs.sym_planes;


figure
hold on
for i = 1:length(conn(:,1))
  idx = find(~isnan(conn(i,:)));
  plot3(pts{1}(conn(i,idx)),pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
end
axis equal

if(sym_planes >= 1)
for i = 1:length(conn(:,1))
  idx = find(~isnan(conn(i,:)));
  plot3(pts{1}(conn(i,idx)),-pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
end  
end

if(sym_planes == 2)
for i = 1:length(conn(:,1))
  idx = find(~isnan(conn(i,:)));
  plot3(-pts{1}(conn(i,idx)),pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
  plot3(-pts{1}(conn(i,idx)),-pts{2}(conn(i,idx)),pts{3}(conn(i,idx)),'b-');
end  
end

clear idx