float* phiInit = readArray_real("phiInit");
float* rhoInit = readArray_real("rhoInit");
float* piInit = readArray_real("piInit");
float* phiInit = readArray_real("phiInit");
float* rhoInit = readArray_real("rhoInit");
float* piInit = readArray_real("piInit");