matlab_class/gcmfaces_maps/m_map_gcmfaces.m

function []=m_map_gcmfaces(fld,varargin);
%object:    gcmfaces front end to m_map
%inputs:    fld is the 2D field to be mapped, or a cell (see below).
%optional:  proj is either the index (integer; 0 by default) of pre-defined
%               projection(s) or parameters to pass to m_proj (cell)
%more so:   other optional paramaters can be provided ONLY AFTER proj,
%           and must take the following form {'name',param1,param2,...}
%           those that are currently active are
%               {'myCaxis',myCaxis} is the color axis ('auto' by default)
%               {'do_m_coast',do_m_coast} adds call to m_coast (1; default) or not (0).
%               {'myCmap',myCmap} is the colormap name ('jet' by default)
%               {'doHold',1} indicates to 'hold on' and superimpose e.g. a contour
%               {'doCbar',1} indicates to include the colorbar (1 by default)
%               {'doLabl',1} indicates to include the axes labels (1 by default)
%               {'doFit',1} indicates to exclude white space (0 by default)
%               {'myShading','flat'} indicates to use flat shading instead of interp
%
%notes:     - for proj==0 (i.e. the default) three panels will be plotted :
%           a mercator projection over mid-latitude, and two polar
%           stereographic projections. The other predefined proj are
%               -1  mollweide or cylindrical projection
%               1   mercator projection only
%               2   arctic stereographic projection only
%               3   antarctic stereographic projection only
%               1.1 and 1.2 are atlantic mercator projections
%           - myTitle is currently not used; it will take the form
%              {'myTitle',myTitle} is the title (none by default)
%           - myCaxis can be specified with more than 2 values, in which
%             case gcmfaces_cb will be used.
%           - if fld is a 2D field we use pcolor to display it
%           - if fld is a cell one can specify the plotting tool.
%             The cell  must then have the form {plotType,FLD,varargin}
%             where plotType is e.g. 'pcolor' or 'contour', FLD is
%             the 2D field to be plotted, and varargin are options
%             to pass over to the contour command.
%           - Hence .e.g. fld={'contour',mygrid.Depth,'r'} will draw
%             red contours, while fld=mygrid.Depth will color shade.

%check that m_map is in the path
aa=which('m_proj'); if isempty(aa); error('this function requires m_map that is missing'); end;

global mygrid;

%get optional parameters
if nargin>1; proj=varargin{1}; else; proj=0; end;
if iscell(proj);
    error('not yet implemented');
else;
    choicePlot=proj;
