clear,
start
close all


disp(' ')
disp('~~ Plotting your results for a NEW Region of Interest ~~')
disp('~~ PhysaliaSIM tool  - Mediterranean Sea ~~') 
disp(' ')
disp(' ')
disp('REMINDER: Your prediction runs are saved in the files')
disp('allmaps_nameofyourrun_yearsincluded.mat:')
disp(' ')

eval(['dir ',myresultsdir,'allmaps*.mat'])

disp('  ')
cual=input('Please, write the name of your run: ','s');

disp(' ')
eval(['load auxmat_',cual,'.mat'])
disp(' ')
disp('Now, write the Lats & Lons of your NEW Region of Interest:')
disp(' ')
minlonreg=input('Min longitude = ');
maxlonreg=input('Max longitude = ');
disp(' ')
minlatreg=input('Min latitude = ');
maxlatreg=input('Max latitude = ');
disp(' ')
disp('The NEW plots will be ready in a moment ...')
disp(' ')

eval(['save ',myresultsdir,'auxmat_replot.mat'])

%---2D maps of MedSea for the year prescribed by USER


PhysaliaSIM_dib2D_maps_MedSea_replot
%--- TIME plots of Probab and Encounter risk in USER'S domain

PhysaliaSIM_time1D_cal_userdomain_replot

%---- 2D maps on specific days prescribed by USER
PhysaliaSIM_dib2D_onspecificdays_replot

%
