readSurferData

The function was written to read surfer 6 binary GRD data file specifically suited for reading out models and doing calculations in matlab.

% Function prototype
function [xRange,yRange,matData] = readSurferData(fName,varargin)
% this function was written to read surfer 6 binary GRD data file
% specifically suited for reading out models obtained from
% other softwares
% details of the surfer 6 binary GRD file can be found:
% https://grapherhelp.goldensoftware.com/subsys/surfer_6_grid_file_format.htm

%%
% S. Koley
% Gran Sasso Science Institute
% Dept. of Physics
% L'Aquila, Italy
% 2022

%%
% inputs
% first arg = filename
% second arg = 0, if not output plot, else set to 1, default = 0
% third arg: values in the data matrix to be set to NaN, this corresponds
% to blank values, default to 10^18

%%
% output format
% matData: is the matrix containing the grid values with the topmost
% row corresponding to the minimum y, and first column corresponding
% to the mimimum x
% xRange  = vector x values (xLow->xHigh)
% yRange = vector of y values (yLow->yHigh)