diff --git a/fslselectvols.cc b/fslselectvols.cc
index 6b1d723fe9fa375fb1e0ebbf918ef84f3b5bc185..6befbf43bf68401ca2dd59b53940b2edbc92b438 100644
--- a/fslselectvols.cc
+++ b/fslselectvols.cc
@@ -146,7 +146,7 @@ int main(int argc,char *argv[]){
 	    v=0;v2=0;
 	    for(int i=1;i<=list.Ncols();i++){	  
 	      v += data[list(1,i)](x,y,z);
-	      v2 += v*v;
+	      v2 += data[list(1,i)](x,y,z)*data[list(1,i)](x,y,z);
 	      if(!meancalc && !varcalc)
 		subdata[i-1](x,y,z)= data[list(1,i)](x,y,z);
 	    }