Changeset 5033


Ignore:
Timestamp:
08/06/10 14:54:55 (15 years ago)
Author:
Mathieu Morlighem
Message:

Added CM drag square

Location:
issm/trunk/test
Files:
18 added
7 edited

Legend:

Unmodified
Added
Removed

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

    r5026 r5033  
    109109elseif (id==335), name='SquareSheetConstrainedSteaH3dSerial';
    110110elseif (id==336), name='SquareSheetConstrainedSteaH3dParallel';
     111elseif (id==337), name='SquareSheetConstrainedCMDragM2dSerial';
     112elseif (id==338), name='SquareSheetConstrainedCMDragM2dParallel';
     113elseif (id==339), name='SquareSheetConstrainedCMDragM3dSerial';
     114elseif (id==340), name='SquareSheetConstrainedCMDragM3dParallel';
     115elseif (id==341), name='SquareSheetConstrainedCMDragP3dSerial';
     116elseif (id==342), name='SquareSheetConstrainedCMDragP3dParallel';
     117elseif (id==343), name='SquareSheetConstrainedCMDragS3dSerial';
     118elseif (id==344), name='SquareSheetConstrainedCMDragS3dParallel';
    111119elseif (id==401), name='SquareSheetShelfDiagM2dSerial';
    112120elseif (id==402), name='SquareSheetShelfDiagM2dParallel';

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

    r4996 r5033  
    1 md=mesh(model,'../Exp/Square.exp',180000);
     1md=mesh(model,'../Exp/Square.exp',200000);
    22md=geography(md,'','');
    33md=parameterize(md,'../Par/SquareSheetConstrained.par');
    4 md=extrude(md,4,1);
    5 md=setelementstype(md,'hutter','all');
     4md=extrude(md,3,1);
     5md=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.vx_obs=md.vx; md.vy_obs=md.vy;
     20
    621md.cluster='none';
    7 md.dt=0;
    8 md=solve(md,'analysis_type',SteadystateSolutionEnum);
     22md=solve(md,'analysis_type',DiagnosticSolutionEnum);

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

    r4996 r5033  
    1 field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
    2 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
     1field_names     ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy'};
     2field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     3field_values={...
     4        PatchToVec(md.results.DiagnosticSolution.Gradient),...
     5        md.results.DiagnosticSolution.J,...
     6        PatchToVec(md.results.DiagnosticSolution.DragCoefficient),...
     7        PatchToVec(md.results.DiagnosticSolution.Pressure),...
     8        PatchToVec(md.results.DiagnosticSolution.Vel),...
     9        PatchToVec(md.results.DiagnosticSolution.Vx),...
     10        PatchToVec(md.results.DiagnosticSolution.Vy)
    311};
    4 field_values={...
    5         PatchToVec(md.results.SteadystateSolution.Vx),...
    6         PatchToVec(md.results.SteadystateSolution.Vy),...
    7         PatchToVec(md.results.SteadystateSolution.Vz),...
    8         PatchToVec(md.results.SteadystateSolution.Vel),...
    9         PatchToVec(md.results.SteadystateSolution.Pressure),...
    10         PatchToVec(md.results.SteadystateSolution.Temperature),...
    11         PatchToVec(md.results.SteadystateSolution.MeltingRate),...
    12         };

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

    r4996 r5033  
    1 md=mesh(model,'../Exp/Square.exp',180000);
     1md=mesh(model,'../Exp/Square.exp',200000);
    22md=geography(md,'','');
    33md=parameterize(md,'../Par/SquareSheetConstrained.par');
    4 md=extrude(md,4,1);
    5 md=setelementstype(md,'hutter','all');
     4md=extrude(md,3,1);
     5md=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.vx_obs=md.vx; md.vy_obs=md.vy;
     20
    621md.cluster=oshostname();
    7 md.dt=0;
    8 md=solve(md,'analysis_type',SteadystateSolutionEnum);
     22md=solve(md,'analysis_type',DiagnosticSolutionEnum);

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

    r4996 r5033  
    1 field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
    2 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
     1field_names     ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy'};
     2field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     3field_values={...
     4        PatchToVec(md.results.DiagnosticSolution.Gradient),...
     5        md.results.DiagnosticSolution.J,...
     6        PatchToVec(md.results.DiagnosticSolution.DragCoefficient),...
     7        PatchToVec(md.results.DiagnosticSolution.Pressure),...
     8        PatchToVec(md.results.DiagnosticSolution.Vel),...
     9        PatchToVec(md.results.DiagnosticSolution.Vx),...
     10        PatchToVec(md.results.DiagnosticSolution.Vy)
    311};
    4 field_values={...
    5         PatchToVec(md.results.SteadystateSolution.Vx),...
    6         PatchToVec(md.results.SteadystateSolution.Vy),...
    7         PatchToVec(md.results.SteadystateSolution.Vz),...
    8         PatchToVec(md.results.SteadystateSolution.Vel),...
    9         PatchToVec(md.results.SteadystateSolution.Pressure),...
    10         PatchToVec(md.results.SteadystateSolution.Temperature),...
    11         PatchToVec(md.results.SteadystateSolution.MeltingRate),...
    12         };
Note: See TracChangeset for help on using the changeset viewer.