clear all;fclose all;close all;clc

cd f:\mctools\mc_toolbox; mcsettings
addpath(fileparts(mfilename('fullpath')));

% edit getCalls getIsCalledBy getInputSize getInputVariables roundoff double2evalstr getFileName getVolume
writemessage('init');

%% getInputVariables / getInputSize
clear;fclose all;close all;clc
quickstart internet
d = readTransectData('Jarkus Data','Noord-Holland','03000','2006');
xInitial    = d.xe(~isnan(d.ze)); %keep only the points with non-NaN z-values
zInitial    = d.ze(~isnan(d.ze)); %keep only the points with non-NaN z-values

dbstop in getVolume at 42
dbstopcurrent
[Volume, result, Boundaries] = getVolume(xInitial, zInitial);
dbclear in getVolume at 42
dbstopcurrent
writemessage(3,'Example of getInputVariables completed');

%% roundoff
clear;fclose all;close all;clc
X = 87643.7852;
for n = -2:2:4
    Xround = roundoff(X, n);
    disp(num2str(Xround))
end

%% getdefaults
dbstop in roundoff at 40
Xround = roundoff(X);
dbclear in roundoff at 40

%% double2evalstr
clear;fclose all;close all;clc

dbstop in double2evalstr_test at 9
run double2evalstr_test
dbstopcurrent

dbstop getHsig_t_test at 30
dbstop getHsig_t_test at 41
run getHsig_t_test
dbstopcurrent
writemessage(4,'Example of double2evalstr_test completed');

%% getFileName
clear;fclose all;close all;clc
directory = 'D:\heijer\My Documents\Work\TU Delft\PhD project\Prob2B\ParametricStudy\DenHeijer2008\Case_10';
extension = 'txt';
exception = 'DP_Case_10_maxRD=068.7.txt';
fileid = [];
dbstopcurrent
fname = getFileName(directory, extension, exception, fileid)
writemessage(5,'Example of getFileName completed');

%% writemessage
writemessage(7,'Staring example of writemessage');
dbstopcurrent
SavedMessages = writemessage('get')

%% IsCalls
fun = 'getDuneErosion_VTV2006';
directory = 'D:\heijer\My Documents\Program Files\McTools\mc_toolbox\mc_applications\mc_ucit\getDuneErosion\';
IsCalls = getCalls(fun, directory);
IsCalls = getCalls(fun, directory, 'quiet');
dbstopcurrent
writemessage(1,'Example of IsCalls completed');

%% IsCalledBy
clear;fclose all;close all;clc
fun = 'getVolume';
directory = 'D:\heijer\My Documents\Program Files\McTools\mc_toolbox\mc_applications\mc_ucit\getDuneErosion\';
[IsCalledBy, Line, Column] = getIsCalledBy(fun, directory);
[IsCalledBy, Line, Column] = getIsCalledBy(fun, directory, 'quiet');
dbstopcurrent
writemessage(2,'Example of IsCalledBy completed');

%% FD_GUI_figure
cd( directory);
run FD_GUI_figure
dbstopcurrent
writemessage(6,'Example of FD_GUI completed');
clc

%% guidisp
guidisp('set','FD_GUI_message_panel');
guidisp('test message panel');
guidisp('set',[]);
guidisp('Test Command line');
dbstopcurrent
writemessage(7,'Example of guidisp completed');

%% dbstate
dbstop in getDuneErosion_test at 6
dbstop in getDuneErosion_VTV2006 at 12
getDuneErosion_test('VTV2006',1);
dbstopcurrent
writemessage(7,'Example of dbstate completed');

%% dbstopcurrent
cd('D:\heijer\My Documents\Work\Deltares\Projects\McTools');
dbstopcurrent;
pause(0.1);

%% writemessage
SavedMessages = writemessage('get')
dbclear all

INTEGER*4 FUNCTION fortomat(val1,val2)
!DEC$ ATTRIBUTES DLLEXPORT,ALIAS : "fortomat" :: fortomat
INTEGER*4 val1,val2    
FORTOMAT=val1*val2
END FUNCTION

real*4 FUNCTION fortomatreal(val1,val2)
!DEC$ ATTRIBUTES DLLEXPORT,ALIAS : "fortomatreal" :: fortomatreal
real*4 val1,val2    
fortomatreal=val1*val2
END FUNCTION

function out = fortomat(var1,var2)
%FORTOMAT
%
% fortomat(val1,val2) 
% where var* can be integer or float

   OPT.libname      = 'ekufbqegjhbgoebgo';
   OPT.dllname      = 'endecdll.dll';;
   OPT.hdrname      = 'endecdll.h'
   OPT.debug        = 1;

   %% Load library
   %% ------------
   if libisloaded(OPT.libaliasname)~=1
   loadlibrary(OPT.dllname,OPT.hdrname,'alias',OPT.libname)
   end;

   %% Inquire library
   %% ------------
   if OPT.debug
   libfunctions(OPT.libaliasname,'-full')
   end

   %% Call library function
   %% ------------
   if isa(var1,'integer') && isa(var1,'integer')
   out = calllib(OPT.libname,'fortomat',var1,var2);
   elseif isa(var1,'float') && isa(var1,'float')
   out = calllib(OPT.libname,'fortomat',var1,var2);
   else
   error(['No matching datatype in ',])
   end

   %% Unload library
   %% ------------
   if libisloaded(OPT.libname)==1
    unloadlibrary(OPT.libname);
   end;
%% EOF

C 
C CHECK FOR PROPER NUMBER OF ARGUMENTS 
C       
IF (NRHS .NE. 2) THEN         
CALL MEXERRMSGTXT('YPRIME requires two input arguments')       
ELSEIF (NLHS .GT. 1) THEN         
CALL MEXERRMSGTXT('YPRIME requires one output argument')       
ENDIF

C 
C CHECK THE DIMENSIONS OF Y.  IT CAN BE 4 X 1 OR 1 X 4. 
C
M = MXGETM(PRHS(2))
N = MXGETN(PRHS(2)) 
C       
IF ((MAX(M,N) .NE. 4) .OR. (MIN(M,N) .NE. 1)) THEN         
CALL MEXERRMSGTXT('YPRIME requires that Y be a 4 x 1 vector')       
ENDIF