open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Neil,

Can you send your steering file please? I should have the other input files to run your model now otherwise.

Chi-Tuan

Thanks again Chi-Tuan,

I have attached my WATERSOURCES file just including the first few time steps; it is quite large so I was concerned this may be an issue so I have cut it right down for testing.

I have also attached my cas file - the only part I have modified is the sources section.

Kind regards,

Neil

Hello Neil,

The problem is because the number of tracers times the number of sources lead to exceed the MAXKEYWORDS value which is the maximum number of keywords which is limited at 300 in the official version.

One way to solve is to change in the declarations_special module (in sources/utils/special folder) but directly where your release of TELEMAC is installed not in a FORTRAN file:
INTEGER, PARAMETER :: MAXKEYWORD = 300
to
INTEGER, PARAMETER :: MAXKEYWORD = 400 or more if needed
and recompile the whole system.

That needs that you can do the last actions.

Hope this helps,

Chi-Tuan

Hi Chi-Tuan,

Sorry for the slow response. I finally have Telemac v8 installed so I am just doing some testing with this.

Regarding the MAXKEYWORD issue - that sounds familiar now, I remember reading about this somewhere. I did ask my IT services to make the modification that you suggested on V8 and then to recompile. Unfortunately he said it failed to recompile; not sure why this might be. I have asked him to have another go - I'll let you know how it goes!

Many thanks,

Neil

Oh dear - it seems I'm back to square one again by moving to v8!

I have done some testing with V8p1r1 and whilst I can get my model to work once again with 5 tracers. If I try to increase the number of tracers back to 10 it won't run. I have kept everything else the same (I note that there are no updates to streamline.f, org_charac_type1.F or declerations_parallel.F).

I have attached my error and output files which gives the error. I have also attached my steering file but this is as it was before.

Sorry to have to go backwards on this but it would be nice to move my model to the most recent version of Telemac which was one of the first pieces of advice given!

Once again many thanks for your help and patience!

Hi,

Just wondering if anyone else has had success with running 10 or more tracers on V8 (In parallel).

I still haven't found a way to run 10 tracers on v8 so looks like I am going to have to go back to V7 for now which works for my model.

Thanks,

Neil

Hello Neil,

Can you summarize your installation + setup?
- installation of release v8p1r1 + change of MAXKEYWORDS from 300 to 400 e.g. or not? This suggesion, only if using too many tracers and too many source points. Did you also change the few lines with the value 14 rather than 10 (in your installation and/or FORTRAN FILE), as suggested in a previous post to be able to run Thompson Boundary Conditions, for the following 3 subroutines:
- declarations_parallel.F with DOUBLE PRECISION :: BASKET(14) ! VARIABLES INTERPOLATED AT THE FOOT,
- org_charac_type1.F with INTEGER, PARAMETER :: MAX_BASKET_SIZE=14,
- streamline.f with INTEGER, PARAMETER :: MAX_BASKET_SIZE=14 ! LARGE,
- no sources points in your steering file, but 10 tracers, LIPS scheme for every advection of tracers (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) + Thompson conditions for one liquid boundary, the last one = 6th (OPTION FOR LIQUID BOUNDARIES = 1;1;1;1;1;2)?

For SCHEME FOR ADVECTION OF TRACERS, you can also use schemes 13 or 14 with SCHEME FOR ADVECTION OF TRACERS = 13;13;13;13;13;13;13;13;13;13 or 14;14;14;14;14;14;14;14;14;14 + SCHEME OPTION FOR ADVECTION OF TRACERS = 1;1;1;1;1;1;1;1;1;1.
They are more diffusive, but may be less sensitive to time step.


Hope this helps,

Chi-Tuan

PS: post www.opentelemac.org/index.php/kunena/16-...dd-salt-temp-tracers
shows running 10 tracers is possible with release v8p1r1.
But if you still want to use Thompson BC, you have to do the changes suggested in the previous post (value 10 -> 14).

Hi Pham,

Firstly, many thanks for getting back to me again.

I think most of the details you are asking about are in the message above. To summarise:

Version: V8p1r1

MAXKEYWORDS: 300 (As you say I found in V7 it wasn’t necessary to change this to run with 10 tracers and 0 watersources, once I got this working I had to get a colleague who works on the cluster to make the changes on my behalf but I asked them to change on V7 as V8 wasn’t working – I can ask them to change and recompile on V8 but I don’t think this is the issue here).

Changes to subroutines: Yes all 3 as you suggested.

Keywords for tracers and liquid boundaries; Yes, as you state (6th liquid boundary Is an ocean boundary hence the use of Thompson-conditions)
I have tried the changes that you propose for the scheme for advection of tracers but they haven’t.

Although it isn’t always detailed in the output or error files usually I come back the error related to the ‘MAX_BASKET_SIZE’:

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

It seems that the changes to the subroutines that worked for v7 are not working for v8.

Changing the MAX_BASKET_SIZE is not mentioned in the link to the post you provided, do you know if @shivanesh.rao implemented these changes to the 3 subroutines also?

Many thanks,

Neil

Hi,

I have just noticed that in the latest cas file I attached the FORTRAN file containing the 3 modified subroutines was commented out.

Just for clarity this is just something I did as part of some testing and I do not usually have this commented and did not when getting the error messages linked.

Sorry if this caused any confusion.

Thanks,

Neil

I have been experimenting with different options for SCHEME FOR ADVECTION OF TRACERS (this is still on V7 as I can't get v8 to work).

Originally I went with your suggestion Chi-Tuan of 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 as I believed this was part of the solution to the issue to the MAX_BASKET_SIZE size I was having. As it turns out after some days of simulations my model is unexpectedly crashing at a certain time point. Looking at the output I notice an issue with my salinity tracer which is set to a value of 35 at my ocean boundary at some points in the model (tidal flats and coastline) it goes above this up to around 40. Although there is no specific error in the output I believe the error is linked to the tracer as when tracers are removed the model runs fine. I tried to attach the results file up to just before the crash but it is to large to upload. If anyone wants to take a look I can share via onedrive.

I have since tried SCHEME FOR ADVECTION OF TRACERS = 13;13;13;13;13;13;13;13;13;13 or 14;14;14;14;14;14;14;14;14;14 + SCHEME OPTION FOR ADVECTION OF TRACERS = 1;1;1;1;1;1;1;1;1;1 as per your suggestions Chi-Tuan but my model won't start running with these at all. I have read the manual in detail and experimented with multiple options with tidal flats as detailed in the manual to try to get these to work but with no success. I also tried SCHEME FOR ADVECTION OF TRACER = 15 but again the model won't run. I have attached a typical output/error file but it doesn't suggest much to me.

The only thing that seems to work is reverting back to SCHEME FOR ADVECTION OF TRACERS = 1;1;1;1;1;1;1;1;1;1. I understand that this is not recommended for tidal flats but the results at least look OK.

Any comments very welcome,

Many thanks,

Neil