function output = getSolid(path, filePrefix,resinparamobj,xAxis,slice,res)
cd(path)
a = dir([path '\dose*']);
nLayers = size(a,1)-1;

for cLayer =1:nLayers
    doseFile = strcat(filePrefix,num2str(cLayer-1,'%02.0f'),'.csv');
    doseLayer = readmatrix(doseFile);
    if(xAxis)
        doseSlice(cLayer,:) = doseLayer(slice,:);
    else
        doseSlice(cLayer,:) = doseLayer(:,slice);
    end
end

doseSlice(nLayers+1,:)=zeros(1,size(doseSlice,2));

phi = 1-exp(-resinparamobj.K*doseSlice);
csvwrite('phi2DMat.csv',phi);
xPos = res:res:res*size(doseSlice,2);
zPos = res:res:res*size(doseSlice,1);
h = contour(xPos,zPos,phi,[resinparamobj.phi_c resinparamobj.phi_c]);

if size(h,2)>2
    xyContour(:,1)=h(1,2:end);
    xyContour(:,2)=h(2,2:end);
    csvwrite('contour2D.csv',xyContour);
end

output = true;
end