load output\t.out
i=1;
while i == 1
   r=input('\n Time series: \n \n 1)  Precipitation \n 2)  Omega \n 3)  v \n 4)  u \n 5)  Surface T \n 6)  Surface SW \n \n 0) Exit \n');
if r == 1
   load output\precip.out
   pt=precip(:,2);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Precipitation, Sea')
   xlabel('time (days)')
   ylabel('Precipitation rate (mm/day)')
   %
   pause
   %
   pt=precip(:,3);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Precipitation, Land')
   xlabel('time (days)')
   ylabel('Precipitation rate (mm/day)')
elseif r == 2
   load output\omegat.out
   pt=omegat(:,2);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Omega, Sea')
   xlabel('time (days)')
   ylabel('Level 9 Omega (mb/hour)')
   %
   pause
   %
   pt=omegat(:,3);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Omega, Land')
   xlabel('time (days)')
   ylabel('Level 9 Omega (mb/hour)')
elseif r == 3
   load output\vt.out
   pt=vt(:,2);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Meridional velocity')
   xlabel('time (days)')
   ylabel('Surface v (m/s)')
elseif r == 4
   load output\ut.out
   pt=ut(:,2);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('U, Sea')
   xlabel('time (days)')
   ylabel('Surface u (m/s)')
   %
   pause
   %
   pt=ut(:,3);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('U, Land')
   xlabel('time (days)')
   ylabel('Surface u (m/s)')
elseif r == 5
   load output\sst.out
   pt=sst(:,2);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Surface Temperature, Sea')
   xlabel('time (days)')
   ylabel('Surface T (^o C)')
   %
   pause
   %
   pt=sst(:,3);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Surface Temperature, Land')
   xlabel('time (days)')
   ylabel('Surface T (^o C)')
elseif r == 6
   load output\sw.out
   pt=sw(:,2);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Surface Absorbed Shortwave, Sea')
   xlabel('time (days)')
   ylabel('Surface SW (W/m^2)')
   %
   pause
   %
   pt=sw(:,3);
   h=plot(t,pt);
   set(h,'linewidth',2);
   set(gca,'fontweight','bold','fontsize',12)
   title('Surface Absorbed Shortwave, Land')
   xlabel('time (days)')
   ylabel('Surface SW (W/m^2)')
elseif r == 0
   break
end
end

   