% url = 'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/jarkus/profiles/transect_r2013_filled.nc';
url = 'C:\Users\bieman\Downloads\transect_r2013_filled.nc';
SaveDir = 'd:\checkouts\openearthmodels\deltares\profile_fluctuation_dune_erosion\matlab\results_report_20140730';

% Name    = 'Vlieland';
% ID      = 5004898;
% station1 = 'Steunpunt Waddenzee';
% station2 = 'Eierlandse Gat';
% SaveDir = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\normaal';
% SaveDirPV = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\profiel_variatie';
% TargetPfID = 2;
% Years   = [1991:2001 2003:2013];
% ReturnFreqs = [5e-4; 5e-5; 1e-5];
% TRDAmarker = 'm^';
% VolumeMean  = -2.7;
% VolumeSTD   = 57.1;


Name    = 'Texel';
ID      = 6002011;
station1 = 'Den Helder';
station2 = 'Eierlandse Gat';
SaveDir = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\normaal';
SaveDirPV = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\profiel_variatie';
TargetPfID = 2;
Years   = [1991:2001 2003:2013];
ReturnFreqs = [2.5e-4; 2.5e-5; 1e-5];
TRDAmarker = 'm^';
VolumeMean  = -6.3;
VolumeSTD   = 72.3;
% 
% Name    = 'Egmond';
% ID      = 7003975;
% station1 = 'Den Helder';
% station2 = 'IJmuiden';
% SaveDir = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\normaal';
% SaveDirPV = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\profiel_variatie';
% TargetPfID = 3;
% Years   = [1991:2001 2003:2010 2012:2013];
% ReturnFreqs = [1e-3; 1e-4; 1e-5];
% VolumeMean  = -3.4;
% VolumeSTD   = 82.5;

% Name    = 'Katwijk';
% ID      = 8008950;
% station1 = 'Hoek van Holland';
% station2 = 'IJmuiden';
% SaveDir = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\normaal';
% SaveDirPV = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\profiel_variatie';
% TargetPfID = 3;
% Years   = [1999:2001 2003:2013];%[1991:2001 2003:2013];
% ReturnFreqs = [1e-3; 1e-4; 1e-5];
% VolumeMean  = 7.7;
% VolumeSTD   = 89.1;
% 
% Name    = 'sGravenzande_VVarNew';
% ID      = 9011450;
% station1 = 'Hoek van Holland';
% station2 = 'IJmuiden';
% SaveDir = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\normaal';
% SaveDirPV = 'n:\Projects\1209000\1209436\F. Other information\Temporele Effecten\Figuren\Rt_diagrammen\profiel_variatie';
% TargetPfID = 3;
% Years   = [1999:2001 2003:2013];%[1991:2001 2003:2013];
% ReturnFreqs = [1e-3; 1e-4; 1e-5];
% VolumeMean  = 1.0;
% VolumeSTD   = 66.9;

[MCresult, FORMresult] = profile_fluctuation_probcalc(...
    'jarkus_url', url, ...
    'year', Years, ...
    'ReturnFreqs', ReturnFreqs, ...
    'id', ID, ...
    'ProfileVariation', true, ...
    'VariationVolumeMean', VolumeMean, ...
    'VariationVolumeSTD', VolumeSTD, ...
    'station1', station1, 'station2', station2);