gmaze_pv/visu/eg_view_Timeserie_pl2.m

%DEF All components of the relative vorticity

% Map projection:
m_proj('mercator','long',subdomain.limlon,'lat',subdomain.limlat);
%m_proj('mercator','long',subdomain2.limlon,'lat',subdomain2.limlat);
%m_proj('mercator','long',subdomain.limlon,'lat',[25 40]);
%m_proj('mercator','long',[subdomain.limlon(1) 360-24],'lat',[25 50]);

% Which variables to plot:
wvar = [1:3];
iz = 2; % Surface
%iz = 6; % Core of the Gulf Stream
%iz = 22; % Under the Gulf Stream

figure(12);clf;hold on;iw=1;jw=length(wvar);
  
for ip = 1 : length(wvar)

% Variables loop:
  % Default:
  CBAR = 'h'; % Colorbar orientation
  Tcontour = [17 19]; Tcolor = [0 0 0]; % Theta contours
  Hcontour = -[0:200:600]; Hcolor = [0 0 0]; % MLD contours
  unit = ''; % Default unit
  load mapanom2 ; N = 256; c = [0 0]; cmap = jet; % Colormaping
  showT = 1; % Iso-Theta contours
  showW = 0; % Windstress arrows
  showH = 0; colorW = 'w'; % Mixed Layer Depth
  showCLIM = 0; % Show CLIMODE region box
  CONT  = 0; % Contours instead of pcolor
  CONTshlab = 0; % Show label for contours plot
  colorCOAST = [0 0 0]; % Land color
  SHADE = 'flat'; % shading option
  
  N     = 32;
  c     = [12 30];
  cx    = [-(1+2*c(1)) 1+2*c(2)]*6e-5; cmap = mapanom;
  titf  = 'OMEGA';
  
switch wvar(ip)
 case 1
  C     = -squeeze(OX(iz,:,:));
  Clon  = OXlon; Clat = OXlat;
  tit   = strcat('\omega_x = - \partial v / \partial z');
  showW = 0; % Windstress
  showH = 0;
  showCLIM = 1;
  unit  = '1/s';
 case 2
  C     = -squeeze(OY(iz,:,:));
  Clon  = OYlon; Clat = OYlat;
  tit   = strcat('\omega_y =  \partial u / \partial z');
  showW = 0; % Windstress
  showH = 0;
  showCLIM = 1;
  %N     = 256;
  %c     = [10 10];
  %cx    = [-(1+2*c(1)) 1+2*c(2)]*6e-5; cmap = mapanom;
  unit  = '1/s';
 case 3
  C     = squeeze(ZETA(iz,:,:));
  Clon  = ZETAlon; Clat = ZETAlat;
  tit   = strcat('\zeta = \partial v / \partial x - \partial u / \partial y');
  showW = 0; % Windstress
  showH = 0;
  showCLIM = 1;
  %N     = 256;
  %c     = [10 10];
  %cx    = [-(1+2*c(1)) 1+2*c(2)]*6e-5; cmap = mapanom;
  unit  = '1/s';
end %switch what to plot


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Draw variable:
sp=subplot(iw,jw,ip);hold on
if CONT ~= 1
  m_pcolor(Clon,Clat,C);
  shading(SHADE);
  colormap(logcolormap(N,c(1),c(2),cmap));
else 
  [cs,h] = m_contourf(Clon,Clat,C,CONTv);
  colormap(mycolormap(logcolormap(N,c(1),c(2),cmap),length(CONTv)));
  if CONTshlab
    clabel(cs,h,'labelspacing',200,'fontsize',8)
  end
end
caxis(cx);
if ip == 2
  ccol = colorbar(CBAR,'fontsize',10);
  ctitle(ccol,unit);
  posiC = get(ccol,'position');
  set(ccol,'position',[.2 posiC(2) 1-2*.2 .02]);
