Changeset 5047


Ignore:
Timestamp:
08/06/10 16:02:13 (15 years ago)
Author:
Mathieu Morlighem
Message:

Added nightly run steadystate

Location:
issm/trunk/test
Files:
13 edited

Legend:

Unmodified
Added
Removed

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

    r5033 r5047  
    145145elseif (id==427), name='SquareSheetShelfDiadM3dDakotaMassFluxSerial';
    146146elseif (id==428), name='SquareSheetShelfDiadM3dDakotaMassFluxParallel';
     147elseif (id==429), name='SquareSheetShelfCMDragSteaM3dSerial';
     148elseif (id==430), name='SquareSheetShelfCMDragSteaM3dParallel';
     149elseif (id==431), name='SquareSheetShelfCMDragSteaP3dSerial';
     150elseif (id==432), name='SquareSheetShelfCMDragSteaP3dParallel';
    147151elseif (id==501), name='PigDiagM2dSerial';
    148152elseif (id==502), name='PigDiagM2dParallel';
    … …  
    165169elseif (id==519), name='PigSteaS3dSerial';
    166170elseif (id==520), name='PigSteaS3dParallel';
     171elseif (id==521), name='PigSteaS3dParallel';
    167172elseif (id==601), name='79NorthProg2dSerial';
    168173elseif (id==602), name='79NorthProg2dParallel';

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

    r4998 r5047  
    1 md=mesh(model,'../Exp/Square.exp',150000);
     1md=mesh(model,'../Exp/Square.exp',200000);
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
     4md=extrude(md,3,1);
    45md=setelementstype(md,'macayeal','all');
     6
     7%control parameters
     8md.control_analysis=1;
     9md.control_type='drag_coefficient';
     10md.cm_min=1;
     11md.cm_max=200;
     12md.nsteps=2;
     13md.fit=2*ones(md.nsteps,1);
     14md.weights=ones(md.numberofgrids,1);
     15md.optscal=3*ones(md.nsteps,1);
     16md.maxiter=2*ones(md.nsteps,1);
     17md.cm_jump=0.3*ones(md.nsteps,1);
     18md.cm_noisedmp=2*10^-7;
     19md.dt=0;
     20md.vx_obs=md.vx; md.vy_obs=md.vy;
     21
    522md.cluster='none';
    6 
    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 
    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]);
    15 
    16 md.qmu_method     =dakota_method('nond_l');
    17 
    18 md.qmu_params.direct=true;
    19 md.qmu_params.analysis_driver='diagnostic';
    20 md.qmu_params.evaluation_concurrency=1;
    21 md.qmu_params.interval_type='forward';
    22 md.qmu_params.tabular_graphics_data=true;
    23 
    24 %partitioning
    25 md.npart=20;
    26 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    27 md.part=md.part-1;
    28 md.qmu_analysis=1;
    29 
    30 md.eps_rel=10^-10; %tighten for qmu analyses
    31 
    32 %rm existing results
    33 if exist('./qmu')
    34         system('rm -rf qmu/');
    35 end
    36 
    37 md=solve(md,'analysis_type',DiagnosticSolutionEnum);
     23md=solve(md,'analysis_type',SteadystateSolutionEnum);

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

    r4998 r5047  
    1 md=tres(md,'dakota');
    2 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
    3 field_names     ={'importancefactors'};
    4 field_tolerances={1e-11};
     1field_names     ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'MeltingRate'};
     2field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    53field_values={...
    6          md.results.dakota.importancefactors,...
    7         };
     4        PatchToVec(md.results.SteadystateSolution.Gradient),...
     5        md.results.SteadystateSolution.J,...
     6        PatchToVec(md.results.SteadystateSolution.DragCoefficient),...
     7        PatchToVec(md.results.SteadystateSolution.Pressure),...
     8        PatchToVec(md.results.SteadystateSolution.Vel),...
     9        PatchToVec(md.results.SteadystateSolution.Vx),...
     10        PatchToVec(md.results.SteadystateSolution.Vy),...
     11        PatchToVec(md.results.SteadystateSolution.Vz),...
     12        PatchToVec(md.results.SteadystateSolution.Temperature),...
     13        PatchToVec(md.results.SteadystateSolution.MeltingRate)
     14};

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

    r4998 r5047  
    1 md=mesh(model,'../Exp/Square.exp',150000);
     1md=mesh(model,'../Exp/Square.exp',200000);
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
     4md=extrude(md,3,1);
    45md=setelementstype(md,'macayeal','all');
     6
     7%control parameters
     8md.control_analysis=1;
     9md.control_type='drag_coefficient';
     10md.cm_min=1;
     11md.cm_max=200;
     12md.nsteps=2;
     13md.fit=2*ones(md.nsteps,1);
     14md.weights=ones(md.numberofgrids,1);
     15md.optscal=3*ones(md.nsteps,1);
     16md.maxiter=2*ones(md.nsteps,1);
     17md.cm_jump=0.3*ones(md.nsteps,1);
     18md.cm_noisedmp=2*10^-7;
     19md.dt=0;
     20md.vx_obs=md.vx; md.vy_obs=md.vy;
     21
    522md.cluster=oshostname();
    6 
    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 
    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]);
    15 
    16 md.qmu_method     =dakota_method('nond_l');
    17 
    18 md.qmu_params.direct=true;
    19 md.qmu_params.analysis_driver='diagnostic';
    20 md.qmu_params.evaluation_concurrency=1;
    21 md.qmu_params.interval_type='forward';
    22 md.qmu_params.tabular_graphics_data=true;
    23 
    24 %partitioning
    25 md.npart=20;
    26 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    27 md.part=md.part-1;
    28 md.qmu_analysis=1;
    29 
    30 md.eps_rel=10^-10; %tighten for qmu analyses
    31 
    32 %rm existing results
    33 if exist('./qmu')
    34         system('rm -rf qmu/');
    35 end
    36 
    37 md=solve(md,'analysis_type',DiagnosticSolutionEnum);
     23md=solve(md,'analysis_type',SteadystateSolutionEnum);

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

    r5013 r5047  
    1 md=tres(md,'dakota');
    2 md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
    3 field_names     ={'importancefactors'};
    4 field_tolerances={1e-10};
     1field_names     ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'MeltingRate'};
     2field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    53field_values={...
    6          md.results.dakota.importancefactors,...
    7         };
     4        PatchToVec(md.results.SteadystateSolution.Gradient),...
     5        md.results.SteadystateSolution.J,...
     6        PatchToVec(md.results.SteadystateSolution.DragCoefficient),...
     7        PatchToVec(md.results.SteadystateSolution.Pressure),...
     8        PatchToVec(md.results.SteadystateSolution.Vel),...
     9        PatchToVec(md.results.SteadystateSolution.Vx),...
     10        PatchToVec(md.results.SteadystateSolution.Vy),...
     11        PatchToVec(md.results.SteadystateSolution.Vz),...
     12        PatchToVec(md.results.SteadystateSolution.Temperature),...
     13        PatchToVec(md.results.SteadystateSolution.MeltingRate)
     14};

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

    r4998 r5047  
    1 md=mesh(model,'../Exp/Square.exp',150000);
     1md=mesh(model,'../Exp/Square.exp',200000);
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
    4 md=setelementstype(md,'macayeal','all');
     4md=extrude(md,3,1);
     5md=setelementstype(md,'pattyn','all');
     6
     7%control parameters
     8md.control_analysis=1;
     9md.control_type='drag_coefficient';
     10md.cm_min=1;
     11md.cm_max=200;
     12md.nsteps=2;
     13md.fit=1*ones(md.nsteps,1);
     14md.weights=ones(md.numberofgrids,1);
     15md.optscal=3*ones(md.nsteps,1);
     16md.maxiter=2*ones(md.nsteps,1);
     17md.cm_jump=0.3*ones(md.nsteps,1);
     18md.cm_noisedmp=2*10^-7;
     19md.dt=0;
     20md.vx_obs=md.vx; md.vy_obs=md.vy;
     21
    522md.cluster='none';
    6 
    7 %constrain all velocities to 1 m/yr, in the y-direction
    8 md.spcvelocity(:,1:3)=1;
    9 md.spcvelocity(:,4)=0;
    10 md.spcvelocity(:,5)=1*md.yts;
    11 md.spcvelocity(:,6)=0;
    12 
    13 %Dakota options
    14 md.variables.nuv=normal_uncertain.empty();
    15 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    16 
    17 md.responses=struct();
    18 md.responses.rf =response_function.empty();
    19 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]);
    20 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]);
    21 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]);
    22 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]);
    23 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]);
    24 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]);
    25 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]);
    26 
    27 %mass flux profiles
    28 md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};
    29 md.qmu_mass_flux_profile_directory=pwd;
    30 
    31 md.qmu_method     =dakota_method('nond_l');
    32 
    33 md.qmu_params.direct=true;
    34 md.qmu_params.analysis_driver='diagnostic';
    35 md.qmu_params.evaluation_concurrency=1;
    36 md.qmu_params.interval_type='forward';
    37 md.qmu_params.tabular_graphics_data=true;
    38 
    39 %partitioning
    40 md.npart=20;
    41 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    42 md.part=md.part-1;
    43 md.qmu_analysis=1;
    44 
    45 md.eps_rel=10^-10; %tighten for qmu analyses
    46 
    47 %rm existing results
    48 if exist('./qmu')
    49         system('rm -rf qmu/');
    50 end
    51 
    52 md=solve(md,'analysis_type',DiagnosticSolutionEnum);
     23md=solve(md,'analysis_type',SteadystateSolutionEnum);

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

    r4998 r5047  
    1 md=tres(md,'dakota');
    2 %ok, mass flux of 3 profiles should be -.3, -.3 and the sum, which is -.6
    3 %we recover those mass fluxes through the mean of the response.
    4 %also, we recover the max velo, which should be 1m/s.
    5 %we put all that data in the importancefactors, which we will use to test for success.
    6 %also, check that the stddev are 0.
    7 md.results.dakota.importancefactors=[];
    8 for i=1:7,
    9         md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).mean];
    10 end
    11 for i=1:7,
    12         md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).stddev];
    13 end
    14 field_names     ={'importancefactors'};
    15 field_tolerances={1e-11};
     1field_names     ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'MeltingRate'};
     2field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    163field_values={...
    17          md.results.dakota.importancefactors,...
    18         };
     4        PatchToVec(md.results.SteadystateSolution.Gradient),...
     5        md.results.SteadystateSolution.J,...
     6        PatchToVec(md.results.SteadystateSolution.DragCoefficient),...
     7        PatchToVec(md.results.SteadystateSolution.Pressure),...
     8        PatchToVec(md.results.SteadystateSolution.Vel),...
     9        PatchToVec(md.results.SteadystateSolution.Vx),...
     10        PatchToVec(md.results.SteadystateSolution.Vy),...
     11        PatchToVec(md.results.SteadystateSolution.Vz),...
     12        PatchToVec(md.results.SteadystateSolution.Temperature),...
     13        PatchToVec(md.results.SteadystateSolution.MeltingRate)
     14};

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

    r4998 r5047  
    1 md=mesh(model,'../Exp/Square.exp',150000);
     1md=mesh(model,'../Exp/Square.exp',200000);
    22md=geography(md,'../Exp/SquareShelf.exp','');
    33md=parameterize(md,'../Par/SquareSheetShelf.par');
    4 md=setelementstype(md,'macayeal','all');
     4md=extrude(md,3,1);
     5md=setelementstype(md,'pattyn','all');
     6
     7%control parameters
     8md.control_analysis=1;
     9md.control_type='drag_coefficient';
     10md.cm_min=1;
     11md.cm_max=200;
     12md.nsteps=2;
     13md.fit=1*ones(md.nsteps,1);
     14md.weights=ones(md.numberofgrids,1);
     15md.optscal=3*ones(md.nsteps,1);
     16md.maxiter=2*ones(md.nsteps,1);
     17md.cm_jump=0.3*ones(md.nsteps,1);
     18md.cm_noisedmp=2*10^-7;
     19md.dt=0;
     20md.vx_obs=md.vx; md.vy_obs=md.vy;
     21
    522md.cluster=oshostname();
    6 
    7 %constrain all velocities to 1 m/yr, in the y-direction
    8 md.spcvelocity(:,1:3)=1;
    9 md.spcvelocity(:,4)=0;
    10 md.spcvelocity(:,5)=1*md.yts;
    11 md.spcvelocity(:,6)=0;
    12 
    13 %Dakota options
    14 md.variables.nuv=normal_uncertain.empty();
    15 md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
    16 
    17 md.responses=struct();
    18 md.responses.rf =response_function.empty();
    19 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]);
    20 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]);
    21 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]);
    22 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]);
    23 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]);
    24 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]);
    25 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]);
    26 
    27 %mass flux profiles
    28 md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};
    29 md.qmu_mass_flux_profile_directory=pwd;
    30 
    31 md.qmu_method     =dakota_method('nond_l');
    32 
    33 md.qmu_params.direct=true;
    34 md.qmu_params.analysis_driver='diagnostic';
    35 md.qmu_params.evaluation_concurrency=1;
    36 md.qmu_params.interval_type='forward';
    37 md.qmu_params.tabular_graphics_data=true;
    38 
    39 %partitioning
    40 md.npart=20;
    41 md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
    42 md.part=md.part-1;
    43 md.qmu_analysis=1;
    44 
    45 md.eps_rel=10^-10; %tighten for qmu analyses
    46 
    47 %rm existing results
    48 if exist('./qmu')
    49         system('rm -rf qmu/');
    50 end
    51 
    52 md=solve(md,'analysis_type',DiagnosticSolutionEnum);
     23md=solve(md,'analysis_type',SteadystateSolutionEnum);

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

    r4998 r5047  
    1 md=tres(md,'dakota');
    2 %ok, mass flux of 3 profiles should be -.3, -.3 and the sum, which is -.6
    3 %we recover those mass fluxes through the mean of the response.
    4 %also, we recover the max velo, which should be 1m/s.
    5 %we put all that data in the importancefactors, which we will use to test for success.
    6 %also, check that the stddev are 0.
    7 md.results.dakota.importancefactors=[];
    8 for i=1:7,
    9         md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).mean];
    10 end
    11 for i=1:7,
    12         md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).stddev];
    13 end
    14 field_names     ={'importancefactors'};
    15 field_tolerances={1e-11};
     1field_names     ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'MeltingRate'};
     2field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    163field_values={...
    17          md.results.dakota.importancefactors,...
    18         };
    19 md.results.dakota.importancefactors
     4        PatchToVec(md.results.SteadystateSolution.Gradient),...
     5        md.results.SteadystateSolution.J,...
     6        PatchToVec(md.results.SteadystateSolution.DragCoefficient),...
     7        PatchToVec(md.results.SteadystateSolution.Pressure),...
     8        PatchToVec(md.results.SteadystateSolution.Vel),...
     9        PatchToVec(md.results.SteadystateSolution.Vx),...
     10        PatchToVec(md.results.SteadystateSolution.Vy),...
     11        PatchToVec(md.results.SteadystateSolution.Vz),...
     12        PatchToVec(md.results.SteadystateSolution.Temperature),...
     13        PatchToVec(md.results.SteadystateSolution.MeltingRate)
     14};
Note: See TracChangeset for help on using the changeset viewer.