Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp	(revision 12891)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp	(revision 12892)
@@ -713,5 +713,5 @@
 			input2->GetInputValue(&TemperaturesLgm[iv][month],gauss,month/12.*yts);
 			input3->GetInputValue(&PrecipitationsPresentday[iv][month],gauss,month/12.*yts);
-			monthlyprec[iv][month]=monthlyprec[iv][month]*yts; // convertion to m/yr
+			PrecipitationsPresentday[iv][month]=PrecipitationsPresentday[iv][month]/yts; // converion in m/sec
 		}
 	}
@@ -1568,5 +1568,4 @@
 	}
 	#endif
-
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/m/utils/BC/SetIceSheetBC.m
===================================================================
--- /issm/trunk-jpl/src/m/utils/BC/SetIceSheetBC.m	(revision 12891)
+++ /issm/trunk-jpl/src/m/utils/BC/SetIceSheetBC.m	(revision 12892)
@@ -34,9 +34,9 @@
 
 %Create zeros basal melting rate and surface mass balance if not specified
-if isnan(md.surfaceforcings.precipitation),
+if (isnan(md.surfaceforcings.precipitation) & (md.surfaceforcings.ispdd==1) & (md.surfaceforcings.ispdd==0)),
 	md.surfaceforcings.precipitation=zeros(md.mesh.numberofvertices,1);
 	disp('      no surfaceforcings.precipitation specified: values set as zero');
 end
-if isnan(md.surfaceforcings.mass_balance),
+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');
Index: /issm/trunk-jpl/src/m/utils/BC/SetIceShelfBC.m
===================================================================
--- /issm/trunk-jpl/src/m/utils/BC/SetIceShelfBC.m	(revision 12891)
+++ /issm/trunk-jpl/src/m/utils/BC/SetIceShelfBC.m	(revision 12892)
@@ -65,11 +65,11 @@
 
 %Create zeros basalforcings and surfaceforcings
-if isnan(md.surfaceforcings.precipitation),
+if (isnan(md.surfaceforcings.precipitation) & (md.surfaceforcings.ispdd==1) & (md.surfaceforcings.ispdd==0)),
 	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.mass_balance=zeros(md.mesh.numberofvertices,1);
-	disp('      no surfaceforcings.mass_balance specified: values set as zero');
+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
 if isnan(md.basalforcings.melting_rate),
Index: /issm/trunk-jpl/src/m/utils/BC/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk-jpl/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 12891)
+++ /issm/trunk-jpl/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 12892)
@@ -76,9 +76,9 @@
 
 %Create zeros basalforcings and surfaceforcings
-if isnan(md.surfaceforcings.precipitation),
+if (isnan(md.surfaceforcings.precipitation)& (md.surfaceforcings.ispdd==1) & (md.surfaceforcings.ispdd==0)),
 	md.surfaceforcings.precipitation=zeros(md.mesh.numberofvertices,1);
 	disp('      no surfaceforcings.precipitation specified: values set as zero');
 end
-if isnan(md.surfaceforcings.mass_balance),
+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');