clear all

global MITprof_verbose; %to print notes to screen
if isempty(MITprof_verbose); MITprof_verbose=2; end;

if MITprof_verbose;
fprintf('creating database directory with reference datasets\n');
fprintf('===============================\n\n')
end;
if ~exist('climatology','dir'), status=mkdir('climatology'); end
gcmfaces_global;
eval(['!ln -s ' myenv.gcmfaces_dir 'sample_input/GRIDv4 climatology/']);
eval(['!ln -s ' myenv.gcmfaces_dir 'sample_input/OCCAetcONv4GRID/* climatology/']);

if MITprof_verbose;
fprintf('adding directories to your path\n');
fprintf('===============================\n\n')
end;

global MITprof_dir; MITprof_dir=[pwd '/'];
fid=fopen([MITprof_dir 'MITprof_path.m'],'wt');
aa=sprintf(['global MITprof_dir; MITprof_dir=''' MITprof_dir ''';\n']);
aa=sprintf([aa 'addpath ' MITprof_dir ';\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'profiles_IO_v2/;\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'profiles_IO_external/;\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'profiles_misc/;\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'profiles_process_main_v2/;\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'profiles_stats/;\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'ecco_v4/;\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'profiles_devel/;\n']);
aa=sprintf([aa 'addpath ' MITprof_dir 'climatology/;\n']);
eval(aa);
fid=fopen([MITprof_dir 'MITprof_path.m'],'wt'); fprintf(fid,aa); fclose(fid);

	
if MITprof_verbose;
fprintf('\ngenerating .mat files for reference grid and climatology\n')
fprintf('========================================\n\n')
end;
profiles_prep_load_fields(1);

fprintf('\n\n basic MITprof test: started... \n');

if MITprof_verbose;
fprintf('\nrunning main program on sample data sets\n')
fprintf('========================================\n\n')
end;

%select the data source (and specific params) => wod05 sample
dataset=profiles_prep_select('wod05','90CTD');
%process it
profiles_prep_main(dataset);
%
fprintf(['\n\n wod05 sample -- done -- ' dataset.dirOut dataset.fileOut '.nc was created \n\n']);

%select the data source (and specific params) => argo sample
dataset=profiles_prep_select('argo','indian'); 
%process it
profiles_prep_main(dataset);
%
fprintf(['\n\nargo sample -- done -- ' dataset.dirOut dataset.fileOut '.nc was created \n\n']);

%select the data source (and specific params) => odv sample
dataset=profiles_prep_select('odv','ODVcompact_sample');
%process it
profiles_prep_main(dataset);
%
fprintf(['\n\nodv sample -- done -- ' dataset.dirOut dataset.fileOut '.nc was created \n\n']);

fprintf('\n basic MITprof test: completed. \n');