Refinement from FIner SPs directory is to be same as the MODFLOW 6 version of the problem. Changes include: 

CHD: bottom boundary head is changed from 59.70427 to 59.70180100 to give that slight adjustment to the velocity. 

Use 601 SPs. with temperature (concentration) variation from MF6 model

Solute directory solves for equivalent concentration with solute equation

Heat directory solves directly for temperture with heat equation

Current solves for both heat and solute: NOTE: replace Stallman.pcb (generated from Vistas) with Stallman_2SPECI.PCB which has PCB for both conc and temp. Vistas only generates for 1 species. 