Changeset 9002 for issm/trunk/src/m/utils/BC/SetIceShelfBC.m

Timestamp:

07/15/11 08:49:34 (14 years ago)

Author:

Eric.Larour

Message:

Implemented transient boundary conditions. Applied to thermal core only for now

File:

• issm/trunk/src/m/utils/BC/SetIceShelfBC.m (modified) (1 diff)

issm/trunk/src/m/utils/BC/SetIceShelfBC.m

@@ -93 +93 @@
 
 if (length(md.observed_temperature)==md.numberofnodes),
-        md.spctemperature=[md.nodeonsurface md.observed_temperature]; %impose observed temperature on surface
+        md.spctemperature=NaN*ones(md.numberofnodes,1);
+        pos=find(md.nodeonsurface); md.spctemperature(pos)=md.observed_temperature(pos); %impose observed temperature on surface
         if (length(md.geothermalflux)~=md.numberofnodes),
                 md.geothermalflux=zeros(md.numberofnodes,1);