Changeset 27967 for issm/trunk-jpl/src/m/classes/SMBgradients.m

Timestamp:

10/25/23 12:50:47 (17 months ago)

Author:

badgeley

Message:

CHG: The SMBgradients class now expects smb units of m/yr ice eq.

File:

• issm/trunk-jpl/src/m/classes/SMBgradients.m (modified) (2 diffs)

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

@@ -47 +47 @@
                                 md = checkfield(md,'fieldname','smb.href','timeseries',1,'NaN',1,'Inf',1);
                                 md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1,'Inf',1);
+                                if max(max(abs(md.smb.smbref(1:end-1,:))))<1
+              disp('!!! Warning: SMBgradients now expects smbref to be in m/yr ice eq. instead of mm/yr water eq.');
+            end
                                 md = checkfield(md,'fieldname','smb.b_pos','timeseries',1,'NaN',1,'Inf',1);
                                 md = checkfield(md,'fieldname','smb.b_neg','timeseries',1,'NaN',1,'Inf',1);
@@ -59 +62 @@
                         disp(sprintf('\n   SMB gradients parameters:'));
                         fielddisplay(self,'href',' reference elevation from which deviation is used to calculate SMB adjustment in smb gradients method [m]');
-                        fielddisplay(self,'smbref',' reference smb from which deviation is calculated in smb gradients method [mm/yr water equiv]');
+                        fielddisplay(self,'smbref',' reference smb from which deviation is calculated in smb gradients method [m/yr ice equiv]');
                         fielddisplay(self,'b_pos',' slope of hs - smb regression line for accumulation regime required if smb gradients is activated');
                         fielddisplay(self,'b_neg',' slope of hs - smb regression line for ablation regime required if smb gradients is activated');