close all
clear
latrng = [21.06 70];
lonrng = [-160 -51.62];
fw = 0.24;
fh = fw*((latrng(2)-latrng(1))/180)/((lonrng(2)-lonrng(1))/360);
positions = [0.05 1-1*fh fw fh;
0.05 1-2*fh fw fh;
0.05 1-3*fh fw fh;
0.05 1-4*fh fw fh;
0.05+fw 1-1*fh fw fh;
0.05+fw 1-2*fh fw fh;
0.05+fw 1-3*fh fw fh;
0.05+fw 1-4*fh fw fh;
0.05+2*fw 1-1*fh fw fh;
0.05+2*fw 1-2*fh fw fh;
0.05+2*fw 1-3*fh fw fh;
0.05+2*fw 1-4*fh fw fh;];
Months = {'Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'};
figure('WindowState','maximized')
for Mn = 1:12
subplot('Position',positions(Mn,:))
worldmap(latrng,lonrng)
setm(gca,'Grid','off','MapProjection','miller','parallellabel','off','meridianlabel','off')
text( min(xlim)+0.1*(max(xlim)-min(xlim)), min(ylim)+0.1*(max(ylim)-min(ylim)), Months{Mn},'fontsize',10 )]
end