Changeset 27778
- Timestamp:
- 06/01/23 11:42:37 (22 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/SMBsemic.m
r27773 r27778 105 105 % md.smb = initialize(md.smb,md); 106 106 107 if isnan(self.s0gcm) ,107 if isnan(self.s0gcm) 108 108 if ~isnan(md.geometry.surface) & (numel(md.geometry.surface) == md.mesh.numberofvertices) 109 109 self.s0gcm=md.geometry.surface; … … 114 114 end 115 115 end 116 if isnan(self.mask) ,116 if isnan(self.mask) 117 117 self.mask = 2*ones(md.mesh.numberofvertices,1); 118 118 disp(' no SMBsemic.mask specified: values set as 2 for ice'); … … 120 120 121 121 % update each values. 122 if isnan(self.Tamp) 122 if isnan(self.Tamp) 123 123 self.Tamp= 3*ones(md.mesh.numberofvertices,1); 124 124 disp(' no SMBsemic.Tamp specified: values set as 3.0'); … … 166 166 function md = checkconsistency(self,md,solution,analyses) % {{{ 167 167 168 if ismember('MasstransportAnalysis',analyses) ,168 if ismember('MasstransportAnalysis',analyses) 169 169 md = checkfield(md,'fieldname','smb.desfac','<=',1,'numel',1); 170 170 md = checkfield(md,'fieldname','smb.s0gcm','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]); … … 257 257 end 258 258 % albedo_scheme - 0: none, 1: slater, 2: isba, 3: denby, 4: alex. 259 if self.albedo_scheme == 0 :259 if self.albedo_scheme == 0 260 260 disp(sprintf('\n\tSEMIC snow albedo parameter of None.')); 261 261 disp(sprintf('\t albedo of snow is updated from albedo snow max (alb_smax).')); … … 270 270 fielddisplay(self,'tmin','minimum temperature for which albedo decline become effective. (default: 263.15 K)[unit: K])'); 271 271 fielddisplay(self,'tmax','maxmium temperature for which albedo decline become effective. This value should be fixed. (default: 273.15 K)[unit: K])'); 272 elseif self.albedo_scheme == 2 ,272 elseif self.albedo_scheme == 2 273 273 disp(sprintf('\n\tSEMIC snow albedo parameters for ISBA.? where is citation?')); 274 274 fielddisplay(self,'mcrit','critical melt rate (default: 6e-8) [unit: m/sec]'); … … 276 276 fielddisplay(self,'tau_a','dry albedo decline [unit: 1/day]'); 277 277 fielddisplay(self,'tau_f','wet albedo decline [unit: 1/day]'); 278 elseif self.albedo_scheme == 3 ,278 elseif self.albedo_scheme == 3 279 279 disp(sprintf('\n\tSEMIC snow albedo parameters for Denby et al. (2002 Tellus)')); 280 280 fielddisplay(self,'mcrit','critical melt rate (default: 6e-8) [unit: m/sec]'); 281 elseif self.albedo_scheme == 4 ,281 elseif self.albedo_scheme == 4 282 282 disp(sprintf('\n\tSEMIC snow albedo parameters for Alex.?')); 283 283 fielddisplay(self,'afac','[unit: ?]'); … … 356 356 outputs = self.requested_outputs; 357 357 pos = find(ismember(outputs,'default')); 358 if ~isempty(pos) ,358 if ~isempty(pos) 359 359 outputs(pos) = []; %remove 'default' from outputs 360 360 outputs = [outputs defaultoutputs(self,md)]; %add defaults
Note:
See TracChangeset
for help on using the changeset viewer.