Fixing Vortex in cell

example/Numerics/Vortex_in_cell/Makefile

@@ -4,18 +4,22 @@ CC=mpic++
 
 LDIR =
 
+OPT=
 OBJ_VIC_PETSC = main_vic_petsc.o
 
+vic_petsc_test: OPT += -DTEST_RUN
+vic_petsc_test: vic_petsc
+
 %.o: %.cpp
-	$(CC) -O3 -g -c --std=c++11 -o $@ $< $(INCLUDE_PATH)
+	$(CC) -O3 $(OPT) -g -c --std=c++11 -o $@ $< $(INCLUDE_PATH)
 
 all: vic_petsc
 
 vic_petsc: $(OBJ_VIC_PETSC)
-	$(CC) -o $@ $^ $(CFLAGS) $(LIBS_PATH) $(LIBS_SE2)
+	$(CC) -o $@ $^ $(LIBS_PATH) $(LIBS_SE2)
 	
-run: all
-	mpirun -np 3 ./vic_petsc
+run: vic_petsc_test
+	mpirun -np 4 ./vic_petsc
 	
 .PHONY: clean all
 

example/Numerics/Vortex_in_cell/main_vic_petsc.cpp

@@ -142,6 +142,12 @@ float nu = 1.0/tgtre;
 // float dt = 0.0025 for Re 7500
 float dt = 0.0125;
 
+#ifdef TEST_RUN
+const unsigned int nsteps = 10;
+#else
+const unsigned int nsteps = 10001;
+#endif
+
 // All the properties by index
 constexpr unsigned int vorticity = 0;
 constexpr unsigned int velocity = 0;
@@ -1280,7 +1286,7 @@ int main(int argc, char* argv[])
 	}
 
 	// Time Integration
-	for ( ; i < 10001 ; i++)
+	for ( ; i < nsteps ; i++)
 	{
 		// do step 4-5-6-7
 		do_step(particles,g_vort,g_vel,g_dvort,domain,inte,phi_s,solver);