Changeset 13635

Timestamp:

10/12/12 10:49:05 (12 years ago)

Author:

jschierm

Message:

NEW: Python version of SetIceSheetBC.

Location:

issm/trunk-jpl/src/m/boundaryconditions

Files:

• SetIceSheetBC.m (modified) (1 diff)
• SetIceSheetBC.py (added)

issm/trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.m

@@ -58 +58 @@
         pos=find(md.mesh.vertexonsurface); md.thermal.spctemperature(pos)=md.initialization.temperature(pos); %impose observed temperature on surface
         if (length(md.basalforcings.geothermalflux)~=md.mesh.numberofvertices),
-                md.basalforcings.geothermalflux=50*10^-3*ones(md.mesh.numberofvertices,1); %50 mW/m^2
+                md.basalforcings.geothermalflux=50.*10^-3*ones(md.mesh.numberofvertices,1); %50 mW/m^2
         end
 else