Changeset 28050 for issm/trunk-jpl/src/m/classes/hydrologyglads.m

Timestamp:

01/17/24 17:08:12 (14 months ago)

Author:

sehrenfe

Message:

CHG: GlaDS- new input parameter rheology_B_base and new optional output HydrologySheetDischarge added to hydrology model

File:

• issm/trunk-jpl/src/m/classes/hydrologyglads.m (modified) (4 diffs)

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

@@ -14 +14 @@
                 sheet_alpha               = NaN; 
                 sheet_beta                = NaN; 
+                rheology_B_base           = NaN;
 
                 %Channels
@@ -83 +84 @@
                         md = checkfield(md,'fieldname','hydrology.sheet_alpha', 'numel', [1], '>', 0); 
                         md = checkfield(md,'fieldname','hydrology.sheet_beta', 'numel', [1], '>', 0); 
-
+                        md = checkfield(md,'fieldname','hydrology.rheology_B_base','size',[md.mesh.numberofvertices 1],'>=',0,'NaN',1,'Inf',1);
                         %Channels
                         md = checkfield(md,'fieldname','hydrology.ischannels','numel',[1],'values',[0 1]);
@@ -113 +114 @@
                         fielddisplay(self,'bump_height','typical bump height (h_r) [m]');
                         fielddisplay(self,'omega','transition parameter (omega) []'); 
+                        fielddisplay(self,'rheology_B_base','Ice rheology factor B at base of ice (B) [Pa s^(-1/3)]');
                         disp(sprintf('      CHANNELS'));
                         fielddisplay(self,'ischannels','Do we allow for channels? 1: yes, 0: no');
@@ -143 +145 @@
                         WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','sheet_alpha','format','Double'); 
                         WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','sheet_beta','format','Double'); 
+                        WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','rheology_B_base','format','DoubleMat','mattype',1);
 
                         %Channels