Intriduced qflux and qfluxbc
used coeq to replace 0.001 in heatflux calculation
cleaned set_bnd_vel()
deleted unused variables in set_bnd_vel
deleted unused variables
cleaned set_bnd_pressure
changed the document information
changed the funciton for mass conservation
changed the funciton for mass conservation
add comments
improved the format
add comments
added initial cosimulation code that could be called by modelica