Changeset 22482 for issm/trunk-jpl/src/m/classes/SMBgemb.m

Timestamp:

02/27/18 19:33:05 (7 years ago)

Author:

schlegel

Message:

CHG: add option for forcing albedo only above a specific surface density

File:

• issm/trunk-jpl/src/m/classes/SMBgemb.m (modified) (9 diffs)

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

@@ -60 +60 @@
                                           % 3: density and cloud amount [Greuell & Konzelmann, 1994]
                                           % 4: exponential time decay & wetness [Bougamont & Bamber, 2005]
-                                          % 5: ingest MODIS mode, direct input from aValue parameter applied to surface ice only
 
                 swIdx  = NaN; % apply all SW to top grid cell (0) or allow SW to penetrate surface (1) (default 1)
@@ -89 +88 @@
                 t0dry = NaN; % warm snow timescale (30) 
                 K     = NaN; % time scale temperature coef. (7) 
+                adThresh = NaN; %Apply aIdx method to all areas with densities below this value,
+                                %or else apply direct input value from aValue, allowing albedo to be altered.
+                                                         %Default value is rho water (1023 kg m-3).
 
                 %densities:
@@ -125 +127 @@
                         self.pAir=project3d(md,'vector',self.pAir,'type','node');
 
-                        if (aIdx == 0 | aIdx == 5) & ~isnan(self.aValue)
+                        if (aIdx == 0) & ~isnan(self.aValue)
                                 self.aValue=project3d(md,'vector',self.aValue,'type','node');
                         end
@@ -169 +171 @@
                 self.t0dry = 30;
                 self.K = 7;
+                self.adThresh = 1023;
 
                 self.teValue = ones(mesh.numberofelements,1);
@@ -212 +215 @@
                         md = checkfield(md,'fieldname','smb.teValue','timeseries',1,'NaN',1,'Inf',1,'>=',0,'<=',1);
 
-                        md = checkfield(md,'fieldname','smb.aIdx','NaN',1,'Inf',1,'values',[0,1,2,3,4,5]);
+                        md = checkfield(md,'fieldname','smb.aIdx','NaN',1,'Inf',1,'values',[0,1,2,3,4]);
                         md = checkfield(md,'fieldname','smb.swIdx','NaN',1,'Inf',1,'values',[0,1]);
                         md = checkfield(md,'fieldname','smb.denIdx','NaN',1,'Inf',1,'values',[1,2,3,4,5]);
@@ -223 +226 @@
                         md = checkfield(md,'fieldname','smb.InitDensityScaling','NaN',1,'Inf',1,'>=',0,'<=',1);
                         md = checkfield(md,'fieldname','smb.ThermoDeltaTScaling','NaN',1,'Inf',1,'>=',0,'<=',1);
+                        md = checkfield(md,'fieldname','smb.adThresh','NaN',1,'Inf',1,'>=',0);
 
                         switch self.aIdx,
@@ -278 +282 @@
                         fielddisplay(self,'ThermoDeltaTScaling',{'scaling factor to multiply the thermal diffusion timestep (delta t)'});
                         fielddisplay(self,'outputFreq','output frequency in days (default is monthly, 30)');
+                        fielddisplay(self,'adThresh',{'Apply aIdx method to all areas with densities below this value,','or else apply direct input value from aValue, allowing albedo to be altered.'});
                         fielddisplay(self,'aIdx',{'method for calculating albedo and subsurface absorption (default is 1)',...
                                 '0: direct input from aValue parameter',...
@@ -283 +288 @@
                                 '2: effective grain radius [Brun et al., 2009]',...
                                 '3: density and cloud amount [Greuell & Konzelmann, 1994]',...
-                                '4: exponential time decay & wetness [Bougamont & Bamber, 2005]',...
-                                '5: ingest MODIS mode, direct input from aValue parameter applied to surface ice only'});
+                                '4: exponential time decay & wetness [Bougamont & Bamber, 2005]'})
 
                         fielddisplay(self,'teValue','Outward longwave radiation thermal emissivity forcing at every element (default in code is 1)');
@@ -373 +377 @@
                         WriteData(fid,prefix,'object',self,'class','smb','fieldname','t0dry','format','Double');
                         WriteData(fid,prefix,'object',self,'class','smb','fieldname','K','format','Double');
+                        WriteData(fid,prefix,'object',self,'class','smb','fieldname','adThresh','format','Double');
 
                         WriteData(fid,prefix,'object',self,'class','smb','fieldname','aValue','format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts);