you should be able to compare your imposed salt correction flux with what the model is getting by looking at the following diagnostics (assuming you’ll have a zero sfc_salt_flux_restore when you turn off SSS restoring):
id_pme_correct = register_diag_field('ocean_model','pme_correct', Grd%tracer_axes(1:2),&
Time%model_time, 'precip-evap from flux correction (>0 enters ocean)', &
'(kg/m^3)*(m/sec)', missing_value=missing_value,range=(/-1e6,1e6/))
id_total_ocean_pme_correct = register_diag_field('ocean_model','total_ocean_pme_correct',&
Time%model_time, 'total precip-evap from flux correction (>0 enters ocean)', &
'(kg/sec)/1e15', missing_value=missing_value,range=(/-1e6,1e6/))