Context Navigation

-              r14310
+              r16137
 #md.friction.coefficient=sqrt(md.constants.yts.*(1000.+1000.*sin(md.mesh.x*2.*pi/max(md.mesh.x/2.)).*sin(md.mesh.y*2.*pi/max(md.mesh.x/2.)))./(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.bed)));
 md.friction.coefficient=numpy.sqrt(md.constants.yts*(1000.+1000.*numpy.sin(md.mesh.x.reshape(-1,1)*2.*numpy.pi/numpy.max(md.mesh.x))*numpy.sin(md.mesh.y.reshape(-1,1)*2.*numpy.pi/numpy.max(md.mesh.x))))
+#Take care of iceshelves: no basal drag
+pos=numpy.nonzero(md.mask.elementonfloatingice)[0]
+md.friction.coefficient[md.mesh.elements[pos,:]-1]=0.
+md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
 md.friction.p=numpy.ones((md.mesh.numberofelements,1))
 md.friction.q=numpy.zeros((md.mesh.numberofelements,1))
 …
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 print "      boundary conditions for diagnostic model:"
+print "      boundary conditions for stressbalance model:"
 #Create node on boundary first (because we can not use mesh)
 md=SetIceSheetBC(md)

Note: See TracChangeset for help on using the changeset viewer.

Changeset 16137 for issm/trunk/test/Par/ISMIPC.py