Changeset 5033 for issm/trunk/test/NightlyRun/test340.m

Timestamp:

08/06/10 14:54:55 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added CM drag square

File:

• issm/trunk/test/NightlyRun/test340.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test340.m

@@ -1 @@
-md=mesh(model,'../Exp/Square.exp',180000);
+md=mesh(model,'../Exp/Square.exp',200000);
 md=geography(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=extrude(md,4,1);
-md=setelementstype(md,'hutter','all');
+md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','all');
+
+%control parameters
+md.control_analysis=1;
+md.control_type='drag_coefficient';
+md.cm_min=1;
+md.cm_max=200;
+md.nsteps=2;
+md.fit=2*ones(md.nsteps,1);
+md.weights=ones(md.numberofgrids,1);
+md.optscal=3*ones(md.nsteps,1);
+md.maxiter=2*ones(md.nsteps,1);
+md.cm_jump=0.3*ones(md.nsteps,1);
+md.cm_noisedmp=2*10^-7;
+md.vx_obs=md.vx; md.vy_obs=md.vy;
+
 md.cluster=oshostname();
-md.dt=0;
-md=solve(md,'analysis_type',SteadystateSolutionEnum);
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);