matlab_class/gcmfaces_diags/diags_driver_tex.m

function []=diags_driver_tex(dirMat,setDiags,dirTex,nameTex);
% DIAGS_DRIVER_TEX(dirMat,setDiags,dirTex,nameTex)
%
%    displays multiple sets of diagnostics (setDiags={'profiles',
%    'cost','A','B','C','MLD','D'} by default) from the results
%    stored in dirMat and outputs the plots to tex 
%    ([dirTex nameTex '.tex'] and [dirTex nameTex '*.eps'])
%
%    setDiags is the choice of diagnostics (cell) that may include
%                       'profiles') model to insitu data comparison
%                       'cost') ssh etc. cost functions
%                       'A') trasnports
%                       'B') air-sea fluxes
%                       'C') state variables
%                       'D') global and hemispheric budgets
%                       'MLD') mixed layer depths
%                       'SEAICE') seaice fields
%                       'controls') control vector adjustments

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%determine input/output params:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%directory names:
if iscell(dirMat); dirMatRef=dirMat{2}; dirMat=dirMat{1}; end;
dirMat=[dirMat '/'];
if isempty(who('dirMatRef')); dirMatRef='';
elseif ~isempty(dirMatRef); dirMatRef=[dirMatRef '/'];
end;
if isempty(who('dirTex')); dirTex=''; else; dirTex=[dirTex '/']; end;
if isempty(who('nameTex')); nameTex='myPlots'; end;

%set default setDiags and myTitle
if isempty(setDiags);
  setDiags={'profiles','cost','A','B','C','MLD','D'};
  if strcmp(nameTex,'myPlots'); nameTex='standardAnalysis'; end;
  %
  tmp1=dirMat; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
  myTitle={'gcmfaces','standard analysis of the solution in',tmp1};
else;
  tmp1=dirMat; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
  myTitle={'gcmfaces analysis of the solution in',tmp1};
end;

%set fileTex and create dirTex if needed
if isempty(dirTex); error('dirTex must be specified'); end;
fileTex=[dirTex nameTex '.tex'];
if isempty(dir(dirTex)); mkdir(dirTex); end;

%%%%%%%%%%%%%%%%%%%%%%
%load grid and params:
%%%%%%%%%%%%%%%%%%%%%%

gcmfaces_global; global myparms;
test1=~isempty(dir([dirMat 'basic_diags_ecco_mygrid.mat']));
test2=~isempty(dir([dirMat 'diags_grid_parms.mat']));
if ~test1&~test2;
  error('missing diags_grid_parms.mat')
elseif test2;
  nameGrid='diags_grid_parms.mat';
  suffDiag='diags_set_';
  budgetList='diags_select_budget_list.mat';
else;
  nameGrid='basic_diags_ecco_mygrid.mat';
  suffDiag='basic_diags_ecco_';
  budgetList='basic_diags_ecco_budget_list.mat';
end;

%here we always reload the grid from dirMat to make sure the same one is used throughout
eval(['load ' dirMat nameGrid ';']);

