Changeset 12124

Timestamp:

04/25/12 09:17:55 (13 years ago)

Author:

Eric.Larour

Message:

Almost finished translated to python

File:

• issm/trunk-jpl/test/Par/SquareShelfConstrained.py (modified) (4 diffs)

issm/trunk-jpl/test/Par/SquareShelfConstrained.py

@@ -1 +1 @@
 from numpy import *
+from verbose import *
 import scipy.io as matio
 import InterpFromMeshToMesh2d as im
+from   paterson import  *
 
 #Start defining model parameters here
@@ -14 +16 @@
 md.geometry.surface = md.geometry.bed+md.geometry.thickness
 
-#%Initial velocity 
+#Initial velocity 
 mat=matio.loadmat('../Data/SquareShelfConstrained.data')
 #deal with 'F' oriented matlab matrices!
@@ -24 +26 @@
 md.initialization.vz = zeros(md.mesh.numberofvertices)
 md.initialization.pressure = zeros(md.mesh.numberofvertices)
-#%Materials
+#Materials
 md.initialization.temperature = (273.-20)*ones(md.mesh.numberofvertices)
 md.materials.rheology_B = paterson(md.initialization.temperature)
 md.materials.rheology_n = 3.*ones(md.mesh.numberofelements)
-#%Surface mass balance and basal melting
+#Surface mass balance and basal melting
 md.surfaceforcings.mass_balance = 10.*ones(md.mesh.numberofvertices)
 md.basalforcings.melting_rate = 5.*ones(md.mesh.numberofvertices)
-#%Friction
+#Friction
 pos = nonzero(md.mask.elementonfloatingice)
 md.friction.coefficient = 20.*ones(md.mesh.numberofvertices)
-md.friction.coefficient[(md.mesh.elements[int(pos)-1,:])] = 0.
+md.friction.coefficient[md.mesh.elements[pos,:].astype(int)-1] =0.
 md.friction.p = ones(md.mesh.numberofelements)
 md.friction.q = ones(md.mesh.numberofelements)
-#%Numerical parameters
+#Numerical parameters
 md.diagnostic.viscosity_overshoot = 0.0
 md.prognostic.stabilization = 1.
 md.thermal.stabilization = 1.
-md.verbose = verbose[-1]
+md.verbose = verbose()
 md.settings.waitonlock = 30.
 md.diagnostic.restol = 0.05
@@ -49 +51 @@
 md.timestepping.time_step = 1.
 md.timestepping.final_time = 3.
-#%Deal with boundary conditions:
+#Deal with boundary conditions:
 md = SetIceShelfBC(md)
-#%Change name so that no test have the same name
+#Change name so that no test have the same name
 A = dbstack
 if length(A) > 2.: