Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 27451)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 27452)
@@ -7590,5 +7590,5 @@
 		b=i*nt;
 		c=av*nel;
-		for (int ae=0;ae<loads->nactiveloads;ae++){
+		for(ae=0;ae<loads->nactiveloads;ae++){
 			e=loads->combined_loads_index[ae];
 			a=AlphaIndex[c+e]*viscousnumsteps;
@@ -7695,5 +7695,5 @@
 		//get projection
 		if (spatial_component==1){ //north
-			for (int ae=0;ae<loads->nactiveloads;ae++){
+			for(ae=0;ae<loads->nactiveloads;ae++){
 				e=loads->combined_loads_index[ae];
 				horiz_projection[ae]=cos(2.0*M_PI*reCast<IssmDouble,int>(AzimIndex[av*nel+e])/65535.0); // 65535=2^16-1 = max value of 16 bits unsigned int
@@ -7701,5 +7701,5 @@
 			for(l=0;l<SLGEOM_NUMLOADS;l++){
 				nbar=slgeom->nbar[l];
-				for (ae=0;ae<loads->nactivesubloads[l];ae++){
+				for(ae=0;ae<loads->nactivesubloads[l];ae++){
 					e=loads->combined_subloads_index[l][ae];
 					horiz_projectionsub[l][ae]=cos(2.0*M_PI*reCast<IssmDouble,int>(AzimIndexsub[l][av*nbar+e])/65535.0);
@@ -7708,5 +7708,5 @@
 		}
 		else if (spatial_component==2){ //east
-			for (int ae=0;ae<loads->nactiveloads;ae++){
+			for(ae=0;ae<loads->nactiveloads;ae++){
 				e=loads->combined_loads_index[ae];
 				horiz_projection[ae]=sin(2.0*M_PI*reCast<IssmDouble,int>(AzimIndex[av*nel+e])/65535.0);
@@ -7714,5 +7714,5 @@
 			for(l=0;l<SLGEOM_NUMLOADS;l++){
 				nbar=slgeom->nbar[l];
-				for (ae=0;ae<loads->nactivesubloads[l];ae++){
+				for(ae=0;ae<loads->nactivesubloads[l];ae++){
 					e=loads->combined_subloads_index[l][ae];
 					horiz_projectionsub[l][ae]=sin(2.0*M_PI*reCast<IssmDouble,int>(AzimIndexsub[l][av*nbar+e])/65535.0);
@@ -7722,10 +7722,10 @@
 
 		//project load in the right direction 
-		for (int ae=0;ae<loads->nactiveloads;ae++){
+		for (ae=0;ae<loads->nactiveloads;ae++){
 			projected_loads[ae]=loads->combined_loads[ae]*horiz_projection[ae];
 		}
 		for(l=0;l<SLGEOM_NUMLOADS;l++){
 			nbar=slgeom->nbar[l];
-			for (ae=0;ae<loads->nactivesubloads[l];ae++){
+			for(ae=0;ae<loads->nactivesubloads[l];ae++){
 				projected_subloads[l][ae]=loads->combined_subloads[l][ae]*horiz_projectionsub[l][ae];
 			}
@@ -7734,5 +7734,5 @@
 		//do the convolution
 		c=av*nel;
-		for (int ae=0;ae<loads->nactiveloads;ae++){
+		for(ae=0;ae<loads->nactiveloads;ae++){
 			e=loads->combined_loads_index[ae];
 			a=AlphaIndex[c+e]*viscousnumsteps;
@@ -7744,5 +7744,5 @@
 			nbar=slgeom->nbar[l];
 			c=av*nbar;
-			for (ae=0;ae<loads->nactivesubloads[l];ae++){
+			for(ae=0;ae<loads->nactivesubloads[l];ae++){
 				e=loads->combined_subloads_index[l][ae];
 				a=AlphaIndexsub[l][c+e]*viscousnumsteps;
@@ -7757,5 +7757,4 @@
 
 	//free resources
-
 	xDelete<IssmDouble>(horiz_projection);
 	xDelete<IssmDouble>(projected_loads);