Changeset 5498


Ignore:
Timestamp:
08/23/10 09:23:32 (15 years ago)
Author:
Eric.Larour
Message:

New scaled,indexed,nodal expressions for qmu variables

Location:
issm/trunk/test/NightlyRun
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/test/NightlyRun/test235.m

    r5355 r5498  
    6161%Dakota options
    6262%variables
    63 md.variables.rheology_B=normal_uncertain('RheologyBbar',1,.5);
     63md.variables.rheology_B=normal_uncertain('scaled_RheologyBbar',1,.5);
    6464
    6565%responses
    … …  
    7474md.qmu_params.evaluation_concurrency=1;
    7575md.qmu_params.interval_type='forward';
    76 md.qmu_params.tabular_graphics_data=true;
    7776
    7877%imperative!
    … …  
    8584%Fields and tolerances to track changes
    8685md=tres(md,'dakota');
    87 md.results.dakota.importancefactors=importancefactors(md,'RheologyBbar','MaxVel')';
     86md.results.dakota.importancefactors=importancefactors(md,'scaled_RheologyBbar','MaxVel')';
    8887field_names     ={'importancefactors'};
    8988field_tolerances={1e-11};

  • issm/trunk/test/NightlyRun/test236.m

    r5355 r5498  
    6161%Dakota options
    6262%variables
    63 md.variables.rheology_B=normal_uncertain('RheologyBbar',1,.5);
     63md.variables.rheology_B=normal_uncertain('scaled_RheologyBbar',1,.05);
    6464
    6565%responses
    … …  
    7474md.qmu_params.evaluation_concurrency=1;
    7575md.qmu_params.interval_type='forward';
    76 md.qmu_params.tabular_graphics_data=true;
    7776
    7877%imperative!
    … …  
    8584%Fields and tolerances to track changes
    8685md=tres(md,'dakota');
    87 md.results.dakota.importancefactors=importancefactors(md,'RheologyBbar','MaxVel')';
     86md.results.dakota.importancefactors=importancefactors(md,'scaled_RheologyBbar','MaxVel')';
    8887field_names     ={'importancefactors'};
    8988field_tolerances={1e-13};

  • issm/trunk/test/NightlyRun/test423.m

    r5214 r5498  
    1414%variables
    1515md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
    16 
    17 %md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,.01);
    18 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,.01*ones(md.npart,1));
     16md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,.01);
     17%md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,.01*ones(md.npart,1));
    1918
    2019%responses
    … …  
    2928md.qmu_params.evaluation_concurrency=1;
    3029md.qmu_params.interval_type='forward';
    31 md.qmu_params.tabular_graphics_data=true;
    3230
    3331
    … …  
    4038%Fields and tolerances to track changes
    4139md=tres(md,'dakota');
    42 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     40md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
    4341field_names     ={'importancefactors'};
    4442field_tolerances={1e-11};

  • issm/trunk/test/NightlyRun/test424.m

    r5214 r5498  
    1515%variables
    1616md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
    17 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     17md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
    1818
    1919%responses
    … …  
    2828md.qmu_params.evaluation_concurrency=1;
    2929md.qmu_params.interval_type='forward';
    30 md.qmu_params.tabular_graphics_data=true;
    3130
    3231%imperative!
    … …  
    3837%Fields and tolerances to track changes
    3938md=tres(md,'dakota');
    40 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     39md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
    4140field_names     ={'importancefactors'};
    4241field_tolerances={1e-10};

  • issm/trunk/test/NightlyRun/test425.m

    r5193 r5498  
    1414%variables
    1515md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
    16 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     16md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
    1717
    1818%responses
    … …  
    2727md.qmu_params.evaluation_concurrency=1;
    2828md.qmu_params.interval_type='forward';
    29 md.qmu_params.tabular_graphics_data=true;
    3029
    3130
    … …  
    3938%Fields and tolerances to track changes
    4039md=tres(md,'dakota');
    41 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     40md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
    4241field_names     ={'importancefactors'};
    4342field_tolerances={1e-11};

  • issm/trunk/test/NightlyRun/test426.m

    r5303 r5498  
    1414%variables
    1515md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
    16 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     16md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
    1717
    1818%responses
    … …  
    2727md.qmu_params.evaluation_concurrency=1;
    2828md.qmu_params.interval_type='forward';
    29 md.qmu_params.tabular_graphics_data=true;
    3029
    3130
    … …  
    3938%Fields and tolerances to track changes
    4039md=tres(md,'dakota');
    41 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
     40md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
    4241field_names     ={'importancefactors'};
    4342field_tolerances={1e-10};

  • issm/trunk/test/NightlyRun/test427.m

    r5193 r5498  
    2121
    2222%variables
    23 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     23md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
    2424
    2525%responses
    2626md.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]);
     27md.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     28md.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     29md.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     30md.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     31md.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     32md.responses.massFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    3333
    3434%mass flux profiles
    … …  
    4444md.qmu_params.evaluation_concurrency=1;
    4545md.qmu_params.interval_type='forward';
    46 md.qmu_params.tabular_graphics_data=true;
    4746md.qmu_analysis=1;
    4847md.eps_rel=10^-5; %tighten for qmu analyses

  • issm/trunk/test/NightlyRun/test428.m

    r5193 r5498  
    2121
    2222%variables
    23 md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
     23md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
    2424
    2525%responses
    2626md.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]);
     27md.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     28md.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     29md.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     30md.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     31md.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     32md.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    3333
    3434%mass flux profiles
    … …  
    4444md.qmu_params.evaluation_concurrency=1;
    4545md.qmu_params.interval_type='forward';
    46 md.qmu_params.tabular_graphics_data=true;
    4746md.qmu_analysis=1;
    4847md.eps_rel=10^-5; %tighten for qmu analyses

  • issm/trunk/test/NightlyRun/test434.m

    r5118 r5498  
    1414
    1515%Dakota options
    16 md.variables.nuv=normal_uncertain.empty();
    17 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
     16md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
    1817
    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]);
     18md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     19md.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     20md.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     21md.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     22md.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     23md.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     24md.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    2825
    2926%mass flux profiles
    3027md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};
    3128md.qmu_mass_flux_profile_directory=pwd;
    32 
    3329
    3430
Note: See TracChangeset for help on using the changeset viewer.