Changeset 13665
- Timestamp:
- 10/15/12 11:41:07 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test512.py
r13664 r13665 8 8 from solve import * 9 9 from MatlabFuncs import * 10 from triangle import *11 from setmask import *12 from parameterize import *13 from setflowequation import *14 from solve import *15 10 16 md=triangle(model ,'../Exp/Pig.exp',20000.)11 md=triangle(model(),'../Exp/Pig.exp',20000.) 17 12 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 18 13 md=parameterize(md,'../Par/Pig.py') … … 20 15 md=setflowequation(md,'pattyn','all') 21 16 22 23 17 # control parameters 24 25 18 md.inversion.iscontrol=1 26 19 md.inversion.control_parameters=['FrictionCoefficient'] 27 md.inversion.min_parameters=1 *ones(md.mesh.numberofvertices,1)28 md.inversion.max_parameters=200 *ones(md.mesh.numberofvertices,1)20 md.inversion.min_parameters=1.*numpy.ones((md.mesh.numberofvertices,1)) 21 md.inversion.max_parameters=200.*numpy.ones((md.mesh.numberofvertices,1)) 29 22 md.inversion.nsteps=2 30 md.inversion.cost_functions= [103*ones(md.inversion.nsteps,1) 501*ones(md.inversion.nsteps,1)]31 md.inversion.cost_functions_coefficients= ones(md.mesh.numberofvertices,2)32 md.inversion.cost_functions_coefficients(:,2)=2.*10^-733 md.inversion.gradient_scaling=3.* ones(md.inversion.nsteps,1)34 md.inversion.maxiter_per_step=2 *ones(md.inversion.nsteps,1)35 md.inversion.step_threshold= 0.99*ones(md.inversion.nsteps,1)23 md.inversion.cost_functions=numpy.hstack((103*numpy.ones((md.inversion.nsteps,1)),501*numpy.ones((md.inversion.nsteps,1)))) 24 md.inversion.cost_functions_coefficients=numpy.ones((md.mesh.numberofvertices,2)) 25 md.inversion.cost_functions_coefficients[:,1]=2.*10**-7 26 md.inversion.gradient_scaling=3.*numpy.ones((md.inversion.nsteps,1)) 27 md.inversion.maxiter_per_step=2.*numpy.ones((md.inversion.nsteps,1)) 28 md.inversion.step_threshold=2.99*numpy.ones((md.inversion.nsteps,1)) 36 29 md.inversion.vx_obs=md.initialization.vx 37 md.inversion.vy_obs=md.initialization.vy 38 30 md.inversion.vy_obs=md.initialization.vy 39 31 40 32 md.cluster=generic('name',oshostname(),'np',3) 41 33 md=solve(md,DiagnosticSolutionEnum()) 42 34 43 44 35 # Fields and tolerances to track changes 45 46 36 field_names =['Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy'] 47 37 field_tolerances=[1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11]
Note:
See TracChangeset
for help on using the changeset viewer.