Changeset 25309
- Timestamp:
- 07/29/20 02:54:40 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/Par/SquareShelf.py
r25307 r25309 45 45 [md.initialization.vy] = InterpFromMeshToMesh2d(index, x, y, vy, md.mesh.x, md.mesh.y) 46 46 x = 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))47 md.initialization.vz = np.zeros((md.mesh.numberofvertices)) 48 md.initialization.pressure = np.zeros((md.mesh.numberofvertices)) 49 49 50 50 #dbg - begin … … 60 60 61 61 # Materials 62 md.initialization.temperature = (273 - 20) * np.ones((md.mesh.numberofvertices , 1))62 md.initialization.temperature = (273 - 20) * np.ones((md.mesh.numberofvertices)) 63 63 md.materials.rheology_B = paterson(md.initialization.temperature) 64 md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements , 1))64 md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements)) 65 65 66 66 # Friction 67 md.friction.coefficient = 20 * np.ones((md.mesh.numberofvertices , 1))67 md.friction.coefficient = 20 * np.ones((md.mesh.numberofvertices)) 68 68 md.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))69 md.friction.p = np.ones((md.mesh.numberofelements)) 70 md.friction.q = np.ones((md.mesh.numberofelements)) 71 71 72 72 # Numerical parameters
Note:
See TracChangeset
for help on using the changeset viewer.