Changeset 8318

Timestamp:

05/16/11 16:56:15 (14 years ago)

Author:

seroussi

Message:

extrude only if not nan

File:

• issm/trunk/src/m/model/extrude.m (modified) (1 diff)

issm/trunk/src/m/model/extrude.m

@@ -153 +153 @@
 md.dhdt=project3d(md,md.dhdt,'node');
 md.firn_layer=project3d(md,md.firn_layer,'node',md.numlayers);
-md.forcing_accumulation=[project3d(md,md.forcing_accumulation(1:end-1,:),'node');md.forcing_accumulation(end,:)];
+if ~isnan(md.forcing_accumulation), md.forcing_accumulation=[project3d(md,md.forcing_accumulation(1:end-1,:),'node');md.forcing_accumulation(end,:)]; end
 
 %results