Changeset 750
- Timestamp:
- 06/03/09 16:36:40 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/classes/public/isresultconsistent.m
r716 r750 146 146 else 147 147 fields1={['results.transient(' num2str(iter) ').vx'],['results.transient(' num2str(iter) ').vy'],... 148 ['results.transient(' num2str(iter) ').v z'],['results.transient(' num2str(iter) ').vel'],...148 ['results.transient(' num2str(iter) ').vel'],... 149 149 ['results.transient(' num2str(iter) ').bed'],['results.transient(' num2str(iter) ').surface'],... 150 150 ['results.transient(' num2str(iter) ').thickness']}; … … 173 173 174 174 %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 178 180 end 179 181 end
Note:
See TracChangeset
for help on using the changeset viewer.