Changeset 27967 for issm/trunk-jpl/src

Timestamp:

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

Author:

badgeley

Message:

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

Location:

issm/trunk-jpl/src

Files:

• c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp (modified) (1 diff)
• m/classes/SMBgradients.m (modified) (2 diffs)

issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp

@@ -67 +67 @@
                         }
 
-                        smb[v]=smb[v]/1000*rho_water/rho_ice;      // SMB in m/y ice
                 }  //end of the loop over the vertices
 

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');