source: issm/branches/trunk-larour-NatGeoScience2016/test/Par/ISMIPD.py@ 21759

Last change on this file since 21759 was 21759, checked in by Eric.Larour, 8 years ago

CHG: merged branch back to trunk-jpl 21754.
  • Property svn:executable set to *
File size: 975 bytes
Line 
1import numpy
2from SetIceSheetBC import SetIceSheetBC
3
4#Ok, start defining model parameters here
5
6print " creating thickness"
7md.geometry.surface=2000.-md.mesh.x*numpy.tan(0.1*numpy.pi/180.) #to have z>0
8md.geometry.base=md.geometry.surface-1000.
9md.geometry.thickness=md.geometry.surface-md.geometry.base
10
11print " creating drag"
12md.friction.coefficient=numpy.sqrt(md.constants.yts*(1000.+1000.*numpy.sin(md.mesh.x*2.*numpy.pi/numpy.max(md.mesh.x))))
13md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
14md.friction.p=numpy.ones((md.mesh.numberofelements))
15md.friction.q=numpy.zeros((md.mesh.numberofelements))
16
17print " creating flow law parameter"
18md.materials.rheology_B=6.8067*10**7*numpy.ones((md.mesh.numberofvertices))
19md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements))
20
21print " boundary conditions for stressbalance model:"
22#Create node on boundary first (because we can not use mesh)
23md=SetIceSheetBC(md)
Note: See TracBrowser for help on using the repository browser.