Changeset 24214 for issm/trunk-jpl/test/NightlyRun/test430.py
- Timestamp:
- 10/11/19 00:27:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test430.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.) … … 24 24 25 25 #Boundary conditions: 26 md.mask.ice_levelset = - np.ones((md.mesh.numberofvertices,))26 md.mask.ice_levelset = - np.ones((md.mesh.numberofvertices, )) 27 27 md.mask.ice_levelset[np.where(md.mesh.x == max(md.mesh.x))] = 0. 28 28 md.stressbalance.spcvx[:] = float('NaN') … … 30 30 md.stressbalance.spcvz[:] = float('NaN') 31 31 posA = np.intersect1d(np.array(np.where(md.mesh.y < 1000000.1)), np.array(np.where(md.mesh.y > 999999.9))) 32 posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > - 0.1)))32 posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > - 0.1))) 33 33 pos = np.unique(np.concatenate((posA, posB))) 34 34 md.stressbalance.spcvy[pos] = 0. 35 pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > - 0.1)))35 pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > - 0.1))) 36 36 md.stressbalance.spcvx[pos2] = 0. 37 37 md.stressbalance.spcvy[pos2] = 0. 38 38 39 md.materials.rheology_B = 1. / ((10**- 25)**(1. / 3.)) * np.ones((md.mesh.numberofvertices,))39 md.materials.rheology_B = 1. / ((10**- 25)**(1. / 3.)) * np.ones((md.mesh.numberofvertices, )) 40 40 md.materials.rheology_law = 'None' 41 md.friction.coefficient[:] = np.sqrt(10**7) * np.ones((md.mesh.numberofvertices, ))42 md.friction.p = 3. * np.ones((md.mesh.numberofelements, ))41 md.friction.coefficient[:] = np.sqrt(10**7) * np.ones((md.mesh.numberofvertices, )) 42 md.friction.p = 3. * np.ones((md.mesh.numberofelements, )) 43 43 md.smb.mass_balance[:] = 1. 44 44 md.basalforcings.groundedice_melting_rate[:] = 0.
Note:
See TracChangeset
for help on using the changeset viewer.