Changeset 24711 for issm/trunk-jpl/test/NightlyRun/test701.py

Timestamp:

04/16/20 02:44:38 (5 years ago)

Author:

bdef

Message:

NEW:Finishing implementation of substeps fro 3D

File:

• issm/trunk-jpl/test/NightlyRun/test701.py (modified) (4 diffs)

issm/trunk-jpl/test/NightlyRun/test701.py

@@ -13 +13 @@
 md = bamgflowband(model(), x, b + h, b, 'hmax', 80.)
 
-print(isinstance(md, model))
-
 #Geometry  #interp1d returns a function to be called on md.mesh.x
 md.geometry.surface = np.interp(md.mesh.x, x, b + h)
@@ -21 +19 @@
 
 #mask
-md.mask.ice_levelset = -np.ones((md.mesh.numberofvertices, ))
+md.mask.ice_levelset = -np.ones((md.mesh.numberofvertices))
 md.mask.ice_levelset[np.where(md.mesh.vertexflags(2))] = 0.
-md.mask.groundedice_levelset = np.zeros((md.mesh.numberofvertices, )) - 0.5
+md.mask.groundedice_levelset = np.zeros((md.mesh.numberofvertices)) - 0.5
 
 #materials
-md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices, ))
+md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices))
 md.materials.rheology_B = paterson(md.initialization.temperature)
-md.materials.rheology_n = 3. * np.ones((md.mesh.numberofelements, ))
+md.materials.rheology_n = 3. * np.ones((md.mesh.numberofelements))
 
 #friction
-md.friction.coefficient = np.zeros((md.mesh.numberofvertices, ))
+md.friction.coefficient = np.zeros((md.mesh.numberofvertices))
 md.friction.coefficient[np.where(md.mesh.vertexflags(1))] = 20.
-md.friction.p = np.ones((md.mesh.numberofelements, ))
-md.friction.q = np.ones((md.mesh.numberofelements, ))
+md.friction.p = np.ones((md.mesh.numberofelements))
+md.friction.q = np.ones((md.mesh.numberofelements))
 
 #Boundary conditions
@@ -47 +45 @@
 
 #Misc
-print(isinstance(md, model))
-print(type(md))
 md = setflowequation(md, 'FS', 'all')
 md.stressbalance.abstol = np.nan
@@ -67 +63 @@
 for i in ['MINI', 'MINIcondensed', 'TaylorHood', 'LATaylorHood', 'CrouzeixRaviart', 'LACrouzeixRaviart']:
     print(' ')
-    print(' == == == Testing ' + i + ' Full - Stokes Finite element == == = ')
+    print('=====Testing ' + i + ' Full-Stokes Finite element=====')
     md.flowequation.fe_FS = i
     md = solve(md, 'Stressbalance')