Changeset 17022 for issm/trunk-jpl/src/m/classes/hydrologydc.m

Timestamp:

12/16/13 10:31:22 (11 years ago)

Author:

bdef

Message:

CHG:the sediment transmitivity is now a mat and no more a constant

File:

• issm/trunk-jpl/src/m/classes/hydrologydc.m (modified) (3 diffs)

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

@@ -17 +17 @@
 
                 spcsediment_head         = NaN;
+                sediment_transmitivity   = NaN;
                 sediment_compressibility = 0;
                 sediment_porosity        = 0;
                 sediment_thickness       = 0;
-                sediment_transmitivity   = 0;
+
 
                 spcepl_head              = NaN;
@@ -89 +90 @@
                         md = checkfield(md,'fieldname','hydrology.sediment_porosity','>',0,'numel',1);
                         md = checkfield(md,'fieldname','hydrology.sediment_thickness','>',0,'numel',1);
-                        md = checkfield(md,'fieldname','hydrology.sediment_transmitivity','>',0,'numel',1);
+                        md = checkfield(md,'fieldname','hydrology.sediment_transmitivity','>=',0,'size',[md.mesh.numberofvertices 1]);
 
                         if obj.isefficientlayer==1,
@@ -164 +165 @@
                         WriteData(fid,'object',obj,'fieldname','sediment_porosity','format','Double');                  
                         WriteData(fid,'object',obj,'fieldname','sediment_thickness','format','Double');
-                        WriteData(fid,'object',obj,'fieldname','sediment_transmitivity','format','Double');             
+                        WriteData(fid,'object',obj,'fieldname','sediment_transmitivity','format','DoubleMat','mattype',1');             
 
                         if obj.isefficientlayer==1,