Changeset 13466
- Timestamp:
- 09/27/12 10:51:08 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/Par/SquareShelf.py
r13465 r13466 15 15 ymax=max(md.mesh.y) 16 16 17 md.geometry.thickness=hmax+(hmin-hmax)*(md.mesh.y -ymin)/(ymax-ymin)17 md.geometry.thickness=hmax+(hmin-hmax)*(md.mesh.y.reshape(-1,1)-ymin)/(ymax-ymin) 18 18 md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness 19 19 md.geometry.surface=md.geometry.bed+md.geometry.thickness … … 39 39 [md.initialization.vx]=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y) 40 40 [md.initialization.vy]=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y) 41 md.initialization.vz=zeros( md.mesh.numberofvertices)42 md.initialization.pressure=zeros( md.mesh.numberofvertices)41 md.initialization.vz=zeros((md.mesh.numberofvertices,1)) 42 md.initialization.pressure=zeros((md.mesh.numberofvertices,1)) 43 43 44 44 #dbg - begin … … 55 55 56 56 #Materials 57 md.initialization.temperature = (273.-20.)*ones( md.mesh.numberofvertices)57 md.initialization.temperature = (273.-20.)*ones((md.mesh.numberofvertices,1)) 58 58 md.materials.rheology_B = paterson(md.initialization.temperature) 59 md.materials.rheology_n = 3.*ones( md.mesh.numberofelements)59 md.materials.rheology_n = 3.*ones((md.mesh.numberofelements,1)) 60 60 61 61 #Friction 62 62 pos=nonzero(md.mask.elementonfloatingice) 63 md.friction.coefficient = 20.*ones( md.mesh.numberofvertices)63 md.friction.coefficient = 20.*ones((md.mesh.numberofvertices,1)) 64 64 md.friction.coefficient[md.mesh.elements[pos,:].astype(int)-1] = 0. 65 md.friction.p = ones( md.mesh.numberofelements)66 md.friction.q = ones( md.mesh.numberofelements)65 md.friction.p = ones((md.mesh.numberofelements,1)) 66 md.friction.q = ones((md.mesh.numberofelements,1)) 67 67 68 68 #Numerical parameters
Note:
See TracChangeset
for help on using the changeset viewer.