Changeset 750


Ignore:
Timestamp:
06/03/09 16:36:40 (15 years ago)
Author:
Mathieu Morlighem
Message:

minor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/classes/public/isresultconsistent.m

    r716 r750  
    146146                else
    147147                        fields1={['results.transient(' num2str(iter) ').vx'],['results.transient(' num2str(iter) ').vy'],...
    148                                 ['results.transient(' num2str(iter) ').vz'],['results.transient(' num2str(iter) ').vel'],...
     148                                ['results.transient(' num2str(iter) ').vel'],...
    149149                                ['results.transient(' num2str(iter) ').bed'],['results.transient(' num2str(iter) ').surface'],...
    150150                                ['results.transient(' num2str(iter) ').thickness']};
     
    173173
    174174                %check melting (<=0 via penalties)
    175                 if any(abs(md.results.transient(iter).melting(md.numberofgrids2d+1:end))>tolerance)
    176                         disp(['''thermal'' result not consistent: increase penalty_melting (negative melting)']);
    177                         bool=0; return;
     175                if strcmpi(md.type,'3d'),
     176                        if any(abs(md.results.transient(iter).melting(md.numberofgrids2d+1:end))>tolerance)
     177                                disp(['''thermal'' result not consistent: increase penalty_melting (negative melting)']);
     178                                bool=0; return;
     179                        end
    178180                end
    179181        end
Note: See TracChangeset for help on using the changeset viewer.