Changeset 21927
- Timestamp:
- 08/09/17 17:55:25 (8 years ago)
- Location:
- issm/trunk-jpl/test
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test444.m
r21820 r21927 20 20 in=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,'../Exp/SquareHalfRight.exp','node',1); 21 21 regionalmask(find(in))=1; 22 md.transient.requested_outputs={'default','GroundedArea 1','FloatingArea1','TotalFloatingBmb1','TotalGroundedBmb1','TotalSmb1',...22 md.transient.requested_outputs={'default','GroundedArea','FloatingArea','IceVolumeAboveFloatation','GroundedArea1','FloatingArea1','TotalFloatingBmb1','TotalGroundedBmb1','TotalSmb1',... 23 23 'IceMass1','IceVolume1','IceVolumeAboveFloatation1','IceVolumeAboveFloatation'}; 24 24 md.outputdefinition.definitions{1}=regionaloutput('name','GroundedArea1','outputnamestring','GroundedArea','mask',regionalmask,... … … 38 38 md.outputdefinition.definitions{8}=regionaloutput('name','TotalSmb1','outputnamestring','TotalSmb','mask',regionalmask,... 39 39 'definitionstring','Outputdefinition8'); 40 md.outputdefinition.definitions{9}=regionaloutput('name','TotalSmb2','outputnamestring','TotalSmb','mask',regionalmask,... 41 'definitionstring','Outputdefinition9'); 40 42 41 43 md=extrude(md,3,1); … … 58 60 md.qmu.responses.TotalGroundedBmb1=response_function('Outputdefinition4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 59 61 md.qmu.responses.TotalSmb1=response_function('Outputdefinition8',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 62 md.qmu.responses.TotalSmb2=response_function('Outputdefinition9',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 63 md.qmu.responses.FloatingArea=response_function('FloatingArea',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 60 64 61 65 md.qmu.method =dakota_method('nond_samp'); … … 94 98 %we put all the mean and stdev data in the montecarlo field, which we will use to test for success. 95 99 md.results.dakota.montecarlo=[]; 96 for i=1: 8,100 for i=1:11, 97 101 md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).mean]; 98 102 end 99 for i=1: 8,103 for i=1:11, 100 104 md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).stddev]; 101 105 end
Note:
See TracChangeset
for help on using the changeset viewer.