0001 function plot_results(Action,varargin)
0002
0003
0004
0005
0006
0007
0008
0009 global index x y u v vel
0010
0011 if nargin==0,
0012 Action='initialize';
0013 series=1;
0014 end
0015
0016 if nargin==1,
0017 if length(Action)==1,
0018 series=Action;
0019 Action='initialize';
0020 end
0021 end
0022
0023
0024
0025 if nargin==2,
0026 series=varargin{1};
0027 end
0028
0029 switch Action
0030
0031
0032 case 'initialize'
0033 LocalOpenFig(series);
0034 case 'ToggleButton'
0035
0036
0037 value=findobj(gcf,'Tag','Choice');
0038 value=popupstr(value);
0039
0040
0041 year=findobj(gcf,'Tag','Version');
0042 year_value=popupstr(year);
0043
0044 zoom=findobj(gcf,'Tag','Zoom');
0045 zoom_value=popupstr(zoom);
0046
0047 ser=findobj(gcf,'Tag','Series');
0048 series=popupstr(ser);
0049 series=str2num(series);
0050
0051 diff=findobj(gcf,'Tag','Diff');
0052 diff_value=popupstr(diff);
0053
0054 diff_choice=findobj(gcf,'Tag','Diff_choice');
0055 diff_choice=popupstr(diff_choice);
0056
0057 range_choice=findobj(gcf,'Tag','Range_choice');
0058 range_choice=popupstr(range_choice);
0059 if strcmp(range_choice,'auto'),
0060 range_choice=0;
0061 else
0062 range_choice=1;
0063 end
0064
0065 range=findobj(gcf,'Tag','Range');
0066 range=deblank(get(range,'String'));
0067
0068
0069
0070
0071
0072 switch zoom_value,
0073 case 'RONNE'
0074 zo=0;
0075 case 'HIR'
0076 zo=1;
0077 case 'LASSITER'
0078 zo=2;
0079 case 'MIDDLE'
0080 zo=3;
0081
0082 end
0083
0084
0085 switch value
0086 case 'obs_ronne'
0087 plot_obs_ronne(year_value,zo,series,range_choice,range);
0088 case 'obs_hir_vx'
0089 plot_obs_hir_vx(year_value,zo,series,range_choice,range);
0090 case 'obs_hir_vy'
0091 plot_obs_hir_vy(year_value,zo,series,range_choice,range);
0092 case 'obs_hir_vel'
0093 plot_obs_hir_vel(year_value,zo,series,range_choice,range);
0094 case 'u'
0095 plot_u(year_value,zo,series,range_choice,range);
0096 case 'diff_u_hir'
0097 plot_diff_u_hir(year_value,zo,series,range_choice,range);
0098 case 'v'
0099 plot_v(year_value,zo,series,range_choice,range);
0100 case 'v_prism'
0101 plot_v_prism(year_value,zo,series,range_choice,range);
0102 case 'diff_v_hir'
0103 plot_diff_v_hir(year_value,zo,series,range_choice,range);
0104 case 'vel'
0105 plot_vel(year_value,zo,series,range_choice,range);
0106 case 'vel_prism'
0107 plot_vel_prism(year_value,zo,series,range_choice,range);
0108 case 'diff_vel_ronne'
0109 plot_diff_vel_ronne(year_value,zo,series,range_choice,range);
0110 case 'diff_vel_hir'
0111 plot_diff_vel_hir(year_value,zo,series,range_choice,range);
0112 case 'rotation_hir'
0113 plot_rotation_hir(year_value,zo,series,range_choice,range);
0114 case 'rotation_model'
0115 plot_rotation_model(year_value,zo,series,range_choice,range);
0116 case 'diff_rotation'
0117 plot_diff_rotation(year_value,zo,series,range_choice,range);
0118 case 'quiver_hir'
0119 plot_quiver_hir(year_value,zo,series,range_choice,range);
0120 case 'quiver_las'
0121 plot_quiver_las(year_value,zo,series,range_choice,range);
0122 case 'quiver_model'
0123 plot_quiver_model(year_value,zo,series,range_choice,range);
0124 case 'mesh'
0125 plot_mesh(year_value,zo,series,range_choice,range);
0126 case 'obs_LAS_vx'
0127 plot_obs_LAS_vx(year_value,zo,series,range_choice,range);
0128 case 'obs_LAS_vy'
0129 plot_obs_LAS_vy(year_value,zo,series,range_choice,range);
0130 case 'obs_LAS_vel'
0131 plot_obs_LAS_vel(year_value,zo,series,range_choice,range);
0132 case 'diff_u_LAS'
0133 plot_diff_u_LAS(year_value,zo,series,range_choice,range);
0134 case 'diff_v_LAS'
0135 plot_diff_v_LAS(year_value,zo,series,range_choice,range);
0136 case 'diff_vel_LAS'
0137 plot_diff_vel_LAS(year_value,zo,series,range_choice,range);
0138
0139
0140
0141
0142
0143
0144
0145 case 'tau1'
0146 plot_tau1(year_value,zo,series,range_choice,range);
0147 case 'tau2'
0148 plot_tau2(year_value,zo,series,range_choice,range);
0149 case 'teta1'
0150 plot_teta1(year_value,zo,series,range_choice,range);
0151 case 'teta2'
0152 plot_teta2(year_value,zo,series,range_choice,range);
0153
0154
0155 case 'thick'
0156 plot_thick(year_value,zo,series,range_choice,range);
0157 case 'thickbar'
0158 plot_thickbar(year_value,zo,series,range_choice,range);
0159
0160
0161 otherwise
0162 plot_results('Refresh');
0163 end
0164
0165
0166
0167 case 'Refresh'
0168 year=findobj(gcf,'Tag','Version');
0169 year_value=popupstr(year);
0170 quantity=findobj(gcf,'Tag','Choice');
0171 quantity_value=popupstr(quantity);
0172 zoom=findobj(gcf,'Tag','Zoom');
0173 zoom_value=popupstr(zoom);
0174 ser=findobj(gcf,'Tag','Series');
0175 series=popupstr(ser);
0176 diff=findobj(gcf,'Tag','Diff');
0177 diff_value=popupstr(diff);
0178 range_choice=findobj(gcf,'Tag','Range_choice');
0179 range_choice=popupstr(range_choice);
0180 if strcmp(range_choice,'auto'),
0181 range_choice=0;
0182 else
0183 range_choice=1;
0184 end
0185
0186 range=findobj(gcf,'Tag','Range');
0187 range=deblank(get(range,'String'));
0188
0189
0190
0191 switch zoom_value,
0192 case 'RONNE'
0193 zo=0;
0194 case 'HIR'
0195 zo=1;
0196 case 'LASSITER'
0197 zo=2;
0198 case 'MIDDLE'
0199 zo=3;
0200 end
0201
0202 eval(['plot_' quantity_value '(' year_value ',' num2str(zo) ',' series ',' range_choice ',' range ' )']);
0203
0204
0205 end
0206
0207
0208
0209
0210
0211
0212
0213
0214 function LocalOpenFig(series);
0215 StdUnit = 'point';
0216 StdColor = get(0,'DefaultUIcontrolBackgroundColor');
0217 PointsPerPixel = 72/get(0,'ScreenPixelsPerInch');
0218
0219 a = figure('Color',[0.8 0.8 0.8], ...
0220 'Position',[25 45 560 420], ...
0221 'IntegerHandle','on', ...
0222 'Number','off',...
0223 'Name','Results',...
0224 'Tag','Results Visualisation');
0225
0226 b = axes('Parent',a, ...
0227 'Units',StdUnit, ...
0228 'Box','on', ...
0229 'Color',[1 1 1], ...
0230 'Position',PointsPerPixel*[40 45 360 330], ...
0231 'Units','Normalized',...
0232 'XColor',[0 0 0], ...
0233 'XGrid','off', ...
0234 'XTick',[],...
0235 'YColor',[0 0 0], ...
0236 'YGrid','off', ...
0237 'YTick',[], ...
0238 'ZColor',[0 0 0], ...
0239 'ZGrid','off', ...
0240 'ZTick',[]);
0241
0242 set(get(b,'XLabel'),'visible','off');
0243 set(get(b,'YLabel'),'visible','off');
0244 set(get(b,'ZLabel'),'visible','off');
0245
0246
0247 b = axes('Parent',a,'Units',StdUnit,'Color',[1 1 1],'Position',...
0248 PointsPerPixel*[40 45 360 330],'Units','Normalized', 'Tag','ActiveAxes','Box','on','XColor',[0 0 0], ...
0249 'XGrid','off', ...
0250 'YColor',[0 0 0], ...
0251 'YGrid','off');
0252 names_function(1,series);
0253 load ../data/matlab/commondata/names.mat
0254 eval(['load ' root directory '\data']);
0255 aa=['load ' root directory '/' results '_results'];
0256 eval(aa)
0257 aa=['load ' root directory '/' maillage '_datamesh'];
0258 eval(aa)
0259 aa=['load ' root directory '/' maillage '_fem_parameters'];
0260 eval(aa)
0261
0262 u=u*yts;
0263 v=v*yts;
0264
0265
0266
0267
0268 trisurf(index,x/1e3,y/1e3,u),shading interp,view(2);
0269 caxis([-500 100]);
0270 colorbar
0271 set(gca,'Tag','ActiveAxes');
0272
0273
0274
0275 b = uicontrol('Parent',a, ...
0276 'Unit',StdUnit, ...
0277 'BackgroundColor',StdColor, ...
0278 'Position',PointsPerPixel*[425 295 100 50], ...
0279 'Units','Normalized',...
0280 'Style','frame');
0281
0282 b = uicontrol('Parent',a, ...
0283 'Unit',StdUnit, ...
0284 'BackgroundColor',StdColor, ...
0285 'HorizontalAlignment','left', ...
0286 'Position',PointsPerPixel*[430 320 80 20], ...
0287 'Units','Normalized',...
0288 'String','Series', ...
0289 'Style','text');
0290
0291
0292 str='';
0293 i=1;
0294 while exist([root 'RONNE_w' num2str(i)],'file')==7,
0295 str=[str num2str(i) '|'];
0296 i=i+1;
0297 end
0298 str=str(1:length(str)-1);
0299
0300
0301 b = uicontrol('Parent',a, ...
0302 'Unit',StdUnit, ...
0303 'BackgroundColor',[1 1 1], ...
0304 'Style','popup', ...
0305 'Position',PointsPerPixel*[430 300 90 20], ...
0306 'Units','Normalized',...
0307 'String',str, ...
0308 'callback','plot_results(''ToggleButton'');', ...
0309 'Tag','Series',...
0310 'Value',series);
0311
0312
0313
0314
0315
0316 b = uicontrol('Parent',a, ...
0317 'Unit',StdUnit, ...
0318 'BackgroundColor',StdColor, ...
0319 'Position',PointsPerPixel*[425 245 100 50], ...
0320 'Units','Normalized',...
0321 'Style','frame');
0322
0323
0324 b = uicontrol('Parent',a, ...
0325 'Unit',StdUnit, ...
0326 'BackgroundColor',StdColor, ...
0327 'HorizontalAlignment','left', ...
0328 'Position',PointsPerPixel*[430 270 80 20], ...
0329 'Units','Normalized',...
0330 'String','Version', ...
0331 'Style','text');
0332
0333
0334
0335 str='';
0336 i=1;
0337 names_function(1,series);
0338 load ../data/matlab/commondata/names.mat
0339 while exist_by_names([root directory '\' results '_results.mat']),
0340 str=[str num2str(i) '|'];
0341 i=i+1;
0342 names_function(i,series);
0343 load ../data/matlab/commondata/names.mat
0344 end
0345 str=str(1:length(str)-1);
0346
0347 b = uicontrol('Parent',a, ...
0348 'Unit',StdUnit, ...
0349 'BackgroundColor',[1 1 1], ...
0350 'Style','popup', ...
0351 'Position',PointsPerPixel*[430 250 90 20], ...
0352 'Units','Normalized',...
0353 'String',str, ...
0354 'callback','plot_results(''ToggleButton'');', ...
0355 'Tag','Version');
0356
0357
0358
0359
0360 b = uicontrol('Parent',a, ...
0361 'Unit',StdUnit, ...
0362 'BackgroundColor',StdColor, ...
0363 'Position',PointsPerPixel*[425 195 100 50], ...
0364 'Units','Normalized',...
0365 'Style','frame');
0366
0367 b = uicontrol('Parent',a, ...
0368 'Unit',StdUnit,'BackgroundColor',StdColor, ...
0369 'HorizontalAlignment','left', ...
0370 'Position',PointsPerPixel*[430 220 80 20], ...
0371 'Units','Normalized',...
0372 'String','Quantity to view', ...
0373 'Style','text');
0374
0375
0376
0377 b = uicontrol('Parent',a, ...
0378 'Unit',StdUnit, ...
0379 'BackgroundColor',[1 1 1], ...
0380 'Style','popup', ...
0381 'Position',PointsPerPixel*[430 200 90 20], ...
0382 'Units','Normalized',...
0383 'String','u|v|v_prism|vel|vel_prism|obs_hir_vx|obs_hir_vy|obs_hir_vel|diff_u_hir|diff_v_hir|diff_vel_hir|obs_ronne|diff_vel_ronne|obs_LAS_vx|obs_LAS_vy|obs_LAS_vel|diff_u_LAS|diff_v_LAS|diff_vel_LAS|rotation_model|rotation_hir|diff_rotation|quiver_hir|quiver_las|quiver_model|tau1|tau2|teta1|teta2|thick|thickbar|mesh', ...
0384 'callback','plot_results(''ToggleButton'');', ...
0385 'Tag','Choice');
0386
0387
0388
0389
0390
0391 b = uicontrol('Parent',a, ...
0392 'Unit',StdUnit, ...
0393 'BackgroundColor',StdColor, ...
0394 'Position',PointsPerPixel*[425 145 100 50], ...
0395 'Units','Normalized',...
0396 'Style','frame');
0397
0398 b = uicontrol('Parent',a, ...
0399 'Unit',StdUnit, ...
0400 'BackgroundColor',StdColor, ...
0401 'HorizontalAlignment','left', ...
0402 'Position',PointsPerPixel*[430 170 80 20], ...
0403 'Units','Normalized',...
0404 'String','Zoom', ...
0405 'Style','text');
0406
0407
0408
0409 b = uicontrol('Parent',a, ...
0410 'Unit',StdUnit, ...
0411 'BackgroundColor',[1 1 1], ...
0412 'Style','popup', ...
0413 'Position',PointsPerPixel*[430 150 90 20], ...
0414 'Units','Normalized',...
0415 'String','RONNE|HIR|LASSITER|MIDDLE', ...
0416 'callback','plot_results(''ToggleButton'');', ...
0417 'Tag','Zoom');
0418
0419
0420
0421 set(a,'visible','on','HandleVisibility','callback');
0422
0423
0424 b = uicontrol('Parent',a, ...
0425 'Unit',StdUnit, ...
0426 'BackgroundColor',StdColor, ...
0427 'Position',PointsPerPixel*[425 115 50 30], ...
0428 'Units','Normalized',...
0429 'Style','frame');
0430
0431
0432 b = uicontrol('Parent',a, ...
0433 'Unit',StdUnit, ...
0434 'BackgroundColor',StdColor, ...
0435 'HorizontalAlignment','left', ...
0436 'Position',PointsPerPixel*[430 130 40 10], ...
0437 'Units','Normalized',...
0438 'String','Diff Version', ...
0439 'Style','text');
0440
0441
0442
0443 str='';
0444 i=1;
0445 names_function(1,series);
0446 load ../data/matlab/commondata/names.mat
0447 while exist_by_names([root directory '\' results '_results.mat']),
0448 str=[str num2str(i) '|'];
0449 i=i+1;
0450 names_function(i,series);
0451 load ../data/matlab/commondata/names.mat
0452 end
0453 str=str(1:length(str)-1);
0454
0455 b = uicontrol('Parent',a, ...
0456 'Unit',StdUnit, ...
0457 'BackgroundColor',[1 1 1], ...
0458 'Style','popup', ...
0459 'Position',PointsPerPixel*[430 120 40 10], ...
0460 'Units','Normalized',...
0461 'String',str, ...
0462 'callback','plot_results(''ToggleButton'');', ...
0463 'Tag','Diff');
0464
0465
0466
0467
0468
0469 b = uicontrol('Parent',a, ...
0470 'Unit',StdUnit, ...
0471 'BackgroundColor',StdColor, ...
0472 'Position',PointsPerPixel*[475 115 50 30], ...
0473 'Units','Normalized',...
0474 'Style','frame');
0475
0476 b = uicontrol('Parent',a, ...
0477 'Unit',StdUnit, ...
0478 'BackgroundColor',StdColor, ...
0479 'HorizontalAlignment','left', ...
0480 'Position',PointsPerPixel*[480 130 40 10], ...
0481 'Units','Normalized',...
0482 'String','Diff mesh', ...
0483 'Style','text');
0484
0485
0486
0487 b = uicontrol('Parent',a, ...
0488 'Unit',StdUnit, ...
0489 'BackgroundColor',[1 1 1], ...
0490 'Style','popup', ...
0491 'Position',PointsPerPixel*[480 120 40 10], ...
0492 'Units','Normalized',...
0493 'String','Version|Diff_version', ...
0494 'callback','plot_results(''ToggleButton'');', ...
0495 'Tag','Mesh_choice');
0496
0497
0498
0499 set(a,'visible','on','HandleVisibility','callback');
0500
0501
0502
0503 b = uicontrol('Parent',a, ...
0504 'Unit',StdUnit, ...
0505 'BackgroundColor',StdColor, ...
0506 'Position',PointsPerPixel*[425 85 50 30], ...
0507 'Units','Normalized',...
0508 'Style','frame');
0509
0510
0511 b = uicontrol('Parent',a, ...
0512 'Unit',StdUnit, ...
0513 'BackgroundColor',StdColor, ...
0514 'HorizontalAlignment','left', ...
0515 'Position',PointsPerPixel*[430 100 40 10], ...
0516 'Units','Normalized',...
0517 'String','Plot differentials', ...
0518 'Style','text');
0519
0520
0521
0522 b = uicontrol('Parent',a, ...
0523 'Unit',StdUnit, ...
0524 'BackgroundColor',[1 1 1], ...
0525 'Style','popup', ...
0526 'Position',PointsPerPixel*[430 90 40 10], ...
0527 'Units','Normalized',...
0528 'String','no|yes', ...
0529 'callback','plot_results(''ToggleButton'');', ...
0530 'Tag','Diff_choice');
0531
0532
0533
0534 set(a,'visible','on','HandleVisibility','callback');
0535
0536 b = uicontrol('Parent',a, ...
0537 'Unit',StdUnit, ...
0538 'BackgroundColor',StdColor, ...
0539 'Position',PointsPerPixel*[475 85 50 30], ...
0540 'Units','Normalized',...
0541 'Style','frame');
0542
0543 b = uicontrol('Parent',a, ...
0544 'Unit',StdUnit, ...
0545 'BackgroundColor',StdColor, ...
0546 'HorizontalAlignment','left', ...
0547 'Position',PointsPerPixel*[485 105 30 5], ...
0548 'Units','Normalized',...
0549 'String','Color Range', ...
0550 'Style','text');
0551
0552 b = uicontrol('Parent',a, ...
0553 'Unit',StdUnit, ...
0554 'BackgroundColor',[1 1 1], ...
0555 'Position',PointsPerPixel*[480 95 40 10], ...
0556 'Units','Normalized',...
0557 'String','[ , ]',...
0558 'callback','plot_results(''ToggleButton'');', ...
0559 'Style','edit',...
0560 'Tag','Range');
0561
0562
0563 b = uicontrol('Parent',a, ...
0564 'Unit',StdUnit, ...
0565 'BackgroundColor',[1 1 1], ...
0566 'Style','popup', ...
0567 'Position',PointsPerPixel*[480 90 40 5], ...
0568 'Units','Normalized',...
0569 'String','auto|manual', ...
0570 'callback','plot_results(''ToggleButton'');', ...
0571 'Tag','Range_choice');
0572
0573
0574
0575 set(a,'visible','on','HandleVisibility','callback');
0576
0577
0578 b = uicontrol('Parent',a, ...
0579 'Unit',StdUnit, ...
0580 'BackgroundColor',StdColor, ...
0581 'Position',PointsPerPixel*[475 55 50 30], ...
0582 'Units','Normalized',...
0583 'Style','frame');
0584
0585 b = uicontrol('Parent',a, ...
0586 'Unit',StdUnit, ...
0587 'BackgroundColor',StdColor, ...
0588 'HorizontalAlignment','left', ...
0589 'Position',PointsPerPixel*[485 70 30 10], ...
0590 'Units','Normalized',...
0591 'String','Quiver sparsity', ...
0592 'Style','text');
0593
0594 b = uicontrol('Parent',a, ...
0595 'Unit',StdUnit, ...
0596 'BackgroundColor',[1 1 1], ...
0597 'Position',PointsPerPixel*[480 60 40 10], ...
0598 'Units','Normalized',...
0599 'String','3',...
0600 'callback','plot_results(''ToggleButton'');', ...
0601 'Style','edit',...
0602 'Tag','Quiver_Sparsity');
0603
0604
0605 set(a,'visible','on','HandleVisibility','callback');
0606
0607
0608
0609 function plot_obs_ronne(num,zo,series,range_choice,range);
0610
0611 names_function(num,series)
0612 load ../data/matlab/commondata/names.mat
0613 load ../data/matlab/commondata/data
0614 exist=exist_by_names([root directory '\' results '_results.mat']);
0615
0616 if exist==1,
0617
0618 aa=['load ' root directory '/' results '_results'];
0619 eval(aa)
0620 aa=['load ' root directory '/' maillage '_datamesh'];
0621 eval(aa)
0622 aa=['load ' root directory '/' maillage '_fem_parameters'];
0623 eval(aa)
0624 if exist_by_names([root directory '\' maillage '_observations_Ian.mat']),
0625 eval(['load ' root directory '/' maillage '_observations_Ian'])
0626 v_input=sqrt(u_RONNE.^2+v_RONNE.^2)*yts;
0627 end
0628
0629
0630
0631 trisurf(index,x/1e3,y/1e3,v_input),shading interp,view(2);
0632 xlabel('X (km)','fontname','arial','fontsize',16)
0633 ylabel('Y (km)','fontname','arial','fontsize',16)
0634 title(texlabel('Observed velocity magnitude (m.a^(-1)). BAS data.'),'fontname','arial','fontsize',18);
0635 set(gca,'fontname','arial','fontsize',16)
0636
0637 zoom_choice(zo);
0638 colormap(jet)
0639 caxis([0 1500]);
0640 colorbar
0641 set(gca,'Tag','ActiveAxes')
0642 end
0643
0644
0645
0646
0647
0648 function plot_obs_hir_vx(num,zo,series,range_choice,range);
0649
0650 names_function(num,series)
0651 load ../data/matlab/commondata/names.mat
0652 load ../data/matlab/commondata/data
0653 exist=exist_by_names([root directory '\' results '_results.mat']);
0654
0655 if exist==1,
0656
0657 aa=['load ' root directory '/' results '_results'];
0658 eval(aa)
0659 aa=['load ' root directory '/' maillage '_datamesh'];
0660 eval(aa)
0661 aa=['load ' root directory '/' maillage '_fem_parameters'];
0662 eval(aa)
0663 aa=['load ' root directory '/' maillage '_observations'];
0664 eval(aa)
0665
0666
0667 trisurf(index,x/1e3,y/1e3,u_HIR),shading interp,view(2);
0668 xlabel('X (km)','fontname','arial','fontsize',16)
0669 ylabel('Y (km)','fontname','arial','fontsize',16)
0670 title(texlabel('Observed velocity Vx (m.a^(-1)) around HIR. From RSAT1'),'fontname','arial','fontsize',18);
0671 set(gca,'fontname','arial','fontsize',16)
0672
0673
0674 zoom_choice(zo);
0675 colormap(jet)
0676 caxis([-500 100]);
0677 colorbar
0678 set(gca,'Tag','ActiveAxes')
0679 end
0680
0681
0682
0683
0684
0685
0686
0687
0688
0689 function plot_obs_hir_vy(num,zo,series,range_choice,range);
0690
0691 names_function(num,series)
0692 load ../data/matlab/commondata/names.mat
0693 load ../data/matlab/commondata/data
0694 exist=exist_by_names([root directory '\' results '_results.mat']);
0695
0696 if exist==1,
0697
0698 aa=['load ' root directory '/' results '_results'];
0699 eval(aa)
0700 aa=['load ' root directory '/' maillage '_datamesh'];
0701 eval(aa)
0702 aa=['load ' root directory '/' maillage '_fem_parameters'];
0703 eval(aa)
0704 aa=['load ' root directory '/' maillage '_observations'];
0705 eval(aa)
0706
0707 trisurf(index,x/1e3,y/1e3,v_HIR),shading interp,view(2);
0708 xlabel('X (km)','fontname','arial','fontsize',16)
0709 ylabel('Y (km)','fontname','arial','fontsize',16)
0710 title(texlabel('Observed velocity Vy (m.a^(-1)) around HIR. From RSAT1'),'fontname','arial','fontsize',18);
0711 set(gca,'fontname','arial','fontsize',16)
0712
0713 zoom_choice(zo);
0714 colormap(jet)
0715 caxis([0 1500]);
0716 colorbar
0717 set(gca,'Tag','ActiveAxes')
0718 end
0719
0720
0721
0722
0723
0724
0725
0726
0727
0728 function plot_obs_hir_vel(num,zo,series,range_choice,range);
0729
0730 names_function(num,series)
0731 load ../data/matlab/commondata/names.mat
0732 load ../data/matlab/commondata/data
0733 exist=exist_by_names([root directory '\' results '_results.mat']);
0734
0735 if exist==1,
0736
0737 aa=['load ' root directory '/' results '_results'];
0738 eval(aa)
0739 aa=['load ' root directory '/' maillage '_datamesh'];
0740 eval(aa)
0741 aa=['load ' root directory '/' maillage '_fem_parameters'];
0742 eval(aa)
0743 aa=['load ' root directory '/' maillage '_observations'];
0744 eval(aa)
0745
0746 vel_HIR=sqrt(u_HIR.^2+v_HIR.^2);
0747
0748 trisurf(index,x/1e3,y/1e3,vel_HIR),shading interp,view(2);
0749 xlabel('X (km)','fontname','arial','fontsize',16)
0750 ylabel('Y (km)','fontname','arial','fontsize',16)
0751 title(texlabel('Observed velocity magnitude (m.a^(-1)) around HIR. From RSAT1'),'fontname','arial','fontsize',18);
0752 set(gca,'fontname','arial','fontsize',16)
0753
0754 zoom_choice(zo);
0755
0756 colormap(jet)
0757 caxis([0 1500]);
0758 colorbar
0759 set(gca,'Tag','ActiveAxes')
0760 end
0761
0762
0763
0764
0765
0766
0767
0768 function plot_u(num,zo,series,range_choice,range);
0769 diff_choice=findobj(gcf,'Tag','Diff_choice');
0770 diff_choice=popupstr(diff_choice);
0771 if strcmp(diff_choice,'no')==1,
0772
0773
0774 names_function(num,series)
0775 load ../data/matlab/commondata/names.mat
0776 load ../data/matlab/commondata/data
0777 exist=exist_by_names([root directory '\' results '_results.mat']);
0778
0779 if exist==1,
0780
0781 aa=['load ' root directory '/' results '_results'];
0782 eval(aa)
0783 aa=['load ' root directory '/' maillage '_datamesh'];
0784 eval(aa)
0785 aa=['load ' root directory '/' maillage '_fem_parameters'];
0786 eval(aa)
0787 u=u*yts;
0788 v=v*yts;
0789 trisurf(index,x/1e3,y/1e3,u),shading interp,view(2);
0790 xlabel('X (km)','fontname','arial','fontsize',16)
0791 ylabel('Y (km)','fontname','arial','fontsize',16)
0792 title(texlabel('Computed Vx (m.a^(-1)) around HIR.'),'fontname','arial','fontsize',18);
0793 set(gca,'fontname','arial','fontsize',16)
0794
0795
0796 zoom_choice(zo);
0797 colormap(jet)
0798 if range_choice==0,
0799 caxis([-500 500]);
0800 else
0801 eval(['caxis(' range ');']);
0802 end
0803
0804 colorbar
0805 set(gca,'Tag','ActiveAxes')
0806 end
0807
0808 else
0809 diff=findobj(gcf,'Tag','Diff');
0810 diff_value=popupstr(diff);
0811
0812 mesh=findobj(gcf,'Tag','Mesh_choice');
0813 mesh=popupstr(mesh)
0814 if strcmp(mesh,'Version')==1,
0815 mesh=1;
0816 else
0817 mesh=0;
0818 end
0819
0820
0821
0822 names_function(num,series)
0823 load ../data/matlab/commondata/names.mat
0824 load ../data/matlab/commondata/data
0825 exist=exist_by_names([root directory '\' results '_results.mat']);
0826
0827 if exist==1,
0828
0829 aa=['load ' root directory '/' results '_results'];
0830 eval(aa)
0831 aa=['load ' root directory '/' maillage '_datamesh'];
0832 eval(aa)
0833 aa=['load ' root directory '/' maillage '_fem_parameters'];
0834 eval(aa)
0835 aa=['load ' root directory '/' maillage '_observations'];
0836 eval(aa)
0837 vel1=u*yts;
0838
0839 x1=x;
0840 y1=y;
0841 index1=index;
0842 end
0843 names_function(diff_value,series)
0844 load ../data/matlab/commondata/names.mat
0845 load ../data/matlab/commondata/data
0846 exist=exist_by_names([root directory '\' results '_results.mat']);
0847
0848 if exist==1,
0849
0850 aa=['load ' root directory '/' results '_results'];
0851 eval(aa)
0852 aa=['load ' root directory '/' maillage '_datamesh'];
0853 eval(aa)
0854 aa=['load ' root directory '/' maillage '_fem_parameters'];
0855 eval(aa)
0856 aa=['load ' root directory '/' maillage '_observations'];
0857 eval(aa)
0858 vel2=u*yts;
0859 x2=x;
0860 y2=y;
0861 index2=index;
0862 end
0863
0864 if mesh==1,
0865 vel2=griddata(x2,y2,vel2,x1,y1);
0866 trisurf(index1,x1/1e3,y1/1e3,vel1-vel2),shading interp,view(2);
0867 title(texlabel(['Difference vx model #' num2str(num) '- vx model #' num2str(diff_value) ' on mesh ' num2str(num)]),'fontname','arial','fontsize',18);
0868
0869 else
0870 vel1=griddata(x1,y1,vel1,x2,y2);
0871 trisurf(index2,x2/1e3,y2/1e3,vel1-vel2),shading interp,view(2);
0872 title(texlabel(['Difference vx model #' num2str(num) '- vx model #' num2str(diff_value) ' on mesh ' num2str(diff_value)]),'fontname','arial','fontsize',18);
0873 end
0874
0875 if range_choice==0,
0876 else
0877 eval(['caxis(' range ');']);
0878 end
0879
0880
0881 colorbar;
0882 xlabel('X (km)','fontname','arial','fontsize',16)
0883 ylabel('Y (km)','fontname','arial','fontsize',16)
0884 set(gca,'fontname','arial','fontsize',16)
0885
0886 zoom_choice(zo);
0887 set(gca,'Tag','ActiveAxes')
0888
0889
0890
0891
0892
0893 end
0894
0895
0896
0897 function plot_diff_u_hir(num,zo,series,range_choice,range);
0898
0899 names_function(num,series)
0900 load ../data/matlab/commondata/names.mat
0901 load ../data/matlab/commondata/data
0902 exist=exist_by_names([root directory '\' results '_results.mat']);
0903
0904 if exist==1,
0905
0906 aa=['load ' root directory '/' results '_results'];
0907 eval(aa)
0908 aa=['load ' root directory '/' maillage '_datamesh'];
0909 eval(aa)
0910 aa=['load ' root directory '/' maillage '_fem_parameters'];
0911 eval(aa)
0912 aa=['load ' root directory '/' maillage '_observations'];
0913 eval(aa)
0914
0915 u=u*yts;
0916
0917 trisurf(index,x/1e3,y/1e3,u-u_HIR),shading interp,view(2);
0918 xlabel('X (km)','fontname','arial','fontsize',16)
0919 ylabel('Y (km)','fontname','arial','fontsize',16)
0920 title(texlabel('Vx-Vxobs (m.a^(-1)) around HIR.'),'fontname','arial','fontsize',18);
0921 set(gca,'fontname','arial','fontsize',16)
0922
0923 zoom_choice(zo);
0924 colormap(jet)
0925
0926 colorbar
0927 set(gca,'Tag','ActiveAxes')
0928 end
0929
0930
0931
0932
0933
0934 function plot_v(num,zo,series,range_choice,range);
0935 diff_choice=findobj(gcf,'Tag','Diff_choice');
0936 diff_choice=popupstr(diff_choice);
0937 if strcmp(diff_choice,'no')==1,
0938
0939
0940 names_function(num,series)
0941 load ../data/matlab/commondata/names.mat
0942 load ../data/matlab/commondata/data
0943 exist=exist_by_names([root directory '\' results '_results.mat']);
0944 if exist==1,
0945
0946 aa=['load ' root directory '/' results '_results'];
0947 eval(aa)
0948 aa=['load ' root directory '/' maillage '_datamesh'];
0949 eval(aa)
0950 aa=['load ' root directory '/' maillage '_fem_parameters'];
0951 eval(aa)
0952 u=u*yts;
0953 v=v*yts;
0954 trisurf(index,x/1e3,y/1e3,v),shading interp,view(2);
0955 xlabel('X (km)','fontname','arial','fontsize',16)
0956 ylabel('Y (km)','fontname','arial','fontsize',16)
0957 title(texlabel('Computed Vy (m.a^(-1)) around HIR.'),'fontname','arial','fontsize',18);
0958 set(gca,'fontname','arial','fontsize',16)
0959
0960 zoom_choice(zo);
0961
0962 colormap(jet)
0963 if range_choice==0,
0964 caxis([0 1500]);
0965 else
0966 eval(['caxis(' range ');']);
0967 end
0968
0969
0970 colorbar
0971 set(gca,'Tag','ActiveAxes')
0972 end
0973
0974 else
0975 diff=findobj(gcf,'Tag','Diff');
0976 diff_value=popupstr(diff);
0977
0978 mesh=findobj(gcf,'Tag','Mesh_choice');
0979 mesh=popupstr(mesh)
0980 if strcmp(mesh,'Version')==1,
0981 mesh=1;
0982 else
0983 mesh=0;
0984 end
0985
0986
0987
0988
0989 names_function(num,series)
0990 load ../data/matlab/commondata/names.mat
0991 load ../data/matlab/commondata/data
0992 exist=exist_by_names([root directory '\' results '_results.mat']);
0993
0994 if exist==1,
0995
0996 aa=['load ' root directory '/' results '_results'];
0997 eval(aa)
0998 aa=['load ' root directory '/' maillage '_datamesh'];
0999 eval(aa)
1000 aa=['load ' root directory '/' maillage '_fem_parameters'];
1001 eval(aa)
1002 aa=['load ' root directory '/' maillage '_observations'];
1003 eval(aa)
1004 vel1=v*yts;
1005
1006 x1=x;
1007 y1=y;
1008 index1=index;
1009 end
1010 names_function(diff_value,series)
1011 load ../data/matlab/commondata/names.mat
1012 load ../data/matlab/commondata/data
1013
1014 exist=exist_by_names([root directory '\' results '_results.mat']);
1015
1016 if exist==1,
1017
1018 aa=['load ' root directory '/' results '_results'];
1019 eval(aa)
1020 aa=['load ' root directory '/' maillage '_datamesh'];
1021 eval(aa)
1022 aa=['load ' root directory '/' maillage '_fem_parameters'];
1023 eval(aa)
1024 aa=['load ' root directory '/' maillage '_observations'];
1025 eval(aa)
1026 vel2=v*yts;
1027 x2=x;
1028 y2=y;
1029 index2=index;
1030 end
1031
1032 if mesh==1,
1033 vel2=griddata(x2,y2,vel2,x1,y1);
1034 trisurf(index1,x1/1e3,y1/1e3,vel1-vel2),shading interp,view(2);
1035 title(texlabel(['Difference vy model #' num2str(num) '- vx model #' num2str(diff_value) ' on mesh ' num2str(num)]),'fontname','arial','fontsize',18);
1036
1037 else
1038 vel1=griddata(x1,y1,vel1,x2,y2);
1039 trisurf(index2,x2/1e3,y2/1e3,vel1-vel2),shading interp,view(2);
1040 title(texlabel(['Difference vy model #' num2str(num) '- vx model #' num2str(diff_value) ' on mesh ' num2str(diff_value)]),'fontname','arial','fontsize',18);
1041 end
1042
1043 if range_choice==0,
1044 else
1045 eval(['caxis(' range ');']);
1046 end
1047
1048
1049 colorbar;
1050 xlabel('X (km)','fontname','arial','fontsize',16)
1051 ylabel('Y (km)','fontname','arial','fontsize',16)
1052 set(gca,'fontname','arial','fontsize',16)
1053
1054 zoom_choice(zo);
1055 set(gca,'Tag','ActiveAxes')
1056
1057
1058
1059
1060
1061
1062 end
1063
1064
1065
1066
1067 function plot_v_prism(num,zo,series,range_choice,range);
1068 names_function(num,series)
1069 load ../data/matlab/commondata/names.mat
1070 load ../data/matlab/commondata/data
1071 exist=exist_by_names([root directory '\' results '_results.mat']);
1072 if exist==1,
1073
1074 aa=['load ' root directory '/' results '_results'];
1075 eval(aa)
1076 aa=['load ' root directory '/' maillage '_datamesh'];
1077 eval(aa)
1078 aa=['load ' root directory '/' maillage '_fem_parameters'];
1079 eval(aa)
1080 v=v*yts;
1081
1082 trisurf(index,x/1e3,y/1e3,v),shading interp,view(2);
1083 xlabel('X (km)','fontname','arial','fontsize',16)
1084 ylabel('Y (km)','fontname','arial','fontsize',16)
1085 title(texlabel('Computed velocity magnitude (m.a^(-1)).'),'fontname','arial','fontsize',18);
1086 set(gca,'fontname','arial','fontsize',16)
1087 zoom_choice(zo);
1088 cycle=100;
1089 caxis([0 1600]);
1090 colormap(prism(6*cycle))
1091 colorbar
1092 set(gca,'Tag','ActiveAxes')
1093 end
1094
1095
1096
1097
1098 function plot_diff_v_hir(num,zo,series,range_choice,range);
1099
1100 names_function(num,series)
1101 load ../data/matlab/commondata/names.mat
1102 load ../data/matlab/commondata/data
1103 exist=exist_by_names([root directory '\' results '_results.mat']);
1104
1105 if exist==1,
1106
1107 aa=['load ' root directory '/' results '_results'];
1108 eval(aa)
1109 aa=['load ' root directory '/' maillage '_datamesh'];
1110 eval(aa)
1111 aa=['load ' root directory '/' maillage '_fem_parameters'];
1112 eval(aa)
1113 aa=['load ' root directory '/' maillage '_observations'];
1114 eval(aa)
1115
1116 v=v*yts;
1117
1118 trisurf(index,x/1e3,y/1e3,v-v_HIR),shading interp,view(2);
1119 xlabel('X (km)','fontname','arial','fontsize',16)
1120 ylabel('Y (km)','fontname','arial','fontsize',16)
1121 title(texlabel('Vy-Vyobs (m.a^(-1)) around HIR.'),'fontname','arial','fontsize',18);
1122 set(gca,'fontname','arial','fontsize',16)
1123
1124 zoom_choice(zo);
1125
1126 colormap(jet)
1127 caxis([-500 500]);
1128 colorbar
1129 set(gca,'Tag','ActiveAxes')
1130 end
1131
1132
1133
1134
1135
1136
1137
1138
1139 function plot_vel(num,zo,series,range_choice,range);
1140
1141 diff_choice=findobj(gcf,'Tag','Diff_choice');
1142 diff_choice=popupstr(diff_choice);
1143 if strcmp(diff_choice,'no')==1,
1144
1145 names_function(num,series)
1146 load ../data/matlab/commondata/names.mat
1147 load ../data/matlab/commondata/data
1148 exist=exist_by_names([root directory '\' results '_results.mat']);
1149 if exist==1,
1150
1151 aa=['load ' root directory '/' results '_results'];
1152 eval(aa)
1153 aa=['load ' root directory '/' maillage '_datamesh'];
1154 eval(aa)
1155 aa=['load ' root directory '/' maillage '_fem_parameters'];
1156 eval(aa)
1157 u=u*yts;
1158 v=v*yts;
1159 vel=sqrt(u.^2+v.^2);
1160 trisurf(index,x/1e3,y/1e3,vel),shading interp,view(2);
1161 xlabel('X (km)','fontname','arial','fontsize',16)
1162 ylabel('Y (km)','fontname','arial','fontsize',16)
1163 title(texlabel('Computed velocity magnitude (m.a^(-1)).'),'fontname','arial','fontsize',18);
1164 set(gca,'fontname','arial','fontsize',16)
1165 zoom_choice(zo);
1166 cycle=30;
1167
1168 if range_choice==0,
1169 caxis([0 1500]);
1170
1171 else
1172 eval(['caxis(' range ');']);
1173 end
1174
1175 colormap(jet)
1176 colorbar
1177 set(gca,'Tag','ActiveAxes')
1178
1179 end
1180
1181
1182 else
1183
1184 diff=findobj(gcf,'Tag','Diff');
1185 diff_value=popupstr(diff);
1186
1187
1188 mesh=findobj(gcf,'Tag','Mesh_choice');
1189 mesh=popupstr(mesh)
1190 if strcmp(mesh,'Version')==1,
1191 mesh=1;
1192 else
1193 mesh=0;
1194 end
1195
1196
1197
1198
1199 names_function(num,series)
1200 load ../data/matlab/commondata/names.mat
1201 load ../data/matlab/commondata/data
1202 exist=exist_by_names([root directory '\' results '_results.mat']);
1203
1204 if exist==1,
1205
1206 aa=['load ' root directory '/' results '_results'];
1207 eval(aa)
1208 aa=['load ' root directory '/' maillage '_datamesh'];
1209 eval(aa)
1210 aa=['load ' root directory '/' maillage '_fem_parameters'];
1211 eval(aa)
1212 aa=['load ' root directory '/' maillage '_observations'];
1213 eval(aa)
1214 vel1=sqrt(u.^2+v.^2)*yts;
1215
1216 x1=x;
1217 y1=y;
1218 index1=index;
1219 end
1220 names_function(diff_value,series)
1221 load ../data/matlab/commondata/names.mat
1222
1223 load ../data/matlab/commondata/data
1224 exist=exist_by_names([root directory '\' results '_results.mat']);
1225
1226 if exist==1,
1227
1228 aa=['load ' root directory '/' results '_results'];
1229 eval(aa)
1230 aa=['load ' root directory '/' maillage '_datamesh'];
1231 eval(aa)
1232 aa=['load ' root directory '/' maillage '_fem_parameters'];
1233 eval(aa)
1234 aa=['load ' root directory '/' maillage '_observations'];
1235 eval(aa)
1236 vel2=sqrt(u.^2+v.^2)*yts;
1237 x2=x;
1238 y2=y;
1239 index2=index;
1240 end
1241
1242
1243 if mesh==1,
1244 vel2=griddata(x2,y2,vel2,x1,y1);
1245 trisurf(index1,x1/1e3,y1/1e3,vel1-vel2),shading interp,view(2);
1246 title(texlabel(['Difference vel model #' num2str(num) '- vx model #' num2str(diff_value) ' on mesh ' num2str(num)]),'fontname','arial','fontsize',18);
1247
1248 else
1249 vel1=griddata(x1,y1,vel1,x2,y2);
1250 trisurf(index2,x2/1e3,y2/1e3,vel1-vel2),shading interp,view(2);
1251 title(texlabel(['Difference vel model #' num2str(num) '- vx model #' num2str(diff_value) ' on mesh ' num2str(diff_value)]),'fontname','arial','fontsize',18);
1252 end
1253
1254
1255 if range_choice==0,
1256 else
1257 eval(['caxis(' range ');']);
1258 end
1259
1260 colorbar;
1261 xlabel('X (km)','fontname','arial','fontsize',16)
1262 ylabel('Y (km)','fontname','arial','fontsize',16)
1263 set(gca,'fontname','arial','fontsize',16)
1264
1265 zoom_choice(zo);
1266 set(gca,'Tag','ActiveAxes')
1267
1268
1269
1270
1271
1272
1273
1274
1275 end
1276
1277
1278
1279
1280 function plot_vel_prism(num,zo,series,range_choice,range);
1281 names_function(num,series)
1282 load ../data/matlab/commondata/names.mat
1283 load ../data/matlab/commondata/data
1284 exist=exist_by_names([root directory '\' results '_results.mat']);
1285 if exist==1,
1286
1287 aa=['load ' root directory '/' results '_results'];
1288 eval(aa)
1289 aa=['load ' root directory '/' maillage '_datamesh'];
1290 eval(aa)
1291 aa=['load ' root directory '/' maillage '_fem_parameters'];
1292 eval(aa)
1293 u=u*yts;
1294 v=v*yts;
1295 vel=sqrt(u.^2+v.^2);
1296
1297 trisurf(index,x/1e3,y/1e3,vel),shading interp,view(2);
1298 xlabel('X (km)','fontname','arial','fontsize',16)
1299 ylabel('Y (km)','fontname','arial','fontsize',16)
1300 title(texlabel('Computed velocity magnitude (m.a^(-1)).'),'fontname','arial','fontsize',18);
1301 set(gca,'fontname','arial','fontsize',16)
1302 zoom_choice(zo);
1303 cycle=32;
1304 caxis([0 1600]);
1305 colormap(prism(6*cycle))
1306 colorbar
1307 set(gca,'Tag','ActiveAxes')
1308 end
1309
1310
1311
1312
1313
1314 function plot_diff_vel_hir(num,zo,series,range_choice,range);
1315
1316 names_function(num,series)
1317 load ../data/matlab/commondata/names.mat
1318 load ../data/matlab/commondata/data
1319 exist=exist_by_names([root directory '\' results '_results.mat']);
1320
1321 if exist==1,
1322
1323 aa=['load ' root directory '/' results '_results'];
1324 eval(aa)
1325 aa=['load ' root directory '/' maillage '_datamesh'];
1326 eval(aa)
1327 aa=['load ' root directory '/' maillage '_fem_parameters'];
1328 eval(aa)
1329 aa=['load ' root directory '/' maillage '_observations'];
1330 eval(aa)
1331
1332 v=v*yts;
1333 u=u*yts;
1334 vel=sqrt(u.^2+v.^2);
1335 vel_HIR=sqrt(u_HIR.^2+v_HIR.^2);
1336
1337 xmin=-1.09*10^6;
1338 xmax=-960000;
1339 ymin=8.1*10^5;
1340 ymax=9.4*10^5;
1341 diff=vel-vel_HIR;
1342 diff2=abs(vel-vel_HIR);
1343 pos=find(diff2<=50);
1344 diff(pos)=NaN;
1345 trisurf(index,x/1e3,y/1e3,diff),shading interp,view(2);
1346 xlabel('X (km)','fontname','arial','fontsize',16)
1347 ylabel('Y (km)','fontname','arial','fontsize',16)
1348 title(texlabel('V-Vobs (m.a^(-1)) around HIR.'),'fontname','arial','fontsize',18);
1349 set(gca,'fontname','arial','fontsize',16)
1350
1351 zoom_choice(zo);
1352 colormap(jet)
1353 caxis([-500 500]);
1354 colorbar
1355 set(gca,'Tag','ActiveAxes')
1356 end
1357
1358
1359
1360
1361
1362
1363 function plot_diff_vel_ronne(num,zo,series,range_choice,range);
1364
1365 names_function(num,series)
1366 load ../data/matlab/commondata/names.mat
1367 load ../data/matlab/commondata/data
1368 exist=exist_by_names([root directory '\' results '_results.mat']);
1369
1370 if exist==1,
1371
1372 aa=['load ' root directory '/' results '_results'];
1373 eval(aa)
1374 aa=['load ' root directory '/' maillage '_datamesh'];
1375 eval(aa)
1376 aa=['load ' root directory '/' maillage '_fem_parameters'];
1377 eval(aa)
1378 aa=['load ' root directory '/' maillage '_observations'];
1379 eval(aa)
1380
1381 if exist_by_names([root directory '\' maillage '_observations_Ian.mat']),
1382 eval(['load ' root directory '/' maillage '_observations_Ian'])
1383 v_input=sqrt(u_RONNE.^2+v_RONNE.^2)*yts;
1384 end
1385
1386
1387 v=v*yts;
1388 u=u*yts;
1389 vel=sqrt(u.^2+v.^2);
1390
1391 trisurf(index,x/1e3,y/1e3,vel-v_input),shading interp,view(2);
1392 xlabel('X (km)','fontname','arial','fontsize',16)
1393 ylabel('Y (km)','fontname','arial','fontsize',16)
1394 title(texlabel('V-Vobs (m.a^(-1)).'),'fontname','arial','fontsize',18);
1395 set(gca,'fontname','arial','fontsize',16)
1396
1397 zoom_choice(zo);
1398
1399 colormap(jet)
1400 caxis([-500 500]);
1401 colorbar
1402 set(gca,'Tag','ActiveAxes')
1403 end
1404
1405
1406
1407
1408
1409 function plot_obs_LAS_vx(num,zo,series,range_choice,range);
1410
1411 names_function(num,series)
1412 load ../data/matlab/commondata/names.mat
1413 load ../data/matlab/commondata/data
1414 exist=exist_by_names([root directory '\' results '_results.mat']);
1415
1416 if exist==1,
1417
1418 aa=['load ' root directory '/' results '_results'];
1419 eval(aa)
1420 aa=['load ' root directory '/' maillage '_datamesh'];
1421 eval(aa)
1422 aa=['load ' root directory '/' maillage '_fem_parameters'];
1423 eval(aa)
1424 aa=['load ' root directory '/' maillage '_observations_LAS'];
1425 eval(aa)
1426
1427
1428 trisurf(index,x/1e3,y/1e3,u_LAS),shading interp,view(2);
1429 xlabel('X (km)','fontname','arial','fontsize',16)
1430 ylabel('Y (km)','fontname','arial','fontsize',16)
1431 title(texlabel('Observed vx (m.a^(-1)) around LAS. From RSAT1'),'fontname','arial','fontsize',18);
1432 set(gca,'fontname','arial','fontsize',16)
1433
1434 zoom_choice(zo);
1435
1436 colormap(jet)
1437
1438 colorbar
1439 set(gca,'Tag','ActiveAxes')
1440 end
1441
1442
1443
1444 function plot_obs_LAS_vy(num,zo,series,range_choice,range);
1445
1446 names_function(num,series)
1447 load ../data/matlab/commondata/names.mat
1448 load ../data/matlab/commondata/data
1449 exist=exist_by_names([root directory '\' results '_results.mat']);
1450
1451 if exist==1,
1452
1453 aa=['load ' root directory '/' results '_results'];
1454 eval(aa)
1455 aa=['load ' root directory '/' maillage '_datamesh'];
1456 eval(aa)
1457 aa=['load ' root directory '/' maillage '_fem_parameters'];
1458 eval(aa)
1459 aa=['load ' root directory '/' maillage '_observations_LAS'];
1460 eval(aa)
1461
1462
1463 trisurf(index,x/1e3,y/1e3,v_LAS),shading interp,view(2);
1464 xlabel('X (km)','fontname','arial','fontsize',16)
1465 ylabel('Y (km)','fontname','arial','fontsize',16)
1466 title(texlabel('Observed vy around LASSITER (m.a^(-1)) around LAS. From RSAT1'),'fontname','arial','fontsize',18);
1467 set(gca,'fontname','arial','fontsize',16)
1468
1469 zoom_choice(zo);
1470
1471 colormap(jet)
1472
1473 colorbar
1474 set(gca,'Tag','ActiveAxes')
1475 end
1476
1477
1478
1479 function plot_obs_LAS_vel(num,zo,series,range_choice,range);
1480
1481 names_function(num,series)
1482 load ../data/matlab/commondata/names.mat
1483 load ../data/matlab/commondata/data
1484 exist=exist_by_names([root directory '\' results '_results.mat']);
1485
1486 if exist==1,
1487
1488 aa=['load ' root directory '/' results '_results'];
1489 eval(aa)
1490 aa=['load ' root directory '/' maillage '_datamesh'];
1491 eval(aa)
1492 aa=['load ' root directory '/' maillage '_fem_parameters'];
1493 eval(aa)
1494 aa=['load ' root directory '/' maillage '_observations_LAS'];
1495 eval(aa)
1496 vel_LAS=sqrt(u_LAS.^2+v_LAS.^2);
1497
1498 trisurf(index,x/1e3,y/1e3,vel_LAS),shading interp,view(2);
1499 xlabel('X (km)','fontname','arial','fontsize',16)
1500 ylabel('Y (km)','fontname','arial','fontsize',16)
1501 title(texlabel('Observed vel_LAS (m.a^(-1)) around LAS. From RSAT1'),'fontname','arial','fontsize',18);
1502 set(gca,'fontname','arial','fontsize',16)
1503
1504 zoom_choice(zo);
1505
1506 colormap(jet)
1507
1508 colorbar
1509 set(gca,'Tag','ActiveAxes')
1510 end
1511
1512
1513
1514 function plot_diff_u_LAS(num,zo,series,range_choice,range);
1515
1516 names_function(num,series)
1517 load ../data/matlab/commondata/names.mat
1518 load ../data/matlab/commondata/data
1519 exist=exist_by_names([root directory '\' results '_results.mat']);
1520
1521 if exist==1,
1522
1523 aa=['load ' root directory '/' results '_results'];
1524 eval(aa)
1525 aa=['load ' root directory '/' maillage '_datamesh'];
1526 eval(aa)
1527 aa=['load ' root directory '/' maillage '_fem_parameters'];
1528 eval(aa)
1529 aa=['load ' root directory '/' maillage '_observations_LAS'];
1530 eval(aa)
1531
1532 u=u*yts;
1533 pos=find(u_LAS==0);
1534 u_LAS(pos)=NaN;
1535 xmin=-1.09*10^6;
1536 xmax=-960000;
1537 ymin=8.1*10^5;
1538 ymax=9.4*10^5;
1539 diff=u-u_LAS;
1540 trisurf(index,x/1e3,y/1e3,diff),shading interp,view(2);
1541 xlabel('X (km)','fontname','arial','fontsize',16)
1542 ylabel('Y (km)','fontname','arial','fontsize',16)
1543 title(texlabel('u-u_LAS (m.a^(-1)) around LASSITER.'),'fontname','arial','fontsize',18);
1544 set(gca,'fontname','arial','fontsize',16)
1545
1546 zoom_choice(zo);
1547 colormap(jet)
1548 caxis([-500 500]);
1549 colorbar
1550 set(gca,'Tag','ActiveAxes')
1551 end
1552
1553
1554 function plot_diff_v_LAS(num,zo,series,range_choice,range);
1555
1556 names_function(num,series)
1557 load ../data/matlab/commondata/names.mat
1558 load ../data/matlab/commondata/data
1559 exist=exist_by_names([root directory '\' results '_results.mat']);
1560
1561 if exist==1,
1562
1563 aa=['load ' root directory '/' results '_results'];
1564 eval(aa)
1565 aa=['load ' root directory '/' maillage '_datamesh'];
1566 eval(aa)
1567 aa=['load ' root directory '/' maillage '_fem_parameters'];
1568 eval(aa)
1569 aa=['load ' root directory '/' maillage '_observations_LAS'];
1570 eval(aa)
1571
1572 v=v*yts;
1573 pos=find(v_LAS==0);
1574 v_LAS(pos)=NaN;
1575 xmin=-1.09*10^6;
1576 xmax=-960000;
1577 ymin=8.1*10^5;
1578 ymax=9.4*10^5;
1579 diff=v-v_LAS;
1580 trisurf(index,x/1e3,y/1e3,diff),shading interp,view(2);
1581 xlabel('X (km)','fontname','arial','fontsize',16)
1582 ylabel('Y (km)','fontname','arial','fontsize',16)
1583 title(texlabel('v-v_LAS (m.a^(-1)) around LASSITER.'),'fontname','arial','fontsize',18);
1584 set(gca,'fontname','arial','fontsize',16)
1585
1586 zoom_choice(zo);
1587 colormap(jet)
1588 caxis([-500 500]);
1589 colorbar
1590 set(gca,'Tag','ActiveAxes')
1591 end
1592
1593 function plot_diff_vel_LAS(num,zo,series,range_choice,range);
1594
1595 names_function(num,series)
1596 load ../data/matlab/commondata/names.mat
1597 load ../data/matlab/commondata/data
1598 exist=exist_by_names([root directory '\' results '_results.mat']);
1599
1600 if exist==1,
1601
1602 aa=['load ' root directory '/' results '_results'];
1603 eval(aa)
1604 aa=['load ' root directory '/' maillage '_datamesh'];
1605 eval(aa)
1606 aa=['load ' root directory '/' maillage '_fem_parameters'];
1607 eval(aa)
1608 aa=['load ' root directory '/' maillage '_observations_LAS'];
1609 eval(aa)
1610
1611 u=u*yts;
1612 v=v*yts;
1613 vel=sqrt(u.^2+v.^2);
1614 vel_LAS=sqrt(u_LAS.^2+v_LAS.^2);
1615
1616
1617
1618
1619 diff2=abs(vel-vel_LAS);
1620 pos=find(diff2<=50);
1621
1622 xmin=-1.09*10^6;
1623 xmax=-960000;
1624 ymin=8.1*10^5;
1625 ymax=9.4*10^5;
1626 diff=vel-vel_LAS;
1627 diff(pos)=NaN;
1628 trisurf(index,x/1e3,y/1e3,diff),shading interp,view(2);
1629 xlabel('X (km)','fontname','arial','fontsize',16)
1630 ylabel('Y (km)','fontname','arial','fontsize',16)
1631 title(texlabel('vel-vel_LAS (m.a^(-1)) around LASSITER.'),'fontname','arial','fontsize',18);
1632 set(gca,'fontname','arial','fontsize',16)
1633
1634 zoom_choice(zo);
1635 colormap(jet)
1636 caxis([-500 500]);
1637 colorbar
1638 set(gca,'Tag','ActiveAxes')
1639 end
1640
1641
1642
1643
1644
1645 function plot_tau1(num,zo,series,range_choice,range);
1646 names_function(num,series)
1647 load ../data/matlab/commondata/names.mat
1648 load ../data/matlab/commondata/data
1649 exist=exist_by_names([root directory '\' results '_results.mat']);
1650 if exist==1,
1651
1652
1653 aa=['load ' root directory '/' results '_results'];
1654 eval(aa)
1655 aa=['load ' root directory '/' maillage '_datamesh'];
1656 eval(aa)
1657 aa=['load ' root directory '/' maillage '_fem_parameters'];
1658 eval(aa)
1659 [tau_1 tau_2 teta1 teta2]=tau_principal(tau_xx,tau_yy,tau_xy);
1660 summer=[1;1;1];
1661 colormap(jet)
1662 cla
1663 patch(x(index)'/1e3,y(index)'/1e3,[1; 1; 1]*tau_1'),shading interp,view(2)
1664 xlabel('X (km)','fontname','arial','fontsize',16)
1665 ylabel('Y (km)','fontname','arial','fontsize',16)
1666 title(texlabel('First component of stress (10^5 Pa.s).'),'fontname','arial','fontsize',18);
1667 set(gca,'fontname','arial','fontsize',16)
1668 caxis([-6 6]);
1669 colorbar
1670 set(gca,'Tag','ActiveAxes')
1671
1672 end
1673
1674
1675
1676
1677
1678 function plot_tau2(num,zo,series,range_choice,range);
1679 names_function(num,series)
1680 load ../data/matlab/commondata/names.mat
1681 load ../data/matlab/commondata/data
1682 exist=exist_by_names([root directory '\' results '_results.mat']);
1683 if exist==1,
1684
1685 aa=['load ' root directory '/' results '_results'];
1686 eval(aa)
1687 aa=['load ' root directory '/' maillage '_datamesh'];
1688 eval(aa)
1689 aa=['load ' root directory '/' maillage '_fem_parameters'];
1690 eval(aa)
1691 [tau_1 tau_2 teta1 teta2]=tau_principal(tau_xx,tau_yy,tau_xy);
1692 colormap(jet)
1693 cla
1694 patch(x(index)'/1e3,y(index)'/1e3,[1; 1; 1]*tau_2'),shading interp,view(2)
1695 xlabel('X (km)','fontname','arial','fontsize',16)
1696 ylabel('Y (km)','fontname','arial','fontsize',16)
1697 title(texlabel('Second component of stress (10^5 Pa.s).'),'fontname','arial','fontsize',18);
1698 set(gca,'fontname','arial','fontsize',16)
1699 caxis([-6 6]);
1700 colorbar
1701 set(gca,'Tag','ActiveAxes')
1702 end
1703
1704
1705
1706
1707
1708
1709
1710 function plot_teta1(num,zo,series,range_choice,range);
1711 names_function(num,series)
1712 load ../data/matlab/commondata/names.mat
1713 load ../data/matlab/commondata/data
1714 exist=exist_by_names([root directory '\' results '_results.mat']);
1715 if exist==1,
1716
1717 aa=['load ' root directory '/' results '_results'];
1718
1719 eval(aa)
1720 aa=['load ' root directory '/' maillage '_datamesh'];
1721 eval(aa)
1722 aa=['load ' root directory '/' maillage '_fem_parameters'];
1723 eval(aa)
1724 [tau_1 tau_2 teta1 teta2]=tau_principal(tau_xx,tau_yy,tau_xy);
1725 colormap(jet)
1726 cla
1727 patch(x(index)'/1e3,y(index)'/1e3,[1; 1; 1]*teta1'),shading interp,view(2)
1728 xlabel('X (km)','fontname','arial','fontsize',16)
1729 ylabel('Y (km)','fontname','arial','fontsize',16)
1730 title(texlabel('Angle of first component of stress (Pa.s).'),'fontname','arial','fontsize',18);
1731 set(gca,'fontname','arial','fontsize',16)
1732
1733 caxis([-pi pi])
1734 colorbar
1735 set(gca,'Tag','ActiveAxes')
1736 end
1737
1738
1739
1740
1741
1742
1743
1744 function plot_teta2(num,zo,series,range_choice,range);
1745 names_function(num,series)
1746 load ../data/matlab/commondata/names.mat
1747 load ../data/matlab/commondata/data
1748 exist=exist_by_names([root directory '\' results '_results.mat']);
1749 if exist==1,
1750
1751 aa=['load ' root directory '/' results '_results'];
1752 eval(aa)
1753 aa=['load ' root directory '/' maillage '_datamesh'];
1754 eval(aa)
1755 aa=['load ' root directory '/' maillage '_fem_parameters'];
1756 eval(aa)
1757 [tau_1 tau_2 teta1 teta2]=tau_principal(tau_xx,tau_yy,tau_xy);
1758 colormap(jet)
1759 cla
1760 patch(x(index)'/1e3,y(index)'/1e3,[1; 1; 1]*teta2'),shading interp,view(2)
1761 xlabel('X (km)','fontname','arial','fontsize',16)
1762 ylabel('Y (km)','fontname','arial','fontsize',16)
1763 title(texlabel('Angle of second component of stress (Pa.s).'),'fontname','arial','fontsize',18);
1764 set(gca,'fontname','arial','fontsize',16)
1765
1766 caxis([-pi pi])
1767 colorbar
1768 set(gca,'Tag','ActiveAxes')
1769 end
1770
1771
1772
1773 function plot_thickbar(num,zo,series,range_choice,range);
1774 names_function(num,series)
1775 load ../data/matlab/commondata/names.mat
1776 load ../data/matlab/commondata/data
1777 exist=exist_by_names([root directory '\' results '_results.mat']);
1778 if exist==1,
1779
1780 aa=['load ' root directory '/' results '_results'];
1781 eval(aa)
1782 aa=['load ' root directory '/' maillage '_datamesh'];
1783 eval(aa)
1784 aa=['load ' root directory '/' maillage '_fem_parameters'];
1785 eval(aa)
1786 colormap(jet)
1787 cla
1788 patch(x(index)',y(index)',[1; 1; 1]*z_thick_bar'),shading interp,view(2)
1789 xlabel('X (km)','fontname','arial','fontsize',16)
1790 ylabel('Y (km)','fontname','arial','fontsize',16)
1791 title(texlabel('Thickness interpolated on the elements (m).'),'fontname','arial','fontsize',18);
1792 set(gca,'fontname','arial','fontsize',16)
1793
1794 if range_choice==0,
1795 caxis([0 500]);
1796 else
1797 eval(['caxis(' range ');']);
1798 end
1799
1800
1801
1802 colorbar
1803 set(gca,'Tag','ActiveAxes')
1804 end
1805
1806
1807
1808
1809
1810 function plot_thick(num,zo,series,range_choice,range);
1811
1812 diff_choice=findobj(gcf,'Tag','Diff_choice');
1813 diff_choice=popupstr(diff_choice);
1814 if strcmp(diff_choice,'no')==1,
1815
1816
1817 names_function(num,series)
1818 load ../data/matlab/commondata/names.mat
1819 load ../data/matlab/commondata/data
1820 exist=exist_by_names([root directory '\' results '_results.mat']);
1821
1822 if exist==1,
1823
1824 aa=['load ' root directory '/' results '_results'];
1825 eval(aa)
1826 aa=['load ' root directory '/' maillage '_datamesh'];
1827 eval(aa)
1828 aa=['load ' root directory '/' maillage '_fem_parameters'];
1829 eval(aa)
1830
1831 xmin=-1.09*10^6;
1832 xmax=-960000;
1833 ymin=8.1*10^5;
1834 ymax=9.4*10^5;
1835 trisurf(index,x/1e3,y/1e3,z_thick),shading interp,view(2);
1836 xlabel('X (km)','fontname','arial','fontsize',16)
1837 ylabel('Y (km)','fontname','arial','fontsize',16)
1838 title(texlabel('Thickness (m)'),'fontname','arial','fontsize',18);
1839 set(gca,'fontname','arial','fontsize',16)
1840
1841 zoom_choice(zo);
1842 colormap(jet)
1843 if range_choice==0,
1844 caxis([0 500]);
1845 else
1846 eval(['caxis(' range ');']);
1847 end
1848
1849 colorbar
1850 set(gca,'Tag','ActiveAxes')
1851 end
1852
1853
1854 else
1855
1856 diff=findobj(gcf,'Tag','Diff');
1857
1858 diff_value=popupstr(diff);
1859
1860
1861 mesh=findobj(gcf,'Tag','Mesh_choice');
1862 mesh=popupstr(mesh)
1863 if strcmp(mesh,'Version')==1,
1864 mesh=1;
1865 else
1866 mesh=0;
1867 end
1868
1869
1870
1871
1872 names_function(num,series)
1873 load ../data/matlab/commondata/names.mat
1874 load ../data/matlab/commondata/data
1875 exist=exist_by_names([root directory '\' results '_results.mat']);
1876
1877 if exist==1,
1878
1879 aa=['load ' root directory '/' results '_results'];
1880 eval(aa)
1881 aa=['load ' root directory '/' maillage '_datamesh'];
1882 eval(aa)
1883 aa=['load ' root directory '/' maillage '_fem_parameters'];
1884 eval(aa)
1885 aa=['load ' root directory '/' maillage '_observations'];
1886 eval(aa)
1887 z_thick1=z_thick;
1888 x1=x;
1889 y1=y;
1890 index1=index;
1891 end
1892 names_function(diff_value,series)
1893 load ../data/matlab/commondata/names.mat
1894 load ../data/matlab/commondata/data
1895 exist=exist_by_names([root directory '\' results '_results.mat']);
1896
1897 if exist==1,
1898
1899 aa=['load ' root directory '/' results '_results'];
1900 eval(aa)
1901 aa=['load ' root directory '/' maillage '_datamesh'];
1902 eval(aa)
1903 aa=['load ' root directory '/' maillage '_fem_parameters'];
1904 eval(aa)
1905 aa=['load ' root directory '/' maillage '_observations'];
1906 eval(aa)
1907 z_thick2=z_thick;
1908 x2=x;
1909 y2=y;
1910 index2=index;
1911 end
1912
1913
1914 if mesh==1,
1915 z_thick2=griddata(x2,y2,z_thick2,x1,y1);
1916 trisurf(index1,x1/1e3,y1/1e3,z_thick1-z_thick2),shading interp,view(2);
1917 title(texlabel(['Difference z_thick model #' num2str(num) '- z_thick model #' num2str(diff_value) ' on mesh ' num2str(num)]),'fontname','arial','fontsize',18);
1918
1919 else
1920 z_thick1=griddata(x1,y1,z_thick1,x2,y2);
1921 trisurf(index2,x2/1e3,y2/1e3,z_thick1-z_thick2),shading interp,view(2);
1922 title(texlabel(['Difference z_thick model #' num2str(num) '- z_thick model #' num2str(diff_value) ' on mesh ' num2str(diff_value)]),'fontname','arial','fontsize',18);
1923 end
1924
1925
1926 if range_choice==0,
1927 else
1928 eval(['caxis(' range ');']);
1929 end
1930
1931 colorbar;
1932 xlabel('X (km)','fontname','arial','fontsize',16)
1933 ylabel('Y (km)','fontname','arial','fontsize',16)
1934 set(gca,'fontname','arial','fontsize',16)
1935
1936 zoom_choice(zo);
1937 set(gca,'Tag','ActiveAxes')
1938
1939 end
1940
1941
1942
1943
1944 function plot_rotation_hir(num,zo,series,range_choice,range);
1945
1946 names_function(num,series)
1947 load ../data/matlab/commondata/names.mat
1948 load ../data/matlab/commondata/data
1949 exist=exist_by_names([root directory '\' results '_results.mat']);
1950
1951 if exist==1,
1952
1953 aa=['load ' root directory '/' results '_results'];
1954 eval(aa)
1955 aa=['load ' root directory '/' maillage '_datamesh'];
1956 eval(aa)
1957 aa=['load ' root directory '/' maillage '_fem_parameters'];
1958 eval(aa)
1959 aa=['load ' root directory '/' maillage '_observations'];
1960 eval(aa)
1961 angle=atan2(v_HIR,u_HIR)/pi*180;
1962 trisurf(index,x/1e3,y/1e3,angle),shading interp,view(2);
1963 xlabel('X (km)','fontname','arial','fontsize',16)
1964 ylabel('Y (km)','fontname','arial','fontsize',16)
1965 title(texlabel('Observed velocity Vx (m.a^(-1)) around HIR. From RSAT1'),'fontname','arial','fontsize',18);
1966 set(gca,'fontname','arial','fontsize',16)
1967
1968
1969 zoom_choice(zo);
1970
1971 colormap(jet)
1972
1973 colorbar
1974 set(gca,'Tag','ActiveAxes')
1975 end
1976
1977
1978
1979
1980
1981 function plot_quiver_hir(num,zo,series,range_choice,range);
1982
1983 names_function(num,series)
1984 load ../data/matlab/commondata/names.mat
1985 load ../data/matlab/commondata/data
1986 exist=exist_by_names([root directory '\' results '_results.mat']);
1987
1988 if exist==1,
1989
1990 aa=['load ' root directory '/' results '_results'];
1991 eval(aa)
1992 aa=['load ' root directory '/' maillage '_datamesh'];
1993 eval(aa)
1994 aa=['load ' root directory '/' maillage '_fem_parameters'];
1995 eval(aa)
1996 aa=['load ' root directory '/' maillage '_observations'];
1997 eval(aa)
1998
1999 delta=findobj(gcf,'Tag','Quiver_Sparsity');
2000 delta=str2num(deblank(get(delta,'String')));
2001
2002
2003 flag=1:delta:nods;
2004 quiver(x(flag)/1e3,y(flag)/1e3,u_HIR(flag),v_HIR(flag));
2005
2006
2007
2008
2009 xlabel('X (km)','fontname','arial','fontsize',16)
2010 ylabel('Y (km)','fontname','arial','fontsize',16)
2011 title(texlabel('Observed velocity Vx (m.a^(-1)) around HIR. From RSAT1'),'fontname','arial','fontsize',18);
2012 set(gca,'fontname','arial','fontsize',16)
2013 zoom_choice(zo);
2014 colormap(jet)
2015 set(gca,'Tag','ActiveAxes')
2016 end
2017
2018
2019
2020 function plot_quiver_las(num,zo,series,range_choice,range);
2021
2022 names_function(num,series)
2023 load ../data/matlab/commondata/names.mat
2024 load ../data/matlab/commondata/data
2025 exist=exist_by_names([root directory '\' results '_results.mat']);
2026
2027
2028 if exist==1,
2029
2030 aa=['load ' root directory '/' results '_results'];
2031 eval(aa)
2032 aa=['load ' root directory '/' maillage '_datamesh'];
2033 eval(aa)
2034 aa=['load ' root directory '/' maillage '_fem_parameters'];
2035 eval(aa)
2036 eval(['load ' root directory '/' maillage '_observations_LAS']);
2037
2038
2039
2040
2041 delta=findobj(gcf,'Tag','Quiver_Sparsity');
2042 delta=str2num(deblank(get(delta,'String')));
2043
2044
2045 flag=1:delta:nods;
2046 quiver(x(flag)/1e3,y(flag)/1e3,u_LAS(flag),v_LAS(flag));
2047
2048 xlabel('X (km)','fontname','arial','fontsize',16)
2049 ylabel('Y (km)','fontname','arial','fontsize',16)
2050 title(texlabel('Observed velocity Vx (m.a^(-1)) around LASSITER. From RSAT1'),'fontname','arial','fontsize',18);
2051 set(gca,'fontname','arial','fontsize',16)
2052
2053
2054 zoom_choice(zo);
2055 colormap(jet)
2056 colorbar
2057 set(gca,'Tag','ActiveAxes')
2058 end
2059
2060
2061
2062
2063
2064
2065
2066 function plot_quiver_model(num,zo,series,range_choice,range);
2067
2068 names_function(num,series)
2069 load ../data/matlab/commondata/names.mat
2070 load ../data/matlab/commondata/data
2071 exist=exist_by_names([root directory '\' results '_results.mat']);
2072
2073 if exist==1,
2074
2075 aa=['load ' root directory '/' results '_results'];
2076 eval(aa)
2077 aa=['load ' root directory '/' maillage '_datamesh'];
2078 eval(aa)
2079 aa=['load ' root directory '/' maillage '_fem_parameters'];
2080 eval(aa)
2081
2082
2083
2084 delta=findobj(gcf,'Tag','Quiver_Sparsity');
2085 delta=str2num(deblank(get(delta,'String')));
2086
2087
2088 flag=1:delta:nods;
2089 quiver(x(flag)/1e3,y(flag)/1e3,u(flag),v(flag));
2090
2091 xlabel('X (km)','fontname','arial','fontsize',16)
2092 ylabel('Y (km)','fontname','arial','fontsize',16)
2093 title(texlabel('Observed velocity Vx (m.a^(-1)) around HIR. From RSAT1'),'fontname','arial','fontsize',18);
2094 set(gca,'fontname','arial','fontsize',16)
2095
2096
2097 zoom_choice(zo);
2098 colormap(jet)
2099 colorbar
2100 set(gca,'Tag','ActiveAxes')
2101 end
2102
2103
2104
2105
2106
2107
2108
2109
2110 function plot_rotation_model(num,zo,series,range_choice,range);
2111
2112 names_function(num,series)
2113 load ../data/matlab/commondata/names.mat
2114 load ../data/matlab/commondata/data
2115 exist=exist_by_names([root directory '\' results '_results.mat']);
2116
2117 if exist==1,
2118
2119 aa=['load ' root directory '/' results '_results'];
2120 eval(aa)
2121 aa=['load ' root directory '/' maillage '_datamesh'];
2122 eval(aa)
2123 aa=['load ' root directory '/' maillage '_fem_parameters'];
2124 eval(aa)
2125 angle=atan2(v,u)/pi*180;
2126 trisurf(index,x/1e3,y/1e3,angle),shading interp,view(2);
2127 xlabel('X (km)','fontname','arial','fontsize',16)
2128 ylabel('Y (km)','fontname','arial','fontsize',16)
2129 title(texlabel('Observed velocity Vx (m.a^(-1)) around HIR. From RSAT1'),'fontname','arial','fontsize',18);
2130 set(gca,'fontname','arial','fontsize',16)
2131
2132
2133 zoom_choice(zo);
2134 colormap(jet)
2135
2136 colorbar
2137 set(gca,'Tag','ActiveAxes')
2138 end
2139
2140
2141
2142
2143
2144
2145
2146
2147
2148 function plot_diff_rotation(num,zo,series,range_choice,range);
2149
2150 names_function(num,series)
2151 load ../data/matlab/commondata/names.mat
2152 load ../data/matlab/commondata/data
2153 exist=exist_by_names([root directory '\' results '_results.mat']);
2154
2155 if exist==1,
2156
2157 aa=['load ' root directory '/' results '_results'];
2158 eval(aa)
2159 aa=['load ' root directory '/' maillage '_datamesh'];
2160 eval(aa)
2161 aa=['load ' root directory '/' maillage '_fem_parameters'];
2162 eval(aa)
2163 aa=['load ' root directory '/' maillage '_observations'];
2164 eval(aa)
2165 angle_model=atan2(v,u)/pi*180;
2166 angle_obs=atan2(v_HIR,u_HIR)/pi*180;
2167 trisurf(index,x/1e3,y/1e3,angle_model-angle_obs),shading interp,view(2);
2168 xlabel('X (km)','fontname','arial','fontsize',16)
2169 ylabel('Y (km)','fontname','arial','fontsize',16)
2170 title(texlabel('Rotation angle differential. alpha_model-alpha_obs'),'fontname','arial','fontsize',18);
2171 set(gca,'fontname','arial','fontsize',16)
2172
2173
2174 zoom_choice(zo);
2175 colormap(jet)
2176 caxis([-50 50]);
2177 colorbar
2178
2179 set(gca,'Tag','ActiveAxes')
2180 end
2181
2182
2183
2184 function plot_mesh(num,zo,series,range_choice,range);
2185
2186 names_function(num,series)
2187 load ../data/matlab/commondata/names.mat
2188 load ../data/matlab/commondata/data
2189 exist=exist_by_names([root directory '\' results '_results.mat']);
2190
2191 if exist==1,
2192
2193 aa=['load ' root directory '/' results '_results'];
2194 eval(aa)
2195 aa=['load ' root directory '/' maillage '_datamesh'];
2196 eval(aa)
2197 aa=['load ' root directory '/' maillage '_fem_parameters'];
2198 eval(aa)
2199 aa=['load ' root directory '/' maillage '_observations'];
2200 eval(aa)
2201 u=zeros(nods,1);
2202
2203 trisurf(index,x/1e3,y/1e3,u),shading faceted,view(2);
2204 xlabel('X (km)','fontname','arial','fontsize',16)
2205 ylabel('Y (km)','fontname','arial','fontsize',16)
2206 title(texlabel('Mesh'),'fontname','arial','fontsize',18);
2207 set(gca,'fontname','arial','fontsize',16)
2208
2209
2210 zoom_choice(zo);
2211 set(gca,'Tag','ActiveAxes')
2212 end
2213
2214
2215
2216
2217
2218
2219
2220
2221
2222
2223
2224 function names_function(num,series,range_choice,range)
2225 location=find(pwd=='\');
2226 location=location(length(location));
2227 present_directory=pwd;
2228 present_directory=present_directory(location+1:length(present_directory));
2229
2230 root=['..\data\' present_directory '\'];
2231
2232 directory=['RONNE_w' num2str(series)];
2233 rad_maillage='RONNE_w';
2234 rad_results='w';
2235
2236 maillage=[rad_maillage num2str(num) '_mesh'];
2237 results=[rad_results num2str(num)];
2238
2239 save ../data/matlab/commondata/names.mat maillage results directory root
2240
2241 function f=exist_by_names(str);
2242 fid=fopen(str,'r');
2243 if fid==-1,
2244 f=0;
2245 else
2246 f=1;
2247 end
2248
2249
2250
2251 function zoom_choice(zo);
2252 xmin=-1.09*10^6;
2253 xmax=-960000;
2254
2255 ymin=8.1*10^5;
2256 ymax=9.4*10^5;
2257 switch zo,
2258 case 0,
2259 xmin=-1.5*10^6;
2260 xmax=-.5*10^6;
2261 ymin=1*10^5;
2262 ymax=11*10^5;
2263 xlim([xmin/1e3 xmax/1e3]);
2264 ylim([ymin/1e3 ymax/1e3]);
2265 case 1,
2266 xmin=-1.09*10^6;
2267 xmax=-960000;
2268 ymin=8.1*10^5;
2269 ymax=9.4*10^5;
2270 xlim([xmin/1e3 xmax/1e3]);
2271 ylim([ymin/1e3 ymax/1e3]);
2272 case 2,
2273 xmin=-1.46*10^6;
2274 xmax=-1.32*10^6;
2275 ymin=5.5*10^5;
2276 ymax=8.42*10^5;
2277 xlim([xmin/1e3 xmax/1e3]);
2278
2279 ylim([ymin/1e3 ymax/1e3]);
2280 case 3,
2281 xmin=-1.26*10^6;
2282 xmax=-9.7*10^5;
2283 ymin=8.20*10^5;
2284 ymax=9.40*10^5;
2285 xlim([xmin/1e3 xmax/1e3]);
2286 ylim([ymin/1e3 ymax/1e3]);
2287
2288
2289 end