Changeset 15021

Timestamp:

05/14/13 10:07:08 (12 years ago)

Author:

seroussi

Message:

BUG: fixed problem with gia in collapse

File:

• issm/trunk-jpl/src/m/classes/model/model.m (modified) (2 diffs)

issm/trunk-jpl/src/m/classes/model/model.m

@@ -158 +158 @@
                         if ~isnan(md.initialization.temperature),md.initialization.temperature=DepthAverage(md,md.initialization.temperature);end;
 
+                        %gia
+                        if ~isnan(md.gia.mantle_viscosity), md.gia.mantle_viscosity=project2d(md,md.gia.mantle_viscosity,1); end
+                        if ~isnan(md.gia.lithosphere_thickness), md.gia.lithosphere_thickness=project2d(md,md.gia.lithosphere_thickness,1); end
+
                         %bedinfo and surface info
                         md.mesh.elementonbed=ones(md.mesh.numberofelements2d,1);
@@ -209 +213 @@
                         md.geometry.surface=project2d(md,md.geometry.surface,1);
                         md.geometry.thickness=project2d(md,md.geometry.thickness,1);
-                        md.gia.mantle_viscosity=project2d(md,md.gia.mantle_viscosity,1);
-                        md.gia.lithosphere_thickness=project2d(md,md.gia.lithosphere_thickness,1);
                         md.geometry.bed=project2d(md,md.geometry.bed,1);
                         md.geometry.bathymetry=project2d(md,md.geometry.bathymetry,1);