I’m not sure why your wavefront data is all zeros.  Perhaps you can share more of your Matlab code?

When I use the same sample file, and run the following code, it seems to work okay:

% Create new wavefront analysis
New_WF_map = TheSystem.Analyses.New_WavefrontMap(); 

% Run Analysis & gets results
New_WF_map.ApplyAndWaitForCompletion();
WF_Map_Data = New_WF_map.Results.DataGrids(1).Values.double;

% plot results
figure, imagesc(WF_Map_Data), axis image 
colorbar

Hi Jeff, thanks for your reply. Here are the parts of the code which take the wavefront analysis:

Creating the analyses and wavefront analysis are done through these lines

% Create analysis
    TheAnalyses = TheSystem.Analyses;
    New_WF_map = TheAnalyses.New_Analysis(ZOSAPI.Analysis.AnalysisIDM.WavefrontMap);

After creating the analysis the field and wavelength number are set to 1

    %Wavefront settings
    New_WF_Settings = New_WF_map.GetSettings();
    New_WF_Settings.Wavelength.SetWavelengthNumber(1);
    New_WF_Settings.Field.SetFieldNumber(1);

I then apply the analysis and grab the results

    % Run Analysis & gets results
    New_WF_map.ApplyAndWaitForCompletion();
    New_WF_map_results = New_WF_map.GetResults();

Then the results are plotted with the following code

    %Plot wavefront
    WF_Map_Data = New_WF_map_results.DataGrids(1).Values.double;
    figure
    imagesc(WF_Map_Data);
    colorbar;
    title('Wavefront Map 2')
    colormap(jet);
    xlabel('Px');
    ylabel('Py');

Looks fine to me.  When I run your code with the “Double Gauss 5 degree field.zos” lens I get the expected result: