Changeset 23793 for issm/trunk-jpl/test/NightlyRun/test321.py

Timestamp:

03/13/19 03:17:46 (6 years ago)

Author:

bdef

Message:

pep8 compliance of NTs

File:

• issm/trunk-jpl/test/NightlyRun/test321.py (modified) (1 diff)

issm/trunk-jpl/test/NightlyRun/test321.py

@@ -10 +10 @@
 
 
-md=triangle(model(),'../Exp/Square.exp',200000.)
-md=setmask(md,'','')
-md=parameterize(md,'../Par/SquareSheetConstrained.py')
-md.extrude(3,1.)
-md=setflowequation(md,'HO','all')
+md = triangle(model(), '../Exp/Square.exp', 200000.)
+md = setmask(md, '', '')
+md = parameterize(md, '../Par/SquareSheetConstrained.py')
+md.extrude(3, 1.)
+md = setflowequation(md, 'HO', 'all')
 
 #control parameters
-md.inversion.iscontrol=1
-md.inversion.control_parameters=['FrictionCoefficient']
-md.inversion.min_parameters=1.*np.ones((md.mesh.numberofvertices,len(md.inversion.control_parameters)))
-md.inversion.max_parameters=200.*np.ones((md.mesh.numberofvertices,len(md.inversion.control_parameters)))
-md.inversion.nsteps=2
-md.inversion.cost_functions=[102,501]
-md.inversion.cost_functions_coefficients=np.ones((md.mesh.numberofvertices,2))
-md.inversion.cost_functions_coefficients[:,1]=2*10**-7
-md.inversion.gradient_scaling=3.*np.ones((md.inversion.nsteps,len(md.inversion.control_parameters)))
-md.inversion.maxiter_per_step=2*np.ones(md.inversion.nsteps)
-md.inversion.step_threshold=0.3*np.ones(md.inversion.nsteps)
-md.inversion.vx_obs=md.initialization.vx
-md.inversion.vy_obs=md.initialization.vy
+md.inversion.iscontrol = 1
+md.inversion.control_parameters = ['FrictionCoefficient']
+md.inversion.min_parameters = 1. * np.ones((md.mesh.numberofvertices, len(md.inversion.control_parameters)))
+md.inversion.max_parameters = 200. * np.ones((md.mesh.numberofvertices, len(md.inversion.control_parameters)))
+md.inversion.nsteps = 2
+md.inversion.cost_functions = [102, 501]
+md.inversion.cost_functions_coefficients = np.ones((md.mesh.numberofvertices, 2))
+md.inversion.cost_functions_coefficients[:, 1] = 2 * 10**-7
+md.inversion.gradient_scaling = 3. * np.ones((md.inversion.nsteps, len(md.inversion.control_parameters)))
+md.inversion.maxiter_per_step = 2 * np.ones(md.inversion.nsteps)
+md.inversion.step_threshold = 0.3 * np.ones(md.inversion.nsteps)
+md.inversion.vx_obs = md.initialization.vx
+md.inversion.vy_obs = md.initialization.vy
 
-md.cluster=generic('name',gethostname(),'np',3)
-md=solve(md,'Stressbalance')
+md.cluster = generic('name', gethostname(), 'np', 3)
+md = solve(md, 'Stressbalance')
 
 #Fields and tolerances to track changes
-field_names     =['Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy']
-field_tolerances=[1e-08,1e-07,1e-10,1e-10,1e-09,1e-09,1e-09]
-field_values=[\
-        md.results.StressbalanceSolution.Gradient1,\
-        md.results.StressbalanceSolution.J,\
-        md.results.StressbalanceSolution.FrictionCoefficient,\
-        md.results.StressbalanceSolution.Pressure,\
-        md.results.StressbalanceSolution.Vel,\
-        md.results.StressbalanceSolution.Vx,\
-        md.results.StressbalanceSolution.Vy,\
-]
+field_names = ['Gradient', 'Misfits', 'FrictionCoefficient', 'Pressure', 'Vel', 'Vx', 'Vy']
+field_tolerances = [1e-08, 1e-07, 1e-10, 1e-10, 1e-09, 1e-09, 1e-09]
+field_values = [md.results.StressbalanceSolution.Gradient1,
+                md.results.StressbalanceSolution.J,
+                md.results.StressbalanceSolution.FrictionCoefficient,
+                md.results.StressbalanceSolution.Pressure,
+                md.results.StressbalanceSolution.Vel,
+                md.results.StressbalanceSolution.Vx,
+                md.results.StressbalanceSolution.Vy]