You should have the files if you have the appropriate version of the Aerospace Toolbox.
NOTE — You have R2018a and the file is for the year 2020. MATLAB is good, however it is not prescient!
The code works here —
model_epoch = '2020';
decimal_year = 2020;
WMMFileName = 'astWMMResults_Epoch_2020_decyear_2020.mat';
load(WMMFileName);
% Read in continent land areas for plot overlay using Mapping Toolbox function, shaperead.
landAreas = shaperead('landareas.shp','UseGeoCoords',true);
% Plotting Earth's Magnetic Field Overlaid on Earth Maps
% Load plot formatting data for each of the magnetic parameters.
plotWMM = load('astPlotWMM.mat');
hX = figure;
set(hX,'Position',[0 0 827 620],'Color','white')
astPlotWMMContours( WMMResults, plotWMM, 1, landAreas, geod_lat, geod_lon, decimal_year, model_epoch)
WMMFileName = 'astWMMResults_Epoch_2020_decyear_2020.mat';
load(WMMFileName);
% Read in continent land areas for plot overlay using Mapping Toolbox function, shaperead.
landAreas = shaperead('landareas.shp','UseGeoCoords',true);
% Plotting Earth's Magnetic Field Overlaid on Earth Maps
% Load plot formatting data for each of the magnetic parameters.
plotWMM = load('astPlotWMM.mat');
hX = figure;
set(hX,'Position',[0 0 827 620],'Color','white')
astPlotWMMContours( WMMResults, plotWMM, 1, landAreas, geod_lat, geod_lon, decimal_year, model_epoch)
WMMFileName = 'astWMMResults_Epoch_2020_decyear_2020.mat';
load(WMMFileName);
% Read in continent land areas for plot overlay using Mapping Toolbox function, shaperead.
landAreas = shaperead('landareas.shp','UseGeoCoords',true);
% Plotting Earth's Magnetic Field Overlaid on Earth Maps
% Load plot formatting data for each of the magnetic parameters.
plotWMM = load('astPlotWMM.mat');
hX = figure;
set(hX,'Position',[0 0 827 620],'Color','white')
astPlotWMMContours( WMMResults, plotWMM, 1, landAreas, geod_lat, geod_lon, decimal_year, model_epoch)
WMMFileName = 'astWMMResults_Epoch_2020_decyear_2020.mat';
load(WMMFileName);
% Read in continent land areas for plot overlay using Mapping Toolbox function, shaperead.
landAreas = shaperead('landareas.shp','UseGeoCoords',true);
% Plotting Earth's Magnetic Field Overlaid on Earth Maps
% Load plot formatting data for each of the magnetic parameters.
plotWMM = load('astPlotWMM.mat');
hX = figure;
set(hX,'Position',[0 0 827 620],'Color','white')
astPlotWMMContours( WMMResults, plotWMM, 1, landAreas, geod_lat, geod_lon, decimal_year, model_epoch)
WMMFileName = 'astWMMResults_Epoch_2020_decyear_2020.mat';
load(WMMFileName);
% Read in continent land areas for plot overlay using Mapping Toolbox function, shaperead.
landAreas = shaperead('landareas.shp','UseGeoCoords',true);
% Plotting Earth's Magnetic Field Overlaid on Earth Maps
% Load plot formatting data for each of the magnetic parameters.
plotWMM = load('astPlotWMM.mat');
hX = figure;
set(hX,'Position',[0 0 827 620],'Color','white')
astPlotWMMContours( WMMResults, plotWMM, 1, landAreas, geod_lat, geod_lon, decimal_year, model_epoch)
hDIP = figure;
set(hDIP,'Position',[0 0 827 620],'Color','white')
astPlotWMMContours( WMMResults, plotWMM, 6, landAreas, geod_lat, geod_lon, decimal_year, model_epoch)
hF = figure;
set(hF,'Position',[0 0 827 620],'Color','white')
astPlotWMMContours( WMMResults, plotWMM, 7, landAreas, geod_lat, geod_lon, decimal_year, model_epoch)
.
