Ignore:
Timestamp:
05/06/09 16:05:31 (16 years ago)
Author:
Mathieu Morlighem
Message:

minor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/Validation/ControlMethods/runme.m

    r174 r280  
    11step=[2];
    22division=2;
    3 density=80000;
     3density=50000;
     4startup
    45
    56%step1: initial velocity -> observed velocity
     
    3031        %diagnostic with this B
    3132        md=solve(md,'diagnostic','cielo');
     33        %md.cluster='wilkes'; md=solve(md,'diagnostic_horiz','cielo');
    3234
    3335        %save observed velocities and exact B
     
    5557
    5658        %control method
    57         md.nsteps=3;
     59        md.nsteps=50;
    5860        md.control_type='B';
    5961        md.mincontrolconstraint=0;
     
    6163        md.plot=1;
    6264        md.debug=0;
    63         md.optscal=10^6*ones(md.nsteps,1);
     65        md.optscal=10^8*ones(md.nsteps,1);
    6466        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;
    6769        md.eps_abs=NaN;
    6870
    69         %md.cluster='wilkes';
    70         md.eps_abs=NaN;
     71%       md.cluster='wilkes';
     72%       md=solve(md,'control','ice');
    7173        md=solve(md,'control','cielo');
    7274end
Note: See TracChangeset for help on using the changeset viewer.