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]);