Minimum depth for river

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
feroda

Minimum depth for river

#1 Unread post by feroda »

Hi all,

Hmin in ROMS is the minimum depth be considered(say 20m), while depth for water is not require by the river forcing file.
To add rivers into the model, what the minimum depth should be?

If it is also set to 20m, then most of the river area will be masked by land.

Thanks

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

Re: Minimum depth for river

#2 Unread post by kate »

Are you talking about UV_PSOURCE as your river input? The way that works is that the water comes in along the face of a grid box or several. The grid file determines the depth of that grid box in the model. One of the inputs to the model for the sources is Qshape, determining the vertical structure of that input in the sigma coordinates. For instance, if the river is 10 m deep and your hmin is 20 m, you could have the river come in over the upper half of the grid face. You would need to know the mapping from sigma space to physical space, but that is one of the outputs from ROMS for hmin, hmax, and an intermediate depth. Assuming that the river is at hmin and that your vertical coordinate is linear at hmin, you can apply the input to the upper half of the sigma depths, k=N/2 to k=N.

feroda

Re: Minimum depth for river

#3 Unread post by feroda »

Thank! I followed your advice and get ready for the river configuration.

I did the river case just now, unfortunately, got the errors saying:

GET_CYCLE - starting time for variable: river
is greater than current model time.
TMIN = 1.0000 TDAYS = 0.0000

Quite a few guys here also met this kind of errors after they update the ROMS source code.
I am wondering if there is something wrong in defining the "river_time".

Following the matlab tools from
https://www.myroms.org/index.php?page=Datasets

I defined "river_time" as
double river_time(river_time) ;
river_time:long_name = "river runoff time" ;
river_time:units = "day" ;
river_time:cycle_length = 360. ;
river_time:field = "river_time, scalar, series" ;


Any kind comment on that?

Thanks

Post Reply