Hello,
You need to understand how works the coefficient BETAC in subroutine drsurr.f, so that the density of water including your tracer has the correct value.
The reference tracer value is the value that corresponds to the reference density RHO0. I would say that in your case it would be zero, so that the reference density is the density of pure water (I mean without your tracer, you may have salt also). Then suppose that your tracer concentration is in kg/kg, when it is equal to 1 you must get the density of your tracer RHOtracer.
drsurr.f computes (RHO-RHO0)/RHO0, then you can find BETAC that will give :
RHO = RHOtracer when the concentration equals 1.
With best regards,
Jean-Michel Hervouet