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

Index: ../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.js
===================================================================
--- ../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.js      (revision 19836)
+++ ../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.js      (revision 19837)
@@ -17,7 +17,7 @@
        //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,7 +51,6 @@
        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,7 +60,7 @@
        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;