Changeset 280
- Timestamp:
- 05/06/09 16:05:31 (16 years ago)
- Location:
- issm/trunk/test/Validation/ControlMethods
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/Validation/ControlMethods/Square.par
r161 r280 37 37 38 38 %Parallel options 39 md.np= 8;39 md.np=3; 40 40 md.time=50; 41 41 md.waitonlock=1; -
issm/trunk/test/Validation/ControlMethods/runme.m
r174 r280 1 1 step=[2]; 2 2 division=2; 3 density=80000; 3 density=50000; 4 startup 4 5 5 6 %step1: initial velocity -> observed velocity … … 30 31 %diagnostic with this B 31 32 md=solve(md,'diagnostic','cielo'); 33 %md.cluster='wilkes'; md=solve(md,'diagnostic_horiz','cielo'); 32 34 33 35 %save observed velocities and exact B … … 55 57 56 58 %control method 57 md.nsteps= 3;59 md.nsteps=50; 58 60 md.control_type='B'; 59 61 md.mincontrolconstraint=0; … … 61 63 md.plot=1; 62 64 md.debug=0; 63 md.optscal=10^ 6*ones(md.nsteps,1);65 md.optscal=10^8*ones(md.nsteps,1); 64 66 md.fit=0*ones(md.nsteps,1); 65 md.maxiter= 2*ones(md.nsteps,1);66 md.eps_rel= 0.1;67 md.maxiter=10*ones(md.nsteps,1); 68 md.eps_rel=1; 67 69 md.eps_abs=NaN; 68 70 69 %md.cluster='wilkes';70 md.eps_abs=NaN;71 % md.cluster='wilkes'; 72 % md=solve(md,'control','ice'); 71 73 md=solve(md,'control','cielo'); 72 74 end
Note:
See TracChangeset
for help on using the changeset viewer.