Changeset 6239

Timestamp:

10/11/10 14:32:19 (14 years ago)

Author:

seroussi

Message:

problem with geothermal flux

Location:

issm/trunk/src/m/model

Files:

• extrude.m (modified) (1 diff)
• ismodelselfconsistent.m (modified) (1 diff)

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

@@ -232 +232 @@
 md.melting_rate=project3d(md,md.melting_rate,'node',1); 
 md.observed_temperature=project3d(md,md.observed_temperature,'node'); 
-md.geothermalflux=project3d(md,md.geothermalflux,'node',1); %bedrock only gets geothermal flux
+if ~isnan(md.geothermalflux)
+        md.geothermalflux=project3d(md,md.geothermalflux,'node',1); %bedrock only gets geothermal flux
+end
 
 %increase connectivity if less than 25:

issm/trunk/src/m/model/ismodelselfconsistent.m

@@ -276 +276 @@
 
         %VELOCITIES AND PRESSURE
-        fields={'vx','vy','vz','pressure'};
+        fields={'vx','vy','vz','pressure','geothermalflux'};
         checksize(md,fields,[md.numberofgrids 1]);
         checknan(md,fields);