Changeset 14118


Ignore:
Timestamp:
12/09/12 21:25:20 (12 years ago)
Author:
Eric.Larour
Message:

CHG: minor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/workshop/2012/Talks/10_PIG_UQ/runme.m

    r14029 r14118  
    6969       
    7070        %load model and cross over errors
    71         md = loadmodel('./Models/Pig.Control_drag');
     71        md = loadmodel('../09_PIG_Application/Models/Pig.Control_drag');
    7272        load -mat Models/Pig.CrossOvers
    7373
    7474        %partition the mesh
    7575        md.qmu.numberofpartitions=50;
    76         md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
     76        md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,...
     77        'weighting','on');
    7778        md.qmu.partition=md.qmu.partition-1; %switch partition to c-indexing
    7879
     
    144145       
    145146        %load model
    146         md = loadmodel('./Models/Pig.Control_drag');
     147        md = loadmodel('../09_PIG_Application/Models/Pig.Control_drag');
    147148
    148149        %partition the mesh
     
    155156        %variables
    156157        md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05);
    157         md.qmu.variables.DragCoefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.05);
    158         md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,0.05);
     158        md.qmu.variables.DragCoefficient=normal_uncertain(...
     159        'scaled_FrictionCoefficient',1,0.05);
     160        md.qmu.variables.rheology_B=normal_uncertain(...
     161        'scaled_MaterialsRheologyB',1,0.05);
    159162
    160163        %responses
     
    218221        md = loadmodel('./Models/Pig.Sampling');
    219222       
    220         plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,'linewidth',2,'axis#all','image','unit','km','colorbar','off','title','','meshcolor','b','grid','on');
     223        plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,...
     224        'linewidth',2, 'axis#all','image','unit','km','colorbar','off',...
     225        'title','','meshcolor','b','grid','on');
    221226
    222227end
     
    232237        %retrieve results for this specific profile
    233238        result=md.results.dakota.dresp_dat(md.qmu.numberofpartitions+index);
    234         error;
    235239               
    236240        %plot histogram
    237         plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off','xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor','none','EdgeColor','red');
     241        plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off',...
     242        'xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor',...
     243        'none','EdgeColor','red');
    238244
    239245end
Note: See TracChangeset for help on using the changeset viewer.