Changeset 19294


Ignore:
Timestamp:
04/21/15 14:03:42 (10 years ago)
Author:
Mathieu Morlighem
Message:

CHG: allow models to have no spc, but display warning

Location:
issm/trunk-jpl/src/m/classes
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/classes/flowequation.m

    r19069 r19294  
    134134                end % }}}
    135135                function md = checkconsistency(self,md,solution,analyses) % {{{
    136 
    137136
    138137                        %Early return
  • issm/trunk-jpl/src/m/classes/stressbalance.m

    r19090 r19294  
    151151                        if ((~(any(~isnan(md.stressbalance.spcvx)) | any(~isnan(md.stressbalance.spcvy)))) & ~any(md.mask.groundedice_levelset>0)),
    152152                                md = checkmessage(md,['model is not well posed (singular). You need at least one node with fixed velocity!']);
     153                                disp(sprintf('\n !!! Warning: no spc applied, model might not be well posed if no basal friction is applied, check for solution crash\n'));
    153154                        end
    154155                        %CHECK THAT EACH LINES CONTAINS ONLY NAN VALUES OR NO NAN VALUES
  • issm/trunk-jpl/src/m/classes/stressbalance.py

    r19048 r19294  
    151151#               if ~any((~isnan(md.stressbalance.spcvx)+~isnan(md.stressbalance.spcvy))==2),
    152152                if not numpy.any(numpy.logical_and(numpy.logical_not(numpy.isnan(md.stressbalance.spcvx)),numpy.logical_not(numpy.isnan(md.stressbalance.spcvy)))):
    153                         md.checkmessage("model is not well posed (singular). You need at least one node with fixed velocity!")
     153                        print "\n !!! Warning: no spc applied, model might not be well posed if no basal friction is applied, check for solution crash\n"
    154154                #CHECK THAT EACH LINES CONTAINS ONLY NAN VALUES OR NO NAN VALUES
    155155#               if any(sum(isnan(md.stressbalance.referential),2)~=0 & sum(isnan(md.stressbalance.referential),2)~=6),
Note: See TracChangeset for help on using the changeset viewer.