plot_results

PURPOSE ^

PLOT_RESULTS - plot results

SYNOPSIS ^

function plot_results(Action,varargin)

DESCRIPTION ^

PLOT_RESULTS - plot results

   This routine seems to be obsolete... move to Trash?

   Usage:
      plot_results(Action,varargin)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SUBFUNCTIONS ^

SOURCE CODE ^

0001 function plot_results(Action,varargin)
0002 %PLOT_RESULTS - plot results
0003 %
0004 %   This routine seems to be obsolete... move to Trash?
0005 %
0006 %   Usage:
0007 %      plot_results(Action,varargin)
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    %---Initialize the demonstration
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 %---Add the 3D plot as the default
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 %%%%%%%%%%%%%%%%%%%%%%%%5555
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 %%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %caxis([-200 200]);
0926 colorbar
0927 set(gca,'Tag','ActiveAxes')
0928 end
0929 
0930 
0931 
0932 
0933 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
1066 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
1136 
1137 
1138 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
1279 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
1362 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %caxis([0 1500]);
1438 colorbar
1439 set(gca,'Tag','ActiveAxes')
1440 end
1441 
1442 
1443 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %caxis([0 1500]);
1473 colorbar
1474 set(gca,'Tag','ActiveAxes')
1475 end
1476 
1477 
1478 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %caxis([0 1500]);
1508 colorbar
1509 set(gca,'Tag','ActiveAxes')
1510 end
1511 
1512 
1513 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %pos=find(vel_LAS==0);
1616 %vel_LAS(pos)=NaN;
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
1943 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %caxis([-500 500]);
1973 colorbar
1974 set(gca,'Tag','ActiveAxes')
1975 end
1976 
1977 
1978 
1979 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
1980 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
2019 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
2065 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
2109 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %caxis([-500 500]);
2136 colorbar
2137 set(gca,'Tag','ActiveAxes')
2138 end
2139 
2140 
2141 
2142 
2143 
2144 
2145 
2146 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
2147 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
2217 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
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

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003