With regards to mosfet selection and RDS(on) advice ,
It's true to say having a low RDS(on) does not have any meaning in linear mode,
Though It does play a small role here (e.Load) in that it will determine the minimum resistance that the load can achieve ( minimum Burden Resistance), this then determines the maximum current the load could sink with low input Voltages present. (simple e.g - if we have a Rsense of .1 ohm and a RDS(on) of .1 ohm then the Max current we could possibly sink when input is only 500mV would be .5/.2= 2.5 A ,if you have two of these in parrallel then it would be 5A ).
But also when choosing a Mosfet for linear application I would use RDS(on) a filter search critera and choose the power mosfet with the HIGHEST rds(on) I could accomodate to meet my minimum resistance critera.
The logic here is that the newer generations of mosfet which have been optimised for switching (so optimised for a really Low RDS(on) and high GM which is good for switch efficiency but not meaningfull for linear operation ,and one of the trade offs (no free lunch as usual) for producing a low RDS(on) high transconductance Fet seems to be that the Vgs Zero temp Co 'ZTC point' (which is the point which if we are below we have to have some extra thermal derating at higher Vds) gets even higher so requiring us to have extra derating .
I am not 100% sure why the ZTC point moves up so much with lower rds(on) power mosfets but it will make a large difference to power handling at higher end of Vds .
By selecting an older generation or even newer but unoptimised RDS(on) FET that should be giving you a mosfet with higher ZTC point and so will have much less of this extra thermal derating..
So filtering by this critea and power dissipation should be a good selection method for this application, no need for expensive ones as the real bottle neck here in this application is Case to Sink Thermal resitance , so if theres no need to skimp on numbers you can use as many cases as required in parrallel to get overall low Case to Sink thermal resistance,
also as Mr J.Diddy pointed out opamps are also really cheap so giving each it's own opamp is no big cost.(drive speed is also another reason why it's also beneficial to have each mosfet with it's own opamp for the extra drive current capability you are getting over having one opamp trying to drive multiple large Fets).
(I did A quick fet search which turns up examples like the IRFP250 which is only $1.50 a piece and will be much more suitable than your irfp3306 (notice the much higher RDS(on) 85mO compared to the irfp3306 with 6mO and notice there is much less extra derating required on the soa at the right bottom corner ,this one is used in Maynuo load, someone said ) .
Also to add another reason for choosing a higher RDS(on) power fet for linear mode operation is that there's also a realtionship between Crss (gate to drain capacitance- Cgd) and RDS(on) . Cgd generally gets larger the lower the RDS(On) is made, and having as little Cgd as possible is important for a stable and fast controll loop, in this application it's value is multiplied (miller effect and low Vds effect) and can appear as much as 1000 times larger as the figure given in the mosfet data sheets. This capacitance has to be dis/charged by the driving opamp so it's best if this is kept small.
Regards