Changeset 27778


Ignore:
Timestamp:
06/01/23 11:42:37 (22 months ago)
Author:
jdquinn
Message:

CHG: Clean up

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/classes/SMBsemic.m

    r27773 r27778  
    105105                        %  md.smb = initialize(md.smb,md);
    106106
    107                         if isnan(self.s0gcm),
     107                        if isnan(self.s0gcm)
    108108                                if ~isnan(md.geometry.surface) & (numel(md.geometry.surface) == md.mesh.numberofvertices)
    109109                                        self.s0gcm=md.geometry.surface;
     
    114114                                end
    115115                        end
    116                         if isnan(self.mask),
     116                        if isnan(self.mask)
    117117                                self.mask = 2*ones(md.mesh.numberofvertices,1);
    118118                                disp('      no SMBsemic.mask specified: values set as 2 for ice');
     
    120120
    121121                        % update each values.
    122                         if isnan(self.Tamp) 
     122                        if isnan(self.Tamp)
    123123                                self.Tamp= 3*ones(md.mesh.numberofvertices,1);
    124124                                disp('      no SMBsemic.Tamp specified: values set as 3.0');
     
    166166                function md = checkconsistency(self,md,solution,analyses) % {{{
    167167
    168                         if ismember('MasstransportAnalysis',analyses),
     168                        if ismember('MasstransportAnalysis',analyses)
    169169                                md = checkfield(md,'fieldname','smb.desfac','<=',1,'numel',1);
    170170                                md = checkfield(md,'fieldname','smb.s0gcm','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     
    257257                        end
    258258                        % albedo_scheme - 0: none, 1: slater, 2: isba, 3: denby, 4: alex.
    259          if self.albedo_scheme == 0:
     259         if self.albedo_scheme == 0
    260260            disp(sprintf('\n\tSEMIC snow albedo parameter of None.'));
    261261                                disp(sprintf('\t   albedo of snow is updated from albedo snow max (alb_smax).'));
     
    270270                                fielddisplay(self,'tmin','minimum temperature for which albedo decline become effective. (default: 263.15 K)[unit: K])');
    271271                                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
    273273                                disp(sprintf('\n\tSEMIC snow albedo parameters for ISBA.? where is citation?'));
    274274                                fielddisplay(self,'mcrit','critical melt rate (default: 6e-8) [unit: m/sec]');
     
    276276                                fielddisplay(self,'tau_a','dry albedo decline [unit: 1/day]');
    277277                                fielddisplay(self,'tau_f','wet albedo decline [unit: 1/day]');
    278                         elseif self.albedo_scheme == 3,
     278                        elseif self.albedo_scheme == 3
    279279                                disp(sprintf('\n\tSEMIC snow albedo parameters for Denby et al. (2002 Tellus)'));
    280280                                fielddisplay(self,'mcrit','critical melt rate (default: 6e-8) [unit: m/sec]');
    281                         elseif self.albedo_scheme == 4,
     281                        elseif self.albedo_scheme == 4
    282282                                disp(sprintf('\n\tSEMIC snow albedo parameters for Alex.?'));
    283283                                fielddisplay(self,'afac','[unit: ?]');
     
    356356                        outputs = self.requested_outputs;
    357357                        pos  = find(ismember(outputs,'default'));
    358                         if ~isempty(pos),
     358                        if ~isempty(pos)
    359359                                outputs(pos) = []; %remove 'default' from outputs
    360360                                outputs      = [outputs defaultoutputs(self,md)]; %add defaults
Note: See TracChangeset for help on using the changeset viewer.