hight temp values of my river input

[Entrapmen]

My domain is a rectangular channel and I want to add a river into my domain from one side.
I set the river discharge not as a single source but a series of adjacent sources (31 points - Q=1000 m^3/s for each).

domain.PNG (3.81 KiB) Viewed 4015 times

The run blows up immediately within 15 time steps. I get crazy temperature values on the boundary. Here is a vertical section taken along the mid-channel.

snapshot1.png (6.68 KiB) Viewed 4015 times

What might be the reason for this strange shape and values?

Ozgur

[kate]

You have to specify that you are going to provide the temperature values for the river water - and you have to provide those values. Saying you won't provide values causes the model to use an unconditionally unstable downwind advection right at the river entrance.

[Entrapmen]

I have already specified in my netcdf file river_temp, river_salt and river_flag and checked their array consistencies. They look fine.

[austinctodd]

I believe Kate is referring to whether or not you have "define TS_PSOURCE" in addition to "UV_PSOURCE" in your include file.

[Entrapmen]

I believe Kate is referring to whether or not you have "define TS_PSOURCE" in addition to "UV_PSOURCE" in your include file.

I activate both cpp options. I also add Q_PSOURCE.

[jcwarner]

sed_test1 was set up in a similar manner, but flow E-W. I have not tried that case is in a while, but it could provide guidance for you application.

[austinctodd]

So, do you want to be prescribing a vertical source and a horizontal source? It was my understanding that you should use either

#define UV_PSOURCE
#define TS_PSOURCE

or

#define Q_PSOURCE
#define TS_PSOURCE

This was addressed in earlier posts:
viewtopic.php?f=14&t=1663&p=5984&hilit=uv_psource#p5984

viewtopic.php?f=14&t=1663&p=5984&hilit=uv_psource#p5984

It appears that Q_PSOURCE is meant for a vertical point source, and UV_PSOURCE is meant for horizontal point sources. I would think that, for just a river through one wall, you would only want to use UV_PSOURCE.

[kate]

Also, what do you have for Lsrc and Tsrc?