source: issm/oecreview/Archive/19101-20495/ISSM-19836-19837.diff@ 20498

../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.js

@@ -17 +17 @@
         //node on Dirichlet (boundary and ~icefront)
         if (arguments.length==2){
                 icefront=arguments[1];
-                nodeinsideicefront=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,icefrontfile,'node',2);
+                nodeinsideicefront=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,icefront,'node',2);
                 nodeonicefront=ArrayAnd(md.mesh.vertexonboundary,nodeinsideicefront);
         }
         else if(arguments.length==1){
@@ -51 +51 @@
         var obs=false;
         if((md.inversion.vx_obs.length == md.mesh.numberofvertices) & (md.inversion.vy_obs.length==md.mesh.numberofvertices))obs=true;
 
-
         if(obs==true){
                 console.log('      boundary conditions for stressbalance model: setting spc as observed velocities');
         }
@@ -61 +60 @@
         for(var i=0;i<md.mesh.segments.length;i++){
                 var sum=0;
                 for (var j=0;j<numbernodesfront;j++) sum+=md.mask.ice_levelset[md.mesh.segments[i][j]-1];
-                if(sum!=numbernodesfront){
+                if(sum!=0){
                         for (var j=0;j<numbernodesfront;j++){
                                 if(obs==false){
                                         md.stressbalance.spcvx[md.mesh.segments[i][j]-1]=0;