Index: sm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_a.asv =================================================================== --- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_a.asv (revision 7891) +++ (revision ) @@ -1,389 +1,0 @@ -%----------------------------------------------------------------- -% Author: James Fishbaugh, Laura Perichon & Frank Pattyn -% Date: 06/03/2008 -% Description: Reads experiment data and generates a plot -% of a region of one standard deviation from the -% mean for the different types of models used. -%---------------------------------------------------------------- - -%-------------------------------------------------------- -% THIS IS THE FILE IO SECTION WHERE THE DATA FILES ARE -% READ AND VALUES STORED -%-------------------------------------------------------- - -clear all; -close all; - -% The names of the folders (also prefix for filenames) -dirs = {'aas1' 'aas2' 'ahu1' 'ahu2' 'bds1' 'cma1' 'cma2' 'dpo1' ... - 'fpa1' 'fpa2' 'fsa1' 'jvj1' 'lpe1' 'mbr1' 'mmr1' 'mtk1' ... - 'oga1' 'oso1' 'rhi1' 'rhi2' 'rhi3' 'rhi4' 'rhi5' 'spr1' ... - 'ssu1' 'tpa1' 'yko1'}; -% dirs = {'fpa1'}; -% The entries in this vector correspond to the vector of directory names -% For example: aas2 is full stokes, ahu1 is not, ahu2 is not, etc -% NOTE: THIS NEEDS TO BE CHECKED FOR ACCURACY -fullStokesMask = [1 1 0 0 0 1 0 0 ... - 0 1 0 0 0 0 1 0 ... - 1 0 1 0 1 0 0 1 ... - 1 0 1]; -% fullStokesMask = [0]; -% The location of the data (this is the path to the directory where all the -% sub dirs lie) -% NOTE: CHANGE THIS TO THE LOCATION OF YOUR DATA -location = 'C:\Program Files\MATLAB\R2006a\work\ismip_all\ismip_mat\'; - -% The experiment we're interested in -exp = 'a'; - -% The file we're intersted in -fileNames = {[exp '005'], [exp '010'], [exp '020'], ... - [exp '040'], [exp '080'], [exp '160']}; - -% The desired vector for normalized x. Used for interpolation so every -% data set has the same number of points. -newX = linspace(0.0,1.0,250); -newY = 1/4*ones(250,1)'; - -% Will hold all the full stokes surface velocities -fsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 1))); -% Will hold all the non-full stokes surface velocities -nfsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 0))); - -% Keeps track of the index we are currently on (this is necessary because -% some submittors didn't do all experiments so we can't rely on the loop -% counter to keep track of where we are at) -fsIndex = 1; -nfsIndex = 1; - -% Keeps track of the number of submittors for each domain length (this will -% help us later on in allocating the correct size of data structures) -fsSizes = zeros(1, length(fileNames)); -nfsSizes = zeros(1, length(fileNames)); - -% Loop over all folders -for i=1:length(dirs) - - % Loop over all files - for j=1:length(fileNames) - - % Create the full file path from the base path and current directory - fullFilePath = [location dirs(i) fileNames(j)]; - - % Open the current file for reading - load(cell2mat(fullFilePath)); - - % If that file exists - if (isnan(B(1,1))==0) - - % Store the values - x = B(:,1); - y = B(:,2); - xvel = B(:,3); - yvel = B(:,4); - zvel = B(:,5); - - % Make sure the data is normalized - x = x./max(x); - - % Make sure the data is normalized - y = y./max(y); - - x = inpaint_nans(x, 4); - y = inpaint_nans(y, 4); - xvel = inpaint_nans(xvel, 4); - yvel = inpaint_nans(yvel, 4); - zvel = inpaint_nans(zvel, 4); - - % Calculates the magnitude of the surface velocity - - surfVel = sqrt(xvel.^2 + yvel.^2); - - % If the current data is from a full stokes - if (fullStokesMask(i) == 1) - [i j fsIndex] - % Store the data in the full stokes structure - fsVels{j, fsIndex} = griddata(x, y, surfVel, newX, newY); - fsVels{j, fsIndex} = inpaint_nans(fsVels{j,fsIndex}, 4); - MaxVelfs(i,j)=max(fsVels{j,fsIndex}); - MeanVelfs(i,j)=mean(fsVels{j,fsIndex}); - - % Keep track of the number of entries in the array - if (length(fsVels{j, fsIndex}) > 0) - - % Increment the number who performed this experiment - fsSizes(j) = fsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the full stokes index - fsIndex = fsIndex + 1; - - end - - % Else its non-full stokes - else - [i j nfsIndex] - % Store the data in the non-full stokes structure - nfsVels{j, nfsIndex} = griddata(x, y, surfVel, newX, newY); - nfsVels{j, nfsIndex} = inpaint_nans(nfsVels{j, nfsIndex}, 4); - MaxVelnfs(i,j)=max(nfsVels{j,nfsIndex}); - MeanVelnfs(i,j)=mean(nfsVels{j,nfsIndex}); - - % Keep track of the number of entries in the array - if (length(nfsVels{j, nfsIndex}) > 0) - - % Increment the number who performed this experiment - nfsSizes(j) = nfsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the non-full stokes index - nfsIndex = nfsIndex + 1; - - end - - end - - % We are done with this file, we can close it -% fclose('all'); - - end - - end - -end - -%-------------------------------------------------------- -% THIS IS THE PROCESSING AND DISPLAY SECTION WHERE THE -% MEAN AND STANDARD DEVIATION ARE CALCULATED AND PLOTTED -%-------------------------------------------------------- - -% The titles for the plots -plotTitles = {'5km' '10km' '20km' '40km' '80km' '160km'}; -% The label for the x-axis -xLabel = 'Normalized x'; -% The label for the y-axis -yLabel = 'Velocity (m/yr)'; - -% The font to use for the axis labels -axisLabelFont = 'normal'; -% The font size to use for the axis labels -axisLabelSize = 7; -% The font weight to use for the axis labels -axisFontWeight = 'normal'; - -% The font size to use for tick mark labels -axisTickLabelSize = 7; -% The font weight to use for tick mark labels -axisTickFontWeight = 'normal'; - -% The font size to use for the title -titleFontSize = 10; -% The font weight to use for the title -titleFontWeight = 'bold'; - -% Loop over each domain length -for i=1:length(fileNames) - - % Allocate a matrix for full stokes surface velocities - fsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - fsIndex = 1; - - % Allocate a matrix for non-full stokes surface velocities - nfsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - nfsIndex = 1; - - % Loop over full stokes - for j=1:length(fullStokesMask(fullStokesMask == 1)) - - % If there are data at this entry - if (length(fsVels{i,j}) ~= 0) - - % Add it to the matrix - fsCurDomainVels(:,fsIndex) = fsVels{i,j}; - - % Increment the index - fsIndex = fsIndex + 1; - - end - - end - - % Loop over non-full stokes - for j=1:length(fullStokesMask(fullStokesMask == 0)) - - % If there are data at this entry - if (length(nfsVels{i,j}) ~= 0) - - % Add it to the matrix - nfsCurDomainVels(:,nfsIndex) = nfsVels{i,j}; - - % Increment the index - nfsIndex = nfsIndex + 1; - end - - end - - % Calcuate the mean for each domain length (full stokes) - fsMeanVel = mean(fsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (full stokes) - fsStdVel = std(fsCurDomainVels, 0, 2)'; - - % Calcuate the mean for each domain length (non-full stokes) - nfsMeanVel = mean(nfsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (non-full stokes) - nfsStdVel = std(nfsCurDomainVels, 0, 2)'; - - % Tell MATLAB which subplot we are currently on - subplot(2, 3, i); - - % The fill color (rgb) for full stokes - fsFillColor = [0.254 0.625 0.660]; - % The fill color (rgb) for non-full stokes - nfsFillColor = [0.422 0.688 0.367]; - - % Plot full stokes - [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ... - fsFillColor, fsFillColor, 1, 0.5); - - % Plot non-full stokes - [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ... - nfsFillColor, nfsFillColor, 1, 0.35); - - % Turn the grid on - grid on; - - % Set title properties - th = title(plotTitles(i)); - set(th, 'FontSize', titleFontSize); - set(th, 'FontWeight', titleFontWeight); - - % Set x label properties - xlabel(xLabel); - xh = get(gca, 'xlabel'); - set(xh, 'FontName', axisLabelFont); - set(xh, 'FontSize', axisLabelSize); - set(xh, 'FontWeight', axisFontWeight); - - % Set y label properties - ylabel('Velocity (m/yr)'); - yh = get(gca, 'ylabel'); - set(yh, 'FontName', axisLabelFont); - set(yh, 'FontSize', axisLabelSize); - set(yh, 'FontWeight', axisFontWeight); - - % Set tick mark properties - set(gca, 'FontSize', axisTickLabelSize); - set(gca, 'FontWeight', axisTickFontWeight); - - % This turns the box around the axis on - set(gca, 'Box', 'on' ); - - % These lines plot the mean - hold on; - plot(newX, fsMeanVel, 'k'); - set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2); - plot(newX, nfsMeanVel, 'b'); - set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2); - -end - -% Add the legend to the final subplot (in upper left corner) -legend('FS', 'NFS', 'FS Mean', 'NFS Mean', 'Location','West'); - -set(gcf, 'paperpositionmode','manual','paperorientation', ... - 'portrait','papertype','A4','paperunits', 'centimeters', ... - 'paperposition',[1 4 21 13]); -saveas(gcf, 'expa.eps', 'psc2'); - -% Analysis of maximum velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MaxVelfs) - for j=1:size(MaxVelfs,2) - if MaxVelfs(i,j)==0 - MaxVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - StdDevfs(j)=std(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MaxVelfs(:,j))); -end -for i=1:length(MaxVelnfs) - for j=1:size(MaxVelnfs,2) - if MaxVelnfs(i,j)==0 - MaxVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelnfs,2) - MeanMaxnfs(j)=mean(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - StdDevnfs(j)=std(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MaxVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n'); - -% Analysis of mean velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MaxVelfs) - for j=1:size(MaxVelfs,2) - if MaxVelfs(i,j)==0 - MaxVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - StdDevfs(j)=std(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MaxVelfs(:,j))); -end -for i=1:length(MaxVelnfs) - for j=1:size(MaxVelnfs,2) - if MaxVelnfs(i,j)==0 - MaxVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelnfs,2) - MeanMaxnfs(j)=mean(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - StdDevnfs(j)=std(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MaxVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n'); - Index: sm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_b.asv =================================================================== --- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_b.asv (revision 7891) +++ (revision ) @@ -1,294 +1,0 @@ -%----------------------------------------------------------------- -% Author: James Fishbaugh, Laura Perichon & Frank Pattyn -% Date: 06/03/2008 -% Description: Reads experiment data and generates a plot -% of a region of one standard deviation from the -% mean for the different types of models used. -%---------------------------------------------------------------- - -%-------------------------------------------------------- -% THIS IS THE FILE IO SECTION WHERE THE DATA FILES ARE -% READ AND VALUES STORED -%-------------------------------------------------------- - -clear all; -close all; - -% The names of the folders (also prefix for filenames) -dirs = {'aas1' 'ahu1' 'ahu2' 'bds1' 'cma1' 'cma2' 'dpo1' ... - 'fpa1' 'fpa2' 'fsa1' 'jvj1' 'lpe1' 'mbr1' 'mmr1' 'mtk1' ... - 'oga1' 'rhi1' 'rhi2' 'rhi3' 'rhi4' 'rhi5' 'spr1' ... - 'ssu1' 'tpa1' 'yko1'}; -dirs = {'fsa1'}; -% The entries in this vector correspond to the vector of directory names -% For example: aas2 is full stokes, ahu1 is not, ahu2 is not, etc -% NOTE: THIS NEEDS TO BE CHECKED FOR ACCURACY -fullStokesMask = [1 0 0 0 1 0 0 ... - 0 1 0 0 0 0 1 0 ... - 1 1 0 1 0 0 1 ... - 1 0 1]; -fullStokesMask = [0]; -% The location of the data (this is the path to the directory where all the -% sub dirs lie) -% NOTE: CHANGE THIS TO THE LOCATION OF YOUR DATA -location = 'C:\Program Files\MATLAB\R2006a\work\ismip_all\ismip_mat\'; - -% The experiment we're interested in -exp = 'b'; - -% The file we're intersted in -fileNames = {[exp '005'], [exp '010'], [exp '020'], ... - [exp '040'], [exp '080'], [exp '160']}; - -% The desired vector for normalized x. Used for interpolation so every -% data set has the same number of points. -newX = linspace(0.0,1.0,250); - -% Will hold all the full stokes surface velocities -fsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 1))); -% Will hold all the non-full stokes surface velocities -nfsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 0))); - -% Keeps track of the index we are currently on (this is necessary because -% some submittors didn't do all experiments so we can't rely on the loop -% counter to keep track of where we are at) -fsIndex = 1; -nfsIndex = 1; - -% Keeps track of the number of submittors for each domain length (this will -% help us later on in allocating the correct size of data structures) -fsSizes = zeros(1, length(fileNames)); -nfsSizes = zeros(1, length(fileNames)); - -% Loop over all folders -for i=1:length(dirs) - - % Loop over all files - for j=1:length(fileNames) - - % Create the full file path from the base path and current directory - fullFilePath = [location dirs(i) fileNames(j)]; - - % Open the current file for reading - load(cell2mat(fullFilePath)); - - % If that file exists - if (isnan(B(1,1))==0) - - % Store the values - x = B(:,1); - xvel = B(:,2); - zvel = B(:,3); - - % Make sure the data is normalized - x = x./max(x); - - x = inpaint_nans(x, 4); - xvel = inpaint_nans(xvel, 4); - zvel = inpaint_nans(zvel, 4); - - % Calculates the magnitude of the surface velocity - - surfVel = sqrt(xvel.^2); - - % If the current data is from a full stokes - if (fullStokesMask(i) == 1) - [i j fsIndex] - % Store the data in the full stokes structure - fsVels{j, fsIndex} = interp1(x, surfVel, newX); - fsVels{j, fsIndex} = inpaint_nans(fsVels{j,fsIndex}, 4); - - % Keep track of the number of entries in the array - if (length(fsVels{j, fsIndex}) > 0) - - % Increment the number who performed this experiment - fsSizes(j) = fsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the full stokes index - fsIndex = fsIndex + 1; - - end - - % Else its non-full stokes - else - [i j nfsIndex] - % Store the data in the non-full stokes structure - nfsVels{j, nfsIndex} = interp1(x, surfVel, newX); - nfsVels{j, nfsIndex} = inpaint_nans(nfsVels{j, nfsIndex}, 4); - - % Keep track of the number of entries in the array - if (length(nfsVels{j, nfsIndex}) > 0) - - % Increment the number who performed this experiment - nfsSizes(j) = nfsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the non-full stokes index - nfsIndex = nfsIndex + 1; - - end - - end - - % We are done with this file, we can close it -% fclose('all'); - - end - - end - -end - -%-------------------------------------------------------- -% THIS IS THE PROCESSING AND DISPLAY SECTION WHERE THE -% MEAN AND STANDARD DEVIATION ARE CALCULATED AND PLOTTED -%-------------------------------------------------------- - -% The titles for the plots -plotTitles = {'5km' '10km' '20km' '40km' '80km' '160km'}; -% The label for the x-axis -xLabel = 'Normalized x'; -% The label for the y-axis -yLabel = 'Velocity (m/yr)'; - -% The font to use for the axis labels -axisLabelFont = 'normal'; -% The font size to use for the axis labels -axisLabelSize = 14; -% The font weight to use for the axis labels -axisFontWeight = 'normal'; - -% The font size to use for tick mark labels -axisTickLabelSize = 11; -% The font weight to use for tick mark labels -axisTickFontWeight = 'normal'; - -% The font size to use for the title -titleFontSize = 12; -% The font weight to use for the title -titleFontWeight = 'bold'; - -% Loop over each domain length -for i=1:length(fileNames) - - % Allocate a matrix for full stokes surface velocities - fsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - fsIndex = 1; - - % Allocate a matrix for non-full stokes surface velocities - nfsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - nfsIndex = 1; - - % Loop over full stokes - for j=1:length(fullStokesMask(fullStokesMask == 1)) - - % If there are data at this entry - if (length(fsVels{i,j}) ~= 0) - - % Add it to the matrix - fsCurDomainVels(:,fsIndex) = fsVels{i,j}; - - % Increment the index - fsIndex = fsIndex + 1; - - end - - end - - % Loop over non-full stokes - for j=1:length(fullStokesMask(fullStokesMask == 0)) - - % If there are data at this entry - if (length(nfsVels{i,j}) ~= 0) - - % Add it to the matrix - nfsCurDomainVels(:,nfsIndex) = nfsVels{i,j}; - - % Increment the index - nfsIndex = nfsIndex + 1; - end - - end - - % Calcuate the mean for each domain length (full stokes) - fsMeanVel = mean(fsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (full stokes) - fsStdVel = std(fsCurDomainVels, 0, 2)'; - - % Calcuate the mean for each domain length (non-full stokes) - nfsMeanVel = mean(nfsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (non-full stokes) - nfsStdVel = std(nfsCurDomainVels, 0, 2)'; - - % Tell MATLAB which subplot we are currently on - subplot(2, 3, i); - - % The fill color (rgb) for full stokes - fsFillColor = [0.254 0.625 0.660]; - % The fill color (rgb) for non-full stokes - nfsFillColor = [0.422 0.688 0.367]; - - % Plot full stokes - [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ... - fsFillColor, fsFillColor, 1, 0.5); - - % Plot non-full stokes - [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ... - nfsFillColor, nfsFillColor, 1, 0.35); - - % Turn the grid on - grid on; - - % Set title properties - th = title(plotTitles(i)); - set(th, 'FontSize', titleFontSize); - set(th, 'FontWeight', titleFontWeight); - - % Set x label properties - xlabel(xLabel); - xh = get(gca, 'xlabel'); - set(xh, 'FontName', axisLabelFont); - set(xh, 'FontSize', axisLabelSize); - set(xh, 'FontWeight', axisFontWeight); - - % Set y label properties - ylabel('Velocity (m/yr)'); - yh = get(gca, 'ylabel'); - set(yh, 'FontName', axisLabelFont); - set(yh, 'FontSize', axisLabelSize); - set(yh, 'FontWeight', axisFontWeight); - - % Set tick mark properties - set(gca, 'FontSize', axisTickLabelSize); - set(gca, 'FontWeight', axisTickFontWeight); - - % This turns the box around the axis on - set(gca, 'Box', 'on' ); - - % These lines plot the mean - hold on; - plot(newX, fsMeanVel, 'k'); - set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2); - plot(newX, nfsMeanVel, 'b'); - set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2); - -end - -% Add the legend to the final subplot (in upper left corner) -legend('FS', 'NFS', 'FS Mean', 'NFS Mean', 2); - -set(gcf, 'paperpositionmode','manual','paperorientation', ... - 'landscape','papertype','A4','paperunits', 'centimeters', ... - 'paperposition',[1 4 25 15]); -saveas(gcf, 'expd.eps', 'psc2'); Index: sm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_d.asv =================================================================== --- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_d.asv (revision 7891) +++ (revision ) @@ -1,381 +1,0 @@ -%----------------------------------------------------------------- -% Author: James Fishbaugh, Laura Perichon & Frank Pattyn -% Date: 06/03/2008 -% Description: Reads experiment data and generates a plot -% of a region of one standard deviation from the -% mean for the different types of models used. -%---------------------------------------------------------------- - -%-------------------------------------------------------- -% THIS IS THE FILE IO SECTION WHERE THE DATA FILES ARE -% READ AND VALUES STORED -%-------------------------------------------------------- - -clear all; -close all; - -% The names of the folders (also prefix for filenames) -dirs = {'aas1' 'aas2' 'ahu1' 'ahu2' 'bds1' 'cma1' 'cma2' 'dpo1' ... - 'fpa1' 'fpa2' 'fsa1' 'lpe1' 'mbr1' 'mmr1' 'mtk1' ... - 'oga1' 'oso1' 'rhi1' 'rhi2' 'rhi3' 'rhi4' 'rhi5' 'spr1' ... - 'ssu1' 'tpa1' 'yko1'}; -dirs = {'dpo1' 'lpe1' 'rhi4' 'rhi5'}; -% The entries in this vector correspond to the vector of directory names -% For example: aas2 is full stokes, ahu1 is not, ahu2 is not, etc -% NOTE: THIS NEEDS TO BE CHECKED FOR ACCURACY -fullStokesMask = [1 1 0 0 0 1 0 0 ... - 0 1 0 0 0 1 0 ... - 1 0 1 0 1 0 0 1 ... - 1 0 1]; -fullStokesMask = [0 ]; -% The location of the data (this is the path to the directory where all the -% sub dirs lie) -% NOTE: CHANGE THIS TO THE LOCATION OF YOUR DATA -location = 'C:\Program Files\MATLAB\R2006a\work\ismip_all\ismip_mat\'; - -% The experiment we're interested in -exp = 'd'; - -% The file we're intersted in -fileNames = {[exp '005'], [exp '010'], [exp '020'], ... - [exp '040'], [exp '080'], [exp '160']}; - -% The desired vector for normalized x. Used for interpolation so every -% data set has the same number of points. -newX = linspace(0.0,1.0,250); - -% Will hold all the full stokes surface velocities -fsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 1))); -% Will hold all the non-full stokes surface velocities -nfsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 0))); - -% Keeps track of the index we are currently on (this is necessary because -% some submittors didn't do all experiments so we can't rely on the loop -% counter to keep track of where we are at) -fsIndex = 1; -nfsIndex = 1; - -% Keeps track of the number of submittors for each domain length (this will -% help us later on in allocating the correct size of data structures) -fsSizes = zeros(1, length(fileNames)); -nfsSizes = zeros(1, length(fileNames)); - -% Loop over all folders -for i=1:length(dirs) - - % Loop over all files - for j=1:length(fileNames) - - % Create the full file path from the base path and current directory - fullFilePath = [location dirs(i) fileNames(j)]; - - % Open the current file for reading - load(cell2mat(fullFilePath)); - - % If that file exists - if (isnan(B(1,1))==0) - - % Store the values - x = B(:,1); - xvel = B(:,2); - zvel = B(:,3); - - % Make sure the data is normalized - x = x./max(x); - - x = inpaint_nans(x, 4); - xvel = inpaint_nans(xvel, 4); - zvel = inpaint_nans(zvel, 4); - - % Calculates the magnitude of the surface velocity - - surfVel = sqrt(xvel.^2); - - % If the current data is from a full stokes - if (fullStokesMask(i) == 1) - [i j fsIndex] - % Store the data in the full stokes structure - fsVels{j, fsIndex} = interp1(x, surfVel, newX); - fsVels{j, fsIndex} = inpaint_nans(fsVels{j,fsIndex}, 4); - MaxVelfs(i,j)=max(fsVels{j,fsIndex}); - MeanVelfs(i,j)=mean(fsVels{j,fsIndex}); - - % Keep track of the number of entries in the array - if (length(fsVels{j, fsIndex}) > 0) - - % Increment the number who performed this experiment - fsSizes(j) = fsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the full stokes index - fsIndex = fsIndex + 1; - - end - - % Else its non-full stokes - else - [i j nfsIndex] - % Store the data in the non-full stokes structure - nfsVels{j, nfsIndex} = interp1(x, surfVel, newX); - nfsVels{j, nfsIndex} = inpaint_nans(nfsVels{j, nfsIndex}, 4); - MaxVelnfs(i,j)=max(nfsVels{j,nfsIndex}); - MeanVelnfs(i,j)=mean(nfsVels{j,nfsIndex}); - - % Keep track of the number of entries in the array - if (length(nfsVels{j, nfsIndex}) > 0) - - % Increment the number who performed this experiment - nfsSizes(j) = nfsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the non-full stokes index - nfsIndex = nfsIndex + 1; - - end - - end - - % We are done with this file, we can close it -% fclose('all'); - - end - - end - -end - -%-------------------------------------------------------- -% THIS IS THE PROCESSING AND DISPLAY SECTION WHERE THE -% MEAN AND STANDARD DEVIATION ARE CALCULATED AND PLOTTED -%-------------------------------------------------------- - -% The titles for the plots -plotTitles = {'5km' '10km' '20km' '40km' '80km' '160km'}; -% The label for the x-axis -xLabel = 'Normalized x'; -% The label for the y-axis -yLabel = 'Velocity (m a^{-1})'; - -% The font to use for the axis labels -axisLabelFont = 'normal'; -% The font size to use for the axis labels -axisLabelSize = 7; -% The font weight to use for the axis labels -axisFontWeight = 'normal'; - -% The font size to use for tick mark labels -axisTickLabelSize = 7; -% The font weight to use for tick mark labels -axisTickFontWeight = 'normal'; - -% The font size to use for the title -titleFontSize = 10; -% The font weight to use for the title -titleFontWeight = 'bold'; - -% Loop over each domain length -for i=1:length(fileNames) - - % Allocate a matrix for full stokes surface velocities - fsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - fsIndex = 1; - - % Allocate a matrix for non-full stokes surface velocities - nfsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - nfsIndex = 1; - - % Loop over full stokes - for j=1:length(fullStokesMask(fullStokesMask == 1)) - - % If there are data at this entry - if (length(fsVels{i,j}) ~= 0) - - % Add it to the matrix - fsCurDomainVels(:,fsIndex) = fsVels{i,j}; - - % Increment the index - fsIndex = fsIndex + 1; - - end - - end - - % Loop over non-full stokes - for j=1:length(fullStokesMask(fullStokesMask == 0)) - - % If there are data at this entry - if (length(nfsVels{i,j}) ~= 0) - - % Add it to the matrix - nfsCurDomainVels(:,nfsIndex) = nfsVels{i,j}; - - % Increment the index - nfsIndex = nfsIndex + 1; - end - - end - - % Calcuate the mean for each domain length (full stokes) - fsMeanVel = mean(fsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (full stokes) - fsStdVel = std(fsCurDomainVels, 0, 2)'; - - % Calcuate the mean for each domain length (non-full stokes) - nfsMeanVel = mean(nfsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (non-full stokes) - nfsStdVel = std(nfsCurDomainVels, 0, 2)'; - - % Tell MATLAB which subplot we are currently on - subplot(2, 3, i); - - % The fill color (rgb) for full stokes - fsFillColor = [0.254 0.625 0.660]; - % The fill color (rgb) for non-full stokes - nfsFillColor = [0.422 0.688 0.367]; - - % Plot full stokes - [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ... - fsFillColor, fsFillColor, 1, 0.5); - - % Plot non-full stokes - [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ... - nfsFillColor, nfsFillColor, 1, 0.35); - - % Turn the grid on - grid on; - - % Set title properties - th = title(plotTitles(i)); - set(th, 'FontSize', titleFontSize); - set(th, 'FontWeight', titleFontWeight); - - % Set x label properties - xlabel(xLabel); - xh = get(gca, 'xlabel'); - set(xh, 'FontName', axisLabelFont); - set(xh, 'FontSize', axisLabelSize); - set(xh, 'FontWeight', axisFontWeight); - - % Set y label properties - ylabel('Velocity (m a^{-1})'); - yh = get(gca, 'ylabel'); - set(yh, 'FontName', axisLabelFont); - set(yh, 'FontSize', axisLabelSize); - set(yh, 'FontWeight', axisFontWeight); - - % Set tick mark properties - set(gca, 'FontSize', axisTickLabelSize); - set(gca, 'FontWeight', axisTickFontWeight); - - % This turns the box around the axis on - set(gca, 'Box', 'on' ); - - % These lines plot the mean - hold on; - plot(newX, fsMeanVel, 'k'); - set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2); - plot(newX, nfsMeanVel, 'b'); - set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2); - -end - -% Add the legend to the final subplot (in upper left corner) -legend('FS', 'NFS', 'FS Mean', 'NFS Mean', 'Location','West'); - -set(gcf, 'paperpositionmode','manual','paperorientation', ... - 'portrait','papertype','A4','paperunits', 'centimeters', ... - 'paperposition',[1 4 21 13]); -saveas(gcf, 'expd.eps', 'psc2'); - -% Analysis of maximum velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MaxVelfs) - for j=1:size(MaxVelfs,2) - if MaxVelfs(i,j)==0 - MaxVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - StdDevfs(j)=std(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MaxVelfs(:,j))); -end -for i=1:length(MaxVelnfs) - for j=1:size(MaxVelnfs,2) - if MaxVelnfs(i,j)==0 - MaxVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelnfs,2) - MeanMaxnfs(j)=mean(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - StdDevnfs(j)=std(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MaxVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n'); - -% Analysis of mean surface velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MeanVelfs) - for j=1:size(MeanVelfs,2) - if MeanVelfs(i,j)==0 - MeanVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MeanVelfs(~isnan(MeanVelfs(:,j)),j)); - StdDevfs(j)=std(MeanVelfs(~isnan(MeanVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MeanVelfs(:,j))); -end -for i=1:length(MeanVelnfs) - for j=1:size(MeanVelnfs,2) - if MeanVelnfs(i,j)==0 - MeanVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MeanVelnfs,2) - MeanMaxnfs(j)=mean(MeanVelnfs(~isnan(MeanVelnfs(:,j)),j)); - StdDevnfs(j)=std(MeanVelnfs(~isnan(MeanVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MeanVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n'); - Index: sm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e.asv =================================================================== --- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e.asv (revision 7891) +++ (revision ) @@ -1,380 +1,0 @@ -%----------------------------------------------------------------- -% Author: James Fishbaugh, Laura Perichon & Frank Pattyn -% Date: 06/03/2008 -% Description: Reads experiment data and generates a plot -% of a region of one standard deviation from the -% mean for the different types of models used. -%---------------------------------------------------------------- - -%-------------------------------------------------------- -% THIS IS THE FILE IO SECTION WHERE THE DATA FILES ARE -% READ AND VALUES STORED -%-------------------------------------------------------- - -clear all; -close all; - -% The names of the folders (also prefix for filenames) -dirs = {'aas1' 'aas2' 'ahu1' 'ahu2' 'bds1' 'cma1' 'cma2' 'spr1' 'dpo1' ... - 'fpa1' 'fpa2' 'fsa1' 'lpe1' 'mbr1' 'mmr1' 'mtk1' ... - 'oga1' 'oso1' 'rhi1' 'rhi2' 'rhi3' 'rhi4' 'rhi5' ... - 'ssu1' 'tpa1'}; -dirs = {'jvj1'}; -% The entries in this vector correspond to the vector of directory names -% For example: aas2 is full stokes, ahu1 is not, ahu2 is not, etc -% NOTE: THIS NEEDS TO BE CHECKED FOR ACCURACY -fullStokesMask = [1 1 0 0 0 1 0 0 1 ... - 0 1 0 1 0 0 1 0 ... - 1 0 1 0 1 0 0 ... - 1 0]; -fullStokesMask = [0]; -% The location of the data (this is the path to the directory where all the -% sub dirs lie) -% NOTE: CHANGE THIS TO THE LOCATION OF YOUR DATA -location = 'C:\Program Files\MATLAB\R2006a\work\ismip_all\ismip_mat\'; - -% The experiment we're interested in -exp = 'e'; - -% The file we're intersted in -fileNames = {[exp '000'], [exp '001']}; - -% The desired vector for normalized x. Used for interpolation so every -% data set has the same number of points. -newX = linspace(0.0,1.0,250); - -% Will hold all the full stokes surface velocities -fsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 1))); -% Will hold all the non-full stokes surface velocities -nfsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 0))); - -% Keeps track of the index we are currently on (this is necessary because -% some submittors didn't do all experiments so we can't rely on the loop -% counter to keep track of where we are at) -fsIndex = 1; -nfsIndex = 1; - -% Keeps track of the number of submittors for each domain length (this will -% help us later on in allocating the correct size of data structures) -fsSizes = zeros(1, length(fileNames)); -nfsSizes = zeros(1, length(fileNames)); - -% Loop over all folders -for i=1:length(dirs) - - % Loop over all files - for j=1:length(fileNames) - - % Create the full file path from the base path and current directory - fullFilePath = [location dirs(i) fileNames(j)]; - - % Open the current file for reading - load(cell2mat(fullFilePath)); - - % If that file exists - if (isnan(B(1,1))==0) - - % Store the values - x = B(:,1); - xvel = B(:,2); - zvel = B(:,3); - - % Make sure the data is normalized - x = x./max(x); - - x = inpaint_nans(x, 4); - xvel = inpaint_nans(xvel, 4); - zvel = inpaint_nans(zvel, 4); - - % Calculates the magnitude of the surface velocity - - surfVel = sqrt(xvel.^2); - - % If the current data is from a full stokes - if (fullStokesMask(i) == 1) - [i j fsIndex] - % Store the data in the full stokes structure - fsVels{j, fsIndex} = interp1(x, surfVel, newX); - fsVels{j, fsIndex} = inpaint_nans(fsVels{j,fsIndex}, 4); - MaxVelfs(i,j)=max(fsVels{j,fsIndex}); - MeanVelfs(i,j)=mean(fsVels{j,fsIndex}); - - % Keep track of the number of entries in the array - if (length(fsVels{j, fsIndex}) > 0) - - % Increment the number who performed this experiment - fsSizes(j) = fsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the full stokes index - fsIndex = fsIndex + 1; - - end - - % Else its non-full stokes - else - [i j nfsIndex] - % Store the data in the non-full stokes structure - nfsVels{j, nfsIndex} = interp1(x, surfVel, newX); - nfsVels{j, nfsIndex} = inpaint_nans(nfsVels{j, nfsIndex}, 4); - MaxVelnfs(i,j)=max(nfsVels{j,nfsIndex}); - MeanVelnfs(i,j)=mean(nfsVels{j,nfsIndex}); - - % Keep track of the number of entries in the array - if (length(nfsVels{j, nfsIndex}) > 0) - - % Increment the number who performed this experiment - nfsSizes(j) = nfsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the non-full stokes index - nfsIndex = nfsIndex + 1; - - end - - end - - % We are done with this file, we can close it -% fclose('all'); - - end - - end - -end - -%-------------------------------------------------------- -% THIS IS THE PROCESSING AND DISPLAY SECTION WHERE THE -% MEAN AND STANDARD DEVIATION ARE CALCULATED AND PLOTTED -%-------------------------------------------------------- - -% The titles for the plots -plotTitles = {'No-Slip Bed' 'Slip Bed'}; -% The label for the x-axis -xLabel = 'Normalized x'; -% The label for the y-axis -yLabel = 'Velocity (m a^{-1})'; - -% The font to use for the axis labels -axisLabelFont = 'normal'; -% The font size to use for the axis labels -axisLabelSize = 7; -% The font weight to use for the axis labels -axisFontWeight = 'normal'; - -% The font size to use for tick mark labels -axisTickLabelSize = 7; -% The font weight to use for tick mark labels -axisTickFontWeight = 'normal'; - -% The font size to use for the title -titleFontSize = 10; -% The font weight to use for the title -titleFontWeight = 'bold'; - -% Loop over each domain length -for i=1:length(fileNames) - - % Allocate a matrix for full stokes surface velocities - fsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - fsIndex = 1; - - % Allocate a matrix for non-full stokes surface velocities - nfsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - nfsIndex = 1; - - % Loop over full stokes - for j=1:length(fullStokesMask(fullStokesMask == 1)) - - % If there are data at this entry - if (length(fsVels{i,j}) ~= 0) - - % Add it to the matrix - fsCurDomainVels(:,fsIndex) = fsVels{i,j}; - - % Increment the index - fsIndex = fsIndex + 1; - - end - - end - - % Loop over non-full stokes - for j=1:length(fullStokesMask(fullStokesMask == 0)) - - % If there are data at this entry - if (length(nfsVels{i,j}) ~= 0) - - % Add it to the matrix - nfsCurDomainVels(:,nfsIndex) = nfsVels{i,j}; - - % Increment the index - nfsIndex = nfsIndex + 1; - end - - end - - % Calcuate the mean for each domain length (full stokes) - fsMeanVel = mean(fsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (full stokes) - fsStdVel = std(fsCurDomainVels, 0, 2)'; - - % Calcuate the mean for each domain length (non-full stokes) - nfsMeanVel = mean(nfsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (non-full stokes) - nfsStdVel = std(nfsCurDomainVels, 0, 2)'; - - % Tell MATLAB which subplot we are currently on - subplot(2, 1, i); - - % The fill color (rgb) for full stokes - fsFillColor = [0.254 0.625 0.660]; - % The fill color (rgb) for non-full stokes - nfsFillColor = [0.422 0.688 0.367]; - - % Plot full stokes - [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ... - fsFillColor, fsFillColor, 1, 0.5); - - % Plot non-full stokes - [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ... - nfsFillColor, nfsFillColor, 1, 0.35); - - % Turn the grid on - grid on; - - % Set title properties - th = title(plotTitles(i)); - set(th, 'FontSize', titleFontSize); - set(th, 'FontWeight', titleFontWeight); - - % Set x label properties - xlabel(xLabel); - xh = get(gca, 'xlabel'); - set(xh, 'FontName', axisLabelFont); - set(xh, 'FontSize', axisLabelSize); - set(xh, 'FontWeight', axisFontWeight); - - % Set y label properties - ylabel('Velocity (m a^{-1})'); - yh = get(gca, 'ylabel'); - set(yh, 'FontName', axisLabelFont); - set(yh, 'FontSize', axisLabelSize); - set(yh, 'FontWeight', axisFontWeight); - - % Set tick mark properties - set(gca, 'FontSize', axisTickLabelSize); - set(gca, 'FontWeight', axisTickFontWeight); - - % This turns the box around the axis on - set(gca, 'Box', 'on' ); - - % These lines plot the mean - hold on; - plot(newX, fsMeanVel, 'k'); - set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2); - plot(newX, nfsMeanVel, 'b'); - set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2); - -end - -% Add the legend to the final subplot (in upper left corner) -legend('FS', 'NFS', 'FS Mean', 'NFS Mean', 2); - - -set(gcf, 'paperpositionmode','manual','paperorientation', ... - 'portrait','papertype','A4','paperunits', 'centimeters', ... - 'paperposition',[1 4 11 17]); -saveas(gcf, 'expe.eps', 'psc2'); - -% Analysis of maximum velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MaxVelfs) - for j=1:size(MaxVelfs,2) - if MaxVelfs(i,j)==0 - MaxVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - StdDevfs(j)=std(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MaxVelfs(:,j))); -end -for i=1:length(MaxVelnfs) - for j=1:size(MaxVelnfs,2) - if MaxVelnfs(i,j)==0 - MaxVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelnfs,2) - MeanMaxnfs(j)=mean(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - StdDevnfs(j)=std(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MaxVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n'); - -% Analysis of mean surface velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MeanVelfs) - for j=1:size(MeanVelfs,2) - if MeanVelfs(i,j)==0 - MeanVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MeanVelfs(~isnan(MeanVelfs(:,j)),j)); - StdDevfs(j)=std(MeanVelfs(~isnan(MeanVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MeanVelfs(:,j))); -end -for i=1:length(MeanVelnfs) - for j=1:size(MeanVelnfs,2) - if MeanVelnfs(i,j)==0 - MeanVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MeanVelnfs,2) - MeanMaxnfs(j)=mean(MeanVelnfs(~isnan(MeanVelnfs(:,j)),j)); - StdDevnfs(j)=std(MeanVelnfs(~isnan(MeanVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MeanVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n'); Index: sm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e_txz.asv =================================================================== --- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e_txz.asv (revision 7891) +++ (revision ) @@ -1,298 +1,0 @@ -%----------------------------------------------------------------- -% Author: James Fishbaugh, Laura Perichon & Frank Pattyn -% Date: 06/03/2008 -% Description: Reads experiment data and generates a plot -% of a region of one standard deviation from the -% mean for the different types of models used. -%---------------------------------------------------------------- - -%-------------------------------------------------------- -% THIS IS THE FILE IO SECTION WHERE THE DATA FILES ARE -% READ AND VALUES STORED -%-------------------------------------------------------- - -clear all; -close all; - -% The names of the folders (also prefix for filenames) -dirs = {'aas1' 'aas2' 'ahu1' 'ahu2' 'bds1' 'cma1' 'cma2' 'dpo1' ... - 'fpa1' 'fpa2' 'fsa1' 'jvj1' 'lpe1' 'mbr1' 'mmr1' 'mtk1' ... - 'oga1' 'oso1' 'rhi1' 'rhi2' 'rhi3' 'rhi4' 'rhi5' 'spr1' ... - 'ssu1' 'tpa1' 'yko1'}; -% The entries in this vector correspond to the vector of directory names -% For example: aas2 is full stokes, ahu1 is not, ahu2 is not, etc -% NOTE: THIS NEEDS TO BE CHECKED FOR ACCURACY -fullStokesMask = [1 1 0 0 0 1 0 0 ... - 0 1 0 0 0 0 1 0 ... - 1 0 1 0 1 0 0 1 ... - 1 0 1]; -% The location of the data (this is the path to the directory where all the -% sub dirs lie) -% NOTE: CHANGE THIS TO THE LOCATION OF YOUR DATA -location = 'C:\Program Files\MATLAB\R2006a\work\ismip_all\ismip_mat\'; - -% The experiment we're interested in -exp = 'e'; - -% The file we're intersted in -fileNames = {[exp '001'], [exp '000']}; - -% The desired vector for normalized x. Used for interpolation so every -% data set has the same number of points. -newX = linspace(0.0,1.0,250); - -% Will hold all the full stokes surface velocities -fsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 1))); -% Will hold all the non-full stokes surface velocities -nfsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 0))); - -% Keeps track of the index we are currently on (this is necessary because -% some submittors didn't do all experiments so we can't rely on the loop -% counter to keep track of where we are at) -fsIndex = 1; -nfsIndex = 1; - -% Keeps track of the number of submittors for each domain length (this will -% help us later on in allocating the correct size of data structures) -fsSizes = zeros(1, length(fileNames)); -nfsSizes = zeros(1, length(fileNames)); - -% Loop over all folders -for i=1:length(dirs) - - % Loop over all files - for j=1:length(fileNames) - - % Create the full file path from the base path and current directory - fullFilePath = [location dirs(i) fileNames(j)]; - - % Open the current file for reading - load(cell2mat(fullFilePath)); - - % If that file exists - if (isnan(B(1,1))==0) - - % Store the values - x = B(:,1); - xvel = B(:,2); - zvel = B(:,3); - txz = B(:,4); - - % Make sure the data is normalized - x = x./max(x); - - x = inpaint_nans(x, 4); - xvel = inpaint_nans(xvel, 4); - zvel = inpaint_nans(zvel, 4); - txz = inpaint_nans(txz, 4); - - % Calculates the magnitude of the basal shear - - surfVel = txz; - - % If the current data is from a full stokes - if (fullStokesMask(i) == 1) - [i j fsIndex] - % Store the data in the full stokes structure - fsVels{j, fsIndex} = interp1(x, surfVel, newX); - fsVels{j, fsIndex} = inpaint_nans(fsVels{j,fsIndex}, 4); - - % Keep track of the number of entries in the array - if (length(fsVels{j, fsIndex}) > 0) - - % Increment the number who performed this experiment - fsSizes(j) = fsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the full stokes index - fsIndex = fsIndex + 1; - - end - - % Else its non-full stokes - else - [i j nfsIndex] - % Store the data in the non-full stokes structure - nfsVels{j, nfsIndex} = interp1(x, surfVel, newX); - nfsVels{j, nfsIndex} = inpaint_nans(nfsVels{j, nfsIndex}, 4); - - % Keep track of the number of entries in the array - if (length(nfsVels{j, nfsIndex}) > 0) - - % Increment the number who performed this experiment - nfsSizes(j) = nfsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the non-full stokes index - nfsIndex = nfsIndex + 1; - - end - - end - - % We are done with this file, we can close it -% fclose('all'); - - end - - end - -end - -%-------------------------------------------------------- -% THIS IS THE PROCESSING AND DISPLAY SECTION WHERE THE -% MEAN AND STANDARD DEVIATION ARE CALCULATED AND PLOTTED -%-------------------------------------------------------- - -% The titles for the plots -plotTitles = {'Slip Bed' 'No-Slip Bed'}; -% The label for the x-axis -xLabel = 'Normalized x'; -% The label for the y-axis -yLabel = 'Basal shear stress (kPa)'; - -% The font to use for the axis labels -axisLabelFont = 'normal'; -% The font size to use for the axis labels -axisLabelSize = 7; -% The font weight to use for the axis labels -axisFontWeight = 'normal'; - -% The font size to use for tick mark labels -axisTickLabelSize = 7; -% The font weight to use for tick mark labels -axisTickFontWeight = 'normal'; - -% The font size to use for the title -titleFontSize = 10; -% The font weight to use for the title -titleFontWeight = 'bold'; - -% Loop over each domain length -for i=1:length(fileNames) - - % Allocate a matrix for full stokes surface velocities - fsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - fsIndex = 1; - - % Allocate a matrix for non-full stokes surface velocities - nfsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - nfsIndex = 1; - - % Loop over full stokes - for j=1:length(fullStokesMask(fullStokesMask == 1)) - - % If there are data at this entry - if (length(fsVels{i,j}) ~= 0) - - % Add it to the matrix - fsCurDomainVels(:,fsIndex) = fsVels{i,j}; - - % Increment the index - fsIndex = fsIndex + 1; - - end - - end - - % Loop over non-full stokes - for j=1:length(fullStokesMask(fullStokesMask == 0)) - - % If there are data at this entry - if (length(nfsVels{i,j}) ~= 0) - - % Add it to the matrix - nfsCurDomainVels(:,nfsIndex) = nfsVels{i,j}; - - % Increment the index - nfsIndex = nfsIndex + 1; - end - - end - - % Calcuate the mean for each domain length (full stokes) - fsMeanVel = mean(fsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (full stokes) - fsStdVel = std(fsCurDomainVels, 0, 2)'; - - % Calcuate the mean for each domain length (non-full stokes) - nfsMeanVel = mean(nfsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (non-full stokes) - nfsStdVel = std(nfsCurDomainVels, 0, 2)'; - - % Tell MATLAB which subplot we are currently on - if i==1 - subplot(2, 1, 2); - else - subplot(2, 1, 1); - end - - % The fill color (rgb) for full stokes - fsFillColor = [0.254 0.625 0.660]; - % The fill color (rgb) for non-full stokes - nfsFillColor = [0.422 0.688 0.367]; - - % Plot full stokes - [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ... - fsFillColor, fsFillColor, 1, 0.5); - - % Plot non-full stokes - [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ... - nfsFillColor, nfsFillColor, 1, 0.35); - - % Turn the grid on - grid on; - - % Set title properties - th = title(plotTitles(i)); - set(th, 'FontSize', titleFontSize); - set(th, 'FontWeight', titleFontWeight); - - % Set x label properties - xlabel(xLabel); - xh = get(gca, 'xlabel'); - set(xh, 'FontName', axisLabelFont); - set(xh, 'FontSize', axisLabelSize); - set(xh, 'FontWeight', axisFontWeight); - - % Set y label properties - ylabel('Velocity (m/yr)'); - yh = get(gca, 'ylabel'); - set(yh, 'FontName', axisLabelFont); - set(yh, 'FontSize', axisLabelSize); - set(yh, 'FontWeight', axisFontWeight); - - % Set tick mark properties - set(gca, 'FontSize', axisTickLabelSize); - set(gca, 'FontWeight', axisTickFontWeight); - - % This turns the box around the axis on - set(gca, 'Box', 'on' ); - - % These lines plot the mean - hold on; - plot(newX, fsCurDomainVels, 'k'); - set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2); - plot(newX, nfsCurDomainVels, 'b'); - set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2); - -end - -% Add the legend to the final subplot (in upper left corner) -legend('FS', 'NFS', 'FS Mean', 'NFS Mean', 2); - - -set(gcf, 'paperpositionmode','manual','paperorientation', ... - 'portrait','papertype','A4','paperunits', 'centimeters', ... - 'paperposition',[1 4 11 17]); -saveas(gcf, 'expe.eps', 'psc2'); Index: sm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_f.asv =================================================================== --- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_f.asv (revision 7891) +++ (revision ) @@ -1,414 +1,0 @@ -%----------------------------------------------------------------- -% Author: James Fishbaugh, Laura Perichon & Frank Pattyn -% Date: 06/03/2008 -% Description: Reads experiment data and generates a plot -% of a region of one standard deviation from the -% mean for the different types of models used. -%---------------------------------------------------------------- - -%-------------------------------------------------------- -% THIS IS THE FILE IO SECTION WHERE THE DATA FILES ARE -% READ AND VALUES STORED -%-------------------------------------------------------- - -clear all; -close all; - -% The names of the folders (also prefix for filenames) -dirs = {'aas1' 'aas2' 'ahu1' 'ahu2' 'bds1' 'cma1' 'cma2' 'dpo1' ... - 'fpa1' 'fpa2' 'fsa1' 'jvj1' 'lpe1' 'mbr1' 'mmr1' 'mtk1' ... - 'oga1' 'oso1' 'rhi1' 'rhi2' 'rhi3' 'rhi4' 'rhi5' 'spr1' ... - 'ssu1' 'tpa1' 'yko1'}; -% dirs = {'mtk1'}; -% The entries in this vector correspond to the vector of directory names -% For example: aas2 is full stokes, ahu1 is not, ahu2 is not, etc -% NOTE: THIS NEEDS TO BE CHECKED FOR ACCURACY -fullStokesMask = [1 1 0 0 0 1 0 0 ... - 0 1 0 0 0 0 1 0 ... - 1 0 1 0 1 0 0 1 ... - 1 0 1]; -% fullStokesMask = [0]; -% The location of the data (this is the path to the directory where all the -% sub dirs lie) -% NOTE: CHANGE THIS TO THE LOCATION OF YOUR DATA -location = 'C:\Program Files\MATLAB\R2006a\work\ismip_all\ismip_mat\'; - -% The experiment we're interested in -exp = 'f'; - -% The file we're intersted in -fileNames = {[exp '000'], [exp '001']}; - -% The desired vector for normalized x. Used for interpolation so every -% data set has the same number of points. -newX = linspace(-50.0,50.0,250); -newY = 0*ones(250,1)'; - -% Will hold all the full stokes surface velocities -fsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 1))); -% Will hold all the non-full stokes surface velocities -nfsVels = cell(length(fileNames), length(fullStokesMask(fullStokesMask == 0))); - -% Keeps track of the index we are currently on (this is necessary because -% some submittors didn't do all experiments so we can't rely on the loop -% counter to keep track of where we are at) -fsIndex = 1; -nfsIndex = 1; - -% Keeps track of the number of submittors for each domain length (this will -% help us later on in allocating the correct size of data structures) -fsSizes = zeros(1, length(fileNames)); -nfsSizes = zeros(1, length(fileNames)); - -% Loop over all folders -for i=1:length(dirs) - - % Loop over all files - for j=1:length(fileNames) - - % Create the full file path from the base path and current directory - fullFilePath = [location dirs(i) fileNames(j)]; - - % Open the current file for reading - load(cell2mat(fullFilePath)); - - % If that file exists - if (isnan(B(1,1))==0) - - % Store the values - x = B(:,1); - y = B(:,2); - xvel = B(:,3); - yvel = B(:,4); - zvel = B(:,5); - - % Make sure the data is normalized -% x = x./max(x); - - % Make sure the data is normalized -% y = y./max(y); - - x = inpaint_nans(x, 4); - y = inpaint_nans(y, 4); - xvel = inpaint_nans(xvel, 4); - yvel = inpaint_nans(yvel, 4); - zvel = inpaint_nans(zvel, 4); - - % Calculates the magnitude of the surface velocity - - surfVel = sqrt(yvel.^2 + zvel.^2); - - % If the current data is from a full stokes - if (fullStokesMask(i) == 1) - [i j fsIndex] - % Store the data in the full stokes structure - fsVels{j, fsIndex} = griddata(x, y, surfVel, newX, newY); - fsVels{j, fsIndex} = inpaint_nans(fsVels{j,fsIndex}, 4); - MaxVelfs(i,j)=max(fsVels{j,fsIndex}); - MeanVelfs(i,j)=mean(fsVels{j,fsIndex}); - - % Keep track of the number of entries in the array - if (length(fsVels{j, fsIndex}) > 0) - - % Increment the number who performed this experiment - fsSizes(j) = fsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the full stokes index - fsIndex = fsIndex + 1; - - end - - % Else its non-full stokes - else - [i j nfsIndex] - % Store the data in the non-full stokes structure - nfsVels{j, nfsIndex} = griddata(x, y, surfVel, newX, newY); - nfsVels{j, nfsIndex} = inpaint_nans(nfsVels{j, nfsIndex}, 4); - MaxVelnfs(i,j)=max(nfsVels{j,nfsIndex}); - MeanVelnfs(i,j)=mean(nfsVels{j,nfsIndex}); - - % Keep track of the number of entries in the array - if (length(nfsVels{j, nfsIndex}) > 0) - - % Increment the number who performed this experiment - nfsSizes(j) = nfsSizes(j) + 1; - - end - - % If we are on the last file we can increment the index - if (j == (length(fileNames))) - - % Increment the non-full stokes index - nfsIndex = nfsIndex + 1; - - end - - end - - % We are done with this file, we can close it -% fclose('all'); - - end - - end - -end - -% Analytical solution by Gudmundsson - -% Create the full file path from the base path and current directory -fullFilePath = [location 'ghg1' fileNames(1)]; - -% Open the current file for reading -load(cell2mat(fullFilePath)); - -x = B(:,1); -y = B(:,2); -xvel = B(:,3); -yvel = B(:,4); -zvel = B(:,5); - -x = inpaint_nans(x, 4); -y = inpaint_nans(y, 4); -xvel = inpaint_nans(xvel, 4); -yvel = inpaint_nans(yvel, 4); -zvel = inpaint_nans(zvel, 4); - -anVels{1, 1} = griddata(x, y, surfVel, newX, newY); - -% Calculates the magnitude of the surface velocity - -surfVel = sqrt(yvel.^2 + zvel.^2); - -%-------------------------------------------------------- -% THIS IS THE PROCESSING AND DISPLAY SECTION WHERE THE -% MEAN AND STANDARD DEVIATION ARE CALCULATED AND PLOTTED -%-------------------------------------------------------- - -% The titles for the plots -plotTitles = {'No-Slip Bed' 'Slip Bed'}; -% The label for the x-axis -xLabel = 'Distance for center (km)'; -% The label for the y-axis -yLabel = 'Velocity (m a^{-1})'; - -% The font to use for the axis labels -axisLabelFont = 'normal'; -% The font size to use for the axis labels -axisLabelSize = 7; -% The font weight to use for the axis labels -axisFontWeight = 'normal'; - -% The font size to use for tick mark labels -axisTickLabelSize = 7; -% The font weight to use for tick mark labels -axisTickFontWeight = 'normal'; - -% The font size to use for the title -titleFontSize = 10; -% The font weight to use for the title -titleFontWeight = 'bold'; - -% Loop over each domain length -for i=1:length(fileNames) - - % Allocate a matrix for full stokes surface velocities - fsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - fsIndex = 1; - - % Allocate a matrix for non-full stokes surface velocities - nfsCurDomainVels = zeros(250, fsSizes(i)); - % Again we need a seperate index - nfsIndex = 1; - - % Loop over full stokes - for j=1:length(fullStokesMask(fullStokesMask == 1)) - - % If there are data at this entry - if (length(fsVels{i,j}) ~= 0) - - % Add it to the matrix - fsCurDomainVels(:,fsIndex) = fsVels{i,j}; - - % Increment the index - fsIndex = fsIndex + 1; - - end - - end - - % Loop over non-full stokes - for j=1:length(fullStokesMask(fullStokesMask == 0)) - - % If there are data at this entry - if (length(nfsVels{i,j}) ~= 0) - - % Add it to the matrix - nfsCurDomainVels(:,nfsIndex) = nfsVels{i,j}; - - % Increment the index - nfsIndex = nfsIndex + 1; - end - - end - - % Calcuate the mean for each domain length (full stokes) - fsMeanVel = mean(fsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (full stokes) - fsStdVel = std(fsCurDomainVels, 0, 2)'; - - % Calcuate the mean for each domain length (non-full stokes) - nfsMeanVel = mean(nfsCurDomainVels, 2)'; - % Calculate the std deviation for each domain length (non-full stokes) - nfsStdVel = std(nfsCurDomainVels, 0, 2)'; - - % Tell MATLAB which subplot we are currently on - subplot(2, 1, i); - - % The fill color (rgb) for full stokes - fsFillColor = [0.254 0.625 0.660]; - % The fill color (rgb) for non-full stokes - nfsFillColor = [0.422 0.688 0.367]; - - % Plot full stokes - [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ... - fsFillColor, fsFillColor, 1, 0.5); - - % Plot non-full stokes - [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ... - nfsFillColor, nfsFillColor, 1, 0.35); - - % Turn the grid on - grid on; - - % Set title properties - th = title(plotTitles(i)); - set(th, 'FontSize', titleFontSize); - set(th, 'FontWeight', titleFontWeight); - - % Set x label properties - xlabel(xLabel); - xh = get(gca, 'xlabel'); - set(xh, 'FontName', axisLabelFont); - set(xh, 'FontSize', axisLabelSize); - set(xh, 'FontWeight', axisFontWeight); - - % Set y label properties - ylabel('Velocity (m a^{-1})'); - yh = get(gca, 'ylabel'); - set(yh, 'FontName', axisLabelFont); - set(yh, 'FontSize', axisLabelSize); - set(yh, 'FontWeight', axisFontWeight); - - % Set tick mark properties - set(gca, 'FontSize', axisTickLabelSize); - set(gca, 'FontWeight', axisTickFontWeight); - - % This turns the box around the axis on - set(gca, 'Box', 'on' ); - - % These lines plot the mean - hold on; - plot(newX, fsMeanVel, 'k'); - set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2); - plot(newX, nfsMeanVel, 'b'); - set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2); - -end - -% Add the legend to the final subplot (in upper left corner) -legend('FS', 'NFS', 'FS Mean', 'NFS Mean', 2); - - -set(gcf, 'paperpositionmode','manual','paperorientation', ... - 'portrait','papertype','A4','paperunits', 'centimeters', ... - 'paperposition',[1 4 11 17]); -saveas(gcf, 'expf.eps', 'psc2'); - -% Analysis of maximum velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MaxVelfs) - for j=1:size(MaxVelfs,2) - if MaxVelfs(i,j)==0 - MaxVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - StdDevfs(j)=std(MaxVelfs(~isnan(MaxVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MaxVelfs(:,j))); -end -for i=1:length(MaxVelnfs) - for j=1:size(MaxVelnfs,2) - if MaxVelnfs(i,j)==0 - MaxVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelnfs,2) - MeanMaxnfs(j)=mean(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - StdDevnfs(j)=std(MaxVelnfs(~isnan(MaxVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MaxVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n'); - -% Analysis of mean surface velocities (F. Pattyn) -% Full Stokes stored in MaxVelfs (else MaxVelnfs) -% cols = L, rows = participants, 0 = no result - -for i=1:length(MeanVelfs) - for j=1:size(MeanVelfs,2) - if MeanVelfs(i,j)==0 - MeanVelfs(i,j)=NaN; - end - end -end -for j=1:size(MaxVelfs,2) - MeanMaxfs(j)=mean(MeanVelfs(~isnan(MeanVelfs(:,j)),j)); - StdDevfs(j)=std(MeanVelfs(~isnan(MeanVelfs(:,j)),j)); - Numfs(j)=sum(~isnan(MeanVelfs(:,j))); -end -for i=1:length(MeanVelnfs) - for j=1:size(MeanVelnfs,2) - if MeanVelnfs(i,j)==0 - MeanVelnfs(i,j)=NaN; - end - end -end -for j=1:size(MeanVelnfs,2) - MeanMaxnfs(j)=mean(MeanVelnfs(~isnan(MeanVelnfs(:,j)),j)); - StdDevnfs(j)=std(MeanVelnfs(~isnan(MeanVelnfs(:,j)),j)); - Numnfs(j)=sum(~isnan(MeanVelnfs(:,j))); -end - -% print in format for LaTeX file -fprintf('NFS & '); -for j=1:size(MaxVelnfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxnfs(j), StdDevnfs(j), Numnfs(j)); -end -fprintf('\\\\ \n'); -fprintf('FS & '); -for j=1:size(MaxVelfs,2) - fprintf('%.2f & %.2f & %d & ', MeanMaxfs(j), StdDevfs(j), Numfs(j)); -end -fprintf('\\\\ \n');