Index: /issm/trunk-jpl/src/m/classes/SMBgemb.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBgemb.m	(revision 27195)
+++ /issm/trunk-jpl/src/m/classes/SMBgemb.m	(revision 27196)
@@ -73,5 +73,5 @@
 		% 2: default value of 1, in areas with grain radius below teThresh and areas of dry snow (not bare ice or wet) at the surface
 
-		swIdx  = NaN; % apply all SW to top grid cell (0) or allow SW to penetrate surface (1) (default 1, if aIdx=2 function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002))
+		swIdx  = NaN; %apply all SW to top grid cell (0) or allow SW to penetrate surface (1) (default 0, if swIdx=1 and aIdx=2, function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002))
 
 		denIdx = NaN; %densification model to use (default is 2):
@@ -236,5 +236,5 @@
 			end
 
-			fielddisplay(self,'swIdx','apply all SW to top grid cell (0) or allow SW to penetrate surface (1) [default 1, if aIdx=2 function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002)]');
+			fielddisplay(self,'swIdx','apply all SW to top grid cell (0) or allow SW to penetrate surface (1) [default 0, if swIdx=1 and aIdx=2 function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002)]');
 			fielddisplay(self,'denIdx',{'densification model to use (default is 2):',...
 				'1 = emperical model of Herron and Langway (1980)',...
@@ -343,5 +343,5 @@
 			self.aIdx = 1;
 			self.eIdx = 1;
-			self.swIdx = 1;
+			self.swIdx = 0;
 			self.denIdx = 2;
 			self.dsnowIdx = 1;
Index: /issm/trunk-jpl/src/m/classes/SMBgemb.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBgemb.py	(revision 27195)
+++ /issm/trunk-jpl/src/m/classes/SMBgemb.py	(revision 27196)
@@ -80,5 +80,5 @@
         # 2: default value of 1, in areas with grain radius below teThresh and areas of dry snow (not bare ice or wet) at the surface
 
-        self.swIdx                  = np.nan    # apply all SW to top grid cell (0) or allow SW to penetrate surface (1) (default 1, if aIdx=2 function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002))
+        self.swIdx                  = np.nan    # apply all SW to top grid cell (0) or allow SW to penetrate surface (1) (default 0, if swIdx=1 and aIdx=2, function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002)) 
 
         self.denIdx                 = np.nan    # densification model to use (default is 2):
@@ -238,5 +238,5 @@
             s += '{}\n'.format(fielddisplay(self, 'K', 'time scale temperature coef. (7) [d]'))
 
-        s += '{}\n'.format(fielddisplay(self, 'swIdx', 'apply all SW to top grid cell (0) or allow SW to penetrate surface (1) [default 1, if aIdx=2 function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002)]'))
+        s += '{}\n'.format(fielddisplay(self, 'swIdx', 'apply all SW to top grid cell (0) or allow SW to penetrate surface (1) [default 0, if swIdx=1 and aIdx=2 function of effective radius (Brun et al., 1992) or else dependent on snow density (taken from Bassford, 2002)]'))
         s += '{}\n'.format(fielddisplay(self, 'denIdx', ['densification model to use (default is 2):',
             '1 = emperical model of Herron and Langway (1980)',
@@ -330,5 +330,5 @@
         self.aIdx = 1
         self.eIdx = 1
-        self.swIdx = 1
+        self.swIdx = 0 
         self.denIdx = 2
         self.dsnowIdx = 1
Index: /issm/trunk-jpl/test/NightlyRun/test243.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test243.m	(revision 27195)
+++ /issm/trunk-jpl/test/NightlyRun/test243.m	(revision 27196)
@@ -10,4 +10,5 @@
 md.smb = SMBgemb(md.mesh,md.geometry);
 md.smb.dsnowIdx = 1;
+md.smb.swIdx = 1;
 
 %load hourly surface forcing date from 1979 to 2009:
Index: /issm/trunk-jpl/test/NightlyRun/test243.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test243.py	(revision 27195)
+++ /issm/trunk-jpl/test/NightlyRun/test243.py	(revision 27196)
@@ -21,4 +21,5 @@
 md.smb = SMBgemb(md.mesh, md.geometry)
 md.smb.dsnowIdx = 1
+md.smb.swIdx = 1
 
 #load hourly surface forcing date from 1979 to 2009:
Index: /issm/trunk-jpl/test/NightlyRun/test244.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test244.m	(revision 27195)
+++ /issm/trunk-jpl/test/NightlyRun/test244.m	(revision 27196)
@@ -11,4 +11,5 @@
 md.smb = SMBgemb(md.mesh,md.geometry);
 md.smb.dsnowIdx = 0;
+md.smb.swIdx = 1;
 
 %load hourly surface forcing date from 1979 to 2009:
Index: /issm/trunk-jpl/test/NightlyRun/test244.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test244.py	(revision 27195)
+++ /issm/trunk-jpl/test/NightlyRun/test244.py	(revision 27196)
@@ -30,4 +30,5 @@
 md.smb = SMBgemb(md.mesh, md.geometry)
 md.smb.dsnowIdx = 0
+md.smb.swIdx = 1
 
 #load hourly surface forcing date from 1979 to 2009:
Index: /issm/trunk-jpl/test/NightlyRun/test252.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test252.m	(revision 27195)
+++ /issm/trunk-jpl/test/NightlyRun/test252.m	(revision 27196)
@@ -10,4 +10,5 @@
 md.smb = SMBgemb(md.mesh,md.geometry);
 md.smb.dsnowIdx = 4;
+md.smb.swIdx = 1;
 
 %load hourly surface forcing date from 1979 to 2009:
Index: /issm/trunk-jpl/test/NightlyRun/test252.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test252.py	(revision 27195)
+++ /issm/trunk-jpl/test/NightlyRun/test252.py	(revision 27196)
@@ -23,4 +23,5 @@
 md.smb = SMBgemb(md.mesh, md.geometry)
 md.smb.dsnowIdx = 4
+md.smb.swIdx = 1
 
 #load hourly surface forcing date from 1979 to 2009: