TELEMAC-MASCARET Forum: V8P1r0: Telemac2D+Waqtel EUTRO (can't add salt+temp tracers) (1/1)

Welcome, Guest

TOPIC: V8P1r0: Telemac2D+Waqtel EUTRO (can't add salt+temp tracers)

V8P1r0: Telemac2D+Waqtel EUTRO (can't add salt+temp tracers) 4 years 7 months ago #35795

shivanesh.rao
OFFLINE
Junior Boarder
Posts: 37
Thank you received: 1

Hi all,

I'm having a problem running a parallel model of telemac2d+waqtel eutro model. This is using the latest version.

When I run the model with no tracers (i.e., no salinity or temp...just flow discharges+WAQTEL auto-declared tracers), the full model runs in both serial (--ncsize=1) and parallel (--ncsize=8) modes.

When I run the model with 1 tracer (i.e., EITHER salinity or temp...plus flow discharges+WAQTEL auto-declared tracers), the full model runs in both serial and parallel modes.

However, when I run with 2 tracers (i.e., both salinity and temp...plus flow discharges+WAQTEL auto-declared tracers). The model runs only in serial (1 cpu, --ncsize=1) but fails in parallel (--ncsize=8) with the following error:

USING STREAMLINE VERSION 7.3 FOR CHARACTERISTICS
PARALLEL:: ORG_CHARAC_TYPE1:: NOMB NOT IN RANGE [0..MAX_BASKET_SIZE]
MAX_BASKET_SIZE, NOMB:               10           11

MPI_ABORT was invoked...

If I understand the waqtel (EUTRO), it adds 8 tracers to the existing tracers declared by the user. So when I declare 2 tracers + 8 tracers from EUTRO, it should be 10, thus <=max basket of 10...so should work.

Has anyone come across this? It seems like an additional tracer is added somewhere in the parallel implementation since it works in the --ncsize=1 mode.

Any help would be greatly appreciated.

The administrator has disabled public write access.

V8P1r0: Telemac2D+Waqtel EUTRO (can't add salt+temp tracers) 4 years 7 months ago #35796

pham
OFFLINE
Administrator
Posts: 1559
Thank you received: 602

Hello,

The problem is because you use the characteristics method for advection schemes for both tracers (10 of them) + the 2 components of velocity, thus 12 values.
In parallel, it uses the STREAMLINE module that may call ORG_CHARAC_TYPE1 subroutine where the MAX_BASKET_SIZE variable is hard-coded to 10.
And you try to advect 12 variables.

Anyway, when modelling tracers, it is strongly recommended to use conservative advection schemes, in particular not characteristics.

The WAQTEL examples in 2D and the default value for SCHEME FOR ADVECTION OF TRACERS may have to be changed.

I would advise you to add or change the following options in your case (EUTRO+temperature+salinity):
SCHEME FOR ADVECTION OF TRACERS = 5;5;5;5;5;5;5;5;5;5
SCHEME OPTION FOR ADVECTION OF TRACERS = 4;4;4;4;4;4;4;4;4;4
(with as many values as the number of tracers).
You will then use the LIPS advection scheme which is the default and recommended one in 3D + works on tidal flats.

Hope this helps,

Chi-Tuan

The administrator has disabled public write access.

V8P1r0: Telemac2D+Waqtel EUTRO (can't add salt+temp tracers) 4 years 7 months ago #35797

shivanesh.rao OFFLINE Junior Boarder Posts: 37 Thank you received: 1	Hi Chi-Tuan, Thank you for the speedy reply and explanation. I tried your changes and the model is working now. Thank you for the help. Cheers, Shiv
	The administrator has disabled public write access.