function plot_AIRS_MFx(period,lev,chead)
% period=['20181215-20181222' '20181223-20181231' '20190101-20190108']

cfile=['../data/AIRS_3DST_',period,'.mat'];
load(cfile,'mMFx')
load('../latlonzgrid.mat','long','latg','zg')
[xg,yg]=meshgrid(long,latg);
var=squeeze(mMFx(:,:,lev));

t=tiledlayout(1,5);
t.Padding='compact';
t.TileSpacing='tight';
ax1=nexttile([1 4]);
p1=contourf(ax1,xg,yg,var*1000,'LineStyle','none','LevelStep',0.1);
colorbar('southoutside')
load('../colormap_blue-red.mat')
colormap(CustomColormap1)
clim(ax1,[-5 5])
if lev>18
clim(ax1,[-2 2])
end
yticks(ax1,-90:15:90)
yticklabels(ax1,{'90S' '75S' '60S' '45S' '30S' '15S' 'EQ' '15N' '30N' '45N' '60N' '75N' '90N'})
xticks(ax1,-180:30:180)
xticklabels(ax1,{'180W' '150W' '120W' '90W' '60W' '30W' '0' '30E' '60E' '90E' '120E' '150E' '180E'})
hold on
load topo topo topomap1
topos = circshift(topo,[0 180]);
x=-179:180;
y=-89:90;
p0=contour(ax1,x,y,topos,[0 0],'-k','LineWidth',1);
set(ax1,'fontsize',10)
set(ax1,'fontname','Helvetica')
% set(ax1,'position',[0.1 0.4 0.7 0.7])

ax2=nexttile;
zmean = mean(var,2,'omitnan');
p2=plot(ax2,zmean*1000,squeeze(yg(:,1)));
ylim(ax2,[-90 90])
xlim(ax2,[-1.4 0.6])
xticks(ax2,-1.5:0.25:0.5)
if lev>18
    xlim([-0.4 0.2])
    xticks(ax2,-0.8:0.1:0.4)
end
yticks(ax2,-90:15:90)
yticklabels(ax2,{'90S' '75S' '60S' '45S' '30S' '15S' 'EQ' '15N' '30N' '45N' '60N' '75N' '90N'})
set(ax2,'yaxislocation','right')
set(ax2,'fontsize',10)
set(ax2,'fontname','Helvetica')

title(ax1,['Zonal Momentum Flux (mPa) ',chead],'FontSize',13,'FontWeight','normal','FontName','Helvetica')

% set(ax2,'position',[0.71 0.4 0.19 0.7])

end