Changeset 24214 for issm/trunk-jpl/test/NightlyRun/test442.py
- Timestamp:
- 10/11/19 00:27:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test442.py
r23793 r24214 15 15 md.initialization.vy[:] = 1. 16 16 md.geometry.thickness[:] = 500. - md.mesh.x / 10000. 17 md.geometry.bed = - 100. - md.mesh.x / 1000.18 md.geometry.base = - md.geometry.thickness * md.materials.rho_ice / md.materials.rho_water17 md.geometry.bed = - 100. - md.mesh.x / 1000. 18 md.geometry.base = - md.geometry.thickness * md.materials.rho_ice / md.materials.rho_water 19 19 md.mask.groundedice_levelset = md.geometry.thickness + md.materials.rho_water / md.materials.rho_ice * md.geometry.bed 20 20 pos = np.where(md.mask.groundedice_levelset >= 0.) … … 25 25 26 26 #Boundary conditions: 27 md.mask.ice_levelset = - np.ones((md.mesh.numberofvertices,))27 md.mask.ice_levelset = - np.ones((md.mesh.numberofvertices, )) 28 28 md.mask.ice_levelset[np.where(md.mesh.x == max(md.mesh.x))] = 0. 29 29 md.stressbalance.spcvx[:] = float('Nan') … … 31 31 md.stressbalance.spcvz[:] = float('Nan') 32 32 posA = np.intersect1d(np.array(np.where(md.mesh.y < 1000000.1)), np.array(np.where(md.mesh.y > 999999.9))) 33 posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > - 0.1)))33 posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > - 0.1))) 34 34 pos = np.unique(np.concatenate((posA, posB))) 35 35 md.stressbalance.spcvy[pos] = 0. 36 pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > - 0.1)))36 pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > - 0.1))) 37 37 md.stressbalance.spcvx[pos2] = 0. 38 38 md.stressbalance.spcvy[pos2] = 0. 39 39 40 md.materials.rheology_B = 1. / ((10**- 25)**(1. / 3.)) * np.ones((md.mesh.numberofvertices,))40 md.materials.rheology_B = 1. / ((10**- 25)**(1. / 3.)) * np.ones((md.mesh.numberofvertices, )) 41 41 md.materials.rheology_law = 'None' 42 md.friction.coefficient[:] = np.sqrt(1e7) * np.ones((md.mesh.numberofvertices, ))43 md.friction.p = 3. * np.ones((md.mesh.numberofelements, ))42 md.friction.coefficient[:] = np.sqrt(1e7) * np.ones((md.mesh.numberofvertices, )) 43 md.friction.p = 3. * np.ones((md.mesh.numberofelements, )) 44 44 md.smb.mass_balance[:] = 1. 45 45 md.basalforcings.groundedice_melting_rate[:] = 0.
Note:
See TracChangeset
for help on using the changeset viewer.