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
RevLine 
[14134]1import numpy
[16170]2from SetIceSheetBC import SetIceSheetBC
[14134]3
4#Ok, start defining model parameters here
5
6print " creating thickness"
[21759]7md.geometry.surface=2000.-md.mesh.x*numpy.tan(0.1*numpy.pi/180.) #to have z>0
[17590]8md.geometry.base=md.geometry.surface-1000.
9md.geometry.thickness=md.geometry.surface-md.geometry.base
[14134]10
11print " creating drag"
[21759]12md.friction.coefficient=numpy.sqrt(md.constants.yts*(1000.+1000.*numpy.sin(md.mesh.x*2.*numpy.pi/numpy.max(md.mesh.x))))
[15988]13md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
[21759]14md.friction.p=numpy.ones((md.mesh.numberofelements))
15md.friction.q=numpy.zeros((md.mesh.numberofelements))
[14134]16
17print " creating flow law parameter"
[21759]18md.materials.rheology_B=6.8067*10**7*numpy.ones((md.mesh.numberofvertices))
19md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements))
[14134]20
[15771]21print " boundary conditions for stressbalance model:"
[14134]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.