function [figure1] = Plot_LocalBudgetsandErrors(L,Tobs,THRrel,M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t)


cm = m_colmap('diverging',10);
range = [-3.5:3.5];
TD = (Tobs(:,1)-THRrel(:,1))*1000;

figure1 = figure;

axes1 = axes('Parent',figure1,...
    'Position',[0.3 0.35 0.3 0.3]);
hold all
m_proj('robinson','lon',[-180 180]);
m_coast('patch',[1 1 1],'edgecolor','k');
hold all
[x1,y1] = m_ll2xy(L(:,1),L(:,2));
scatter(x1,y1,10,TD,'filled','MarkerEdgecolor','none')
m_grid('XTick',[],'YTick',[]);
load('roma.mat')
colormap(roma);
caxis([-3.5 3.5])
cb = colorbar('YTick',[-3:1:3],'Fontsize',8,'Location','Eastoutside')

axes2 = axes('Parent',figure1,...
    'Position',[0.1 0.6 0.2 0.15]);
[t] = Local_Budget(M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t,2,1);

axes3 = axes('Parent',figure1,...
    'Position',[0.4 0.65 0.2 0.15]);
[t] = Local_Budget(M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t,1,0);

axes4 = axes('Parent',figure1,...
    'Position',[0.7 0.65 0.2 0.15]);
[t] = Local_Budget(M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t,26,0);

axes5 = axes('Parent',figure1,...
    'Position',[0.1 0.25 0.2 0.15]);
[t] = Local_Budget(M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t,169,0);

axes6 = axes('Parent',figure1,...
    'Position',[0.4 0.2 0.2 0.15]);
[t] = Local_Budget(M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t,390,0);

axes7 = axes('Parent',figure1,...
    'Position',[0.7 0.2 0.2 0.15]);
[t] = Local_Budget(M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t,35,0);

axes8 = axes('Parent',figure1,...
    'Position',[0.775 0.425 0.2 0.15]);
[t] = Local_Budget(M,MHRrel,MSE,MSDSL,MSDSLSE,MBary,MBarySE,MGIA,MGIASE,MIBE,MALTrel,N,t,380,0);


end