Dear All,
I’m analyzing seawater intrusion problem. In FEFLOW modelling the hydraulic head is represented with reference density of freshwater, which means hydraulic head at the seaside (hs) couldn’t be used directly as a boundary condition and should be converted to hydraulic head with reference density of freshwater (h) using h=( (1+ (?s-?0  )/?0 ) hs - ((?s-?0  )/?0 )z). This is easily to be done for the sea side, while for adjusting the initial heads for all the study area as I have different salinity, which means different actual density at every cell of the model, couldn’t be done because I don’t know the specific density corresponds to each salinity.
Please any help is much appreciated how to assign the initial heads considering variable salinity from fresh water to saline water across the study area.

Thanks