Ignore:
Timestamp:
08/09/11 11:50:24 (14 years ago)
Author:
Mathieu Morlighem
Message:

removed balance velocities (not working) for code cleaning

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/model/ismodelselfconsistent.m

    r9183 r9218  
    573573                        if any(sum(abs(md.vx(md.elements)),2)==0 & sum(abs(md.vy(md.elements)),2)==0)
    574574                                message('model not consistent: at least one triangle has all its vertices with a zero velocity');
    575                         end
    576                         %}}}
    577                 case BalancevelocitiesAnalysisEnum,
    578                         % {{{2
    579                         %VELOCITIES MELTING AND ACCUMULATION
    580                         fields={'vx','vy','basal_melting_rate'};
    581                         checksize(md,fields,[md.numberofnodes 1]);
    582                         checknan(md,fields);
    583 
    584                         %SPC
    585                         if any((~isnan(md.spcvx(find(md.nodeonboundary))))~=1),
    586                                 message(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcvx']);
    587                         end
    588 
    589                         if any((~isnan(md.spcvy(find(md.nodeonboundary))))~=1),
    590                                 message(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcvy']);
    591575                        end
    592576                        %}}}
Note: See TracChangeset for help on using the changeset viewer.