Ignore:
Timestamp:
11/22/16 02:31:19 (8 years ago)
Author:
bdef
Message:

CHG: uniformization fix

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test341.py

    r21060 r21408  
    11#Test Name: SquareSheetConstrainedCMm1qn3DragHO
    2 import numpy
     2import numpy as np
    33from model import *
     4from socket import gethostname
    45from m1qn3inversion import *
    56from triangle import *
     
    89from setflowequation import *
    910from solve import *
    10 from MatlabFuncs import *
    1111
    12 md=triangle(model(),'../Exp/Square.exp',200000.);
    13 md=setmask(md,'','');
    14 md=parameterize(md,'../Par/SquareSheetConstrained.py');
     12
     13md=triangle(model(),'../Exp/Square.exp',200000.)
     14md=setmask(md,'','')
     15md=parameterize(md,'../Par/SquareSheetConstrained.py')
    1516md.extrude(3,1.)
    16 md=setflowequation(md,'HO','all');
     17md=setflowequation(md,'HO','all')
    1718
    1819#control parameters
    19 md.inversion=m1qn3inversion(md.inversion);
    20 md.inversion.iscontrol=1;
     20md.inversion=m1qn3inversion(md.inversion)
     21md.inversion.iscontrol=1
    2122md.inversion.control_parameters=['FrictionCoefficient']
    22 md.inversion.min_parameters=1.*numpy.ones((md.mesh.numberofvertices,1))
    23 md.inversion.max_parameters=200.*numpy.ones((md.mesh.numberofvertices,1))
     23md.inversion.min_parameters=1.*np.ones((md.mesh.numberofvertices))
     24md.inversion.max_parameters=200.*np.ones((md.mesh.numberofvertices))
    2425md.inversion.maxsteps=2
    2526md.inversion.maxiter=6
    2627md.inversion.cost_functions=[102,501]
    27 md.inversion.cost_functions_coefficients=numpy.ones((md.mesh.numberofvertices,2))
     28md.inversion.cost_functions_coefficients=np.ones((md.mesh.numberofvertices,2))
    2829md.inversion.cost_functions_coefficients[:,1]=2.*10**-7
    29 md.inversion.vx_obs=md.initialization.vx;
    30 md.inversion.vy_obs=md.initialization.vy;
     30md.inversion.vx_obs=md.initialization.vx
     31md.inversion.vy_obs=md.initialization.vy
    3132
    32 md.cluster=generic('name',oshostname(),'np',3);
    33 md=solve(md,'Stressbalance');
     33md.cluster=generic('name',gethostname(),'np',3)
     34md=solve(md,'Stressbalance')
    3435
    3536#Fields and tolerances to track changes
Note: See TracChangeset for help on using the changeset viewer.