end;
%determine the type of plot
if iscell(fld); myPlot=fld{1}; else; myPlot='pcolor'; fld={'pcolor',fld}; end;
%set more optional paramaters to default values
myCaxis=[]; myTitle=''; myShading='interp'; myCmap='jet';
do_m_coast=1; doHold=0; doCbar=1; doLabl=1; doFit=0;
%set more optional paramaters to user defined values
for ii=2:nargin-1;
    if ~iscell(varargin{ii});
        warning('inputCheck:m_map_gcmfaces_1',...
            ['As of june 2011, m_map_gcmfaces expects \n'...
            '         its optional parameters as cell arrays. \n'...
            '         Argument no. ' num2str(ii+1) ' was ignored \n'...
            '         Type ''help m_map_gcmfaces'' for details.']);
    elseif ~ischar(varargin{ii}{1});
        warning('inputCheck:m_map_gcmfaces_2',...
            ['As of june 2011, m_map_gcmfaces expects \n'...
            '         its optional parameters cell arrays \n'...
            '         to start with character string. \n'...
            '         Argument no. ' num2str(ii+1) ' was ignored \n'...
            '         Type ''help m_map_gcmfaces'' for details.']);
    else;
        if strcmp(varargin{ii}{1},'myCaxis')|...
                strcmp(varargin{ii}{1},'myCmap')|...
                strcmp(varargin{ii}{1},'myShading')|...
                strcmp(varargin{ii}{1},'myTitle')|...
                strcmp(varargin{ii}{1},'doHold')|...
                strcmp(varargin{ii}{1},'doCbar')|...
                strcmp(varargin{ii}{1},'doLabl')|...
                strcmp(varargin{ii}{1},'doFit')|...
                strcmp(varargin{ii}{1},'do_m_coast');
            eval([varargin{ii}{1} '=varargin{ii}{2};']);
        else;
            warning('inputCheck:m_map_gcmfaces_3',...
                ['unknown option ''' varargin{ii}{1} ''' was ignored']);
        end;
    end;
end;

%make parameter inferences
if length(myCaxis)==0;
    plotCBAR=0;
elseif length(myCaxis)==2;
    plotCBAR=1;
else;
    plotCBAR=2;
end;
%
if choicePlot==0&~doHold;
    clf;
elseif ~doHold;
    cla;
else;
    hold on;
end;

%re-group param:
param.plotCBAR=plotCBAR;
param.myCmap=myCmap;
param.myCaxis=myCaxis;
param.do_m_coast=do_m_coast;
param.doHold=doHold;
param.doCbar=doCbar;
param.doLabl=doLabl;
param.doFit=doFit;
param.myPlot=myPlot;

%do the plotting:
if (choicePlot~=0&choicePlot~=1&choicePlot~=2&choicePlot~=3);
    do_my_plot(fld,param,choicePlot,myShading);
end;%if choicePlot==0|choicePlot==1;

if choicePlot==0; subplot(2,1,1); end;
if choicePlot==0|choicePlot==1;
    if mygrid.nFaces~=5&mygrid.nFaces~=6;
       do_my_plot(fld,param,1,myShading);
    else;
       do_my_plot(fld,param,1.1,myShading);
    end;
end;%if choicePlot==0|choicePlot==1;

if choicePlot==0; subplot(2,2,3); end;
if choicePlot==0|choicePlot==2;
    do_my_plot(fld,param,2,myShading);
end;%if choicePlot==0|choicePlot==1;

if choicePlot==0; subplot(2,2,4); end;
if choicePlot==0|choicePlot==3;
    do_my_plot(fld,param,3,myShading);
end;%if choicePlot==0|choicePlot==1;

if plotCBAR==2&strcmp(myPlot,'pcolor')&doCbar;
    cbar=gcmfaces_cmap_cbar(myCaxis,{'myCmap',myCmap});
    if choicePlot==0;
        set(cbar,'Position',[0.88 0.15 0.02 0.75]);
    elseif choicePlot==-1;
        set(cbar,'Position',[0.88 0.35 0.02 0.3]);
    elseif choicePlot==1;
        set(cbar,'Position',[0.88 0.34 0.02 0.35]);
    elseif choicePlot==1.2;
        set(cbar,'Position',[0.88 0.34 0.02 0.35]);
    elseif choicePlot==1.3;
        set(cbar,'Position',[0.88 0.34 0.02 0.35]);
    else;
        set(cbar,'Position',[0.88 0.3 0.02 0.4]);
    end;
end;

if doFit;
    if doLabl&~doCbar; 
        set(gca,'LooseInset',[0.05 0.02 0.03 0.03]);
    else; 
        set(gca,'LooseInset',[0.01 0 0.03 0.03]);
    end;
    tmp1=get(gca,'PlotBoxAspectRatio');
    set(gcf,'PaperPosition',[0 4 8 8/tmp1(1)]);
end;

function []=do_my_plot(fld,param,proj,shad);

gcmfaces_global;

%default m_grid params:
xloc='bottom'; xtic=[-180:60:180]; xticlab=1;
yloc='left';   ytic=[-60:30:60];   yticlab=1;

%choice of projection:
if proj==-1;
    %%m_proj('Miller Cylindrical','lat',[-90 90]);
    %m_proj('Equidistant cylindrical','lat',[-90 90]);
    %m_proj('mollweide','lon',[-180 180],'lat',[-80 80]);
    m_proj('mollweide','lon',[-180 180],'lat',[-88 88]);
    myConv='pcol'; xticlab=0; yticlab=0;
elseif proj==1;
    m_proj('Mercator','lat',[-70 70]);
    myConv='pcol';
elseif proj==1.1;
    m_proj('Mercator','lat',[-70 70],'lon',[0 360]+20);
    myConv='pcol';
    xtic=[-360:60:360]; ytic=[-60:30:60];
elseif proj==1.2;
    m_proj('Equidistant cylindrical','lat',[-90 90],'lon',[0 360]+20);
    myConv='pcol';
    xtic=[-360:60:360]; ytic=[-90:30:90];
elseif proj==1.3;
    m_proj('Equidistant cylindrical','lat',[-90 -12],'lon',[0 360]+20);
    myConv='pcol';
    xtic=[-360:60:360]; ytic=[-90:30:90];
elseif proj==2;
    m_proj('Stereographic','lon',0,'lat',90,'rad',40);
    myConv='convert2arctic';
    yloc='bottom'; ytic=[50:10:90];
elseif proj==2.1;
    m_proj('Stereographic','lon',0,'lat',90,'rad',60);
    myConv='convert2arctic';
    yloc='bottom'; ytic=[30:10:90];
elseif proj==3;
    m_proj('Stereographic','lon',0,'lat',-90,'rad',40);
    myConv='convert2southern';
    xloc='top'; ytic=[-90:10:-50];
elseif proj==3.1;
    m_proj('Stereographic','lon',0,'lat',-90,'rad',60);
    myConv='convert2southern';
    xloc='top'; ytic=[-90:10:-30];
elseif proj==4.1;
    m_proj('mollweide','lat',[25 75],'lon',[-100 30]);
    myConv='pcol';
    xtic=[-100:20:30]; ytic=[30:10:70];
elseif proj==4.2;
    m_proj('mollweide','lat',[-30 30],'lon',[-65 20]);
    myConv='pcol';
    xtic=[-60:10:20]; ytic=[-30:10:30];
elseif proj==4.3;
    m_proj('mollweide','lat',[-75 -25],'lon',[-75 25]);
    myConv='pcol';
    xtic=[-70:20:20]; ytic=[-70:10:-30]; xloc='top';
elseif proj==4.4;
    m_proj('mollweide','lat',[25 75],'lon',[-240 -110]);
    myConv='pcol';
    xtic=[-240:20:-120]; ytic=[30:10:60]; 
elseif proj==4.5;
    m_proj('mollweide','lat',[-30 30],'lon',[-240 -70]);
    myConv='pcol';
    xtic=[-240:20:-70]; ytic=[-30:10:30];
elseif proj==4.6;
    m_proj('mollweide','lat',[-75 -25],'lon',[-215 -60]);
    myConv='pcol';
    xtic=[-210:20:-70]; ytic=[-70:10:-30]; xloc='top'; 
elseif proj==4.7;
    m_proj('mollweide','lat',[-30 30],'lon',[15 155]);
    myConv='pcol';
    xtic=[10:20:160]; ytic=[-30:10:30]; xloc='top'; 
elseif proj==4.8;
    m_proj('mollweide','lat',[-75 -25],'lon',[15 155]);
    myConv='pcol';
    xtic=[10:20:160]; ytic=[-70:10:-30]; xloc='top';
else;
    error('undefined projection');
end;

if ~param.doLabl; xticlab=0; yticlab=0;end; %omit labels

m_grid_opt=['''XaxisLocation'',xloc,''YaxisLocation'',yloc,''xtick'',xtic,''ytick'',ytic'];
if xticlab==0; m_grid_opt=[m_grid_opt ',''xticklabel'',[]']; end;
if yticlab==0; m_grid_opt=[m_grid_opt ',''yticklabel'',[]']; end;

if strcmp(param.myPlot,'pcolor')|strcmp(param.myPlot,'contour');
    x=mygrid.XC;
    %mask out the XC padded zeros
    if isfield(mygrid,'xtrct');
      pt1=mygrid.xtrct.pt1face; pt2=mygrid.xtrct.pt2face;
      if pt1~=pt2;
        for iF=1:6; if iF~=pt1&iF~=pt2; x{iF}(:)=NaN; end; end; 
      end;
      x(x<0)=x(x<0)+360;
    end;
    if strcmp(myConv,'pcol');
        [xx,yy,z]=convert2pcol(x,mygrid.YC,fld{2});
    else;
        eval(['xx=' myConv '(x);']);
        eval(['yy=' myConv '(mygrid.YC);']);
        eval(['z=' myConv '(fld{2});']);
    end;
    [x,y]=m_ll2xy(xx,yy);
    if strcmp(param.myPlot,'pcolor');
        if sum(~isnan(x(:)))>0; pcolor(x,y,z); eval(['shading ' shad ';']); end;
        if param.plotCBAR==0;
            colormap(param.myCmap); if param.doCbar; colorbar; end;
        elseif param.plotCBAR==1;
            caxis(param.myCaxis); colormap(param.myCmap); if param.doCbar; colorbar; end;
        else;
            cbar=gcmfaces_cmap_cbar(param.myCaxis,{'myCmap',param.myCmap}); delete(cbar);
        end;
        if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
        eval(['m_grid(' m_grid_opt ');']);
    elseif strcmp(param.myPlot,'contour');
        if ~param.doHold;
            if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
            eval(['m_grid(' m_grid_opt ');']);
        end;
        if length(fld)==2; fld{3}='k'; end;
        hold on; contour(x,y,z,fld{3:end});
    end;
elseif strcmp(param.myPlot,'plot');
    if ~param.doHold;
        if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
        eval(['m_grid(' m_grid_opt ');']);
    end;
    [x,y]=m_ll2xy(fld{2},fld{3});
    hold on; plot(x,y,fld{4:end});
elseif strcmp(param.myPlot,'scatter');
    if ~param.doHold;
        if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
        eval(['m_grid(' m_grid_opt ');']);
    end;
    [x,y]=m_ll2xy(fld{2},fld{3});
    hold on; scatter(x,y,fld{4:end});
elseif strcmp(param.myPlot,'text');
    if ~param.doHold;
        if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
        eval(['m_grid(' m_grid_opt ');']);
    end;
    [x,y]=m_ll2xy(fld{2},fld{3});
    if length(fld)>4; cc=fld{5}; else; cc='k'; end;
    hold on; hold on; text(x,y,fld{4},'Color',cc,'FontSize',16,'FontWeight','bold');
end;

%add tag spec. to map & proj generated with this routine
set(gca,'Tag',['gfMap' num2str(proj)]);

1	gforget	1.1	function []=m_map_gcmfaces(fld,varargin);
2	gforget	1.3	%object: gcmfaces front end to m_map
3	gforget	1.7	%inputs: fld is the 2D field to be mapped, or a cell (see below).
4	gforget	1.4	%optional: proj is either the index (integer; 0 by default) of pre-defined
5			% projection(s) or parameters to pass to m_proj (cell)
6	gforget	1.7	%more so: other optional paramaters can be provided ONLY AFTER proj,
7	gforget	1.4	% and must take the following form {'name',param1,param2,...}
8	gforget	1.6	% those that are currently active are
9	gforget	1.4	% {'myCaxis',myCaxis} is the color axis ('auto' by default)
10			% {'do_m_coast',do_m_coast} adds call to m_coast (1; default) or not (0).
11	gforget	1.6	% {'myCmap',myCmap} is the colormap name ('jet' by default)
12	gforget	1.7	% {'doHold',1} indicates to 'hold on' and superimpose e.g. a contour
13	gforget	1.8	% {'doCbar',1} indicates to include the colorbar (1 by default)
14			% {'doLabl',1} indicates to include the axes labels (1 by default)
15			% {'doFit',1} indicates to exclude white space (0 by default)
16	gforget	1.9	% {'myShading','flat'} indicates to use flat shading instead of interp
17	gforget	1.3	%
18	gforget	1.4	%notes: - for proj==0 (i.e. the default) three panels will be plotted :
19			% a mercator projection over mid-latitude, and two polar
20	gforget	1.7	% stereographic projections. The other predefined proj are
21	gforget	1.3	% -1 mollweide or cylindrical projection
22			% 1 mercator projection only
23			% 2 arctic stereographic projection only
24			% 3 antarctic stereographic projection only
25	gforget	1.7	% 1.1 and 1.2 are atlantic mercator projections
26	gforget	1.4	% - myTitle is currently not used; it will take the form
27			% {'myTitle',myTitle} is the title (none by default)
28	gforget	1.7	% - myCaxis can be specified with more than 2 values, in which
29	gforget	1.8	% case gcmfaces_cb will be used.
30	gforget	1.7	% - if fld is a 2D field we use pcolor to display it
31			% - if fld is a cell one can specify the plotting tool.
32			% The cell must then have the form {plotType,FLD,varargin}
33	gforget	1.8	% where plotType is e.g. 'pcolor' or 'contour', FLD is
34			% the 2D field to be plotted, and varargin are options
35	gforget	1.7	% to pass over to the contour command.
36	gforget	1.8	% - Hence .e.g. fld={'contour',mygrid.Depth,'r'} will draw
37	gforget	1.7	% red contours, while fld=mygrid.Depth will color shade.
38	gforget	1.1
39	gforget	1.2	%check that m_map is in the path
40			aa=which('m_proj'); if isempty(aa); error('this function requires m_map that is missing'); end;
41
42	gforget	1.1	global mygrid;
43
44	gforget	1.4	%get optional parameters
45			if nargin>1; proj=varargin{1}; else; proj=0; end;
46	gforget	1.7	if iscell(proj);
47			error('not yet implemented');
48	gforget	1.4	else;
49			choicePlot=proj;
50			end;
51	gforget	1.7	%determine the type of plot
52			if iscell(fld); myPlot=fld{1}; else; myPlot='pcolor'; fld={'pcolor',fld}; end;
53	gforget	1.4	%set more optional paramaters to default values
54	gforget	1.10	myCaxis=[]; myTitle=''; myShading='interp'; myCmap='jet';
55	gforget	1.9	do_m_coast=1; doHold=0; doCbar=1; doLabl=1; doFit=0;
56	gforget	1.4	%set more optional paramaters to user defined values
57			for ii=2:nargin-1;
58			if ~iscell(varargin{ii});
59			warning('inputCheck:m_map_gcmfaces_1',...
60			['As of june 2011, m_map_gcmfaces expects \n'...
61			' its optional parameters as cell arrays. \n'...
62			' Argument no. ' num2str(ii+1) ' was ignored \n'...
63			' Type ''help m_map_gcmfaces'' for details.']);
64			elseif ~ischar(varargin{ii}{1});
65			warning('inputCheck:m_map_gcmfaces_2',...
66			['As of june 2011, m_map_gcmfaces expects \n'...
67			' its optional parameters cell arrays \n'...
68			' to start with character string. \n'...
69			' Argument no. ' num2str(ii+1) ' was ignored \n'...
70			' Type ''help m_map_gcmfaces'' for details.']);
71			else;
72			if strcmp(varargin{ii}{1},'myCaxis')\|...
73	gforget	1.6	strcmp(varargin{ii}{1},'myCmap')\|...
74	gforget	1.9	strcmp(varargin{ii}{1},'myShading')\|...
75			strcmp(varargin{ii}{1},'myTitle')\|...
76	gforget	1.7	strcmp(varargin{ii}{1},'doHold')\|...
77	gforget	1.8	strcmp(varargin{ii}{1},'doCbar')\|...
78			strcmp(varargin{ii}{1},'doLabl')\|...
79			strcmp(varargin{ii}{1},'doFit')\|...
80	gforget	1.9	strcmp(varargin{ii}{1},'do_m_coast');
81	gforget	1.4	eval([varargin{ii}{1} '=varargin{ii}{2};']);
82			else;
83			warning('inputCheck:m_map_gcmfaces_3',...
84			['unknown option ''' varargin{ii}{1} ''' was ignored']);
85			end;
86			end;
87			end;
88	gforget	1.7
89	gforget	1.4	%make parameter inferences
90	gforget	1.7	if length(myCaxis)==0;
91			plotCBAR=0;
92			elseif length(myCaxis)==2;
93			plotCBAR=1;
94			else;
95			plotCBAR=2;
96			end;
97			%
98			if choicePlot==0&~doHold;
99			clf;
100			elseif ~doHold;
101			cla;
102			else;
103			hold on;
104			end;
105
106			%re-group param:
107			param.plotCBAR=plotCBAR;
108			param.myCmap=myCmap;
109			param.myCaxis=myCaxis;
110			param.do_m_coast=do_m_coast;
111			param.doHold=doHold;
112	gforget	1.8	param.doCbar=doCbar;
113			param.doLabl=doLabl;
114			param.doFit=doFit;
115	gforget	1.7	param.myPlot=myPlot;
116
117			%do the plotting:
118			if (choicePlot~=0&choicePlot~=1&choicePlot~=2&choicePlot~=3);
119	gforget	1.9	do_my_plot(fld,param,choicePlot,myShading);
120	gforget	1.7	end;%if choicePlot==0\|choicePlot==1;
121	gforget	1.1
122			if choicePlot==0; subplot(2,1,1); end;
123	gforget	1.7	if choicePlot==0\|choicePlot==1;
124	gforget	1.8	if mygrid.nFaces~=5&mygrid.nFaces~=6;
125	gforget	1.9	do_my_plot(fld,param,1,myShading);
126	gforget	1.8	else;
127	gforget	1.9	do_my_plot(fld,param,1.1,myShading);
128	gforget	1.8	end;
129	gforget	1.7	end;%if choicePlot==0\|choicePlot==1;
130	gforget	1.3
131	gforget	1.1	if choicePlot==0; subplot(2,2,3); end;
132	gforget	1.7	if choicePlot==0\|choicePlot==2;
133	gforget	1.9	do_my_plot(fld,param,2,myShading);
134	gforget	1.7	end;%if choicePlot==0\|choicePlot==1;
135	gforget	1.1
136			if choicePlot==0; subplot(2,2,4); end;
137	gforget	1.7	if choicePlot==0\|choicePlot==3;
138	gforget	1.9	do_my_plot(fld,param,3,myShading);
139	gforget	1.7	end;%if choicePlot==0\|choicePlot==1;
140	gforget	1.1
141	gforget	1.8	if plotCBAR==2&strcmp(myPlot,'pcolor')&doCbar;
142	gforget	1.6	cbar=gcmfaces_cmap_cbar(myCaxis,{'myCmap',myCmap});
143	gforget	1.7	if choicePlot==0;
144	gforget	1.14	set(cbar,'Position',[0.88 0.15 0.02 0.75]);
145	gforget	1.3	elseif choicePlot==-1;
146	gforget	1.14	set(cbar,'Position',[0.88 0.35 0.02 0.3]);
147	gforget	1.7	elseif choicePlot==1;
148	gforget	1.14	set(cbar,'Position',[0.88 0.34 0.02 0.35]);
149			elseif choicePlot==1.2;
150			set(cbar,'Position',[0.88 0.34 0.02 0.35]);
151	gforget	1.15	elseif choicePlot==1.3;
152			set(cbar,'Position',[0.88 0.34 0.02 0.35]);
153	gforget	1.14	else;
154			set(cbar,'Position',[0.88 0.3 0.02 0.4]);
155	gforget	1.7	end;
156			end;
157
158	gforget	1.8	if doFit;
159			if doLabl&~doCbar;
160			set(gca,'LooseInset',[0.05 0.02 0.03 0.03]);
161			else;
162			set(gca,'LooseInset',[0.01 0 0.03 0.03]);
163			end;
164			tmp1=get(gca,'PlotBoxAspectRatio');
165			set(gcf,'PaperPosition',[0 4 8 8/tmp1(1)]);
166			end;
167	gforget	1.7
168	gforget	1.9	function []=do_my_plot(fld,param,proj,shad);
169	gforget	1.7
170			gcmfaces_global;
171
172	gforget	1.8	%default m_grid params:
173			xloc='bottom'; xtic=[-180:60:180]; xticlab=1;
174			yloc='left'; ytic=[-60:30:60]; yticlab=1;
175
176			%choice of projection:
177	gforget	1.7	if proj==-1;
178			%%m_proj('Miller Cylindrical','lat',[-90 90]);
179			%m_proj('Equidistant cylindrical','lat',[-90 90]);
180			%m_proj('mollweide','lon',[-180 180],'lat',[-80 80]);
181			m_proj('mollweide','lon',[-180 180],'lat',[-88 88]);
182	gforget	1.8	myConv='pcol'; xticlab=0; yticlab=0;
183	gforget	1.7	elseif proj==1;
184			m_proj('Mercator','lat',[-70 70]);
185	gforget	1.8	myConv='pcol';
186			elseif proj==1.1;
187			m_proj('Mercator','lat',[-70 70],'lon',[0 360]+20);
188			myConv='pcol';
189			xtic=[-360:60:360]; ytic=[-60:30:60];
190	gforget	1.11	elseif proj==1.2;
191			m_proj('Equidistant cylindrical','lat',[-90 90],'lon',[0 360]+20);
192			myConv='pcol';
193			xtic=[-360:60:360]; ytic=[-90:30:90];
194			elseif proj==1.3;
195	gforget	1.15	m_proj('Equidistant cylindrical','lat',[-90 -12],'lon',[0 360]+20);
196	gforget	1.11	myConv='pcol';
197	gforget	1.15	xtic=[-360:60:360]; ytic=[-90:30:90];
198	gforget	1.7	elseif proj==2;
199			m_proj('Stereographic','lon',0,'lat',90,'rad',40);
200	gforget	1.8	myConv='convert2arctic';
201			yloc='bottom'; ytic=[50:10:90];
202			elseif proj==2.1;
203			m_proj('Stereographic','lon',0,'lat',90,'rad',60);
204			myConv='convert2arctic';
205			yloc='bottom'; ytic=[30:10:90];
206	gforget	1.7	elseif proj==3;
207			m_proj('Stereographic','lon',0,'lat',-90,'rad',40);
208	gforget	1.8	myConv='convert2southern';
209			xloc='top'; ytic=[-90:10:-50];
210			elseif proj==3.1;
211			m_proj('Stereographic','lon',0,'lat',-90,'rad',60);
212			myConv='convert2southern';
213			xloc='top'; ytic=[-90:10:-30];
214			elseif proj==4.1;
215			m_proj('mollweide','lat',[25 75],'lon',[-100 30]);
216			myConv='pcol';
217			xtic=[-100:20:30]; ytic=[30:10:70];
218			elseif proj==4.2;
219			m_proj('mollweide','lat',[-30 30],'lon',[-65 20]);
220			myConv='pcol';
221			xtic=[-60:10:20]; ytic=[-30:10:30];
222			elseif proj==4.3;
223			m_proj('mollweide','lat',[-75 -25],'lon',[-75 25]);
224			myConv='pcol';
225			xtic=[-70:20:20]; ytic=[-70:10:-30]; xloc='top';
226			elseif proj==4.4;
227			m_proj('mollweide','lat',[25 75],'lon',[-240 -110]);
228			myConv='pcol';
229			xtic=[-240:20:-120]; ytic=[30:10:60];
230			elseif proj==4.5;
231			m_proj('mollweide','lat',[-30 30],'lon',[-240 -70]);
232			myConv='pcol';
233			xtic=[-240:20:-70]; ytic=[-30:10:30];
234			elseif proj==4.6;
235			m_proj('mollweide','lat',[-75 -25],'lon',[-215 -60]);
236			myConv='pcol';
237			xtic=[-210:20:-70]; ytic=[-70:10:-30]; xloc='top';
238			elseif proj==4.7;
239			m_proj('mollweide','lat',[-30 30],'lon',[15 155]);
240			myConv='pcol';
241			xtic=[10:20:160]; ytic=[-30:10:30]; xloc='top';
242			elseif proj==4.8;
243			m_proj('mollweide','lat',[-75 -25],'lon',[15 155]);
244			myConv='pcol';
245			xtic=[10:20:160]; ytic=[-70:10:-30]; xloc='top';
246			else;
247			error('undefined projection');
248	gforget	1.3	end;
249	gforget	1.1
250	gforget	1.8	if ~param.doLabl; xticlab=0; yticlab=0;end; %omit labels
251
252			m_grid_opt=['''XaxisLocation'',xloc,''YaxisLocation'',yloc,''xtick'',xtic,''ytick'',ytic'];
253			if xticlab==0; m_grid_opt=[m_grid_opt ',''xticklabel'',[]']; end;
254			if yticlab==0; m_grid_opt=[m_grid_opt ',''yticklabel'',[]']; end;
255
256	gforget	1.7	if strcmp(param.myPlot,'pcolor')\|strcmp(param.myPlot,'contour');
257	gforget	1.8	x=mygrid.XC;
258			%mask out the XC padded zeros
259			if isfield(mygrid,'xtrct');
260			pt1=mygrid.xtrct.pt1face; pt2=mygrid.xtrct.pt2face;
261			if pt1~=pt2;
262			for iF=1:6; if iF~=pt1&iF~=pt2; x{iF}(:)=NaN; end; end;
263			end;
264			x(x<0)=x(x<0)+360;
265			end;
266			if strcmp(myConv,'pcol');
267			[xx,yy,z]=convert2pcol(x,mygrid.YC,fld{2});
268	gforget	1.7	else;
269	gforget	1.8	eval(['xx=' myConv '(x);']);
270			eval(['yy=' myConv '(mygrid.YC);']);
271			eval(['z=' myConv '(fld{2});']);
272	gforget	1.7	end;
273			[x,y]=m_ll2xy(xx,yy);
274			if strcmp(param.myPlot,'pcolor');
275	gforget	1.9	if sum(~isnan(x(:)))>0; pcolor(x,y,z); eval(['shading ' shad ';']); end;
276	gforget	1.8	if param.plotCBAR==0;
277			colormap(param.myCmap); if param.doCbar; colorbar; end;
278			elseif param.plotCBAR==1;
279			caxis(param.myCaxis); colormap(param.myCmap); if param.doCbar; colorbar; end;
280			else;
281			cbar=gcmfaces_cmap_cbar(param.myCaxis,{'myCmap',param.myCmap}); delete(cbar);
282	gforget	1.7	end;
283			if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
284	gforget	1.8	eval(['m_grid(' m_grid_opt ');']);
285	gforget	1.7	elseif strcmp(param.myPlot,'contour');
286			if ~param.doHold;
287			if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
288	gforget	1.8	eval(['m_grid(' m_grid_opt ');']);
289	gforget	1.7	end;
290			if length(fld)==2; fld{3}='k'; end;
291			hold on; contour(x,y,z,fld{3:end});
292			end;
293			elseif strcmp(param.myPlot,'plot');
294			if ~param.doHold;
295			if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
296	gforget	1.8	eval(['m_grid(' m_grid_opt ');']);
297	gforget	1.7	end;
298			[x,y]=m_ll2xy(fld{2},fld{3});
299			hold on; plot(x,y,fld{4:end});
300	gforget	1.16	elseif strcmp(param.myPlot,'scatter');
301			if ~param.doHold;
302			if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
303			eval(['m_grid(' m_grid_opt ');']);
304			end;
305			[x,y]=m_ll2xy(fld{2},fld{3});
306			hold on; scatter(x,y,fld{4:end});
307	gforget	1.8	elseif strcmp(param.myPlot,'text');
308			if ~param.doHold;
309			if param.do_m_coast; m_coast('patch',[1 1 1]*.7,'edgecolor','none'); end;
310			eval(['m_grid(' m_grid_opt ');']);
311			end;
312			[x,y]=m_ll2xy(fld{2},fld{3});
313			if length(fld)>4; cc=fld{5}; else; cc='k'; end;
314			hold on; hold on; text(x,y,fld{4},'Color',cc,'FontSize',16,'FontWeight','bold');
315	gforget	1.7	end;
316	gforget	1.1
317	gforget	1.8	%add tag spec. to map & proj generated with this routine
318			set(gca,'Tag',['gfMap' num2str(proj)]);
319