Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m	(revision 13075)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m	(revision 13076)
@@ -75,11 +75,11 @@
 
 
 %Create zeros basalforcings and surfaceforcings
-if (isnan(md.surfaceforcings.precipitation)& (md.surfaceforcings.ispdd==1) & (md.surfaceforcings.ispdd==0)),
+if (isnan(md.surfaceforcings.precipitation) & (md.surfaceforcings.ispdd==1)),
 	md.surfaceforcings.precipitation=zeros(md.mesh.numberofvertices,1);
 	disp('      no surfaceforcings.precipitation specified: values set as zero');
 end
-if isnan(md.surfaceforcings.mass_balance) & (md.surfaceforcings.ispdd==0),
+if (isnan(md.surfaceforcings.mass_balance) & (md.surfaceforcings.ispdd==0)),
 	md.surfaceforcings.mass_balance=zeros(md.mesh.numberofvertices,1);
 	disp('      no surfaceforcings.mass_balance specified: values set as zero');
 end
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.m
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.m	(revision 13075)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.m	(revision 13076)
@@ -33,11 +33,11 @@
 end
 
 %Create zeros basal melting rate and surface mass balance if not specified
-if (isnan(md.surfaceforcings.precipitation) & (md.surfaceforcings.ispdd==1) & (md.surfaceforcings.ispdd==0)),
+if (isnan(md.surfaceforcings.precipitation) & (md.surfaceforcings.ispdd==1)),
 	md.surfaceforcings.precipitation=zeros(md.mesh.numberofvertices,1);
 	disp('      no surfaceforcings.precipitation specified: values set as zero');
 end
-if isnan(md.surfaceforcings.mass_balance) & (md.surfaceforcings.ispdd==0),
+if (isnan(md.surfaceforcings.mass_balance) & (md.surfaceforcings.ispdd==0)),
 	md.surfaceforcings.mass_balance=zeros(md.mesh.numberofvertices,1);
 	disp('      no surfaceforcings.mass_balance specified: values set as zero');
 end
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 13075)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 13076)
@@ -68,10 +68,10 @@
 	md.diagnostic.icefront=pressureload
 
 	#Create zeros basalforcings and surfaceforcings
-	if numpy.all(numpy.isnan(md.surfaceforcings.precipitation)):
+	if numpy.all(numpy.isnan(md.surfaceforcings.precipitation)) and (md.surfaceforcings.ispdd==1):
 		md.surfaceforcings.precipitation=numpy.zeros(md.mesh.numberofvertices)
 		print '      no surfaceforcings.precipitation specified: values set as zero'
-	if numpy.all(numpy.isnan(md.surfaceforcings.mass_balance)):
+	if numpy.all(numpy.isnan(md.surfaceforcings.mass_balance)) and (md.surfaceforcings.ispdd==0):
 		md.surfaceforcings.mass_balance=numpy.zeros(md.mesh.numberofvertices)
 		print '      no surfaceforcings.mass_balance specified: values set as zero'
 	if numpy.all(numpy.isnan(md.basalforcings.melting_rate)):
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.m
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.m	(revision 13075)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.m	(revision 13076)
@@ -64,11 +64,11 @@
 md.diagnostic.icefront=pressureload;
 
 %Create zeros basalforcings and surfaceforcings
-if (isnan(md.surfaceforcings.precipitation) & (md.surfaceforcings.ispdd==1) & (md.surfaceforcings.ispdd==0)),
+if (isnan(md.surfaceforcings.precipitation) & (md.surfaceforcings.ispdd==1)),
 	md.surfaceforcings.precipitation=zeros(md.mesh.numberofvertices,1);
 	disp('      no surfaceforcings.precipitation specified: values set as zero');
 end
-if isnan(md.surfaceforcings.mass_balance) & (md.surfaceforcings.ispdd==0),
+if (isnan(md.surfaceforcings.mass_balance) & (md.surfaceforcings.ispdd==0)),
 	        md.surfaceforcings.mass_balance=zeros(md.mesh.numberofvertices,1);
 		disp('      no surfaceforcings.mass_balance specified: values set as zero');
 end