Hi again! I am modelling close to Rio de la Plata so I am forcing the model with daily river discharges. I have defined by defect SFLX_CORR, that way the fresh water plume is almost non existent and I have no vertical salinity stratification close to Rio de la Plata mouth. I tried with undefined SFLX_CORR, in that case I get the fresh water plume in the oceanic region, but my salinity in regions where RIo the la Plata has no effect get values unrealisticly high (e.g. 36-37 close to the coast)

SO, I am trying now by changing the correction parameter CST in the get_sss.F file SSSFLX_model ~ SSS*(E-P) + CST * (SSS_model - SSS)

As that parameter is supposed to be the same as dqdSST applied in Qcorrection, I changed where it sets scale factor to convert input dQdSST from Watts/m2/Celsius to meter/second:

if defined SFLX_CORR && !defined QCORRECTION

    scldqdt=1./(rho0*Cp)

endif

I changed it with

if defined SFLX_CORR && !defined QCORRECTION

    scldqdt=1./(2.*rho0*Cp)

endif

After 6 month of simulation I can see only about a difference of 0.6 in salinity values. SO I am not sure I am making the change in the right place.

Is there some other way to change dqdt values only for the salinty field? Or should I diminish even more the scldqdt parameter?

Many thanks for any help!

Hi again! again!

I am modelling close to Rio de la Plata so I am forcing the model with daily river discharges. I have had defined by defect SFLX_CORR, and that way the fresh water plume is was almost non existent and I have had no vertical salinity stratification close to Rio de la Plata mouth. mouth.

I tried with undefined SFLX_CORR, in that case I get the fresh water plume in the oceanic region, but my salinity in regions where RIo Rio the la Plata has no effect get values unrealisticly high (e.g. 36-37 close to the coast)

SO, I am trying now So, I tried by changing the correction parameter CST in the get_sss.F file SSSFLX_model ~ SSS*(E-P) + CST * (SSS_model - SSS)

As that parameter value is supposed to be the same as dqdSST applied in Qcorrection, I changed where it sets the scale factor to convert input dQdSST from Watts/m2/Celsius to meter/second:

if defined SFLX_CORR && !defined QCORRECTION

    scldqdt=1./(rho0*Cp)

endif

I changed it with with:

if defined SFLX_CORR && !defined QCORRECTION

    scldqdt=1./(2.*rho0*Cp)

endif

After 6 month of simulation with the new parameter, I can see only about a difference of 0.6 about 0.4 in salinity values. SO values in the area I am suppose to get the water plume. So I am not sure I am making the change in the right place.

Is there some other way to change dqdt values only for the salinty field? Or should I diminish even more the scldqdt parameter? Or, does anyone have another idea about how to solve my issue?

Many thanks for any help!