Print density

residue.cpp

@@ -1072,7 +1072,7 @@ void printSpaceTimeOutput(double t, N_Vector* y, FILE* output, UserData data)
 	}
 
 	for (size_t i = 0; i < data->npts; i++) {
-		setGas(data,ydata,i+1,GAS) ; 
+		setGas(data,ydata,i+1) ; 
 		rho = data->gas->density() ;
 		fprintf(output, "%15.6e\t%15.6e\t",psi(i+1),t);
 		if(i==0){