Changeset 25309


Ignore:
Timestamp:
07/29/20 02:54:40 (5 years ago)
Author:
bdef
Message:

BUG: fixing failures (not errors)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/Par/SquareShelf.py

    r25307 r25309  
    4545[md.initialization.vy] = InterpFromMeshToMesh2d(index, x, y, vy, md.mesh.x, md.mesh.y)
    4646x = y = vx = vy = index = None
    47 md.initialization.vz = np.zeros((md.mesh.numberofvertices, 1))
    48 md.initialization.pressure = np.zeros((md.mesh.numberofvertices, 1))
     47md.initialization.vz = np.zeros((md.mesh.numberofvertices))
     48md.initialization.pressure = np.zeros((md.mesh.numberofvertices))
    4949
    5050#dbg - begin
     
    6060
    6161# Materials
    62 md.initialization.temperature = (273 - 20) * np.ones((md.mesh.numberofvertices, 1))
     62md.initialization.temperature = (273 - 20) * np.ones((md.mesh.numberofvertices))
    6363md.materials.rheology_B = paterson(md.initialization.temperature)
    64 md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements, 1))
     64md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements))
    6565
    6666# Friction
    67 md.friction.coefficient = 20 * np.ones((md.mesh.numberofvertices, 1))
     67md.friction.coefficient = 20 * np.ones((md.mesh.numberofvertices))
    6868md.friction.coefficient[np.where(md.mask.ocean_levelset < 0)[0]] = 0
    69 md.friction.p = np.ones((md.mesh.numberofelements, 1))
    70 md.friction.q = np.ones((md.mesh.numberofelements, 1))
     69md.friction.p = np.ones((md.mesh.numberofelements))
     70md.friction.q = np.ones((md.mesh.numberofelements))
    7171
    7272# Numerical parameters
Note: See TracChangeset for help on using the changeset viewer.