%in case mygrid.memoryLimit=1, load the stuff that was not saved to diags_grid_parms.mat
if mygrid.memoryLimit==1;
        list0={'hFacS','hFacW'};
        for iFld=1:length(list0);
          eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([mygrid.dirGrid ''' list0{iFld} '*'']);']);
        end;
        %
        mygrid.hFacCsurf=mygrid.hFacC;
        for ff=1:mygrid.hFacC.nFaces; mygrid.hFacCsurf{ff}=mygrid.hFacC{ff}(:,:,1); end;
        %
        mskC=mygrid.hFacC; mskC(mskC==0)=NaN; mskC(mskC>0)=1; mygrid.mskC=mskC;
        mskW=mygrid.hFacW; mskW(mskW==0)=NaN; mskW(mskW>0)=1; mygrid.mskW=mskW;
        mskS=mygrid.hFacS; mskS(mskS==0)=NaN; mskS(mskS>0)=1; mygrid.mskS=mskS;
        %
        gcmfaces_lines_zonal;
        mygrid.LATS=[mygrid.LATS_MASKS.lat]';
        [lonPairs,latPairs,names]=gcmfaces_lines_pairs;
        gcmfaces_lines_transp(lonPairs,latPairs,names);
end;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finalize listDiags
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%determined where to display anomalies between runs
doAnomalies=~isempty(dirMatRef);

%check if result is here to plot
doBudget=~isempty(dir([dirMat suffDiag 'D_*']));
doBudget=~isempty(dir([dirMat suffDiag 'D_*']))|...
    ~isempty(dir([dirMat 'diags_set_D/' suffDiag 'D_*']));
doProfiles=~isempty(dir([dirMat 'insitu_cost_all.mat']))|...
    ~isempty(dir([dirMat 'cost/insitu_cost_all.mat']));
doCost=~isempty(dir([dirMat 'cost_altimeter_obs.mat']))|...
    ~isempty(dir([dirMat 'cost/cost_altimeter_obs.mat']));
doCtrl=~isempty(dir([dirMat 'cost_xx_aqh.mat']))|...
    ~isempty(dir([dirMat 'cost/cost_xx_aqh.mat']));

%the following have no code for diff between runs
if doAnomalies;
    doProfiles=0;
    doCost=0;
    doCtrl=0;
end;

%reduce setDiags if needed:
doDiags=ones(1,length(setDiags));
%
for ii=1:length(setDiags);
  if iscell(setDiags{ii});
    if ~doBudget&strcmp(setDiags{ii}{1},'D'); doDiags(ii)=0; end;
  elseif ~doBudget&strcmp(setDiags{ii},'D'); doDiags(ii)=0; 
  elseif ~doProfiles&strcmp(setDiags{ii},'profiles'); doDiags(ii)=0;   
  elseif ~doCost&strcmp(setDiags{ii},'cost'); doDiags(ii)=0;   
  elseif ~doCtrl&strcmp(setDiags{ii},'controls'); doDiags(ii)=0;   
  end;
end;
%
setDiags={setDiags{find(doDiags)}};

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%initialize tex file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if ~isempty(dirMatRef);
    tmp1=dirMat; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
    myTitle={'gcmfaces analysis of the solution in',tmp1};
    tmp1=dirMatRef; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
    myTitle={myTitle{:},' minus ',tmp1};
end;

if isempty(dirMatRef)&~isempty(dir([dirMat '../README']));
  [rdm]=read_readme([dirMat '../README']);
else;
  rdm=[];
end;

write2tex(fileTex,0,myTitle,rdm);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%augment tex file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for ii=1:length(setDiags);

  if iscell(setDiags{ii});
    ttl=input('specify tex file section title:\n');
    write2tex(fileTex,1,ttl,1);
    diags_display({dirMat,dirMatRef},setDiags{ii},dirTex,nameTex);

  elseif strcmp(setDiags{ii},'profiles');
        %in situ profiles fit
        write2tex(fileTex,1,'fit to data',1);
        write2tex(fileTex,1,'fit to in situ data',2);
        insitu_diags(dirMat,0,dirTex,nameTex);
    
  elseif strcmp(setDiags{ii},'cost');
        %altimeter fit
        if isempty(find(strcmp(setDiags,'profiles')));
          write2tex(fileTex,1,'fit to data',1);
        end;
        write2tex(fileTex,1,'fit to altimeter data (RADS)',2);
        cost_altimeter_disp(dirMat,0,'etaglo',dirTex,nameTex);
        cost_altimeter_disp(dirMat,2,'modMobs',dirTex,nameTex);
        cost_altimeter_disp(dirMat,1,'modMobs',dirTex,nameTex);
        cost_altimeter_disp(dirMat,3,'modMobs',dirTex,nameTex);
        cost_altimeter_disp(dirMat,1,'obs',dirTex,nameTex);
        cost_altimeter_disp(dirMat,1,'mod',dirTex,nameTex);
        
        %other cost terms
        write2tex(fileTex,1,'fit to sst data',2);
        cost_sst('',dirMat,0,dirTex,nameTex);
        write2tex(fileTex,1,'fit to grace r4 data',2);
        cost_bp('',dirMat,0,dirTex,nameTex);
        write2tex(fileTex,1,'fit to seaice data',2);
        cost_seaicearea('',dirMat,0,dirTex,nameTex);
    
  elseif strcmp(setDiags{ii},'A');
      write2tex(fileTex,1,'volume, heat and salt transports',1);
      diags_display({dirMat,dirMatRef},'A',dirTex,nameTex);

  elseif strcmp(setDiags{ii},'B');
    write2tex(fileTex,1,'mean and variance maps',1);
    diags_display({dirMat,dirMatRef},'B',dirTex,nameTex);

  elseif strcmp(setDiags{ii},'C');
    write2tex(fileTex,1,'global, zonal, regional averages',1);
    diags_display({dirMat,dirMatRef},'C',dirTex,nameTex);

  elseif strcmp(setDiags{ii},'D');
        budget_list=1;
        if ~isempty(dir([dirMat budgetList]));
            eval(['load ' dirMat budgetList ';']);
        end;
        for kk=budget_list;
            if kk==1;
                tmp1='(top to bottom)';
            else;
                tmp1=sprintf('(%im to bottom)',round(-mygrid.RF(kk)));
            end;
            write2tex(fileTex,1,['budgets : volume, heat and salt ' tmp1],1);
            diags_display({dirMat,dirMatRef},{'D',kk},dirTex,nameTex);
        end;

  elseif strcmp(setDiags{ii},'MLD');
      write2tex(fileTex,1,'mixed layer depth fields',1);
      diags_display({dirMat,dirMatRef},'MLD',dirTex,nameTex);
  elseif strcmp(setDiags{ii},'SEAICE');
      write2tex(fileTex,1,'seaice and snow fields',1);
      diags_display({dirMat,dirMatRef},'SEAICE',dirTex,nameTex);
  elseif strcmp(setDiags{ii},'controls');
        %controls
        write2tex(fileTex,1,'controls',1);
        cost_xx('',dirMat,0,dirTex,nameTex);

  else;
    ttl=input('specify tex file section title:\n');
    write2tex(fileTex,1,ttl,1);
    diags_display({dirMat,dirMatRef},setDiags{ii},dirTex,nameTex);

  end;%if iscell(setDiags{ii});
end;%for ii=1:length(setDiags);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finalize tex file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

write2tex(fileTex,4);
%if isunix;
%write2tex(fileTex,5);
%else;
%fprintf(['\nUsers of pc or mac computers are left to compile the tex file :\n' fileTex '\noutside of matlab \n \n']);
%end;
fprintf(['\n Tex file is ready to be compiled :\n ' fileTex '\n']);

%%%%% get README text information

function [rdm]=read_readme(filReadme);

gcmfaces_global;

rdm=[];

fid=fopen(filReadme,'rt');
while ~feof(fid);
    nn=length(rdm);
    rdm{nn+1} = fgetl(fid);
end;
fclose(fid);


1	gforget	1.1	function []=diags_driver_tex(dirMat,setDiags,dirTex,nameTex);
2	gforget	1.17	% DIAGS_DRIVER_TEX(dirMat,setDiags,dirTex,nameTex)
3			%
4			% displays multiple sets of diagnostics (setDiags={'profiles',
5			% 'cost','A','B','C','MLD','D'} by default) from the results
6			% stored in dirMat and outputs the plots to tex
7			% ([dirTex nameTex '.tex'] and [dirTex nameTex '*.eps'])
8			%
9			% setDiags is the choice of diagnostics (cell) that may include
10	gforget	1.1	% 'profiles') model to insitu data comparison
11			% 'cost') ssh etc. cost functions
12			% 'A') trasnports
13			% 'B') air-sea fluxes
14			% 'C') state variables
15	gforget	1.17	% 'D') global and hemispheric budgets
16			% 'MLD') mixed layer depths
17	gforget	1.5	% 'SEAICE') seaice fields
18			% 'controls') control vector adjustments
19	gforget	1.1
20			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
21			%determine input/output params:
22			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
23
24			%directory names:
25			if iscell(dirMat); dirMatRef=dirMat{2}; dirMat=dirMat{1}; end;
26			dirMat=[dirMat '/'];
27			if isempty(who('dirMatRef')); dirMatRef='';
28			elseif ~isempty(dirMatRef); dirMatRef=[dirMatRef '/'];
29			end;
30			if isempty(who('dirTex')); dirTex=''; else; dirTex=[dirTex '/']; end;
31			if isempty(who('nameTex')); nameTex='myPlots'; end;
32
33	gforget	1.10	%set default setDiags and myTitle
34	gforget	1.1	if isempty(setDiags);
35	gforget	1.14	setDiags={'profiles','cost','A','B','C','MLD','D'};
36	gforget	1.5	if strcmp(nameTex,'myPlots'); nameTex='standardAnalysis'; end;
37	gforget	1.10	%
38			tmp1=dirMat; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
39			myTitle={'gcmfaces','standard analysis of the solution in',tmp1};
40			else;
41			tmp1=dirMat; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
42			myTitle={'gcmfaces analysis of the solution in',tmp1};
43	gforget	1.1	end;
44
45			%set fileTex and create dirTex if needed
46			if isempty(dirTex); error('dirTex must be specified'); end;
47			fileTex=[dirTex nameTex '.tex'];
48	gforget	1.11	if isempty(dir(dirTex)); mkdir(dirTex); end;
49	gforget	1.1
50			%%%%%%%%%%%%%%%%%%%%%%
51			%load grid and params:
52			%%%%%%%%%%%%%%%%%%%%%%
53
54			gcmfaces_global; global myparms;
55			test1=~isempty(dir([dirMat 'basic_diags_ecco_mygrid.mat']));
56			test2=~isempty(dir([dirMat 'diags_grid_parms.mat']));
57			if ~test1&~test2;
58			error('missing diags_grid_parms.mat')
59			elseif test2;
60			nameGrid='diags_grid_parms.mat';
61			suffDiag='diags_set_';
62			budgetList='diags_select_budget_list.mat';
63			else;
64			nameGrid='basic_diags_ecco_mygrid.mat';
65			suffDiag='basic_diags_ecco_';
66			budgetList='basic_diags_ecco_budget_list.mat';
67			end;
68
69			%here we always reload the grid from dirMat to make sure the same one is used throughout
70			eval(['load ' dirMat nameGrid ';']);
71
72			%in case mygrid.memoryLimit=1, load the stuff that was not saved to diags_grid_parms.mat
73			if mygrid.memoryLimit==1;
74			list0={'hFacS','hFacW'};
75			for iFld=1:length(list0);
76			eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([mygrid.dirGrid ''' list0{iFld} '*'']);']);
77			end;
78			%
79			mygrid.hFacCsurf=mygrid.hFacC;
80			for ff=1:mygrid.hFacC.nFaces; mygrid.hFacCsurf{ff}=mygrid.hFacC{ff}(:,:,1); end;
81			%
82			mskC=mygrid.hFacC; mskC(mskC==0)=NaN; mskC(mskC>0)=1; mygrid.mskC=mskC;
83			mskW=mygrid.hFacW; mskW(mskW==0)=NaN; mskW(mskW>0)=1; mygrid.mskW=mskW;
84			mskS=mygrid.hFacS; mskS(mskS==0)=NaN; mskS(mskS>0)=1; mygrid.mskS=mskS;
85			%
86			gcmfaces_lines_zonal;
87			mygrid.LATS=[mygrid.LATS_MASKS.lat]';
88	gforget	1.16	[lonPairs,latPairs,names]=gcmfaces_lines_pairs;
89	gforget	1.1	gcmfaces_lines_transp(lonPairs,latPairs,names);
90			end;
91
92
93			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
94			%finalize listDiags
95			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
96
97			%determined where to display anomalies between runs
98			doAnomalies=~isempty(dirMatRef);
99
100			%check if result is here to plot
101			doBudget=~isempty(dir([dirMat suffDiag 'D_*']));
102	gforget	1.4	doBudget=~isempty(dir([dirMat suffDiag 'D_*']))\|...
103			~isempty(dir([dirMat 'diags_set_D/' suffDiag 'D_*']));
104			doProfiles=~isempty(dir([dirMat 'insitu_cost_all.mat']))\|...
105			~isempty(dir([dirMat 'cost/insitu_cost_all.mat']));
106			doCost=~isempty(dir([dirMat 'cost_altimeter_obs.mat']))\|...
107			~isempty(dir([dirMat 'cost/cost_altimeter_obs.mat']));
108			doCtrl=~isempty(dir([dirMat 'cost_xx_aqh.mat']))\|...
109			~isempty(dir([dirMat 'cost/cost_xx_aqh.mat']));
110	gforget	1.1
111			%the following have no code for diff between runs
112			if doAnomalies;
113			doProfiles=0;
114			doCost=0;
115			doCtrl=0;
116			end;
117
118			%reduce setDiags if needed:
119			doDiags=ones(1,length(setDiags));
120			%
121			for ii=1:length(setDiags);
122			if iscell(setDiags{ii});
123			if ~doBudget&strcmp(setDiags{ii}{1},'D'); doDiags(ii)=0; end;
124			elseif ~doBudget&strcmp(setDiags{ii},'D'); doDiags(ii)=0;
125			elseif ~doProfiles&strcmp(setDiags{ii},'profiles'); doDiags(ii)=0;
126			elseif ~doCost&strcmp(setDiags{ii},'cost'); doDiags(ii)=0;
127			elseif ~doCtrl&strcmp(setDiags{ii},'controls'); doDiags(ii)=0;
128			end;
129			end;
130			%
131			setDiags={setDiags{find(doDiags)}};
132
133			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
134			%initialize tex file
135			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
136
137			if ~isempty(dirMatRef);
138			tmp1=dirMat; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
139	gforget	1.10	myTitle={'gcmfaces analysis of the solution in',tmp1};
140	gforget	1.1	tmp1=dirMatRef; tmp2=strfind(tmp1,'_'); tmp1(tmp2)=' ';
141			myTitle={myTitle{:},' minus ',tmp1};
142			end;
143	gforget	1.10
144			if isempty(dirMatRef)&~isempty(dir([dirMat '../README']));
145			[rdm]=read_readme([dirMat '../README']);
146			else;
147			rdm=[];
148			end;
149
150			write2tex(fileTex,0,myTitle,rdm);
151	gforget	1.1
152			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
153			%augment tex file
154			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
155
156			for ii=1:length(setDiags);
157
158			if iscell(setDiags{ii});
159			ttl=input('specify tex file section title:\n');
160			write2tex(fileTex,1,ttl,1);
161			diags_display({dirMat,dirMatRef},setDiags{ii},dirTex,nameTex);
162
163			elseif strcmp(setDiags{ii},'profiles');
164			%in situ profiles fit
165	gforget	1.5	write2tex(fileTex,1,'fit to data',1);
166			write2tex(fileTex,1,'fit to in situ data',2);
167	gforget	1.2	insitu_diags(dirMat,0,dirTex,nameTex);
168	gforget	1.1
169			elseif strcmp(setDiags{ii},'cost');
170			%altimeter fit
171	gforget	1.9	if isempty(find(strcmp(setDiags,'profiles')));
172			write2tex(fileTex,1,'fit to data',1);
173			end;
174			write2tex(fileTex,1,'fit to altimeter data (RADS)',2);
175	gforget	1.14	cost_altimeter_disp(dirMat,0,'etaglo',dirTex,nameTex);
176	gforget	1.2	cost_altimeter_disp(dirMat,2,'modMobs',dirTex,nameTex);
177			cost_altimeter_disp(dirMat,1,'modMobs',dirTex,nameTex);
178			cost_altimeter_disp(dirMat,3,'modMobs',dirTex,nameTex);
179			cost_altimeter_disp(dirMat,1,'obs',dirTex,nameTex);
180			cost_altimeter_disp(dirMat,1,'mod',dirTex,nameTex);
181	gforget	1.1
182			%other cost terms
183	gforget	1.5	write2tex(fileTex,1,'fit to sst data',2);
184	gforget	1.2	cost_sst('',dirMat,0,dirTex,nameTex);
185	gforget	1.9	write2tex(fileTex,1,'fit to grace r4 data',2);
186	gforget	1.2	cost_bp('',dirMat,0,dirTex,nameTex);
187	gforget	1.9	write2tex(fileTex,1,'fit to seaice data',2);
188	gforget	1.3	cost_seaicearea('',dirMat,0,dirTex,nameTex);
189	gforget	1.1
190			elseif strcmp(setDiags{ii},'A');
191			write2tex(fileTex,1,'volume, heat and salt transports',1);
192			diags_display({dirMat,dirMatRef},'A',dirTex,nameTex);
193
194			elseif strcmp(setDiags{ii},'B');
195			write2tex(fileTex,1,'mean and variance maps',1);
196			diags_display({dirMat,dirMatRef},'B',dirTex,nameTex);
197
198			elseif strcmp(setDiags{ii},'C');
199			write2tex(fileTex,1,'global, zonal, regional averages',1);
200			diags_display({dirMat,dirMatRef},'C',dirTex,nameTex);
201
202			elseif strcmp(setDiags{ii},'D');
203			budget_list=1;
204			if ~isempty(dir([dirMat budgetList]));
205			eval(['load ' dirMat budgetList ';']);
206			end;
207			for kk=budget_list;
208			if kk==1;
209			tmp1='(top to bottom)';
210			else;
211			tmp1=sprintf('(%im to bottom)',round(-mygrid.RF(kk)));
212			end;
213			write2tex(fileTex,1,['budgets : volume, heat and salt ' tmp1],1);
214			diags_display({dirMat,dirMatRef},{'D',kk},dirTex,nameTex);
215			end;
216
217	gforget	1.5	elseif strcmp(setDiags{ii},'MLD');
218	gforget	1.6	write2tex(fileTex,1,'mixed layer depth fields',1);
219	gforget	1.5	diags_display({dirMat,dirMatRef},'MLD',dirTex,nameTex);
220			elseif strcmp(setDiags{ii},'SEAICE');
221			write2tex(fileTex,1,'seaice and snow fields',1);
222			diags_display({dirMat,dirMatRef},'SEAICE',dirTex,nameTex);
223	gforget	1.1	elseif strcmp(setDiags{ii},'controls');
224			%controls
225			write2tex(fileTex,1,'controls',1);
226	gforget	1.2	cost_xx('',dirMat,0,dirTex,nameTex);
227	gforget	1.1
228			else;
229			ttl=input('specify tex file section title:\n');
230			write2tex(fileTex,1,ttl,1);
231			diags_display({dirMat,dirMatRef},setDiags{ii},dirTex,nameTex);
232
233			end;%if iscell(setDiags{ii});
234			end;%for ii=1:length(setDiags);
235
236			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
237			%finalize tex file
238			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
239
240			write2tex(fileTex,4);
241	gforget	1.15	%if isunix;
242			%write2tex(fileTex,5);
243			%else;
244	gforget	1.18	%fprintf(['\nUsers of pc or mac computers are left to compile the tex file :\n' fileTex '\noutside of matlab \n \n']);
245	gforget	1.15	%end;
246	gforget	1.18	fprintf(['\n Tex file is ready to be compiled :\n ' fileTex '\n']);
247	gforget	1.1
248	gforget	1.10	%%%%% get README text information
249
250			function [rdm]=read_readme(filReadme);
251
252			gcmfaces_global;
253
254			rdm=[];
255
256			fid=fopen(filReadme,'rt');
257			while ~feof(fid);
258			nn=length(rdm);
259			rdm{nn+1} = fgetl(fid);
260			end;
261			fclose(fid);
262
263