close all %clear all variable
clear all %clear all variable

% Read in NetCDF file
alpha=ncread('TP0_test_file.nc','alpha');

% Create the face
x=[cos(alpha)];
x1=[x fliplr(x)];
y1=[sqrt(1-(x.^2)) -1.*sqrt(1-(x.^2))] ;

% Create the mouth
x2=[-0.3:0.05:0.3];
y2=2*(x2.^2)-0.4;

% Draw aplo
plot(x1,y1,'b','LineWidth',2 ); %Face
hold on
plot(0.05*x1-0.35,0.05*y1+0.35,'k','LineWidth',2 );  %Eye
plot(0.05*x1+0.35,0.05*y1+0.35,'k','LineWidth',2 ); %Eye
plot(x2,y2,'r','LineWidth',2 );                                    %Mouth
hold off

axis([-1.5 1.5 -1.5 1.5]);