Changeset 1777

Timestamp:

08/19/09 07:35:45 (15 years ago)

Author:

Mathieu Morlighem

Message:

better way to set thermal BC

Location:

issm/trunk/src/m/utils/BC

Files:

• SetIceSheetBC.m (modified) (1 diff)
• SetIceShelfBC.m (modified) (1 diff)
• SetMarineIceSheetBC.m (modified) (1 diff)

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

@@ -38 +38 @@
 if (length(md.observed_temperature)==md.numberofgrids),
         displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
-        md.spctemperature=ones(md.numberofgrids,2); %observed temperauture everywhere
-        md.spctemperature(:,2)=md.observed_temperature; %values of the spcs
+        md.spctemperature=[md.gridonsurface md.observed_temperature]; %impose observed temperature on surface
         if (length(md.geothermalflux)~=md.numberofgrids),
                 md.geothermalflux=50*10^-3*ones(md.numberofgrids,1); %50 mW/m^2

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

@@ -51 +51 @@
 if (length(md.observed_temperature)==md.numberofgrids),
         displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
-        md.spctemperature=ones(md.numberofgrids,2); %surface temperature
-        md.spctemperature(:,2)=md.observed_temperature; %values of the spcs
+        md.spctemperature=[md.gridonsurface md.observed_temperature]; %impose observed temperature on surface
         if (length(md.geothermalflux)~=md.numberofgrids),
                 md.geothermalflux=zeros(md.numberofgrids,1);

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

@@ -63 +63 @@
 if (length(md.observed_temperature)==md.numberofgrids),
         displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
-        md.spctemperature=ones(md.numberofgrids,2); %surface temperature
-        md.spctemperature(:,2)=md.observed_temperature; %values of the spcs
+        md.spctemperature=[md.gridonsurface md.observed_temperature]; %impose observed temperature on surface
         if (length(md.geothermalflux)~=md.numberofgrids),
                 md.geothermalflux=zeros(md.numberofgrids,1);