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

Added CM drag square

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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);
Note: See TracChangeset for help on using the changeset viewer.