Changeset 21408 for issm/trunk-jpl/test/NightlyRun/test341.py
- Timestamp:
- 11/22/16 02:31:19 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test341.py
r21060 r21408 1 1 #Test Name: SquareSheetConstrainedCMm1qn3DragHO 2 import numpy 2 import numpy as np 3 3 from model import * 4 from socket import gethostname 4 5 from m1qn3inversion import * 5 6 from triangle import * … … 8 9 from setflowequation import * 9 10 from solve import * 10 from MatlabFuncs import *11 11 12 md=triangle(model(),'../Exp/Square.exp',200000.); 13 md=setmask(md,'',''); 14 md=parameterize(md,'../Par/SquareSheetConstrained.py'); 12 13 md=triangle(model(),'../Exp/Square.exp',200000.) 14 md=setmask(md,'','') 15 md=parameterize(md,'../Par/SquareSheetConstrained.py') 15 16 md.extrude(3,1.) 16 md=setflowequation(md,'HO','all') ;17 md=setflowequation(md,'HO','all') 17 18 18 19 #control parameters 19 md.inversion=m1qn3inversion(md.inversion) ;20 md.inversion.iscontrol=1 ;20 md.inversion=m1qn3inversion(md.inversion) 21 md.inversion.iscontrol=1 21 22 md.inversion.control_parameters=['FrictionCoefficient'] 22 md.inversion.min_parameters=1.*n umpy.ones((md.mesh.numberofvertices,1))23 md.inversion.max_parameters=200.*n umpy.ones((md.mesh.numberofvertices,1))23 md.inversion.min_parameters=1.*np.ones((md.mesh.numberofvertices)) 24 md.inversion.max_parameters=200.*np.ones((md.mesh.numberofvertices)) 24 25 md.inversion.maxsteps=2 25 26 md.inversion.maxiter=6 26 27 md.inversion.cost_functions=[102,501] 27 md.inversion.cost_functions_coefficients=n umpy.ones((md.mesh.numberofvertices,2))28 md.inversion.cost_functions_coefficients=np.ones((md.mesh.numberofvertices,2)) 28 29 md.inversion.cost_functions_coefficients[:,1]=2.*10**-7 29 md.inversion.vx_obs=md.initialization.vx ;30 md.inversion.vy_obs=md.initialization.vy ;30 md.inversion.vx_obs=md.initialization.vx 31 md.inversion.vy_obs=md.initialization.vy 31 32 32 md.cluster=generic('name', oshostname(),'np',3);33 md=solve(md,'Stressbalance') ;33 md.cluster=generic('name',gethostname(),'np',3) 34 md=solve(md,'Stressbalance') 34 35 35 36 #Fields and tolerances to track changes
Note:
See TracChangeset
for help on using the changeset viewer.