diff --git a/MainScripts/surface3D_combine.m b/MainScripts/surface3D_combine.m
index 595d81782e06fa2743702b99d1e2a07ef1276daa..351d89d3a0c7114e1d2429f595ad1104073ca9c6 100644
--- a/MainScripts/surface3D_combine.m
+++ b/MainScripts/surface3D_combine.m
@@ -226,7 +226,7 @@ end
 function dataCurv = loadElev(PARAMS)
 %% Load a 2D elevation map and format it as a mesh
 fprintf('Loading elevation map\n');
-tiffImage = read(Tiff(PARAMS.curveLoc));
+tiffImage = imread(PARAMS.curveLoc);
 
 % Calculate the scaling factor
 scalingFactor = ceil( size(tiffImage,1)*size(tiffImage,2) / PARAMS.maxTiffImSize);