Changeset 17022

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

Location:

issm/trunk-jpl/src/m/classes

Files:

• hydrologydc.m (modified) (3 diffs)
• model.m (modified) (1 diff)

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,     

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

@@ -726 +726 @@
                                 md.hydrology.spcepl_head=project3d(md,'vector',md.hydrology.spcepl_head,'type','node','layer',1);
                                 md.hydrology.mask_eplactive_node=project3d(md,'vector',md.hydrology.mask_eplactive_node,'type','node','layer',1);
+                                md.hydrology.sediment_transmitivity=project3d(md,'vector',md.hydrology.sediment_transmitivity,'type','node','layer',1);
             end