Changeset 16855
- Timestamp:
- 11/21/13 11:34:49 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/hydrologydc.m
r16764 r16855 23 23 24 24 spcepl_head = NaN; 25 mask_eplactive 25 mask_eplactive_node = NaN; 26 26 epl_compressibility = 0; 27 27 epl_porosity = 0; … … 93 93 if obj.isefficientlayer==1, 94 94 md = checkfield(md,'fieldname','hydrology.spcepl_head','forcing',1); 95 md = checkfield(md,'fieldname','hydrology.mask_eplactive ','size',[md.mesh.numberofvertices 1],'values',[0 1]);95 md = checkfield(md,'fieldname','hydrology.mask_eplactive_node','size',[md.mesh.numberofvertices 1],'values',[0 1]); 96 96 md = checkfield(md,'fieldname','hydrology.epl_compressibility','>',0,'numel',1); 97 97 md = checkfield(md,'fieldname','hydrology.epl_porosity','>',0,'numel',1); … … 134 134 disp(sprintf(' - for the epl layer')); 135 135 fielddisplay(obj,'spcepl_head','epl water head constraints (NaN means no constraint) [m above MSL]'); 136 fielddisplay(obj,'mask_eplactive ','active (1) or not (0) EPL');136 fielddisplay(obj,'mask_eplactive_node','active (1) or not (0) EPL'); 137 137 fielddisplay(obj,'epl_compressibility','epl compressibility [Pa^-1]'); 138 138 fielddisplay(obj,'epl_porosity','epl [dimensionless]'); … … 168 168 if obj.isefficientlayer==1, 169 169 WriteData(fid,'object',obj,'fieldname','spcepl_head','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1); 170 WriteData(fid,'object',obj,'fieldname','mask_eplactive ','format','DoubleMat','mattype',1);170 WriteData(fid,'object',obj,'fieldname','mask_eplactive_node','format','DoubleMat','mattype',1); 171 171 WriteData(fid,'object',obj,'fieldname','epl_compressibility','format','Double'); 172 172 WriteData(fid,'object',obj,'fieldname','epl_porosity','format','Double');
Note:
See TracChangeset
for help on using the changeset viewer.