How to setup the river discharges with only salinity but no temperature?

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
yz966445
Posts: 10
Joined: Fri Mar 10, 2017 11:22 pm
Location: Dalhousie University

How to setup the river discharges with only salinity but no temperature?

#1 Unread post by yz966445 »

Hi everyone,
Could anyone tell me how to setup the river discharge in ROMS model as I need?
I want the added river discharge affect the water salinity but no affecting the ocean temperature.

I am using "define ANA_PSOURCE" and "ana_psource.h" file to add the river discharge to my model.
For my own understanding, I setup the LtracerSrc in .in file as ''LtracerSrc == F T F ! temperature, salinity, inert".
However, I found the added river discharge still has lower temperature than the surrounding water.
Could anyone help me with that problem?


Thanks,
Yi Sui

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: How to setup the river discharges with only salinity but no temperature?

#2 Unread post by kate »

It doesn't make sense numerically to do what you want. As you see, the thing will become unstable. Lacking river temperatures, some have used air temperature instead.

yz966445
Posts: 10
Joined: Fri Mar 10, 2017 11:22 pm
Location: Dalhousie University

Re: How to setup the river discharges with only salinity but no temperature?

#3 Unread post by yz966445 »

Hi Kate,
Thanks for you reply.
Could we do that setup? Let's setup the river water has the salinity = 0.4 psu and has the temperature as the surrounding water.
I think the default setup in ana_psource.h is
SOURCES(ng)%Tsrc(is,k,itemp)=5
SOURCES(ng)%Tsrc(is,k,isalt)=0.4

Do we must setup Tsrc(is,k,itemp) a certain value?

Yi Sui

kate wrote: Wed Feb 10, 2021 6:35 pm It doesn't make sense numerically to do what you want. As you see, the thing will become unstable. Lacking river temperatures, some have used air temperature instead.

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: How to setup the river discharges with only salinity but no temperature?

#4 Unread post by kate »

For numerical stability you must give it something other than the water temperature at that grid point. Using the water temperature there is an unconditionally unstable numerical scheme - downwind advection. In MOM6, the default is to use air temperature instead, something you'd have anyway if you were using bulk fluxes. In ROMS, you'd have to code that in since it isn't the default.

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: How to setup the river discharges with only salinity but no temperature?

#5 Unread post by wilkin »

When I recoded the LwSrc option to introduce point sources as a divergence at the cell centers (instead of the the cell faces with LuvSrc) I did a sequence of tests with a modified RIVERPLUME1 test case and LtracerSrc == F T worked perfectly.

When LTracerSrc is false, tracer mass is added to compensate for the added volume to keep the receiving cell concentration unchanged. In the case of temperature, this is effectively adding heat. The relevant code is in step3d_t.F.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

yz966445
Posts: 10
Joined: Fri Mar 10, 2017 11:22 pm
Location: Dalhousie University

Re: How to setup the river discharges with only salinity but no temperature?

#6 Unread post by yz966445 »

Hi Wilkin,

Could you please tell me more information about how you setup? I want to make the LtracerSrc == F T worked perfectly. I am not sure how to code as you say in proper way. DId you recoded the LwSrc section? I very appreciate your help.

Yi Sui
wilkin wrote: Thu Feb 11, 2021 3:37 pm When I recoded the LwSrc option to introduce point sources as a divergence at the cell centers (instead of the the cell faces with LuvSrc) I did a sequence of tests with a modified RIVERPLUME1 test case and LtracerSrc == F T worked perfectly.

When LTracerSrc is false, tracer mass is added to compensate for the added volume to keep the receiving cell concentration unchanged. In the case of temperature, this is effectively adding heat. The relevant code is in step3d_t.F.

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: How to setup the river discharges with only salinity but no temperature?

#7 Unread post by wilkin »

Could you please tell me more information about how you setup? I want to make the LtracerSrc == F T worked perfectly. I am not sure how to code as you say in proper way. DId you recoded the LwSrc section? I very appreciate your help.
The code update is documented in this ticket from July 2020:
https://www.myroms.org/projects/src/ticket/860
If you are using an earlier version of the code, LwSrc will give you wrong results.

The distinction between LuvSrc and LwSrc is documented in WikiROMS:
https://www.myroms.org/wiki/River_Runoff

All you have to do is make sure you have the correct river_Xposition and river_Eposition (i,j position of the rho-cell into which you want the river to flow, in ROMS FORTRAN index convention counting from 0,0 not 1,1)
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

Post Reply