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

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);