end
title(tit);
m_coast('patch',colorCOAST);
m_grid('xtick',360-[20:5:80],'ytick',[20:2:50]);
set(gcf,'name',titf);


if showT
  [cs,h] = m_contour(Tlon,Tlat,squeeze(T(1,:,:)),Tcontour);
  clabel(cs,h,'fontsize',8,'color',[0 0 0],'labelspacing',200);
  for ih=1:length(h)
    set(h(ih),'edgecolor',Tcolor,'linewidth',1);
  end
end %if show THETA contours

if showW
      dx = 10*diff(Txlon(1:2)); dy = 8*diff(Txlat(1:2));
      dx = 20*diff(Txlon(1:2)); dy = 10*diff(Txlat(1:2));
      lo = [Txlon(1):dx:Txlon(length(Txlon))];
      la = [Txlat(1):dy:Txlat(length(Txlat))];
      [lo la] = meshgrid(lo,la);
      Txn = interp2(Txlat,Txlon,Tx',la,lo);
      Tyn = interp2(Txlat,Txlon,Ty',la,lo);
      s = 2;
      m_quiver(lo,la,Txn,Tyn,s,colorW,'linewidth',1.25);
%      m_quiver(lo,la,-(1+sin(la*pi/180)).*Txn,(1+sin(la*pi/180)).*Tyn,s,'w');
      m_quiver(360-84,47,1,0,s,'w','linewidth',1.25);
        m_text(360-84,48,'1 N/m^2','color','w');
end %if show windstress

if showH
  %[cs,h] = m_contour(MLDlon,MLDlat,MLD,Hcontour);
  cm = flipud(mycolormap(jet,length(Hcontour)));
  cm = mycolormap([linspace(0,0,20); linspace(1,.5,20) ;linspace(0,0,20)]',length(Hcontour));
  cm = [0 1 0 ; 0 .6 0 ; 0 .2 0 ; 0 0 0];
  for ii = 1 : length(Hcontour) 
    [cs,h] = m_contour(MLDlon,MLDlat,MLD,[1 1]*Hcontour(ii)); 
    if ~isempty(cs)
%    clabel(cs,h,'fontsize',8,'color',[0 0 0],'labelspacing',300);
    clabel(cs,h,'fontsize',8,'color',cm(ii,:),'labelspacing',600,'fontweight','bold');
    for ih=1:length(h)
%      set(h(ih),'edgecolor',Hcolor,'linewidth',1);
      set(h(ih),'edgecolor',cm(ii,:),'linewidth',1.2);
    end
    end
  end
end %if show Mixed Layer depth

if showCLIM
  m_line(360-[71 62 62 71 71],[36 36 40.5 40.5 36],'color','r','linewidth',1.5)
end

%suptitle(strcat('Relative vorticity component at depth:',num2str(OYdpt(iz)),'m'));

  
end %for ip
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if 1 % Show the date in big in the upper left corner
  spp=subplot('position',[0 .95 .25 .05]);
  p=patch([0 1 1 0],[0 0 1 1],'w');
  set(spp,'ytick',[],'xtick',[]);
  set(spp,'box','off');
  text(0.1,.5,num2str(TIME(it,:)),'fontsize',16,...
       'fontweight','bold','color','r','verticalalign','middle');
end  


%%%%%%%%%%%%%%%%
drawnow
set(gcf,'position',[4 48 888 430]);
videotimeline(num2str(zeros(size(TIME,1),1)),it,'b')
%videotimeline(TIME,it,'b')
if prtimg
set(gcf,'color','white') 
set(findobj('tag','m_grid_color'),'facecolor','none')
set(gcf,'paperposition',[0.6 6.5 25 14]);
exportj(gcf,1,strcat(outimg,sla,titf,'.',snapshot));
end %if

1	gmaze	1.1	%DEF All components of the relative vorticity
2
3			% Map projection:
4			m_proj('mercator','long',subdomain.limlon,'lat',subdomain.limlat);
5			%m_proj('mercator','long',subdomain2.limlon,'lat',subdomain2.limlat);
6			%m_proj('mercator','long',subdomain.limlon,'lat',[25 40]);
7			%m_proj('mercator','long',[subdomain.limlon(1) 360-24],'lat',[25 50]);
8
9			% Which variables to plot:
10			wvar = [1:3];
11			iz = 2; % Surface
12			%iz = 6; % Core of the Gulf Stream
13			%iz = 22; % Under the Gulf Stream
14
15			figure(12);clf;hold on;iw=1;jw=length(wvar);
16
17			for ip = 1 : length(wvar)
18
19			% Variables loop:
20			% Default:
21			CBAR = 'h'; % Colorbar orientation
22			Tcontour = [17 19]; Tcolor = [0 0 0]; % Theta contours
23			Hcontour = -[0:200:600]; Hcolor = [0 0 0]; % MLD contours
24			unit = ''; % Default unit
25			load mapanom2 ; N = 256; c = [0 0]; cmap = jet; % Colormaping
26			showT = 1; % Iso-Theta contours
27			showW = 0; % Windstress arrows
28			showH = 0; colorW = 'w'; % Mixed Layer Depth
29			showCLIM = 0; % Show CLIMODE region box
30			CONT = 0; % Contours instead of pcolor
31			CONTshlab = 0; % Show label for contours plot
32			colorCOAST = [0 0 0]; % Land color
33			SHADE = 'flat'; % shading option
34
35			N = 32;
36			c = [12 30];
37			cx = [-(1+2c(1)) 1+2c(2)]*6e-5; cmap = mapanom;
38			titf = 'OMEGA';
39
40			switch wvar(ip)
41			case 1
42			C = -squeeze(OX(iz,:,:));
43			Clon = OXlon; Clat = OXlat;
44			tit = strcat('\omega_x = - \partial v / \partial z');
45			showW = 0; % Windstress
46			showH = 0;
47			showCLIM = 1;
48			unit = '1/s';
49			case 2
50			C = -squeeze(OY(iz,:,:));
51			Clon = OYlon; Clat = OYlat;
52			tit = strcat('\omega_y = \partial u / \partial z');
53			showW = 0; % Windstress
54			showH = 0;
55			showCLIM = 1;
56			%N = 256;
57			%c = [10 10];
58			%cx = [-(1+2c(1)) 1+2c(2)]*6e-5; cmap = mapanom;
59			unit = '1/s';
60			case 3
61			C = squeeze(ZETA(iz,:,:));
62			Clon = ZETAlon; Clat = ZETAlat;
63			tit = strcat('\zeta = \partial v / \partial x - \partial u / \partial y');
64			showW = 0; % Windstress
65			showH = 0;
66			showCLIM = 1;
67			%N = 256;
68			%c = [10 10];
69			%cx = [-(1+2c(1)) 1+2c(2)]*6e-5; cmap = mapanom;
70			unit = '1/s';
71			end %switch what to plot
72
73
74			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
75			% Draw variable:
76			sp=subplot(iw,jw,ip);hold on
77			if CONT ~= 1
78			m_pcolor(Clon,Clat,C);
79			shading(SHADE);
80			colormap(logcolormap(N,c(1),c(2),cmap));
81			else
82			[cs,h] = m_contourf(Clon,Clat,C,CONTv);
83			colormap(mycolormap(logcolormap(N,c(1),c(2),cmap),length(CONTv)));
84			if CONTshlab
85			clabel(cs,h,'labelspacing',200,'fontsize',8)
86			end
87			end
88			caxis(cx);
89			if ip == 2
90			ccol = colorbar(CBAR,'fontsize',10);
91			ctitle(ccol,unit);
92			posiC = get(ccol,'position');
93			set(ccol,'position',[.2 posiC(2) 1-2*.2 .02]);
94			end
95			title(tit);
96			m_coast('patch',colorCOAST);
97			m_grid('xtick',360-[20:5:80],'ytick',[20:2:50]);
98			set(gcf,'name',titf);
99
100
101			if showT
102			[cs,h] = m_contour(Tlon,Tlat,squeeze(T(1,:,:)),Tcontour);
103			clabel(cs,h,'fontsize',8,'color',[0 0 0],'labelspacing',200);
104			for ih=1:length(h)
105			set(h(ih),'edgecolor',Tcolor,'linewidth',1);
106			end
107			end %if show THETA contours
108
109			if showW
110			dx = 10diff(Txlon(1:2)); dy = 8diff(Txlat(1:2));
111			dx = 20diff(Txlon(1:2)); dy = 10diff(Txlat(1:2));
112			lo = [Txlon(1):dx:Txlon(length(Txlon))];
113			la = [Txlat(1):dy:Txlat(length(Txlat))];
114			[lo la] = meshgrid(lo,la);
115			Txn = interp2(Txlat,Txlon,Tx',la,lo);
116			Tyn = interp2(Txlat,Txlon,Ty',la,lo);
117			s = 2;
118			m_quiver(lo,la,Txn,Tyn,s,colorW,'linewidth',1.25);
119			% m_quiver(lo,la,-(1+sin(lapi/180)).Txn,(1+sin(lapi/180)).Tyn,s,'w');
120			m_quiver(360-84,47,1,0,s,'w','linewidth',1.25);
121			m_text(360-84,48,'1 N/m^2','color','w');
122			end %if show windstress
123
124			if showH
125			%[cs,h] = m_contour(MLDlon,MLDlat,MLD,Hcontour);
126			cm = flipud(mycolormap(jet,length(Hcontour)));
127			cm = mycolormap([linspace(0,0,20); linspace(1,.5,20) ;linspace(0,0,20)]',length(Hcontour));
128			cm = [0 1 0 ; 0 .6 0 ; 0 .2 0 ; 0 0 0];
129			for ii = 1 : length(Hcontour)
130			[cs,h] = m_contour(MLDlon,MLDlat,MLD,[1 1]*Hcontour(ii));
131			if ~isempty(cs)
132			% clabel(cs,h,'fontsize',8,'color',[0 0 0],'labelspacing',300);
133			clabel(cs,h,'fontsize',8,'color',cm(ii,:),'labelspacing',600,'fontweight','bold');
134			for ih=1:length(h)
135			% set(h(ih),'edgecolor',Hcolor,'linewidth',1);
136			set(h(ih),'edgecolor',cm(ii,:),'linewidth',1.2);
137			end
138			end
139			end
140			end %if show Mixed Layer depth
141
142			if showCLIM
143			m_line(360-[71 62 62 71 71],[36 36 40.5 40.5 36],'color','r','linewidth',1.5)
144			end
145
146			%suptitle(strcat('Relative vorticity component at depth:',num2str(OYdpt(iz)),'m'));
147
148
149			end %for ip
150			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
151
152			if 1 % Show the date in big in the upper left corner
153			spp=subplot('position',[0 .95 .25 .05]);
154			p=patch([0 1 1 0],[0 0 1 1],'w');
155			set(spp,'ytick',[],'xtick',[]);
156			set(spp,'box','off');
157			text(0.1,.5,num2str(TIME(it,:)),'fontsize',16,...
158			'fontweight','bold','color','r','verticalalign','middle');
159			end
160
161
162			%%%%%%%%%%%%%%%%
163			drawnow
164			set(gcf,'position',[4 48 888 430]);
165			videotimeline(num2str(zeros(size(TIME,1),1)),it,'b')
166			%videotimeline(TIME,it,'b')
167			if prtimg
168			set(gcf,'color','white')
169			set(findobj('tag','m_grid_color'),'facecolor','none')
170			set(gcf,'paperposition',[0.6 6.5 25 14]);
171			exportj(gcf,1,strcat(outimg,sla,titf,'.',snapshot));
172			end %if
173