source: issm/oecreview/Archive/14312-15392/ISSM-15020-15021.diff@ 15393

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

@@ -157 +157 @@
                         if ~isnan(md.initialization.vel),md.initialization.vel=DepthAverage(md,md.initialization.vel);end;
                         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);
                         md.mesh.elementonsurface=ones(md.mesh.numberofelements2d,1);
@@ -208 +212 @@
                         %parameters
                         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);
                         md.mesh.vertexonboundary=project2d(md,md.mesh.vertexonboundary,1);