Changeset 25403

Timestamp:

08/16/20 13:49:15 (5 years ago)

Author:

Mathieu Morlighem

Message:

BUG: temporary fix, we have a test that has nbe = nbv

Location:

issm/trunk-jpl/src/m/classes

Files:

• friction.m (modified) (1 diff)
• matice.m (modified) (1 diff)

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

@@ -82 +82 @@
 
                         WriteData(fid,prefix,'name','md.friction.law','data',1,'format','Integer');
-                        if(size(self.coefficient,1)==md.mesh.numberofelements | size(self.coefficient,1)==md.mesh.numberofelements+1),
+                        if(size(self.coefficient,1)==md.mesh.numberofvertices | size(self.coefficient,1)==md.mesh.numberofvertices+1),
+                                mattype=1; tsl = md.mesh.numberofvertices;
+                        else
                                 mattype=2; tsl = md.mesh.numberofelements;
-                        else
-                                mattype=1; tsl = md.mesh.numberofvertices;
                         end
                         WriteData(fid,prefix,'object',self,'fieldname','coefficient','format','DoubleMat','mattype',mattype,'timeserieslength',tsl+1,'yts',md.constants.yts);

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

@@ -175 +175 @@
                         WriteData(fid,prefix,'object',self,'class','materials','fieldname','mixed_layer_capacity','format','Double');
                         WriteData(fid,prefix,'object',self,'class','materials','fieldname','thermal_exchange_velocity','format','Double');
-                        if(size(self.rheology_B,1)==md.mesh.numberofelements | size(self.rheology_B,1)==md.mesh.numberofelements+1),
+                        if(size(self.rheology_B,1)==md.mesh.numberofvertices | size(self.rheology_B,1)==md.mesh.numberofvertices+1),
+                                mattype=1; tsl = md.mesh.numberofvertices;
+                        else
                                 mattype=2; tsl = md.mesh.numberofelements;
-                        else
-                                mattype=1; tsl = md.mesh.numberofvertices;
-                        end
-                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',mattype,'timeserieslength',tsl+mattype,'yts',md.constants.yts);
+                        end
+                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',mattype,'timeserieslength',tsl+1,'yts',md.constants.yts);
                         WriteData(fid,prefix,'object',self,'class','materials','fieldname','rheology_n','format','DoubleMat','mattype',2);
                         WriteData(fid,prefix,'data',self.rheology_law,'name','md.materials.rheology_law','format','String');