Index: ../trunk-jpl/src/m/classes/model/model.m
===================================================================
--- ../trunk-jpl/src/m/classes/model/model.m	(revision 15020)
+++ ../trunk-jpl/src/m/classes/model/model.m	(revision 15021)
@@ -157,6 +157,10 @@
 			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,8 +212,6 @@
 			%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);