'color','w'); 
%  
%  % Plot
%  txt = text(0,0,10,{'Tide Model Driver';'for MATLAB'},'horiz','center',...
%     'vert','middle','fontweight','bold',...
%     'fontangle','italic','fontsize',36); 
%  
%  col = cmocean('phase',length(t)); 
%  txt.Color=col(1,:); 
% 
%  axis(ax); % ensures axis limits match in every frame 
% 
%  % Write the first frame: 
%  gif('tmd_logo_v2.gif','delaytime',1/8,'resolution',200)
%  
%  % Loop through the remaining frames
%  for k = 2:length(t) 
%     
%     % Update ocean colo