Changeset 25400 for issm/trunk-jpl/src/m/classes/matice.m

Timestamp:

08/15/20 13:34:33 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: making friction and rheology B more flexible:

File:

• issm/trunk-jpl/src/m/classes/matice.m (modified) (1 diff)

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');
-                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
+                        if(size(self.rheology_B,1)==md.mesh.numberofelements | size(self.rheology_B,1)==md.mesh.numberofelements+1),
+                                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);
                         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');