Changeset 24764

Timestamp:

04/30/20 12:37:30 (5 years ago)

Author:

Eric.Larour

Message:

CHG: fixed issues with mesh generatioon.

File:

• issm/trunk-jpl/test/SandBox/test2004.m (modified) (8 diffs)

issm/trunk-jpl/test/SandBox/test2004.m

@@ -8 +8 @@
 
 skipmesh=1;
+skipparameterization=0;
 
 %create sealevel model to hold our information: 
@@ -91 +92 @@
 %}}}
 end
+if ~skipparameterization,
 %Parameterize ice sheets : {{{
 
@@ -180 +182 @@
 %}}}
 % ParameterizeContinents {{{
+
+sl.basinindx('continent',{'hemisphereeast','hemispherewest'})
 
 for ind=sl.basinindx('continent',{'hemisphereeast','hemispherewest'}),
@@ -240 +244 @@
         %grounded ice: 
         md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1); 
+        
+        md.mask.glacier_levelset=zeros(md.mesh.numberofvertices,1);
 
         % }}}
@@ -269 +275 @@
 end
 % }}}
-
+end
 %Assemble Earth in 3D {{{
 
@@ -279 +285 @@
 %create earth model by concatenating all the icecaps in 3d: 
 sl.caticecaps('tolerance',tolerance,'loneedgesdetect',loneedgesdetect);
-
-error;
 
 %figure out how each icecap's mesh connects to the larger earth mesh: 
@@ -299 +303 @@
 sl.transfer('mask.land_levelset');
 sl.transfer('mask.glacier_levelset');
-sl.transfer('geometry.thickness');
-sl.transfer('geometry.surface');
-sl.transfer('geometry.base');
 sl.transfer('geometry.bed');
-sl.transfer('initialization.temperature');
-sl.transfer('materials.rheology_B');
 sl.transfer('mesh.lat');
 sl.transfer('mesh.long');
@@ -310 +309 @@
 sl.transfer('slr.Ngia');
 sl.transfer('slr.Ugia');
+error;
 
 %radius: