From cab20ac66fc5637626b6cf6dfbe7fece9c406a04 Mon Sep 17 00:00:00 2001
From: Weiye Wang <weiyewan@usc.edu>
Date: Wed, 19 Oct 2022 15:59:25 -0700
Subject: [PATCH] Print density

---
 residue.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/residue.cpp b/residue.cpp
index de48fde..c211122 100644
--- a/residue.cpp
+++ b/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){