Changeset 13975 for issm/trunk/test/NightlyRun/test513.m
- Timestamp:
- 11/16/12 08:10:16 (12 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 13397-13398,13401,13407-13582,13584-13974
- Property svn:mergeinfo changed
-
issm/trunk/test
- Property svn:mergeinfo changed
-
Property svn:ignore
set to
AdTests
-
issm/trunk/test/NightlyRun/test513.m
r13395 r13975 1 md=triangle(model ,'../Exp/Pig.exp',30000);1 md=triangle(model(),'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,3,1 );4 md=extrude(md,3,1.); 5 5 md=setflowequation(md,'macayeal','all'); 6 6 … … 8 8 md.inversion.iscontrol=1; 9 9 md.inversion.control_parameters={'FrictionCoefficient'}; 10 md.inversion.min_parameters=1 *ones(md.mesh.numberofvertices,1);11 md.inversion.max_parameters=200 *ones(md.mesh.numberofvertices,1);10 md.inversion.min_parameters=1.*ones(md.mesh.numberofvertices,1); 11 md.inversion.max_parameters=200.*ones(md.mesh.numberofvertices,1); 12 12 md.inversion.nsteps=2; 13 13 md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1) 501*ones(md.inversion.nsteps,1)]; 14 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2); md.inversion.cost_functions_coefficients(:,2)=2 *10^-7;15 md.inversion.gradient_scaling=3 *ones(md.inversion.nsteps,1);16 md.inversion.maxiter_per_step=2 *ones(md.inversion.nsteps,1);14 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2); md.inversion.cost_functions_coefficients(:,2)=2.*10^-7; 15 md.inversion.gradient_scaling=3.*ones(md.inversion.nsteps,1); 16 md.inversion.maxiter_per_step=2.*ones(md.inversion.nsteps,1); 17 17 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1); 18 18 md.inversion.vx_obs=md.initialization.vx; md.inversion.vy_obs=md.initialization.vy; 19 md.timestepping.time_step=0 ;19 md.timestepping.time_step=0.; 20 20 21 21 md.thermal.penalty_lock=5; 22 22 md.cluster=generic('name',oshostname(),'np',3); 23 md=solve(md,SteadystateSolutionEnum );23 md=solve(md,SteadystateSolutionEnum()); 24 24 25 25 %Fields and tolerances to track changes 26 field_names ={'Gradient' 'Misfits' 'FrictionCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature''BasalforcingsMeltingRate'};26 field_names ={'Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy','Vz','Temperature','BasalforcingsMeltingRate'}; 27 27 field_tolerances={1e-10,1e-10,1e-10,1e-10,1e-6,1e-6,1e-6,1e-6,1e-6,1e-06}; 28 28 field_values={...
Note:
See TracChangeset
for help on using the changeset viewer.