/[escript]/trunk/downunder/test/python/run_magnetic.py
ViewVC logotype

Diff of /trunk/downunder/test/python/run_magnetic.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 4078 by gross, Thu Nov 15 03:30:59 2012 UTC revision 4079 by gross, Fri Nov 16 07:59:01 2012 UTC
# Line 32  except KeyError: Line 32  except KeyError:
32     WORKDIR='.'     WORKDIR='.'
33    
34  features=[SmoothAnomaly(lx=30*U.km, ly=20*U.km, lz=18.*U.km, \  features=[SmoothAnomaly(lx=30*U.km, ly=20*U.km, lz=18.*U.km, \
35       x=8*U.km, y=3*U.km, depth=2.5*U.km, v_inner=200., v_outer=1e-6),\       x=8*U.km, y=3*U.km, depth=2.5*U.km, v_inner=2., v_outer=1e-6),\
36            SmoothAnomaly(lx=25*U.km, ly=20*U.km, lz=20*U.km,            SmoothAnomaly(lx=25*U.km, ly=20*U.km, lz=20*U.km,
37       x=30*U.km, y=1*U.km, depth=18*U.km, v_inner=-100., v_outer=1e-6),\       x=30*U.km, y=1*U.km, depth=18*U.km, v_inner=-1., v_outer=1e-6),\
38            SmoothAnomaly(lx=30*U.km, ly=20*U.km, lz=18.*U.km, \            SmoothAnomaly(lx=30*U.km, ly=20*U.km, lz=18.*U.km, \
39       x=68*U.km, y=3*U.km, depth=5*U.km, v_inner=200., v_outer=1e-6)]       x=68*U.km, y=3*U.km, depth=5*U.km, v_inner=20., v_outer=1e-6)]
40    
41    
42  logger=logging.getLogger('inv')  logger=logging.getLogger('inv')
# Line 51  domainbuilder.setPadding(10) Line 51  domainbuilder.setPadding(10)
51  domainbuilder.setVerticalExtents(depth=30*U.km, air_layer=10*U.km, num_cells=16)  domainbuilder.setVerticalExtents(depth=30*U.km, air_layer=10*U.km, num_cells=16)
52    
53  inv=MagneticInversion()  inv=MagneticInversion()
54  inv.setSolverTolerance(1e-9)  inv.setSolverTolerance(1e-4)
55  inv.setSolverMaxIterations(100)  inv.setSolverMaxIterations(10)
56  inv.setSolverOptions(initialHessian=100)  inv.setSolverOptions(initialHessian=100)
57  #inv.setWeights(mu_reg=1e-4)  #inv.setWeights(mu_reg=1e-4)
58  inv.setup(domainbuilder)  inv.setup(domainbuilder)

Legend:
Removed from v.4078  
changed lines
  Added in v.4079

  ViewVC Help
Powered by ViewVC 1.1.26