Ignore:
Timestamp:
09/08/11 08:29:34 (14 years ago)
Author:
Mathieu Morlighem
Message:

Added inversion

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/test625.m

    r9664 r9681  
    1111
    1212%control parameters
    13 md.thickness_obs=md.thickness;
    14 md.control_analysis=1;
    15 md.thickness_obs=md.thickness;
    16 md.nsteps=2;
    17 md.control_type={'Vx' 'Vy'};
     13md.inversion.thickness_obs=md.thickness;
     14md.inversion.iscontrol=1;
     15md.inversion.thickness_obs=md.thickness;
     16md.inversion.nsteps=2;
     17md.inversion.control_parameters={'Vx' 'Vy'};
    1818md.prognostic.stabilization=1;
    19 md.optscal=[10/md.constants.yts*ones(md.nsteps,1) 10/md.constants.yts*ones(md.nsteps,1)];
    20 md.cm_min=[-2000*ones(md.numberofnodes,1) -2000*ones(md.numberofnodes,1)];
    21 md.cm_max=[+2000*ones(md.numberofnodes,1) +2000*ones(md.numberofnodes,1)];
    22 md.cm_responses=201*ones(md.nsteps,1);
    23 md.weights=ones(md.numberofnodes,1);
    24 md.maxiter=4*ones(md.nsteps,1);
    25 md.cm_jump=0.99*ones(md.nsteps,1);
     19md.inversion.gradient_scaling=[10/md.constants.yts*ones(md.inversion.nsteps,1) 10/md.constants.yts*ones(md.inversion.nsteps,1)];
     20md.inversion.min_parameters=[-2000*ones(md.numberofnodes,1) -2000*ones(md.numberofnodes,1)];
     21md.inversion.max_parameters=[+2000*ones(md.numberofnodes,1) +2000*ones(md.numberofnodes,1)];
     22md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
     23md.inversion.cost_functions_coefficients=ones(md.numberofnodes,1);
     24md.inversion.maxiter_per_step=4*ones(md.inversion.nsteps,1);
     25md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
    2626
    2727md.cluster=none;
Note: See TracChangeset for help on using the changeset viewer.