00001 if (oCorr == 0) 00002 { 00003 #include <finiteVolume/rhoEqn.H> 00004 } 00005 00006 #include "UEqn.H" 00007 00008 #include "hEqn.H" 00009 00010 // --- PISO loop 00011 for (int corr=0; corr<nCorr; corr++) 00012 { 00013 #include "pEqn.H" 00014 } 00015 00016 turb.correct(); 00017 00018 rho = thermo.rho();