Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 22484)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 22485)
@@ -54,8 +54,8 @@
 				if (issurface[i]==1){
 					spcvector[i*N+j] = heatcapacity*(spcvector[i*N+j]-referencetemperature);
-					spcvectorstatic[i*N+j] = nan("");
+					spcvectorstatic[i*N+j] = NAN;
 				}
 				else{
-					spcvector[i*N+j] = nan("");
+					spcvector[i*N+j] = NAN;
 					spcvectorstatic[i*N+j] = heatcapacity*(spcvectorstatic[i*N+j]-referencetemperature);
 				}
Index: /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 22484)
+++ /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 22485)
@@ -50,6 +50,6 @@
 		for(int j=0;j<N;j++){
 			if (isdynamic){
-				if (issurface[i]==1)spcvectorstatic[i*N+j] = nan("");
-				else spcvector[i*N+j] = nan("");
+				if (issurface[i]==1)spcvectorstatic[i*N+j] = NAN;
+				else spcvector[i*N+j] = NAN;
 			}
 		}