Changeset 5193
- Timestamp:
- 08/12/10 12:07:39 (15 years ago)
- Location:
- issm/trunk/test/NightlyRun
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test231.m
r5119 r5193 7 7 8 8 %solver, to handle penalties 9 md=solversettoasm(md);10 %md=solversettomumps(md);9 %md=solversettoasm(md); 10 md=solversettomumps(md); 11 11 12 12 %rift settings -
issm/trunk/test/NightlyRun/test232.m
r5119 r5193 5 5 md=setelementstype(md,'macayeal','all'); 6 6 md.cluster=oshostname(); 7 md.np= 17 md.np=2 8 8 9 9 %solver, to handle penalties … … 14 14 md.rifts.fill=MelangeEnum(); 15 15 md.rifts.fraction=0; 16 md.penalty_lock=2; 16 md.penalty_lock=10; 17 md.min_mechanical_constraints=0; 17 18 md.rifts.fractionincrement=.1; 18 19 19 md=solve(md,'analysis_type',DiagnosticSolutionEnum); 20 20 -
issm/trunk/test/NightlyRun/test423.m
r5115 r5193 6 6 7 7 %Dakota options 8 md.variables.nuv=normal_uncertain.empty(); 9 md.variables. nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);10 md.variables. nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);8 %variables 9 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01); 10 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01); 11 11 12 md.responses=struct(); 13 md.responses.rf =response_function.empty(); 14 md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 12 %responses 13 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 15 14 15 %method 16 16 md.qmu_method =dakota_method('nond_l'); 17 17 18 %parameters 18 19 md.qmu_params.direct=true; 19 20 md.qmu_params.analysis_driver='diagnostic'; … … 22 23 md.qmu_params.tabular_graphics_data=true; 23 24 25 %partitioning 24 26 md.part=0:1:md.numberofgrids-1; 25 27 md.part=md.part'; … … 27 29 md.qmu_analysis=1; 28 30 31 %imperative! 29 32 md.eps_rel=10^-5; %tighten for qmu analysese 30 33 34 %solve 31 35 md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y'); 32 36 -
issm/trunk/test/NightlyRun/test424.m
r5115 r5193 6 6 7 7 %Dakota options 8 md.variables.nuv=normal_uncertain.empty(); 9 md.variables. nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);10 md.variables. nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);8 %variables 9 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01); 10 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01); 11 11 12 md.responses=struct(); 13 md.responses.rf =response_function.empty(); 14 md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 12 %responses 13 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 15 14 15 %method 16 16 md.qmu_method =dakota_method('nond_l'); 17 17 18 %parameters 18 19 md.qmu_params.direct=true; 19 20 md.qmu_params.analysis_driver='diagnostic'; … … 22 23 md.qmu_params.tabular_graphics_data=true; 23 24 25 %partitioning 24 26 md.part=0:1:md.numberofgrids-1; 25 27 md.part=md.part'; … … 27 29 md.qmu_analysis=1; 28 30 31 %imperative! 29 32 md.eps_rel=10^-5; %tighten for qmu analysese 30 33 34 %solve 31 35 md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y'); 32 36 -
issm/trunk/test/NightlyRun/test425.m
r5115 r5193 6 6 7 7 %Dakota options 8 md.variables.nuv=normal_uncertain.empty();9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);10 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);11 8 12 md.responses=struct(); 13 md.responses.rf =response_function.empty(); 14 md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 9 %partitioning 10 md.npart=20; 11 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); 12 md.part=md.part-1; 15 13 14 %variables 15 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01); 16 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01); 17 18 %responses 19 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 20 21 %method 16 22 md.qmu_method =dakota_method('nond_l'); 17 23 24 %parameters 18 25 md.qmu_params.direct=true; 19 26 md.qmu_params.analysis_driver='diagnostic'; … … 22 29 md.qmu_params.tabular_graphics_data=true; 23 30 24 %partitioning 25 md.npart=20; 26 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); 27 md.part=md.part-1; 31 32 %imperative! 33 md.eps_rel=10^-5; %tighten for qmu analyses 28 34 md.qmu_analysis=1; 29 35 30 md.eps_rel=10^-5; %tighten for qmu analyses 31 36 %solve 32 37 md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y'); 33 38 -
issm/trunk/test/NightlyRun/test426.m
r5115 r5193 6 6 7 7 %Dakota options 8 md.variables.nuv=normal_uncertain.empty();9 md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);10 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);11 8 12 md.responses=struct(); 13 md.responses.rf =response_function.empty(); 14 md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 9 %partitioning 10 md.npart=20; 11 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); 12 md.part=md.part-1; 15 13 14 %variables 15 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01); 16 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01); 17 18 %responses 19 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 20 21 %method 16 22 md.qmu_method =dakota_method('nond_l'); 17 23 24 %parameters 18 25 md.qmu_params.direct=true; 19 26 md.qmu_params.analysis_driver='diagnostic'; … … 22 29 md.qmu_params.tabular_graphics_data=true; 23 30 24 %partitioning 25 md.npart=20; 26 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); 27 md.part=md.part-1; 31 32 %imperative! 33 md.eps_rel=10^-5; %tighten for qmu analyses 28 34 md.qmu_analysis=1; 29 35 30 md.eps_rel=10^-5; %tighten for qmu analyses 31 36 %solve 32 37 md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y'); 33 38 -
issm/trunk/test/NightlyRun/test427.m
r5115 r5193 14 14 15 15 %Dakota options 16 md.variables.nuv=normal_uncertain.empty();17 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);18 16 19 md.responses=struct(); 20 md.responses.rf =response_function.empty(); 21 md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 22 md.responses.rf (end+1)=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 23 md.responses.rf (end+1)=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 24 md.responses.rf (end+1)=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 25 md.responses.rf (end+1)=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 26 md.responses.rf (end+1)=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 27 md.responses.rf (end+1)=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 17 %partitioning 18 md.npart=20; 19 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); 20 md.part=md.part-1; 21 22 %variables 23 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01); 24 25 %responses 26 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 27 md.responses.MassFlux1=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 28 md.responses.MassFlux2=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 29 md.responses.MassFlux3=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 30 md.responses.MassFlux4=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 31 md.responses.MassFlux5=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 32 md.responses.massFlux6=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 28 33 29 34 %mass flux profiles … … 31 36 md.qmu_mass_flux_profile_directory=pwd; 32 37 38 %method 33 39 md.qmu_method =dakota_method('nond_l'); 34 40 41 %parameters 35 42 md.qmu_params.direct=true; 36 43 md.qmu_params.analysis_driver='diagnostic'; … … 38 45 md.qmu_params.interval_type='forward'; 39 46 md.qmu_params.tabular_graphics_data=true; 40 41 %partitioning42 md.npart=20;43 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');44 md.part=md.part-1;45 47 md.qmu_analysis=1; 46 47 48 md.eps_rel=10^-5; %tighten for qmu analyses 48 49 50 %solve 49 51 md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y'); 52 md=tres(md,'dakota'); 50 53 51 54 %Fields and tolerances to track changes 52 md=tres(md,'dakota'); 55 53 56 %ok, mass flux of 3 profiles should be -3 Gt/yr -3 Gt/yr and the sum, which is -6 Gt/yr 54 57 %we recover those mass fluxes through the mean of the response. -
issm/trunk/test/NightlyRun/test428.m
r5115 r5193 14 14 15 15 %Dakota options 16 md.variables.nuv=normal_uncertain.empty();17 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);18 16 19 md.responses=struct(); 20 md.responses.rf =response_function.empty(); 21 md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 22 md.responses.rf (end+1)=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 23 md.responses.rf (end+1)=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 24 md.responses.rf (end+1)=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 25 md.responses.rf (end+1)=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 26 md.responses.rf (end+1)=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 27 md.responses.rf (end+1)=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 17 %partitioning 18 md.npart=20; 19 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on'); 20 md.part=md.part-1; 21 22 %variables 23 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01); 24 25 %responses 26 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 27 md.responses.MassFlux1=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 28 md.responses.MassFlux2=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 29 md.responses.MassFlux3=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 30 md.responses.MassFlux4=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 31 md.responses.MassFlux5=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 32 md.responses.massFlux6=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); 28 33 29 34 %mass flux profiles … … 31 36 md.qmu_mass_flux_profile_directory=pwd; 32 37 38 %method 33 39 md.qmu_method =dakota_method('nond_l'); 34 40 41 %parameters 35 42 md.qmu_params.direct=true; 36 43 md.qmu_params.analysis_driver='diagnostic'; … … 38 45 md.qmu_params.interval_type='forward'; 39 46 md.qmu_params.tabular_graphics_data=true; 40 41 %partitioning42 md.npart=20;43 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');44 md.part=md.part-1;45 47 md.qmu_analysis=1; 46 47 48 md.eps_rel=10^-5; %tighten for qmu analyses 48 49 50 %solve 49 51 md=solve(md,'analysis_type',DiagnosticSolutionEnum,'overwrite','y'); 52 md=tres(md,'dakota'); 50 53 51 54 %Fields and tolerances to track changes 52 md=tres(md,'dakota');53 55 %ok, mass flux of 3 profiles should be -3 Gt/yr -3 Gt/yr and the sum, which is -6 Gt/yr 54 56 %we recover those mass fluxes through the mean of the response.
Note:
See TracChangeset
for help on using the changeset viewer.