Context Navigation

source: trunk/src/uvmat_doc/FUNCTIONS_DOC/plot_field.html @ 37

Last change on this file since 37 was 37, checked in by sommeria, 14 years ago
create_grid.fig , uvmat_doc and all the included files added
File size: 92.4 KB

Line
1	<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
2	"http://www.w3.org/TR/REC-html40/loose.dtd">
3	<html>
4	<head>
5	<title>Description of plot_field</title>
6	<meta name="keywords" content="plot_field">
7	<meta name="description" content="'plot_field': plot any field with the structure defined in the uvmat package">
8	<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
9	<meta name="generator" content="m2html © 2003 Guillaume Flandin">
10	<meta name="robots" content="index, follow">
11	<link type="text/css" rel="stylesheet" href="../m2html.css">
12	</head>
13	<body>
14	<a name="_top"></a>
15	<div><a href="../index.html">Home</a> > <a href="index.html">.</a> > plot_field.m</div>
16
17	<!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png"> Master index</a></td>
18	<td align="right"><a href="index.html">Index for . <img alt=">" border="0" src="../right.png"></a></td></tr></table>-->
19
20	<h1>plot_field
21	</h1>
22
23	<h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
24	<div class="box"><strong>'plot_field': plot any field with the structure defined in the uvmat package</strong></div>
25
26	<h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
27	<div class="box"><strong>function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar) </strong></div>
28
29	<h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
30	<div class="fragment"><pre class="comment">'plot_field': plot any field with the structure defined in the uvmat package
31	------------------------------------------------------------------------
32
33	This function is used by uvmat to plot fields. It automatically chooses the representation
34	appropriate to the input field structure:
35	2D vector fields are represented by arrows, 2D scalar fiedlds by grey scale images or contour plots, 1D fields are represented by usual plot with (abscissa, ordinate).
36	The input field structure is first tested by check_field_structure.m,
37	then split into blocks of related variables by find_field_indices.m.
38	The dimensionality of each block is obtained by this fuction
39	considering the presence of variables with the attribute .Role='coord_x'
40	and/or coord_y and/or coord_z (case of unstructured coordinates), or
41	dimension variables (case of matrices).
42
43	function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)
44
45	OUPUT:
46	PlotType: type of plot: 'text','line'(curve plot),'plane':2D view,'volume'
47	PlotParamOut: structure, representing the updated plotting parameters, in case of automatic scaling
48	haxes: handle of the plotting axis, when a new figure is created.
49
50	INPUT
51	Data: structure describing the field to plot
52	(optional) .ListGlobalAttribute: cell listing the names of the global attributes
53	.Att_1,Att_2... : values of the global attributes
54	(requested) .ListVarName: list of variable names to select (cell array of char strings {'VarName1', 'VarName2',...} )
55	(requested) .VarDimName: list of dimension names for each element of .ListVarName (cell array of string cells)
56	.VarAttribute: cell of attributes for each element of .ListVarName (cell array of structures of the form VarAtt.key=value)
57	(requested) .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName
58
59	Variable attribute .Role :
60	The only variable attribute used for plotting purpose is .Role which can take
61	the values
62	Role = 'scalar': (default) represents a scalar field
63	= 'coord_x', 'coord_y', 'coord_z': represents a separate set of
64	unstructured coordinate x, y or z
65	= 'vector': represents a vector field whose number of components
66	is given by the last dimension (called 'nb_dim')
67	= 'vector_x', 'vector_y', 'vector_z' :represents the x, y or z component of a vector
68	= 'warnflag' : provides a warning flag about the quality of data in a 'Field', default=0, no warning
69	= 'errorflag': provides an error flag marking false data,
70	default=0, no error. Different non zero values can represent different criteria of elimination.</pre></div>
71
72	<!-- crossreference -->
73	<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
74	This function calls:
75	<ul style="list-style-image:url(../matlabicon.gif)">
76	<li><a href="check_field_structure.html" class="code" title="function [DataOut,errormsg]=check_field_structure(Data)">check_field_structure</a> 'check_field_structure': check the representation of fields by a Matlab structure</li><li><a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a> 'find_file_indices': group the variables of a nc-formated Matlab structure into 'fields' with common dimensions</li><li><a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a> 'msgbox_uvmat': associated with GUI msgbox_uvmat.fig to display message boxes, for error, warning or input calls</li><li><a href="plot_text.html" class="code" title="function hdisplay=plot_text(FieldData,hdisplay_in)">plot_text</a> 'plot_text': function for displaying the content of a Matlab structure in a figure</li><li><a href="proj_grid.html" class="code" title="function [A,rangx,rangy]=proj_grid(vec_X,vec_Y,vec_A,rgx_in,rgy_in,npxy_in)">proj_grid</a> 'proj_grid': project fields with unstructured coordinantes on a regular grid</li><li><a href="read_plot_param.html" class="code" title="function PlotParam=read_plot_param(handles)">read_plot_param</a> 'read_plot_param': read plotting parameters from the interface uvmat</li><li><a href="set_col_vec.html" class="code" title="function [colorlist,col_vec,colcode_out]=set_col_vec(colcode,vec_C)">set_col_vec</a> 'set_col_vec': sets the color code for vectors depending on a scalar vec_C and parameters given by the struct colcode</li><li><a href="set_title.html" class="code" title="function TITLE=set_title(Style,ProjMode)">set_title</a> 'set_title': defines the 'TITLE' of a projection object</li><li><a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a> 'warndlg_uvmat': display warning message (error, warning, confirmation) in a given figure</li></ul>
77	This function is called by:
78	<ul style="list-style-image:url(../matlabicon.gif)">
79	<li><a href="geometry_calib.html" class="code" title="function varargout = geometry_calib(varargin)">geometry_calib</a> 'geometry_calib': performs geometric calibration from a set of reference points</li><li><a href="get_field.html" class="code" title="function varargout = get_field(varargin)">get_field</a> 'get_field': display variables and attributes from a Netcdf file, and plot selected fields</li><li><a href="mouse_down.html" class="code" title="function xy=mouse_down(hObject,eventdata)">mouse_down</a> 'mouse_down': function activated when the mouse button is pressed on a figure (callback for 'WindowButtonDownFcn'</li><li><a href="probe_calib.html" class="code" title="function varargout = probe_calib(varargin)">probe_calib</a> 'probe_calib': performs geometric calibration from a set of reference points</li><li><a href="series.html" class="code" title="function varargout = series(varargin)">series</a> 'series': master function associated to the GUI series.m for analysis field series</li><li><a href="update_obj.html" class="code" title="function Object_out=update_obj(UvData,IndexObj,ObjectData,PlotHandles)">update_obj</a> 'update_obj': update the object graph representation and its projection field, record it in the uvmat interface</li><li><a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a> 'uvmat': function associated with the GUI 'uvmat.fig' for images and data field visualization</li></ul>
80	<!-- crossreference -->
81
82	<h2><a name="_subfunctions"></a>SUBFUNCTIONS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
83	<ul style="list-style-image:url(../matlabicon.gif)">
84	<li><a href="#_sub1" class="code">function [AxeData,haxes]=plot_hist(Data,haxes,PlotParam)</a></li><li><a href="#_sub2" class="code">function [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)</a></li><li><a href="#_sub3" class="code">function [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)</a></li><li><a href="#_sub4" class="code">function quiresetn(haxes,x,y,u,v,scale,colorlist,col_vec)</a></li><li><a href="#_sub5" class="code">function YTick=colbartick(MinA,MaxA)</a></li></ul>
85	<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
86	<div class="fragment"><pre>0001 <span class="comment">%'plot_field': plot any field with the structure defined in the uvmat package</span>
87	0002 <span class="comment">%------------------------------------------------------------------------</span>
88	0003 <span class="comment">%</span>
89	0004 <span class="comment">% This function is used by uvmat to plot fields. It automatically chooses the representation</span>
90	0005 <span class="comment">% appropriate to the input field structure:</span>
91	0006 <span class="comment">% 2D vector fields are represented by arrows, 2D scalar fiedlds by grey scale images or contour plots, 1D fields are represented by usual plot with (abscissa, ordinate).</span>
92	0007 <span class="comment">% The input field structure is first tested by check_field_structure.m,</span>
93	0008 <span class="comment">% then split into blocks of related variables by find_field_indices.m.</span>
94	0009 <span class="comment">% The dimensionality of each block is obtained by this fuction</span>
95	0010 <span class="comment">% considering the presence of variables with the attribute .Role='coord_x'</span>
96	0011 <span class="comment">% and/or coord_y and/or coord_z (case of unstructured coordinates), or</span>
97	0012 <span class="comment">% dimension variables (case of matrices).</span>
98	0013 <span class="comment">%</span>
99	0014 <span class="comment">% function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)</span>
100	0015 <span class="comment">%</span>
101	0016 <span class="comment">% OUPUT:</span>
102	0017 <span class="comment">% PlotType: type of plot: 'text','line'(curve plot),'plane':2D view,'volume'</span>
103	0018 <span class="comment">% PlotParamOut: structure, representing the updated plotting parameters, in case of automatic scaling</span>
104	0019 <span class="comment">% haxes: handle of the plotting axis, when a new figure is created.</span>
105	0020 <span class="comment">%</span>
106	0021 <span class="comment">%INPUT</span>
107	0022 <span class="comment">% Data: structure describing the field to plot</span>
108	0023 <span class="comment">% (optional) .ListGlobalAttribute: cell listing the names of the global attributes</span>
109	0024 <span class="comment">% .Att_1,Att_2... : values of the global attributes</span>
110	0025 <span class="comment">% (requested) .ListVarName: list of variable names to select (cell array of char strings {'VarName1', 'VarName2',...} )</span>
111	0026 <span class="comment">% (requested) .VarDimName: list of dimension names for each element of .ListVarName (cell array of string cells)</span>
112	0027 <span class="comment">% .VarAttribute: cell of attributes for each element of .ListVarName (cell array of structures of the form VarAtt.key=value)</span>
113	0028 <span class="comment">% (requested) .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName</span>
114	0029 <span class="comment">%</span>
115	0030 <span class="comment">% Variable attribute .Role :</span>
116	0031 <span class="comment">% The only variable attribute used for plotting purpose is .Role which can take</span>
117	0032 <span class="comment">% the values</span>
118	0033 <span class="comment">% Role = 'scalar': (default) represents a scalar field</span>
119	0034 <span class="comment">% = 'coord_x', 'coord_y', 'coord_z': represents a separate set of</span>
120	0035 <span class="comment">% unstructured coordinate x, y or z</span>
121	0036 <span class="comment">% = 'vector': represents a vector field whose number of components</span>
122	0037 <span class="comment">% is given by the last dimension (called 'nb_dim')</span>
123	0038 <span class="comment">% = 'vector_x', 'vector_y', 'vector_z' :represents the x, y or z component of a vector</span>
124	0039 <span class="comment">% = 'warnflag' : provides a warning flag about the quality of data in a 'Field', default=0, no warning</span>
125	0040 <span class="comment">% = 'errorflag': provides an error flag marking false data,</span>
126	0041 <span class="comment">% default=0, no error. Different non zero values can represent different criteria of elimination.</span>
127	0042
128	0043
129	0044 <span class="comment">% additional elements characterizing the projection object (should not be necessary)--</span>
130	0045 <span class="comment">% Data.Style : style of projection object</span>
131	0046 <span class="comment">% Data.XObject,.YObject: set of coordinates defining the object position;</span>
132	0047 <span class="comment">% Data.ProjMode=type of projection ;</span>
133	0048 <span class="comment">% Data.ProjAngle=angle of projection;</span>
134	0049 <span class="comment">% Data.DX,.DY,.DZ=increments;</span>
135	0050 <span class="comment">% Data.MaxY,MinY: min and max Y</span>
136	0051
137	0052 <span class="comment">% haxes: handle of the plotting axes to update with the new plot. If this input is absent or not a valid axes handle, a new figure is created.</span>
138	0053 <span class="comment">%</span>
139	0054 <span class="comment">% PlotParam: parameters for plotting, as read on the uvmat interface (by function 'read_plot_param.m')</span>
140	0055 <span class="comment">% --scalars--</span>
141	0056 <span class="comment">% .Scalar.MaxA: upper bound (saturation color) for the scalar representation, max(field) by default</span>
142	0057 <span class="comment">% .Scalar.MinA: lower bound (saturation) for the scalar representation, min(field) by default</span>
143	0058 <span class="comment">% .Scalar.AutoScal: =1 (default) lower and upper bounds of the scalar representation set to the min and max of the field</span>
144	0059 <span class="comment">% =0 lower and upper bound imposed by .AMax and .MinA</span>
145	0060 <span class="comment">% .Scalar.BW= 1 black and white representation imposed, =0 by default.</span>
146	0061 <span class="comment">% .Scalar.Contours= 1: represent scalars by contour plots (Matlab function 'contour'); =0 by default</span>
147	0062 <span class="comment">% .IncrA : contour interval</span>
148	0063 <span class="comment">% -- vectors--</span>
149	0064 <span class="comment">% .Vectors.VecScale: scale for the vector representation</span>
150	0065 <span class="comment">% .Vectors.AutoVec: =0 (default) automatic length for vector representation, =1: length set by .VecScale</span>
151	0066 <span class="comment">% .Vectors.HideFalse= 0 (default) false vectors represented in magenta, =1: false vectors not represented;</span>
152	0067 <span class="comment">% .Vectors.HideWarning= 0 (default) vectors marked by warnflag~=0 marked in black, 1: no warning representation;</span>
153	0068 <span class="comment">% .Vectors.decimate4 = 0 (default) all vectors reprtesented, =1: half of the vectors represented along each coordinate</span>
154	0069 <span class="comment">% -- vector color--</span>
155	0070 <span class="comment">% .Vectors.ColorCode= 'black','white': imposed color (default ='blue')</span>
156	0071 <span class="comment">% 'rgb', : three colors red, blue, green depending</span>
157	0072 <span class="comment">% on thresholds .colcode1 and .colcode2 on the input scalar value (C)</span>
158	0073 <span class="comment">% 'brg': like rgb but reversed color order (blue, green, red)</span>
159	0074 <span class="comment">% '64 colors': continuous color from blue to red (multijet)</span>
160	0075 <span class="comment">% .Vectors.colcode1 : first threshold for rgb, first value for'continuous'</span>
161	0076 <span class="comment">% .Vectors.colcode2 : second threshold for rgb, last value (saturation) for 'continuous'</span>
162	0077 <span class="comment">% .Vectors.FixedCbounds; =0 (default): the bounds on C representation are min and max, =1: they are fixed by .Minc and .MaxC</span>
163	0078 <span class="comment">% .Vectors.MinC = imposed minimum of the scalar field used for vector color;</span>
164	0079 <span class="comment">% .Vectors.MaxC = imposed maximum of the scalar field used for vector color;</span>
165	0080 <span class="comment">%</span>
166	0081 <span class="comment">% KeepLim:=0 (default) adjust axes limit to the X,Y data, =1: preserves the previous axes limits</span>
167	0082 <span class="comment">% PosColorbar: if not empty, display a colorbar for B&W images</span>
168	0083 <span class="comment">% imposed position of the colorbar (ex [0.821 0.471 0.019 0.445])</span>
169	0084 <span class="comment">%</span>
170	0085 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span>
171	0086 <span class="comment">% Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.</span>
172	0087 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span>
173	0088 <span class="comment">% This file is part of the toolbox UVMAT.</span>
174	0089 <span class="comment">%</span>
175	0090 <span class="comment">% UVMAT is free software; you can redistribute it and/or modify</span>
176	0091 <span class="comment">% it under the terms of the GNU General Public License as published by</span>
177	0092 <span class="comment">% the Free Software Foundation; either version 2 of the License, or</span>
178	0093 <span class="comment">% (at your option) any later version.</span>
179	0094 <span class="comment">%</span>
180	0095 <span class="comment">% UVMAT is distributed in the hope that it will be useful,</span>
181	0096 <span class="comment">% but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
182	0097 <span class="comment">% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</span>
183	0098 <span class="comment">% GNU General Public License (file UVMAT/COPYING.txt) for more details.</span>
184	0099 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span>
185	0100
186	0101 <a name="_sub0" href="#_subfunctions" class="code">function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)</a>
187	0102 <span class="comment">%default output</span>
188	0103 <span class="keyword">if</span> ~exist(<span class="string">'PlotParam'</span>,<span class="string">'var'</span>),PlotParam=[];<span class="keyword">end</span>;
189	0104 <span class="keyword">if</span> ~exist(<span class="string">'KeepLim'</span>,<span class="string">'var'</span>),KeepLim=0;<span class="keyword">end</span>;
190	0105 <span class="keyword">if</span> ~exist(<span class="string">'PosColorbar'</span>,<span class="string">'var'</span>),PosColorbar=[];<span class="keyword">end</span>;
191	0106 PlotType=<span class="string">'text'</span>; <span class="comment">%default</span>
192	0107 PlotParamOut=PlotParam;<span class="comment">%default</span>
193	0108
194	0109 <span class="comment">% check input structure</span>
195	0110 [Data,errormsg]=<a href="check_field_structure.html" class="code" title="function [DataOut,errormsg]=check_field_structure(Data)">check_field_structure</a>(Data);
196	0111
197	0112 <span class="keyword">if</span> ~isempty(errormsg)
198	0113 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,[<span class="string">'input of plot_field/check_field_structure: '</span> errormsg])
199	0114 display([<span class="string">'input of plot_field/check_field_structure:: '</span> errormsg])
200	0115 <span class="keyword">return</span>
201	0116 <span class="keyword">end</span>
202	0117
203	0118 testnewfig=1;<span class="comment">%test to create a new figure (default)</span>
204	0119 testzoomaxes=0;<span class="comment">%test for the existence of a zoom secondary figure attached to the plotting axes</span>
205	0120 <span class="keyword">if</span> exist(<span class="string">'haxes'</span>,<span class="string">'var'</span>)
206	0121 <span class="keyword">if</span> ishandle(haxes)
207	0122 <span class="keyword">if</span> isequal(get(haxes,<span class="string">'Type'</span>),<span class="string">'axes'</span>)
208	0123 <span class="comment">% hfig=get(haxes,'Parent');</span>
209	0124 axes(haxes)
210	0125 testnewfig=0;
211	0126 AxeData=get(haxes,<span class="string">'UserData'</span>);
212	0127 <span class="keyword">if</span> isfield(AxeData,<span class="string">'ZoomAxes'</span>)&& ishandle(AxeData.ZoomAxes)
213	0128 <span class="keyword">if</span> isequal(get(AxeData.ZoomAxes,<span class="string">'Type'</span>),<span class="string">'axes'</span>)
214	0129 testzoomaxes=1;
215	0130 zoomaxes=AxeData.ZoomAxes;
216	0131 <span class="keyword">end</span>
217	0132 <span class="keyword">end</span>
218	0133 <span class="keyword">end</span>
219	0134 <span class="keyword">end</span>
220	0135 <span class="keyword">end</span>
221	0136 <span class="keyword">if</span> testnewfig<span class="comment">% create a new figure and axes if the plotting axes does not exist</span>
222	0137 hfig=figure;
223	0138 <span class="keyword">if</span> isfield(Data,<span class="string">'IndexObj'</span>)&isfield(Data,<span class="string">'Style'</span>)&isfield(Data,<span class="string">'ProjMode'</span>)
224	0139 figname=[num2str(Data.IndexObj) <span class="string">'-'</span> <a href="set_title.html" class="code" title="function TITLE=set_title(Style,ProjMode)">set_title</a>(Data.Style,Data.ProjMode)];
225	0140 set(hfig,<span class="string">'Name'</span>,figname)
226	0141 <span class="keyword">end</span>
227	0142 testhandle=0;
228	0143 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'text_display_1'</span>)& ishandle(PlotParam.text_display_1)
229	0144 set(hfig,<span class="string">'UserData'</span>,PlotParam)
230	0145 testhandle=1;
231	0146 <span class="keyword">end</span>
232	0147 set(hfig,<span class="string">'Units'</span>,<span class="string">'normalized'</span>)
233	0148 set(hfig,<span class="string">'WindowButtonDownFcn'</span>,<span class="string">'mouse_down'</span>)
234	0149 <span class="comment">%set(hfig,'WindowButtonMotionFcn',{'mouse_motion',PlotParam})%set mouse action function</span>
235	0150 set(hfig,<span class="string">'WindowButtonMotionFcn'</span>,<span class="string">'mouse_motion'</span>)<span class="comment">%set mouse action function</span>
236	0151 set(hfig,<span class="string">'WindowButtonUpFcn'</span>,<span class="string">'mouse_up'</span>)<span class="comment">%set mouse action function</span>
237	0152 haxes=axes;
238	0153 set(haxes,<span class="string">'position'</span>,[0.13,0.2,0.775,0.73])
239	0154 <span class="keyword">end</span>
240	0155 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'text_display_1'</span>)& ishandle(PlotParam.text_display_1)
241	0156 PlotParam=<a href="read_plot_param.html" class="code" title="function PlotParam=read_plot_param(handles)">read_plot_param</a>(PlotParam);
242	0157 <span class="keyword">end</span>
243	0158 <span class="keyword">if</span> testnewfig
244	0159 PlotParam.NextPlot=<span class="string">'add'</span>; <span class="comment">%parameter for plot_profile and plot_hist</span>
245	0160 <span class="keyword">end</span>
246	0161 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'Auto_xy'</span>) && isequal(PlotParam.Auto_xy,1)
247	0162 set(haxes,<span class="string">'DataAspectRatioMode'</span>,<span class="string">'auto'</span>)<span class="comment">%automatic aspect ratio</span>
248	0163 <span class="keyword">end</span>
249	0164
250	0165 <span class="comment">% check the cells of fields :</span>
251	0166 testnbdim=1;
252	0167 [CellVarIndex,NbDim,VarType,errormsg]=<a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a>(Data);
253	0168
254	0169 <span class="keyword">if</span> ~isempty(errormsg)
255	0170 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,[<span class="string">'input of plot_field/find_field_indices: '</span> errormsg])
256	0171 display([<span class="string">'input of plot_field: '</span> errormsg])
257	0172 <span class="keyword">return</span>
258	0173 <span class="keyword">end</span>
259	0174 <span class="keyword">if</span> ~isfield(Data,<span class="string">'NbDim'</span>) <span class="comment">%& ~isfield(Data,'Style')%determine the space dimensionb if not defined: choose the kind of plot</span>
260	0175 [Data.NbDim,imax]=max(NbDim);
261	0176 <span class="keyword">end</span>
262	0177 <span class="keyword">if</span> isequal(Data.NbDim,0) <span class="comment">% TODO: chech whether this function is still used, replace by plot_profile ?</span>
263	0178 <span class="keyword">if</span> isfield(Data,<span class="string">'Style'</span>) & isequal(Data.Style,<span class="string">'points'</span>)
264	0179 AxeData=<a href="plot_text.html" class="code" title="function hdisplay=plot_text(FieldData,hdisplay_in)">plot_text</a>(Data,haxes);
265	0180 PlotType=<span class="string">'text'</span>;
266	0181 <span class="keyword">else</span>
267	0182 [AxeData,haxes]=<a href="#_sub1" class="code" title="subfunction [AxeData,haxes]=plot_hist(Data,haxes,PlotParam)">plot_hist</a>(Data,haxes,PlotParam);
268	0183 <span class="keyword">end</span>
269	0184 <span class="keyword">elseif</span> isequal(Data.NbDim,1)
270	0185 [AxeData,haxes]=<a href="#_sub2" class="code" title="subfunction [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)">plot_profile</a>(Data,CellVarIndex,VarType,haxes,PlotParam);<span class="comment">%</span>
271	0186 <span class="keyword">if</span> testzoomaxes
272	0187 [AxeData,zoomaxes,PlotParamOut]=<a href="#_sub2" class="code" title="subfunction [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)">plot_profile</a>(Data,CellVarIndex,VarType,zoomaxes,PlotParam);
273	0188 AxeData.ZoomAxes=zoomaxes;
274	0189 <span class="keyword">end</span>
275	0190 PlotType=<span class="string">'line'</span>;
276	0191 <span class="keyword">elseif</span> isequal(Data.NbDim,2)
277	0192 ind_select=find(NbDim>=2);
278	0193 <span class="keyword">if</span> numel(ind_select)>2
279	0194 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,[<span class="string">'more than two fields to map'</span>])
280	0195 display([<span class="string">'more than two fields to map'</span>])
281	0196 <span class="keyword">return</span>
282	0197 <span class="keyword">end</span>
283	0198 [AxeData,haxes,PlotParamOut,PlotType]=<a href="#_sub3" class="code" title="subfunction [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)">plot_plane</a>(Data,CellVarIndex(ind_select),VarType(ind_select),haxes,PlotParam,KeepLim,PosColorbar);
284	0199 <span class="keyword">if</span> testzoomaxes
285	0200 [AxeData,zoomaxes,PlotParamOut]=<a href="#_sub3" class="code" title="subfunction [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)">plot_plane</a>(Data,CellVarIndex(ind_select),VarType(ind_select),zoomaxes,PlotParam,1,PosColorbar);
286	0201 AxeData.ZoomAxes=zoomaxes;
287	0202 <span class="keyword">end</span>
288	0203 <span class="keyword">elseif</span> isequal(Data.NbDim,3)
289	0204 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'volume plot not implemented yet'</span>)
290	0205 <span class="keyword">return</span>
291	0206 <span class="comment">%plot_volume(haxes,Data,PlotParam)% A FAIRE</span>
292	0207 <span class="comment">%PlotType='volume';</span>
293	0208 <span class="keyword">else</span>
294	0209 testnbdim=0;
295	0210 <span class="keyword">end</span>
296	0211
297	0212 <span class="comment">%display (or delete) error message</span>
298	0213 htext=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'hTxt'</span>);
299	0214 <span class="keyword">if</span> isfield(Data,<span class="string">'Txt'</span>)
300	0215 <span class="keyword">if</span> isempty(htext)
301	0216 Xlim=get(haxes,<span class="string">'XLim'</span>);
302	0217 Ylim=get(haxes,<span class="string">'YLim'</span>);
303	0218 htext=text(Xlim(1),(Ylim(1)+Ylim(2))/2,Data.Txt,<span class="string">'Tag'</span>,<span class="string">'hTxt'</span>,<span class="string">'Color'</span>,<span class="string">'r'</span>);
304	0219 set(htext,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)
305	0220 <span class="keyword">else</span>
306	0221 set(htext,<span class="string">'String'</span>,Data.Txt)
307	0222 <span class="keyword">end</span>
308	0223 <span class="keyword">elseif</span> ~isempty(htext)
309	0224 delete(htext)
310	0225 <span class="keyword">end</span>
311	0226
312	0227 <span class="comment">% set graph aspect ratio</span>
313	0228
314	0229
315	0230 set(haxes,<span class="string">'UserData'</span>,AxeData)
316	0231 <span class="comment">%set(haxes,'Tag','uvmat');</span>
317	0232
318	0233
319	0234 <span class="comment">%-------------------------------------------</span>
320	0235 <a name="_sub1" href="#_subfunctions" class="code">function [AxeData,haxes]=plot_hist(Data,haxes,PlotParam)</a><span class="comment">% TODO: chech whether this function is still used, replace by plot_profile ?</span>
321	0236 <span class="comment">%------------------------------------------</span>
322	0237 AxeData=get(haxes,<span class="string">'UserData'</span>); <span class="comment">%defau</span>
323	0238 hfig=get(haxes,<span class="string">'parent'</span>);
324	0239 <span class="keyword">if</span> ~isfield(Data,<span class="string">'ListVarName'</span>)
325	0240 <span class="keyword">return</span>
326	0241 <span class="keyword">end</span>
327	0242 ColorOrder=[1 0 0;0 0.5 0;0 0 1;0 0.75 0.75;0.75 0 0.75;0.75 0.75 0;0.25 0.25 0.25];
328	0243 set(haxes,<span class="string">'ColorOrder'</span>,ColorOrder)
329	0244 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'NextPlot'</span>)
330	0245 set(haxes,<span class="string">'NextPlot'</span>,PlotParam.NextPlot)
331	0246 <span class="keyword">end</span>
332	0247 charplot=<span class="string">'''-'''</span>;
333	0248 iplot=0;
334	0249 legend_str={};
335	0250 label_str=<span class="string">''</span>;
336	0251 textmean={};
337	0252 plotstr=<span class="string">'plot('</span>;
338	0253 <span class="keyword">for</span> ilist=1:length(Data.ListVarName)
339	0254 VarName=Data.ListVarName{ilist};
340	0255 eval([<span class="string">'['</span> VarName <span class="string">'hist,'</span> VarName <span class="string">'val]=hist(double(Data.'</span> VarName <span class="string">'),100);'</span>]);<span class="comment">%coordinate variable set as c</span>
341	0256 plotstr=[plotstr VarName <span class="string">'val,'</span> VarName <span class="string">'hist,'</span> charplot <span class="string">','</span>];
342	0257 eval([<span class="string">'nbcomponent2=size(Data.'</span> VarName <span class="string">',2);'</span>]);
343	0258 eval([<span class="string">'nbcomponent1=size(Data.'</span> VarName <span class="string">',1);'</span>]);
344	0259 eval([<span class="string">'varmean=mean(double(Data.'</span> VarName <span class="string">'));'</span>]);<span class="comment">%mean value</span>
345	0260 textmean=[textmean; {[VarName <span class="string">'mean= '</span> num2str(varmean,4)]}];
346	0261 <span class="keyword">if</span> nbcomponent1==1\| nbcomponent2==1
347	0262 legend_str=[legend_str {VarName}]; <span class="comment">%variable with one component</span>
348	0263 <span class="keyword">else</span>
349	0264 <span class="keyword">for</span> ic=1:min(nbcomponent1,nbcomponent2)
350	0265 legend_str=[legend_str [VarName <span class="string">'_'</span> num2str(ic)]]; <span class="comment">%variable with severals components</span>
351	0266 <span class="comment">% labeled by their index (e.g. color component)</span>
352	0267 <span class="keyword">end</span>
353	0268 <span class="keyword">end</span>
354	0269 label_str=[label_str <span class="string">' '</span> VarName];
355	0270 <span class="keyword">end</span>
356	0271 <span class="keyword">if</span> ~isequal(plotstr,<span class="string">'plot('</span>)
357	0272 plotstr(end)=<span class="string">')'</span>;
358	0273 eval(plotstr)
359	0274 hlegend=findobj(hfig,<span class="string">'Tag'</span>,<span class="string">'legend'</span>);<span class="comment">%find existing legend on the plot</span>
360	0275 <span class="keyword">if</span> ~isempty(hlegend)
361	0276 legend_old=get(hlegend,<span class="string">'String'</span>);
362	0277 <span class="keyword">if</span> isequal(size(legend_old,1),size(legend_str,1))
363	0278 legend_str=[legend_old legend_str];
364	0279 <span class="keyword">end</span>
365	0280 <span class="keyword">end</span>
366	0281 legend(legend_str)
367	0282 xlabel(label_str)
368	0283 ylabel(<span class="string">'nb values'</span>)
369	0284 grid on
370	0285 title_str=<span class="string">''</span>;
371	0286 <span class="keyword">if</span> isfield(Data,<span class="string">'filename'</span>)
372	0287 [Path, title_str, ext]=fileparts(Data.filename);
373	0288 title_str=[title_str ext];
374	0289 <span class="keyword">end</span>
375	0290 <span class="keyword">if</span> isfield(Data,<span class="string">'Action'</span>)
376	0291 <span class="keyword">if</span> ~isequal(title_str,<span class="string">''</span>)
377	0292 title_str=[title_str <span class="string">', '</span>];
378	0293 <span class="keyword">end</span>
379	0294 title_str=[title_str Data.Action];
380	0295 <span class="keyword">end</span>
381	0296 htitle=title(title_str);
382	0297 set(htitle,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span>
383	0298 hlist=findobj(gcf,<span class="string">'Style'</span>,<span class="string">'listbox'</span>);
384	0299 <span class="keyword">if</span> isempty(hlist)
385	0300 uicontrol(<span class="string">'Style'</span>,<span class="string">'popupmenu'</span>,<span class="string">'Position'</span>,[20 20 200 20],<span class="string">'String'</span>,textmean);
386	0301 <span class="keyword">else</span>
387	0302 set(hlist(1),<span class="string">'String'</span>,textmean)
388	0303 <span class="keyword">end</span>
389	0304 <span class="keyword">end</span>
390	0305 AxeData=Data;
391	0306
392	0307
393	0308 <span class="comment">%----------------------------------------------------------</span>
394	0309 <a name="_sub2" href="#_subfunctions" class="code">function [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)</a>
395	0310 <span class="comment">%-----------------------------------------------------------</span>
396	0311 axes(haxes)
397	0312 hfig=get(haxes,<span class="string">'parent'</span>);
398	0313 AxeData=data;
399	0314 ColorOrder=[1 0 0;0 0.5 0;0 0 1;0 0.75 0.75;0.75 0 0.75;0.75 0.75 0;0.25 0.25 0.25];
400	0315 set(haxes,<span class="string">'ColorOrder'</span>,ColorOrder)
401	0316 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'NextPlot'</span>)
402	0317 set(haxes,<span class="string">'NextPlot'</span>,PlotParam.NextPlot)
403	0318 <span class="keyword">end</span>
404	0319 legend_str={};
405	0320
406	0321 <span class="comment">%initiates string of the plot command</span>
407	0322 plotstr=<span class="string">'plot('</span>;
408	0323 textmean={};
409	0324 abscissa_name=<span class="string">''</span>;
410	0325 coord_x_index=[];
411	0326 <span class="keyword">for</span> icell=1:length(CellVarIndex)
412	0327 testfalse=0;
413	0328 VarIndex=CellVarIndex{icell};<span class="comment">% indices of the selected variables in the list data.ListVarName</span>
414	0329 DimCell=data.VarDimName{VarIndex(1)};
415	0330 <span class="keyword">if</span> ischar(DimCell)
416	0331 DimCell={DimCell};
417	0332 <span class="keyword">end</span>
418	0333 XName=DimCell{1}; <span class="comment">%first dimension considered as abscissa</span>
419	0334 coord_x_index_cell=VarType{icell}.coord(1);
420	0335 <span class="keyword">if</span> isequal(coord_x_index_cell,0)
421	0336 <span class="keyword">continue</span> <span class="comment">% the cell has no abscissa, skip it</span>
422	0337 <span class="keyword">end</span>
423	0338 <span class="keyword">if</span> ~isempty(coord_x_index)&&~isequal(coord_x_index_cell,coord_x_index)
424	0339 <span class="keyword">continue</span> <span class="comment">%all the selected variables must have the same first dimension</span>
425	0340 <span class="keyword">else</span>
426	0341 coord_x_index=coord_x_index_cell;
427	0342 <span class="keyword">end</span>
428	0343 testplot=ones(size(data.ListVarName));<span class="comment">%default test for plotted variables</span>
429	0344 testcoordvar=0;
430	0345 charplot_0=<span class="string">'''-'''</span>;<span class="comment">%default</span>
431	0346 <span class="keyword">if</span> isfield(data,<span class="string">'ObjectProjMode'</span>)& isequal(data.ObjectProjMode,<span class="string">'projection'</span>)
432	0347 charplot_0=<span class="string">'''+'''</span>;
433	0348 <span class="keyword">end</span>
434	0349 xtitle=<span class="string">''</span>;
435	0350
436	0351 xtitle=data.ListVarName{coord_x_index};
437	0352 eval([<span class="string">'coord_x{icell}=data.'</span> data.ListVarName{coord_x_index} <span class="string">';'</span>]);<span class="comment">%coordinate variable set as coord_x</span>
438	0353 <span class="keyword">if</span> isfield(data,<span class="string">'VarAttribute'</span>)&& numel(data.VarAttribute)>=coord_x_index && isfield(data.VarAttribute{coord_x_index},<span class="string">'units'</span>)
439	0354 xtitle=[xtitle <span class="string">'('</span> data.VarAttribute{coord_x_index}.units <span class="string">')'</span>];
440	0355 <span class="keyword">end</span>
441	0356 eval([<span class="string">'coord_x{icell}=data.'</span> data.ListVarName{coord_x_index} <span class="string">';'</span>]);<span class="comment">%coordinate variable set as coord_x</span>
442	0357 testcoordvar=1;
443	0358 testplot(coord_x_index)=0;
444	0359 <span class="keyword">if</span> ~isempty(VarType{icell}.ancillary')
445	0360 testplot(VarType{icell}.ancillary)=0;
446	0361 <span class="keyword">end</span>
447	0362 <span class="keyword">if</span> ~isempty(VarType{icell}.warnflag')
448	0363 testplot(VarType{icell}.warnflag)=0;
449	0364 <span class="keyword">end</span>
450	0365 <span class="keyword">if</span> ~isempty(VarType{icell}.discrete')
451	0366 charplot_0=<span class="string">'''+'''</span>;
452	0367 <span class="keyword">else</span>
453	0368 charplot_0=<span class="string">'''-'''</span>;
454	0369 <span class="keyword">end</span>
455	0370 <span class="keyword">if</span> testcoordvar==0
456	0371 coord_x{icell}=[1:data.DimValue(DimIndices(1))];<span class="comment">%abscissa by default if no coordinate variable</span>
457	0372 <span class="comment">% charplot_0='''-''';</span>
458	0373 <span class="keyword">end</span>
459	0374 <span class="keyword">if</span> isfield(data,<span class="string">'VarAttribute'</span>)
460	0375 VarAttribute=data.VarAttribute;
461	0376 <span class="keyword">for</span> ivar=1:length(VarIndex)
462	0377 <span class="keyword">if</span> length(VarAttribute)>=VarIndex(ivar) & isfield(VarAttribute{VarIndex(ivar)},<span class="string">'long_name'</span>)
463	0378 plotname{VarIndex(ivar)}=VarAttribute{VarIndex(ivar)}.long_name;
464	0379 <span class="keyword">else</span>
465	0380 plotname{VarIndex(ivar)}=data.ListVarName{VarIndex(ivar)};<span class="comment">%name for display in plot A METTRE</span>
466	0381 <span class="keyword">end</span>
467	0382 <span class="keyword">end</span>
468	0383 <span class="keyword">end</span>
469	0384 <span class="keyword">for</span> ivar=1:length(VarIndex)
470	0385 <span class="keyword">if</span> testplot(VarIndex(ivar))
471	0386 VarName=data.ListVarName{VarIndex(ivar)};
472	0387 eval([<span class="string">'data.'</span> VarName <span class="string">'=squeeze(data.'</span> VarName <span class="string">');'</span>])
473	0388 <span class="keyword">if</span> isequal(VarName,<span class="string">'A'</span>)
474	0389 charplot=<span class="string">'''-'''</span>;
475	0390 <span class="keyword">else</span>
476	0391 charplot=charplot_0;
477	0392 <span class="keyword">end</span>
478	0393 plotstr=[plotstr <span class="string">'coord_x{'</span> num2str(icell) <span class="string">'},data.'</span> VarName <span class="string">','</span> charplot <span class="string">','</span>];
479	0394 eval([<span class="string">'nbcomponent2=size(data.'</span> VarName <span class="string">',2);'</span>]);
480	0395 eval([<span class="string">'nbcomponent1=size(data.'</span> VarName <span class="string">',1);'</span>]);
481	0396 eval([<span class="string">'varmean=mean(double(data.'</span> VarName <span class="string">'));'</span>]);<span class="comment">%mean value</span>
482	0397 textmean=[textmean; {[VarName <span class="string">'mean= '</span> num2str(varmean,4)]}];
483	0398 <span class="keyword">if</span> nbcomponent1==1\| nbcomponent2==1
484	0399 legend_str=[legend_str {VarName}]; <span class="comment">%variable with one component</span>
485	0400 <span class="keyword">else</span> <span class="comment">%variable with severals components</span>
486	0401 <span class="keyword">for</span> ic=1:min(nbcomponent1,nbcomponent2)
487	0402 legend_str=[legend_str [VarName <span class="string">'_'</span> num2str(ic)]]; <span class="comment">%variable with severals components</span>
488	0403 <span class="keyword">end</span> <span class="comment">% labeled by their index (e.g. color component)</span>
489	0404 <span class="keyword">end</span>
490	0405 <span class="keyword">end</span>
491	0406 <span class="keyword">end</span>
492	0407 <span class="keyword">end</span>
493	0408 <span class="keyword">if</span> ~isequal(plotstr,<span class="string">'plot('</span>)
494	0409 plotstr(end)=<span class="string">')'</span>;
495	0410 <span class="comment">%execute plot (instruction plotstr)</span>
496	0411 eval(plotstr)
497	0412 <span class="comment">%%%%%</span>
498	0413 grid on
499	0414 hxlabel=xlabel(xtitle);
500	0415 set(hxlabel,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span>
501	0416 <span class="keyword">if</span> length(legend_str)>=1
502	0417 hylabel=ylabel(legend_str{end});
503	0418 set(hylabel,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span>
504	0419 <span class="keyword">end</span>
505	0420 <span class="keyword">if</span> ~isempty(legend_str)
506	0421 hlegend=findobj(hfig,<span class="string">'Tag'</span>,<span class="string">'legend'</span>);
507	0422 <span class="keyword">if</span> isempty(hlegend)
508	0423 hlegend=legend(legend_str);
509	0424 txt=ver;
510	0425 Release=txt(1).Release;
511	0426 relnumb=str2num(Release(3:4));
512	0427 <span class="keyword">if</span> relnumb >= 14
513	0428 set(hlegend,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span>
514	0429 <span class="keyword">end</span>
515	0430 <span class="keyword">else</span>
516	0431 legend_old=get(hlegend,<span class="string">'String'</span>);
517	0432 <span class="keyword">if</span> isequal(size(legend_old,1),size(legend_str,1))&~isequal(legend_old,legend_str)
518	0433 set(hlegend,<span class="string">'String'</span>,[legend_old legend_str]);
519	0434 <span class="keyword">end</span>
520	0435 <span class="keyword">end</span>
521	0436 <span class="keyword">end</span>
522	0437 title_str=<span class="string">''</span>;
523	0438 <span class="keyword">if</span> isfield(data,<span class="string">'filename'</span>)
524	0439 [Path, title_str, ext]=fileparts(data.filename);
525	0440 title_str=[title_str ext];
526	0441 <span class="keyword">end</span>
527	0442 <span class="keyword">if</span> isfield(data,<span class="string">'Action'</span>)
528	0443 <span class="keyword">if</span> ~isequal(title_str,<span class="string">''</span>)
529	0444 title_str=[title_str <span class="string">', '</span>];
530	0445 <span class="keyword">end</span>
531	0446 title_str=[title_str data.Action];
532	0447 <span class="keyword">end</span>
533	0448 htitle=title(title_str);
534	0449 txt=ver;
535	0450 Release=txt(1).Release;
536	0451 relnumb=str2num(Release(3:4));
537	0452 <span class="keyword">if</span> relnumb >= 14
538	0453 set(htitle,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span>
539	0454 <span class="keyword">end</span>
540	0455 <span class="comment">% A REPRENDRE Mean</span>
541	0456 <span class="comment">% hlist=findobj(gcf,'Style','listbox','Tag','liststat');</span>
542	0457 <span class="comment">% if isempty(hlist)</span>
543	0458 <span class="comment">% 'text'</span>
544	0459 <span class="comment">% textmean</span>
545	0460 <span class="comment">% set(gca,'position',[0.13,0.2,0.775,0.73])</span>
546	0461 <span class="comment">% uicontrol('Style','popupmenu','Position',[20 20 200 20],'String',textmean,'Tag','liststat');</span>
547	0462 <span class="comment">% else</span>
548	0463 <span class="comment">% set(hlist(1),'String',textmean)</span>
549	0464 <span class="comment">% end</span>
550	0465 <span class="keyword">end</span>
551	0466
552	0467
553	0468 <span class="comment">%---------------------------------------</span>
554	0469 <span class="comment">% plot_plane</span>
555	0470 <span class="comment">%----------------------------------------</span>
556	0471 <a name="_sub3" href="#_subfunctions" class="code">function [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)</a>
557	0472
558	0473 <span class="comment">%default plotting parameters</span>
559	0474 PlotType=<span class="string">'plane'</span>;<span class="comment">%default</span>
560	0475 <span class="keyword">if</span> ~exist(<span class="string">'PlotParam'</span>,<span class="string">'var'</span>)
561	0476 PlotParam=[];
562	0477 <span class="keyword">end</span>
563	0478 <span class="keyword">if</span> ~isfield(PlotParam,<span class="string">'Scalar'</span>)
564	0479 PlotParam.Scalar=[];
565	0480 <span class="keyword">end</span>
566	0481 <span class="keyword">if</span> ~isfield(PlotParam,<span class="string">'Vectors'</span>)
567	0482 PlotParam.Vectors=[];
568	0483 <span class="keyword">end</span>
569	0484 PlotParamOut=PlotParam;<span class="comment">%default</span>
570	0485
571	0486 <span class="comment">%plotting axes</span>
572	0487 hfig=get(haxes,<span class="string">'parent'</span>);
573	0488 hcol=findobj(hfig,<span class="string">'Tag'</span>,<span class="string">'Colorbar'</span>); <span class="comment">%look for colorbar axes</span>
574	0489 hima=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'ima'</span>);<span class="comment">% search existing image in the current axes</span>
575	0490 AxeData=get(haxes,<span class="string">'UserData'</span>); <span class="comment">%default</span>
576	0491 <span class="keyword">if</span> ~isstruct(AxeData)<span class="comment">% AxeData must be a structure</span>
577	0492 AxeData=[];
578	0493 <span class="keyword">end</span>
579	0494 AxeData.NbDim=2;
580	0495 <span class="keyword">if</span> isfield(Data,<span class="string">'ObjectCoord'</span>)
581	0496 AxeData.ObjectCoord=Data.ObjectCoord;
582	0497 <span class="keyword">end</span>
583	0498
584	0499 test_ima=0; <span class="comment">%default: test for image or map plot</span>
585	0500 test_vec=0; <span class="comment">%default: test for vector plots</span>
586	0501 test_black=0;
587	0502 test_false=0;
588	0503 test_C=0;
589	0504 <span class="keyword">for</span> icell=1:length(CellVarIndex) <span class="comment">% length(CellVarIndex) =1 or 2 (from the calling function)</span>
590	0505 <span class="comment">% VarIndex=CellVarIndex{icell};</span>
591	0506 VarType=VarTypeCell{icell};
592	0507 ivar_X=VarType.coord_x; <span class="comment">% defines (unique) index for the variable representing unstructured x coordinate (default =[])</span>
593	0508 ivar_Y=VarType.coord_y; <span class="comment">% defines (unique)index for the variable representing unstructured y coordinate (default =[])</span>
594	0509 ivar_U=VarType.vector_x; <span class="comment">% defines (unique) index for the variable representing x vector component (default =[])</span>
595	0510 ivar_V=VarType.vector_y; <span class="comment">% defines (unique) index for the variable representing y vector component (default =[])</span>
596	0511 ivar_C=[VarType.scalar VarType.image VarType.color VarType.ancillary]; <span class="comment">%defines index (indices) for the scalar or ancillary fields</span>
597	0512 <span class="keyword">if</span> numel(ivar_C)>1
598	0513 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'too many scalar inputs in plot_field.m'</span>)
599	0514 <span class="keyword">return</span>
600	0515 <span class="keyword">end</span>
601	0516 ivar_F=VarType.warnflag; <span class="comment">%defines index (unique) for warning flag variable</span>
602	0517 ivar_FF=VarType.errorflag; <span class="comment">%defines index (unique) for error flag variable</span>
603	0518 idim_Y=[];
604	0519 test_grid=0;
605	0520 <span class="keyword">if</span> ~isempty(ivar_U) && ~isempty(ivar_V)<span class="comment">% vector components detected</span>
606	0521 <span class="keyword">if</span> test_vec
607	0522 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'attempt to plot two vector fields'</span>)
608	0523 <span class="keyword">return</span>
609	0524 <span class="keyword">else</span>
610	0525 test_vec=1;
611	0526 eval([<span class="string">'vec_U=Data.'</span> Data.ListVarName{ivar_U} <span class="string">';'</span>])
612	0527 eval([<span class="string">'vec_V=Data.'</span> Data.ListVarName{ivar_V} <span class="string">';'</span>])
613	0528 <span class="keyword">if</span> ~isempty(ivar_X) && ~isempty(ivar_Y)<span class="comment">% 2D field (with unstructured coordinates or structured ones (then ivar_X and ivar_Y empty)</span>
614	0529 eval([<span class="string">'vec_X=Data.'</span> Data.ListVarName{ivar_X} <span class="string">';'</span>])
615	0530 eval([<span class="string">'vec_Y=Data.'</span> Data.ListVarName{ivar_Y} <span class="string">';'</span>])
616	0531 <span class="keyword">elseif</span> numel(VarType.coord)==2 & VarType.coord~=[0 0];<span class="comment">%coordinates defines by dimension variables</span>
617	0532 eval([<span class="string">'y=Data.'</span> Data.ListVarName{VarType.coord(1)} <span class="string">';'</span>])
618	0533 eval([<span class="string">'x=Data.'</span> Data.ListVarName{VarType.coord(2)} <span class="string">';'</span>])
619	0534 <span class="keyword">if</span> numel(y)==2 <span class="comment">% y defined by first and last values on aregular mesh</span>
620	0535 y=linspace(y(1),y(2),size(vec_U,1));
621	0536 <span class="keyword">end</span>
622	0537 <span class="keyword">if</span> numel(x)==2 <span class="comment">% y defined by first and last values on aregular mesh</span>
623	0538 x=linspace(x(1),x(2),size(vec_U,2));
624	0539 <span class="keyword">end</span>
625	0540 [vec_X,vec_Y]=meshgrid(x,y);
626	0541 <span class="keyword">else</span>
627	0542 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'error in plot_field: bad coordinate definition for vector field'</span>)
628	0543 <span class="keyword">return</span>
629	0544 <span class="keyword">end</span>
630	0545 <span class="keyword">if</span> ~isempty(ivar_C)
631	0546 eval([<span class="string">'vec_C=Data.'</span> Data.ListVarName{ivar_C} <span class="string">';'</span>]) ;
632	0547 vec_C=reshape(vec_C,1,numel(vec_C));
633	0548 test_C=1;
634	0549 <span class="keyword">end</span>
635	0550 <span class="keyword">if</span> ~(isfield(PlotParam.Vectors,<span class="string">'HideWarning'</span>)&& isequal(PlotParam.Vectors.HideWarning,1))&&~isempty(ivar_F) &&~test_black
636	0551 <span class="keyword">if</span> numel(ivar_F)>1
637	0552 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple warning flags in plot_field.m, make a selection'</span>,<span class="string">'ERROR'</span>)
638	0553 <span class="keyword">return</span>
639	0554 <span class="keyword">elseif</span> test_vec
640	0555 eval([<span class="string">'AxeData.F=Data.'</span> Data.ListVarName{ivar_F} <span class="string">';'</span>]) <span class="comment">% flags for false vectors</span>
641	0556 test_black=1;
642	0557 <span class="keyword">end</span>
643	0558 <span class="keyword">end</span>
644	0559 <span class="keyword">if</span> ~isempty(ivar_FF) && ~test_false
645	0560 <span class="keyword">if</span> numel(ivar_FF)>1
646	0561 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple error flags in plot_field.m, make a selection'</span>,<span class="string">'ERROR'</span>)
647	0562 <span class="keyword">return</span>
648	0563 <span class="keyword">else</span> <span class="comment">%A VOIR CAS IMAGES</span>
649	0564 eval([<span class="string">'AxeData.FF=Data.'</span> Data.ListVarName{ivar_FF} <span class="string">';'</span>]) <span class="comment">% flags for false vectors</span>
650	0565 test_false=1;
651	0566 <span class="keyword">end</span>
652	0567 <span class="keyword">end</span>
653	0568 <span class="keyword">end</span>
654	0569 <span class="keyword">elseif</span> ~isempty(ivar_C) <span class="comment">%scalar or image</span>
655	0570 <span class="keyword">if</span> test_ima
656	0571 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'attempt to plot two scalar fields or images'</span>)
657	0572 <span class="keyword">return</span>
658	0573 <span class="keyword">end</span>
659	0574 eval([<span class="string">'A=Data.'</span> Data.ListVarName{ivar_C} <span class="string">';'</span>]) ;<span class="comment">% scalar represented as color image</span>
660	0575 test_ima=1;
661	0576 <span class="keyword">if</span> ~isempty(ivar_X) && ~isempty(ivar_Y)<span class="comment">% 2D field (with unstructured coordinates or structured ones (then ivar_X and ivar_Y empty)</span>
662	0577 eval([<span class="string">'AX=Data.'</span> Data.ListVarName{ivar_X} <span class="string">';'</span>])
663	0578 eval([<span class="string">'AY=Data.'</span> Data.ListVarName{ivar_Y} <span class="string">';'</span>])
664	0579 [A,AX,AY]=<a href="proj_grid.html" class="code" title="function [A,rangx,rangy]=proj_grid(vec_X,vec_Y,vec_A,rgx_in,rgy_in,npxy_in)">proj_grid</a>(AX',AY',A',[],[],<span class="string">'np>256'</span>); <span class="comment">% interpolate on a grid</span>
665	0580 <span class="keyword">elseif</span> numel(VarType.coord)>=2 & VarType.coord(1:2)~=[0 0];<span class="comment">%structured coordinates</span>
666	0581 eval([<span class="string">'AY=Data.'</span> Data.ListVarName{VarType.coord(1)} <span class="string">';'</span>])
667	0582 eval([<span class="string">'AX=Data.'</span> Data.ListVarName{VarType.coord(2)} <span class="string">';'</span>])
668	0583 test_interp_X=0; <span class="comment">%default, regularly meshed X coordinate</span>
669	0584 test_interp_Y=0; <span class="comment">%default, regularly meshed Y coordinate</span>
670	0585 <span class="keyword">if</span> numel(AY)>2
671	0586 DAY=diff(AY);
672	0587 DAY_min=min(DAY);
673	0588 DAY_max=max(DAY);
674	0589 <span class="keyword">if</span> sign(DAY_min)~=sign(DAY_max);<span class="comment">% =1 for increasing values, 0 otherwise</span>
675	0590 errormsg=[<span class="string">'errror in plot_field.m: non monotonic dimension variable # '</span> ListVarName{VarType.coord(1)} ];
676	0591 <span class="keyword">return</span>
677	0592 <span class="keyword">end</span>
678	0593 test_interp_Y=(DAY_max-DAY_min)> 0.0001*abs(DAY_max);
679	0594 <span class="keyword">end</span>
680	0595 <span class="keyword">if</span> numel(AX)>2
681	0596 DAX=diff(AX);
682	0597 DAX_min=min(DAX);
683	0598 DAX_max=max(DAX);
684	0599 <span class="keyword">if</span> sign(DAX_min)~=sign(DAX_max);<span class="comment">% =1 for increasing values, 0 otherwise</span>
685	0600 errormsg=[<span class="string">'errror in plot_field.m: non monotonic dimension variable # '</span> ListVarName{VarType.coord(2)} ];
686	0601 <span class="keyword">return</span>
687	0602 <span class="keyword">end</span>
688	0603 test_interp_X=(DAX_max-DAX_min)> 0.0001*abs(DAX_max);
689	0604 <span class="keyword">end</span>
690	0605 <span class="keyword">if</span> test_interp_Y
691	0606 npxy(1)=max([256 floor((AY(end)-AY(1))/DAY_min) floor((AY(end)-AY(1))/DAY_max)]);
692	0607 yI=linspace(AY(1),AY(end),npxy(1));
693	0608 <span class="keyword">if</span> ~test_interp_X
694	0609 xI=linspace(AX(1),AX(2),size(A,2));<span class="comment">%default</span>
695	0610 AX=xI;
696	0611 <span class="keyword">end</span>
697	0612 <span class="keyword">end</span>
698	0613 <span class="keyword">if</span> test_interp_X
699	0614 npxy(1)=max([256 floor((AX(end)-AX(1))/DAX_min) floor((AX(end)-AX(1))/DAX_max)]);
700	0615 xI=linspace(AX(1),AX(end),npxy(2));
701	0616 <span class="keyword">if</span> ~test_interp_Y
702	0617 yI=linspace(AY(1),AY(end),size(A,1));
703	0618 AY=yI;
704	0619 <span class="keyword">end</span>
705	0620 <span class="keyword">end</span>
706	0621 <span class="keyword">if</span> test_interp_X \|\| test_interp_Y
707	0622 [AX2D,AY2D]=meshgrid(AX,AY);
708	0623 A=interp2(AX2D,AY2D,double(A),xI,yI');
709	0624 <span class="keyword">end</span>
710	0625 AX=[AX(1) AX(end)];<span class="comment">% keep only the lower and upper bounds for image represnetation</span>
711	0626 AY=[AY(1) AY(end)];
712	0627 <span class="keyword">else</span>
713	0628 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'bad coordinate definition for vector field'</span>)
714	0629 <span class="keyword">return</span>
715	0630 <span class="keyword">end</span>
716	0631 <span class="keyword">end</span>
717	0632 <span class="keyword">if</span> isfield(Data,<span class="string">'VarAttribute'</span>)
718	0633 VarAttribute=Data.VarAttribute;
719	0634 <span class="keyword">end</span>
720	0635 <span class="keyword">end</span>
721	0636
722	0637 <span class="comment">%%%%%%%%%%%%%%%%%%%%% image or scalar plot %%%%%%%%%%%%%%%%%%%%%%%%%%</span>
723	0638
724	0639 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'Contours'</span>)
725	0640 PlotParam.Scalar.Contours=0; <span class="comment">%default</span>
726	0641 <span class="keyword">end</span>
727	0642 PlotParamOut=PlotParam; <span class="comment">%default</span>
728	0643 <span class="keyword">if</span> test_ima
729	0644 <span class="comment">% distinguish B/W and color images</span>
730	0645 np=size(A);<span class="comment">%size of image</span>
731	0646 siz=size(np);
732	0647 <span class="comment">%set the color map</span>
733	0648 <span class="keyword">if</span> siz(2)==2 <span class="comment">%for black and white images</span>
734	0649 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'AutoScal'</span>)
735	0650 PlotParam.Scalar.AutoScal=0;<span class="comment">%default</span>
736	0651 <span class="keyword">end</span>
737	0652 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'MinA'</span>)
738	0653 PlotParam.Scalar.MinA=[];<span class="comment">%default</span>
739	0654 <span class="keyword">end</span>
740	0655 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'MaxA'</span>)
741	0656 PlotParam.Scalar.MaxA=[];<span class="comment">%default</span>
742	0657 <span class="keyword">end</span>
743	0658 <span class="keyword">if</span> isequal(PlotParam.Scalar.AutoScal,0)\|isempty(PlotParam.Scalar.MinA)\|~isa(PlotParam.Scalar.MinA,<span class="string">'double'</span>) <span class="comment">%correct if there is no numerical data in edit box</span>
744	0659 MinA=double(min(min(A)));
745	0660 <span class="keyword">else</span>
746	0661 MinA=PlotParam.Scalar.MinA;
747	0662 <span class="keyword">end</span>;
748	0663 <span class="keyword">if</span> isequal(PlotParam.Scalar.AutoScal,0)\|isempty(PlotParam.Scalar.MaxA)\|~isa(PlotParam.Scalar.MaxA,<span class="string">'double'</span>) <span class="comment">%correct if there is no numerical data in edit box</span>
749	0664 MaxA=double(max(max(A)));
750	0665 <span class="keyword">else</span>
751	0666 MaxA=PlotParam.Scalar.MaxA;
752	0667 <span class="keyword">end</span>;
753	0668 PlotParamOut.Scalar.MinA=MinA;
754	0669 PlotParamOut.Scalar.MaxA=MaxA;
755	0670 axes(haxes)
756	0671 <span class="keyword">if</span> isequal(PlotParam.Scalar.Contours,1)
757	0672 <span class="keyword">if</span> ~isempty(hima) & ishandle(hima)
758	0673 delete(hima)
759	0674 <span class="keyword">end</span>
760	0675 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'IncrA'</span>)
761	0676 PlotParam.Scalar.IncrA=[];
762	0677 <span class="keyword">end</span>
763	0678 <span class="keyword">if</span> isempty(PlotParam.Scalar.IncrA)<span class="comment">% \| PlotParam.Scalar.AutoScal==0</span>
764	0679 cont=<a href="#_sub5" class="code" title="subfunction YTick=colbartick(MinA,MaxA)">colbartick</a>(MinA,MaxA);
765	0680 intercont=cont(2)-cont(1);<span class="comment">%default</span>
766	0681 PlotParamOut.Scalar.IncrA=intercont;
767	0682 <span class="keyword">else</span>
768	0683 intercont=PlotParam.Scalar.IncrA;
769	0684 <span class="keyword">end</span>
770	0685 B=A;
771	0686 abscontmin=intercont*floor(MinA/intercont);
772	0687 abscontmax=intercont*ceil(MaxA/intercont);
773	0688 contmin=intercont*floor(min(min(B))/intercont);
774	0689 contmax=intercont*ceil(max(max(B))/intercont);
775	0690 cont_pos_plus=[0:intercont:contmax];
776	0691 cont_pos_min=[double(contmin):intercont:-intercont];
777	0692 cont_pos=[cont_pos_min cont_pos_plus];
778	0693 sizpx=(AX(end)-AX(1))/(np(2)-1);
779	0694 sizpy=(AY(1)-AY(end))/(np(1)-1);
780	0695 x_cont=[AX(1):sizpx:AX(end)]; <span class="comment">% pixel x coordinates for image display</span>
781	0696 y_cont=[AY(1):-sizpy:AY(end)]; <span class="comment">% pixel x coordinates for image display</span>
782	0697 txt=ver;<span class="comment">%version of Matlab</span>
783	0698 Release=txt(1).Release;
784	0699 relnumb=str2num(Release(3:4));
785	0700 <span class="keyword">if</span> relnumb >= 14
786	0701 vec=linspace(0,1,(abscontmax-abscontmin)/intercont);<span class="comment">%define a greyscale colormap with steps intercont</span>
787	0702 map=[vec' vec' vec'];
788	0703 colormap(map);
789	0704 [var,hcontour]=contour(x_cont,y_cont,B,cont_pos);
790	0705 set(hcontour,<span class="string">'Fill'</span>,<span class="string">'on'</span>)
791	0706 set(hcontour,<span class="string">'LineStyle'</span>,<span class="string">'none'</span>)
792	0707 hold on
793	0708 <span class="keyword">end</span>
794	0709 [var_p,hcontour_p]=contour(x_cont,y_cont,B,cont_pos_plus,<span class="string">'k-'</span>);
795	0710 hold on
796	0711 [var_m,hcontour_m]=contour(x_cont,y_cont,B,cont_pos_min,<span class="string">':'</span>);
797	0712 set(hcontour_m,<span class="string">'LineColor'</span>,[1 1 1])
798	0713 hold off
799	0714 caxis([abscontmin abscontmax])
800	0715 colormap(map);
801	0716 <span class="keyword">end</span>
802	0717 <span class="keyword">if</span> ~isequal(PlotParam.Scalar.Contours,1)
803	0718 <span class="comment">% rescale the grey levels with min and max, put a grey scale colorbar</span>
804	0719 <span class="keyword">if</span> (isfield(PlotParam.Scalar,<span class="string">'BW'</span>)& (isequal(PlotParam.Scalar.BW,1))\|isa(A,<span class="string">'uint8'</span>)\| isa(A,<span class="string">'uint16'</span>))<span class="comment">%images</span>
805	0720 B=A;
806	0721 vec=linspace(0,1,255);<span class="comment">%define a linear greyscale colormap</span>
807	0722 map=[vec' vec' vec'];
808	0723 colormap(map); <span class="comment">%grey scale color map</span>
809	0724 <span class="keyword">else</span>
810	0725 B=A;
811	0726 colormap(<span class="string">'default'</span>); <span class="comment">% standard faulse colors for div, vort , scalar fields</span>
812	0727 <span class="keyword">end</span>
813	0728 <span class="keyword">end</span>
814	0729 <span class="keyword">elseif</span> siz(2)==3 <span class="comment">%color images</span>
815	0730 axes(haxes)
816	0731 B=uint8(A);
817	0732 MinA=0;
818	0733 MaxA=255;
819	0734 <span class="keyword">end</span>
820	0735 <span class="keyword">if</span> ~isequal(PlotParam.Scalar.Contours,1)
821	0736 <span class="comment">%interpolate to increase resolution</span>
822	0737 test_interp=1;
823	0738 <span class="keyword">if</span> max(np) <= 64
824	0739 npxy=8*np;<span class="comment">% increase the resolution 8 times</span>
825	0740 <span class="keyword">elseif</span> max(np) <= 128
826	0741 npxy=4*np;<span class="comment">% increase the resolution 4 times</span>
827	0742 <span class="keyword">elseif</span> max(np) <= 256
828	0743 npxy=2*np;<span class="comment">% increase the resolution 2 times</span>
829	0744 <span class="keyword">else</span>
830	0745 npxy=np;
831	0746 test_interp=0; <span class="comment">% no interpolation done</span>
832	0747 <span class="keyword">end</span>
833	0748 <span class="keyword">if</span> test_interp==1<span class="comment">%if we interpolate</span>
834	0749 x=[AX(1):(AX(2)-AX(1))/(np(2)-1):AX(2)];
835	0750 y=[AY(1):(AY(2)-AY(1))/(np(1)-1):AY(2)];
836	0751 [X,Y]=meshgrid(x,y);
837	0752 <span class="comment">% xi=[AX(1):(AX(2)-AX(1))/(npxy(2)-1):AX(2)];</span>
838	0753 <span class="comment">% yi=[AY(1):(AY(2)-AY(1))/(npxy(1)-1):AY(2)];</span>
839	0754
840	0755 xi=linspace(AX(1),AX(2),npxy(2));
841	0756 yi=linspace(AY(1),AY(2),npxy(1));
842	0757 B = interp2(X,Y,double(B),xi,yi');
843	0758 <span class="keyword">end</span>
844	0759 <span class="keyword">if</span> isempty(hima)
845	0760 tag=get(haxes,<span class="string">'Tag'</span>);
846	0761 hima=imagesc(AX,AY,B,[MinA MaxA]);
847	0762 set(hima,<span class="string">'Tag'</span>,<span class="string">'ima'</span>,<span class="string">'HitTest'</span>,<span class="string">'off'</span>)
848	0763 set(haxes,<span class="string">'Tag'</span>,tag);<span class="comment">%preserve the axes tag (removed by image fct !!!)</span>
849	0764 <span class="comment">%axis(haxes,'image'); % preserve the image aspect ratio %TESTMOD</span>
850	0765 <span class="keyword">else</span>
851	0766 set(hima,<span class="string">'CData'</span>,B);
852	0767 <span class="keyword">if</span> MinA<MaxA
853	0768 caxis([MinA MaxA])
854	0769 <span class="keyword">else</span>
855	0770 caxis([MinA MinA+1])
856	0771 <span class="keyword">end</span>
857	0772 set(hima,<span class="string">'XData'</span>,AX);
858	0773 set(hima,<span class="string">'YData'</span>,AY);
859	0774 <span class="keyword">end</span>
860	0775 <span class="keyword">end</span>
861	0776 <span class="keyword">if</span> ~isstruct(AxeData)
862	0777 AxeData=[];
863	0778 <span class="keyword">end</span>
864	0779 AxeData.A=A;
865	0780 AxeData.AX=[AX(1) AX(end)];
866	0781 AxeData.AY=[AY(1) AY(end)];
867	0782 test_ima=1;
868	0783 <span class="comment">%display the colorbar code for B/W images if Poscolorbar not empty</span>
869	0784 <span class="keyword">if</span> siz(2)==2 & exist(<span class="string">'PosColorbar'</span>,<span class="string">'var'</span>)& ~isempty(PosColorbar)
870	0785 <span class="keyword">if</span> isempty(hcol)\|~ishandle(hcol)
871	0786 hcol=colorbar;<span class="comment">%create new colorbar</span>
872	0787 <span class="keyword">end</span>
873	0788 <span class="keyword">if</span> length(PosColorbar)==4
874	0789 set(hcol,<span class="string">'Position'</span>,PosColorbar)
875	0790 <span class="keyword">end</span>
876	0791 YTick=0;<span class="comment">%default</span>
877	0792 <span class="keyword">if</span> MaxA>MinA
878	0793 <span class="keyword">if</span> isequal(PlotParam.Scalar.Contours,1)
879	0794 colbarlim=get(hcol,<span class="string">'YLim'</span>);
880	0795 scale_bar=(colbarlim(2)-colbarlim(1))/(abscontmax-abscontmin);
881	0796 YTick=cont_pos(2:end-1);
882	0797 YTick_scaled=colbarlim(1)+scale_bar*(YTick-abscontmin);
883	0798 set(hcol,<span class="string">'YTick'</span>,YTick_scaled);
884	0799 <span class="keyword">elseif</span> (isfield(PlotParam.Scalar,<span class="string">'BW'</span>) & isequal(PlotParam.Scalar.BW,1))\|isa(A,<span class="string">'uint8'</span>)\| isa(A,<span class="string">'uint16'</span>)<span class="comment">%images</span>
885	0800 hi=get(hcol,<span class="string">'children'</span>);
886	0801 <span class="keyword">if</span> iscell(hi)<span class="comment">%multiple images in colorbar</span>
887	0802 hi=hi{1};
888	0803 <span class="keyword">end</span>
889	0804 set(hi,<span class="string">'YData'</span>,[MinA MaxA])
890	0805 set(hi,<span class="string">'CData'</span>,[1:256]')
891	0806 set(hcol,<span class="string">'YLim'</span>,[MinA MaxA])
892	0807 YTick=<a href="#_sub5" class="code" title="subfunction YTick=colbartick(MinA,MaxA)">colbartick</a>(MinA,MaxA);
893	0808 set(hcol,<span class="string">'YTick'</span>,YTick)
894	0809 <span class="keyword">else</span>
895	0810 hi=get(hcol,<span class="string">'children'</span>);
896	0811 <span class="keyword">if</span> iscell(hi)<span class="comment">%multiple images in colorbar</span>
897	0812 hi=hi{1};
898	0813 <span class="keyword">end</span>
899	0814 set(hi,<span class="string">'YData'</span>,[MinA MaxA])
900	0815 set(hi,<span class="string">'CData'</span>,[1:64]')
901	0816 YTick=<a href="#_sub5" class="code" title="subfunction YTick=colbartick(MinA,MaxA)">colbartick</a>(MinA,MaxA);
902	0817 set(hcol,<span class="string">'YLim'</span>,[MinA MaxA])
903	0818 set(hcol,<span class="string">'YTick'</span>,YTick)
904	0819 <span class="keyword">end</span>
905	0820 set(hcol,<span class="string">'Yticklabel'</span>,num2str(YTick'));
906	0821 <span class="keyword">end</span>
907	0822 <span class="keyword">elseif</span> ishandle(hcol)
908	0823 delete(hcol); <span class="comment">%erase existing colorbar if not needed</span>
909	0824 <span class="keyword">end</span>
910	0825 <span class="keyword">else</span><span class="comment">%no scalar plot</span>
911	0826 <span class="keyword">if</span> ~isempty(hima) && ishandle(hima)
912	0827 delete(hima)
913	0828 <span class="keyword">end</span>
914	0829 <span class="keyword">if</span> ~isempty(hcol)&& ishandle(hcol)
915	0830 delete(hcol)
916	0831 <span class="keyword">end</span>
917	0832 AxeData.A=[];
918	0833 AxeData.AX=[];
919	0834 AxeData.AY=[];
920	0835 PlotParamOut=rmfield(PlotParamOut,<span class="string">'Scalar'</span>);
921	0836 <span class="keyword">end</span>
922	0837
923	0838 <span class="comment">%%%%%%%%%%%%%%%%%%%%% vector plot %%%%%%%%%%%%%%%%%%%%%%%%%%</span>
924	0839 <span class="keyword">if</span> test_vec
925	0840 <span class="comment">%vector scale representation</span>
926	0841 <span class="keyword">if</span> size(vec_U,1)==numel(vec_Y) && size(vec_U,2)==numel(vec_X); <span class="comment">% x, y coordinate variables</span>
927	0842 [vec_X,vec_Y]=meshgrid(vec_X,vec_Y);
928	0843 <span class="keyword">end</span>
929	0844 vec_X=reshape(vec_X,1,numel(vec_X));<span class="comment">%reshape in matlab vectors</span>
930	0845 vec_Y=reshape(vec_Y,1,numel(vec_Y));
931	0846 vec_U=reshape(vec_U,1,numel(vec_U));
932	0847 vec_V=reshape(vec_V,1,numel(vec_V));
933	0848 MinMaxX=max(vec_X)-min(vec_X);
934	0849 MinMaxY=max(vec_Y)-min(vec_Y);
935	0850 AxeData.Mesh=sqrt((MinMaxX*MinMaxY)/length(vec_X));
936	0851 <span class="keyword">if</span> ~isfield(PlotParam.Vectors,<span class="string">'AutoVec'</span>) \|\| isequal(PlotParam.Vectors.AutoVec,0)\|\| ~isfield(PlotParam.Vectors,<span class="string">'VecScale'</span>)<span class="keyword">...</span>
937	0852 \|\|isempty(PlotParam.Vectors.VecScale)\|\|~isa(PlotParam.Vectors.VecScale,<span class="string">'double'</span>) <span class="comment">%automatic vector scale</span>
938	0853 scale=[];
939	0854 <span class="keyword">if</span> test_false <span class="comment">%remove false vectors</span>
940	0855 indsel=find(AxeData.FF==0);<span class="comment">%indsel =indices of good vectors</span>
941	0856 <span class="keyword">else</span>
942	0857 indsel=[1:numel(vec_X)];<span class="comment">%</span>
943	0858 <span class="keyword">end</span>
944	0859 <span class="keyword">if</span> isempty(vec_U)
945	0860 scale=1;
946	0861 <span class="keyword">else</span>
947	0862 <span class="keyword">if</span> isempty(indsel)
948	0863 MaxU=max(abs(vec_U));
949	0864 MaxV=max(abs(vec_V));
950	0865 <span class="keyword">else</span>
951	0866 MaxU=max(abs(vec_U(indsel)));
952	0867 MaxV=max(abs(vec_V(indsel)));
953	0868 <span class="keyword">end</span>
954	0869 scale=MinMaxX/(max(MaxU,MaxV)*50);
955	0870 PlotParam.Vectors.VecScale=scale;<span class="comment">%update the 'scale' display</span>
956	0871 <span class="keyword">end</span>
957	0872 <span class="keyword">else</span>
958	0873 scale=PlotParam.Vectors.VecScale; <span class="comment">%impose the length of vector representation</span>
959	0874 <span class="keyword">end</span>;
960	0875
961	0876 <span class="comment">%record vectors on the plotting axes</span>
962	0877 <span class="keyword">if</span> test_C==0
963	0878 vec_C=ones(1,numel(vec_X));
964	0879 <span class="keyword">end</span>
965	0880 AxeData.X=vec_X';
966	0881 AxeData.Y=vec_Y';
967	0882 AxeData.U=vec_U';
968	0883 AxeData.V=vec_V';
969	0884 AxeData.C=vec_C';
970	0885 <span class="keyword">if</span> isfield(Data,<span class="string">'W'</span>)
971	0886 AxeData.W=Data.W;
972	0887 <span class="keyword">end</span>
973	0888
974	0889 <span class="comment">%decimate by a factor 2 in vector mesh(4 in nbre of vectors)</span>
975	0890 <span class="keyword">if</span> isfield(PlotParam.Vectors,<span class="string">'decimate4'</span>)&isequal(PlotParam.Vectors.decimate4,1)
976	0891 diffy=diff(vec_Y); <span class="comment">%difference dy=vec_Y(i+1)-vec_Y(i)</span>
977	0892 dy_thresh=max(abs(diffy))/2;
978	0893 ind_jump=find(abs(diffy) > dy_thresh); <span class="comment">%indices with diff(vec_Y)> max/2, detect change of line</span>
979	0894 ind_sel=[1:ind_jump(1)];<span class="comment">%select the first line</span>
980	0895 <span class="keyword">for</span> i=2:2:length(ind_jump)-1
981	0896 ind_sel=[ind_sel [ind_jump(i)+1:ind_jump(i+1)]];<span class="comment">% select the odd lines</span>
982	0897 <span class="keyword">end</span>
983	0898 nb_sel=length(ind_sel);
984	0899 ind_sel=ind_sel([1:2:nb_sel]);<span class="comment">% take half the points on a line</span>
985	0900 vec_X=vec_X(ind_sel);
986	0901 vec_Y=vec_Y(ind_sel);
987	0902 vec_U=vec_U(ind_sel);
988	0903 vec_V=vec_V(ind_sel);
989	0904 vec_C=vec_C(ind_sel);
990	0905 <span class="keyword">if</span> isfield(Data,<span class="string">'F'</span>)&& isequal(size(Data.F),size(vec_X))sdfsdf
991	0906 Data.F=Data.F(ind_sel);
992	0907 <span class="keyword">end</span>
993	0908 <span class="keyword">if</span> isfield(Data,<span class="string">'FF'</span>)&& isequal(size(Data.FF),size(vec_X))
994	0909 Data.FF=Data.FF(ind_sel);
995	0910 <span class="keyword">end</span>
996	0911 <span class="keyword">end</span>
997	0912
998	0913 <span class="comment">%get main level color code</span>
999	0914 [colorlist,col_vec,PlotParamOut.Vectors]=<a href="set_col_vec.html" class="code" title="function [colorlist,col_vec,colcode_out]=set_col_vec(colcode,vec_C)">set_col_vec</a>(PlotParam.Vectors,vec_C);
1000	0915 <span class="comment">% take flags into account: add flag colors to the list of colors</span>
1001	0916 sizlist=size(colorlist);
1002	0917 nbcolor=sizlist(1);
1003	0918 <span class="keyword">if</span> test_black && isfield (Data,<span class="string">'F'</span>)
1004	0919 nbcolor=nbcolor+1;
1005	0920 colorlist(nbcolor,:)=[0 0 0]; <span class="comment">%add black to the list of colors</span>
1006	0921 <span class="keyword">if</span> isfield(Data,<span class="string">'FF'</span>)
1007	0922 ind_flag=find(Data.F~=1 & Data.FF==0); <span class="comment">%flag warning but not faulse</span>
1008	0923 <span class="keyword">else</span>
1009	0924 ind_flag=find(Data.F~=1);
1010	0925 <span class="keyword">end</span>
1011	0926 col_vec(ind_flag)=nbcolor;
1012	0927 <span class="keyword">end</span>
1013	0928 nbcolor=nbcolor+1;
1014	0929 <span class="keyword">if</span> test_false && isfield (Data,<span class="string">'FF'</span>)
1015	0930 ind_flag=find(Data.FF~=0);
1016	0931 <span class="keyword">if</span> isfield(PlotParam.Vectors,<span class="string">'HideFalse'</span>) & PlotParam.Vectors.HideFalse==1
1017	0932 colorlist(nbcolor,:)=[NaN NaN NaN];<span class="comment">% no plot of false vectors</span>
1018	0933 <span class="keyword">else</span>
1019	0934 colorlist(nbcolor,:)=[1 0 1];<span class="comment">% magenta color</span>
1020	0935 <span class="keyword">end</span>
1021	0936 col_vec(ind_flag)=nbcolor;
1022	0937 <span class="keyword">end</span>
1023	0938 <span class="comment">%plot vectors:</span>
1024	0939 <a href="#_sub4" class="code" title="subfunction quiresetn(haxes,x,y,u,v,scale,colorlist,col_vec)">quiresetn</a>(haxes,vec_X,vec_Y,vec_U,vec_V,scale,colorlist,col_vec);
1025	0940 <span class="keyword">else</span>
1026	0941 hvec=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'vel'</span>);
1027	0942 <span class="keyword">if</span> ~isempty(hvec)
1028	0943 delete(hvec);
1029	0944 <span class="keyword">end</span>
1030	0945 AxeData.X=[];
1031	0946 AxeData.Y=[];
1032	0947 AxeData.U=[];
1033	0948 AxeData.V=[];
1034	0949 AxeData.C=[];
1035	0950 AxeData.W=[];
1036	0951 AxeData.Mesh=[];
1037	0952 PlotParamOut=rmfield(PlotParamOut,<span class="string">'Vectors'</span>);
1038	0953 <span class="keyword">end</span>
1039	0954 <span class="keyword">if</span> isfield(Data,<span class="string">'Z'</span>)
1040	0955 AxeData.Z=Data.Z;<span class="comment">% A REVOIR</span>
1041	0956 <span class="keyword">end</span>
1042	0957 listfields={<span class="string">'AY'</span>,<span class="string">'AX'</span>,<span class="string">'A'</span>,<span class="string">'X'</span>,<span class="string">'Y'</span>,<span class="string">'U'</span>,<span class="string">'V'</span>,<span class="string">'C'</span>,<span class="string">'W'</span>,<span class="string">'F'</span>,<span class="string">'FF'</span>};
1043	0958 listdim={<span class="string">'AY'</span>,<span class="string">'AX'</span>,{<span class="string">'AY'</span>,<span class="string">'AX'</span>},<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>};
1044	0959 Role={<span class="string">'coord_y'</span>,<span class="string">'coord_x'</span>,<span class="string">'scalar'</span>,<span class="string">'coord_x'</span>,<span class="string">'coord_y'</span>,<span class="string">'vector_x'</span>,<span class="string">'vector_y'</span>,<span class="string">'scalar'</span>,<span class="string">'vector_z'</span>,<span class="string">'warnflag'</span>,<span class="string">'falseflag'</span>};
1045	0960 ind_select=[];
1046	0961 nbvar=0;
1047	0962
1048	0963 <span class="keyword">for</span> ilist=1:numel(listfields)
1049	0964 eval([<span class="string">'testvar=isfield(AxeData,listfields{ilist}) && ~isempty(AxeData.'</span> listfields{ilist} <span class="string">');'</span>])
1050	0965 <span class="keyword">if</span> testvar
1051	0966 nbvar=nbvar+1;
1052	0967 AxeData.ListVarName{nbvar}=listfields{ilist};
1053	0968 AxeData.VarDimName{nbvar}=listdim{ilist};
1054	0969 AxeData.VarAttribute{nbvar}.Role=Role{ilist};
1055	0970 <span class="keyword">end</span>
1056	0971 <span class="keyword">end</span>
1057	0972
1058	0973 <span class="comment">% adjust the size of the plot to include the whole field, except if KeepLim=1</span>
1059	0974 <span class="keyword">if</span> ~(exist(<span class="string">'KeepLim'</span>,<span class="string">'var'</span>) && isequal(KeepLim,1)) <span class="comment">%adjust the graph limits*</span>
1060	0975 test_lim=0;
1061	0976 <span class="keyword">if</span> test_vec
1062	0977 Xlim=[min(vec_X) max(vec_X)];
1063	0978 Ylim=[min(vec_Y) max(vec_Y)];
1064	0979 test_lim=1;
1065	0980 <span class="keyword">if</span> test_ima<span class="comment">%both background image and vectors coexist, take the wider bound</span>
1066	0981 Xlim(1)=min(AX(1),Xlim(1));
1067	0982 Xlim(2)=max(AX(end),Xlim(2));
1068	0983 Ylim(1)=min(AY(end),Ylim(1));
1069	0984 Ylim(2)=max(AY(1),Ylim(2));
1070	0985 <span class="keyword">end</span>
1071	0986 <span class="keyword">elseif</span> test_ima <span class="comment">%only image plot</span>
1072	0987 Xlim(1)=min(AX(1),AX(end));
1073	0988 Xlim(2)=max(AX(1),AX(end));
1074	0989 Ylim(1)=min(AY(1),AY(end));
1075	0990 Ylim(2)=max(AY(1),AY(end));
1076	0991 test_lim=1;
1077	0992 <span class="keyword">end</span>
1078	0993 <span class="keyword">if</span> test_lim
1079	0994 set(haxes,<span class="string">'XLim'</span>,Xlim);<span class="comment">% set x limits of frame in axes coordinates</span>
1080	0995 set(haxes,<span class="string">'YLim'</span>,Ylim);<span class="comment">% set y limits of frame in axes coordinate</span>
1081	0996 <span class="keyword">end</span>
1082	0997 <span class="keyword">end</span>
1083	0998 <span class="keyword">if</span> ~(isfield(PlotParam,<span class="string">'Auto_xy'</span>) && isequal(PlotParam.Auto_xy,1))
1084	0999 set(haxes,<span class="string">'DataAspectRatio'</span>,[1 1 1])
1085	1000 <span class="keyword">end</span>
1086	1001 set(haxes,<span class="string">'YDir'</span>,<span class="string">'normal'</span>)
1087	1002
1088	1003 <span class="comment">%---------------------------------------------</span>
1089	1004 <span class="comment">%function for plotting vectors</span>
1090	1005 <span class="comment">%------------------------------------------------</span>
1091	1006 <span class="comment">%INPUT:</span>
1092	1007 <span class="comment">% haxes: handles of the plotting axes</span>
1093	1008 <span class="comment">%x,y,u,v: vectors coordinates and vector components to plot, arrays withb the same dimension</span>
1094	1009 <span class="comment">% scale: scaling factor for vector length representation</span>
1095	1010 <span class="comment">%colorlist(icolor,:): list of vector colors, dim (nbcolor,3), depending on color #i</span>
1096	1011 <span class="comment">%col_vec: matlab vector setting the color number #i for each velocity vector</span>
1097	1012 <a name="_sub4" href="#_subfunctions" class="code">function quiresetn(haxes,x,y,u,v,scale,colorlist,col_vec)</a>
1098	1013
1099	1014 <span class="comment">%define arrows</span>
1100	1015 theta=0.5 ;<span class="comment">%angle arrow</span>
1101	1016 alpha=0.3 ;<span class="comment">%length arrow</span>
1102	1017 rot=alpha*[cos(theta) -sin(theta); sin(theta) cos(theta)]';
1103	1018 <span class="comment">%find the existing lines</span>
1104	1019 <span class="comment">%h=findobj(gca,'Type','Line');% search existing lines in the current axes</span>
1105	1020 h=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'vel'</span>);<span class="comment">% search existing lines in the current axes</span>
1106	1021 sizh=size(h);
1107	1022 set(h,<span class="string">'EraseMode'</span>,<span class="string">'xor'</span>);
1108	1023 set(haxes,<span class="string">'NextPlot'</span>,<span class="string">'replacechildren'</span>);
1109	1024 <span class="comment">%htext=findobj('Tag','scalevec');</span>
1110	1025
1111	1026 <span class="comment">%drawnow</span>
1112	1027 <span class="comment">%create lines (if no lines) or modify them</span>
1113	1028 <span class="keyword">if</span> ~isequal(size(col_vec),size(x))
1114	1029 col_vec=ones(size(x));<span class="comment">% case of error in col_vec input</span>
1115	1030 <span class="keyword">end</span>
1116	1031 sizlist=size(colorlist);
1117	1032 ncolor=sizlist(1);
1118	1033
1119	1034 <span class="keyword">for</span> icolor=1:ncolor
1120	1035 <span class="comment">%determine the line positions for each color icolor</span>
1121	1036 ind=find(col_vec==icolor);
1122	1037 xc=x(ind);
1123	1038 yc=y(ind);
1124	1039 uc=u(ind)*scale;
1125	1040 vc=v(ind)*scale;
1126	1041 n=size(xc);
1127	1042 xN=NaN*ones(size(xc));
1128	1043 matx=[xc(:) xc(:)+uc(:) xN(:)]';
1129	1044 matx=reshape(matx,1,3*n(2));
1130	1045 maty=[yc(:) yc(:)+vc(:) xN(:)]';
1131	1046 maty=reshape(maty,1,3*n(2));
1132	1047
1133	1048 <span class="comment">%determine arrow heads</span>
1134	1049 arrowplus=rot*[uc;vc];
1135	1050 arrowmoins=rot'*[uc;vc];
1136	1051 x1=xc+uc-arrowplus(1,:);
1137	1052 x2=xc+uc;
1138	1053 x3=xc+uc-arrowmoins(1,:);
1139	1054 y1=yc+vc-arrowplus(2,:);
1140	1055 y2=yc+vc;
1141	1056 y3=yc+vc-arrowmoins(2,:);
1142	1057 matxar=[x1(:) x2(:) x3(:) xN(:)]';
1143	1058 matxar=reshape(matxar,1,4*n(2));
1144	1059 matyar=[y1(:) y2(:) y3(:) xN(:)]';
1145	1060 matyar=reshape(matyar,1,4*n(2));
1146	1061 <span class="comment">%draw the line or modify the existing ones</span>
1147	1062 isn=isnan(colorlist(icolor,:));<span class="comment">%test if color NaN</span>
1148	1063 <span class="keyword">if</span> 2*icolor > sizh(1) <span class="comment">%if icolor exceeds the number of existing ones</span>
1149	1064 axes(haxes)
1150	1065 <span class="keyword">if</span> ~isn(1) <span class="comment">%if the vectors are visible color not nan</span>
1151	1066 <span class="keyword">if</span> n(2)>0
1152	1067 hold on
1153	1068 line(matx,maty,<span class="string">'Color'</span>,colorlist(icolor,:),<span class="string">'Tag'</span>,<span class="string">'vel'</span>);<span class="comment">% plot new lines</span>
1154	1069 line(matxar,matyar,<span class="string">'Color'</span>,colorlist(icolor,:),<span class="string">'Tag'</span>,<span class="string">'vel'</span>);<span class="comment">% plot arrows</span>
1155	1070 <span class="keyword">end</span>
1156	1071 <span class="keyword">end</span>
1157	1072 <span class="keyword">else</span>
1158	1073 <span class="keyword">if</span> isn(1)
1159	1074 delete(h(2*icolor-1))
1160	1075 delete(h(2*icolor))
1161	1076 <span class="keyword">else</span>
1162	1077 set(h(2*icolor-1),<span class="string">'Xdata'</span>,matx,<span class="string">'Ydata'</span>,maty);
1163	1078 set(h(2*icolor-1),<span class="string">'Color'</span>,colorlist(icolor,:));
1164	1079 set(h(2*icolor-1),<span class="string">'EraseMode'</span>,<span class="string">'xor'</span>);
1165	1080 set(h(2*icolor),<span class="string">'Xdata'</span>,matxar,<span class="string">'Ydata'</span>,matyar);
1166	1081 set(h(2*icolor),<span class="string">'Color'</span>,colorlist(icolor,:));
1167	1082 set(h(2*icolor),<span class="string">'EraseMode'</span>,<span class="string">'xor'</span>);
1168	1083 <span class="keyword">end</span>
1169	1084 <span class="keyword">end</span>
1170	1085 <span class="keyword">end</span>
1171	1086 <span class="keyword">if</span> sizh(1) > 2*ncolor
1172	1087 <span class="keyword">for</span> icolor=ncolor+1 : sizh(1)/2<span class="comment">%delete additional objects</span>
1173	1088 delete(h(2*icolor-1))
1174	1089 delete(h(2*icolor))
1175	1090 <span class="keyword">end</span>
1176	1091 <span class="keyword">end</span>
1177	1092
1178	1093 <span class="comment">%---------------------------------------</span>
1179	1094 <span class="comment">%determine tick positions for colorbar</span>
1180	1095 <span class="comment">%------------------------------------</span>
1181	1096 <a name="_sub5" href="#_subfunctions" class="code">function YTick=colbartick(MinA,MaxA)</a>
1182	1097 <span class="comment">%determine tick positions with "simple" values between MinA and MaxA</span>
1183	1098 YTick=0;<span class="comment">%default</span>
1184	1099 maxabs=max([abs(MinA) abs(MaxA)]);
1185	1100 <span class="keyword">if</span> maxabs>0
1186	1101 ord=10^(floor(log10(maxabs)));<span class="comment">%order of magnitude</span>
1187	1102 div=1;
1188	1103 siz2=1;
1189	1104 <span class="keyword">while</span> siz2<2
1190	1105 <span class="comment">% values=[-9:div:9];</span>
1191	1106 values=-10:div:10;
1192	1107 ind=find((ordvalues-MaxA)<0 & (ordvalues-MinA)>0);<span class="comment">%indices of 'values' such that MinA<ord*values<MaxA</span>
1193	1108 siz=size(ind);
1194	1109 <span class="keyword">if</span> siz(2)<4<span class="comment">%if there are less than 4 selected values (4 levels)</span>
1195	1110 values=[-9:0.5*div:9];
1196	1111 ind=find((ordvalues-MaxA)<0 & (ordvalues-MinA)>0);
1197	1112 <span class="keyword">end</span>
1198	1113 siz2=size(ind,2);
1199	1114 <span class="comment">% siz2=siz(2)</span>
1200	1115 div=div/10;
1201	1116 <span class="keyword">end</span>
1202	1117 YTick=ord*values(ind);
1203	1118 <span class="keyword">end</span></pre></div>
1204	<hr><address>Generated on Fri 13-Nov-2009 11:17:03 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> © 2003</address>
1205	</body>
1206	</html>

Note: See TracBrowser for help on using the repository browser.

Download in other formats: