Commits

Peter Lichtner committed 5a25619

Fixed bug in mphase boundary conditions.

Comments (0)

Files changed (1)

src/pflotran/mphase.F90

 ! Reaction terms here
 ! Note if iireac > 0, then it is the node global index
 
+!#if 0
   if (option%ntrandof > 0) then
     if (iireac > 0) then
       !H2O
 !       option%flow_dt*1D-3
     endif
   endif
-  
+!#endif
 ! if (option%use_isothermal)then
 !   Res(1:option%nflowdof) = mol(:)
 ! else
   fluxm = 0.d0
   fluxe = 0.d0
   v_darcy = 0.d0
+  vv_darcy = 0.d0
   density_ave = 0.d0
   q = 0.d0
   ukvr = 0.d0