Changeset 5026


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

Added Square CM B tests

Location:
issm/trunk/test
Files:
17 added
5 edited

Legend:

Unmodified
Added
Removed

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

    r5025 r5026  
    6363elseif (id==221), name='SquareShelfTranS3dSerial';
    6464elseif (id==222), name='SquareShelfTranS3dParallel';
    65 elseif (id==223), name='SquareShelfCMB2dSerial';
    66 elseif (id==224), name='SquareShelfCMB2dParallel';
    67 elseif (id==225), name='SquareShelfCMB3dSerial';
    68 elseif (id==226), name='SquareShelfCMB3dParallel';
     65elseif (id==223), name='SquareShelfCMBM2dSerial';
     66elseif (id==224), name='SquareShelfCMBM2dParallel';
     67elseif (id==225), name='SquareShelfCMBM3dSerial';
     68elseif (id==226), name='SquareShelfCMBM3dParallel';
     69elseif (id==227), name='SquareShelfCMBP3dSerial';
     70elseif (id==228), name='SquareShelfCMBP3dParallel';
     71elseif (id==229), name='SquareShelfCMBS3dSerial';
     72elseif (id==230), name='SquareShelfCMBS3dParallel';
    6973elseif (id==231), name='SquareShelfDiagM2dRiftSerial';
    7074elseif (id==232), name='SquareShelfDiagM2dRiftParallel';

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

    r4998 r5026  
    1 md=mesh(model,'../Exp/Square.exp','../Exp/Rifts.exp',50000);
    2 md=meshprocessrifts(md,'../Exp/Square.exp');
     1md=mesh(model,'../Exp/Square.exp',200000);
    32md=geography(md,'all','');
    43md=parameterize(md,'../Par/SquareShelf.par');
    5 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='rheology_B';
     10md.cm_min=10^6;
     11md.cm_max=2*10^9;
     12md.nsteps=2;
     13md.fit=0*ones(md.nsteps,1);
     14md.weights=ones(md.numberofgrids,1);
     15md.optscal=10^7*ones(md.nsteps,1);
     16md.maxiter=2*ones(md.nsteps,1);
     17md.cm_jump=0.3*ones(md.nsteps,1);
     18md.vx_obs=md.vx; md.vy_obs=md.vy;
     19
    620md.cluster='none';
    7 md.solverstring=' -mat_type aijmumps -ksp_type preonly -pc_type lu -mat_mumps_icntl_14 80 ';
    8 md=solversettoasm(md);
    9 error('rifts not working for now')
    1021md=solve(md,'analysis_type',DiagnosticSolutionEnum);

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

    r4998 r5026  
    1 field_names     ={'Vx','Vy','Vel','Pressure'};
    2 field_tolerances={1e-13,1e-13,1e-13,1e-13};
     1field_names     ={'Gradient' 'Misfits' 'RheologyB2d' '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};
    33field_values={...
     4        PatchToVec(md.results.DiagnosticSolution.Gradient),...
     5        md.results.DiagnosticSolution.J,...
     6        PatchToVec(md.results.DiagnosticSolution.RheologyB),...
     7        PatchToVec(md.results.DiagnosticSolution.Pressure),...
     8        PatchToVec(md.results.DiagnosticSolution.Vel),...
    49        PatchToVec(md.results.DiagnosticSolution.Vx),...
    5         PatchToVec(md.results.DiagnosticSolution.Vy),...
    6         PatchToVec(md.results.DiagnosticSolution.Vel),...
    7         PatchToVec(md.results.DiagnosticSolution.Pressure),...
    8         };
     10        PatchToVec(md.results.DiagnosticSolution.Vy)
     11};

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

    r4998 r5026  
    1 md=mesh(model,'../Exp/Square.exp','../Exp/Rifts.exp',50000);
    2 md=meshprocessrifts(md,'../Exp/Square.exp');
     1md=mesh(model,'../Exp/Square.exp',200000);
    32md=geography(md,'all','');
    43md=parameterize(md,'../Par/SquareShelf.par');
    5 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='rheology_B';
     10md.cm_min=10^6;
     11md.cm_max=2*10^9;
     12md.nsteps=2;
     13md.fit=0*ones(md.nsteps,1);
     14md.weights=ones(md.numberofgrids,1);
     15md.optscal=10^7*ones(md.nsteps,1);
     16md.maxiter=2*ones(md.nsteps,1);
     17md.cm_jump=0.3*ones(md.nsteps,1);
     18md.vx_obs=md.vx; md.vy_obs=md.vy;
     19
    620md.cluster=oshostname();
    7 md.solverstring=' -mat_type aijmumps -ksp_type preonly -pc_type lu -mat_mumps_icntl_14 80 ';
    8 md=solversettomumps(md);
    9 error('rifts not working for now')
    1021md=solve(md,'analysis_type',DiagnosticSolutionEnum);

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

    r4998 r5026  
    1 field_names     ={'Vx','Vy','Vel','Pressure'};
    2 field_tolerances={1e-13,1e-13,1e-13,1e-13};
     1field_names     ={'Gradient' 'Misfits' 'RheologyB2d' '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};
    33field_values={...
     4        PatchToVec(md.results.DiagnosticSolution.Gradient),...
     5        md.results.DiagnosticSolution.J,...
     6        PatchToVec(md.results.DiagnosticSolution.RheologyB),...
     7        PatchToVec(md.results.DiagnosticSolution.Pressure),...
     8        PatchToVec(md.results.DiagnosticSolution.Vel),...
    49        PatchToVec(md.results.DiagnosticSolution.Vx),...
    5         PatchToVec(md.results.DiagnosticSolution.Vy),...
    6         PatchToVec(md.results.DiagnosticSolution.Vel),...
    7         PatchToVec(md.results.DiagnosticSolution.Pressure),...
    8         };
     10        PatchToVec(md.results.DiagnosticSolution.Vy)
     11};
Note: See TracChangeset for help on using the changeset viewer.