open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi everyone,

I've searched on the forums but haven't found a solution to what I'm trying to do. It may be out there, apologies if this is a re-post.

I'm in the process of re-running some simulations of tides in the North Sea (Through the Pentland Firth and surrounding area). My initial experiment used the TPXO database, which has over a dozen tidal constituents. The model will then be used to determine available energy output from tidal turbines (Effectively just an area of increased bottom friction).

We're trying to get this work published, but have been strongly recommended to re-run the model with only the M2 (And possibly S2) tides.

Does anyone have a solution on how to force single tidal constituents in this area? I've so far tried turning off all other constituents in TPXO.f by setting frequency to zero, but this simply leaves me with the rest of the components as a free surface gradient.

Kind regards,

Ryan

Hi,

To do such specific treatments, you should change the building of the indices array of given constituents and the number of constituents used in the module TPXO (CCIND and NCON) when calling TPXO_BORD. You can look at what is done for schematic tides (TIDALTYPE between 2 and 6) with the building of CCIND2 and NCON2 variables and replace CCIND and NCON by these two variables. In my opinion, what you want to do is between what is usually done with option OPTION FOR TIDAL BOUNDARY CONDITIONS=1 and what is done with option OPTION FOR TIDAL BOUNDARY CONDITIONS=4.

Hope this helps,

Chi-Tuan

Hi,

I do not know if you have tried to change the TPXO module to deal with only M2 constituent. You can try the enclosed file and tell me if it is OK for you (or other users). If it is OK, we can integrate this in the next release.

Hope this helps,

Chi-Tuan

Pham,

Many thanks for your response. I'm in the process of looking into the tpxo.f code. I've had a quick read-through of the code you sent, may I ask what you've changed in the file you've shared?

Ryan

Hi,

If you use the tkdiff program (available for Linux or Windows), you can see the differences between 2 ASCII files. The original file is tpxo.f of version 7.1 release 1 (tag v7p1r1).
The differences are the extra option with TIDALTYPE = 8 (please use the keyword OPTION FOR TIDAL BOUNDARY CONDITIONS with the value 8 in your steering file).

Chi-Tuan

Hi,

I have the same problem as the thread. I would to use only M2 and M4 tides. i am using the tpxo database. I tried to dig into tpxo.f but didn't get what's going on in it. I tried the tpxo_m2.f. But I get this error :

t2dfort.f:2151:9:

USE BIEF
2
t2dfort.f:2154:17:

INTEGER LNG,LU
1
Error: Symbol 'lng' at (1) conflicts with symbol from module 'declarations_special', use-associated at (2)
t2dfort.f:3022:17:

t2dfort.f:3019:9:

USE BIEF
2
t2dfort.f:3022:17:

INTEGER LNG,LU
1
Error: Symbol 'lng' at (1) conflicts with symbol from module 'declarations_special', use-associated at (2)
t2dfort.f:3886:17:

t2dfort.f:3883:9:

USE BIEF
2
t2dfort.f:3886:17:

INTEGER LNG,LU
1
Error: Symbol 'lng' at (1) conflicts with symbol from module 'declarations_special', use-associated at (2)
runCAS: could not compile: tpxo.exe


I've been told that it is possible to just change the parameters in the subroutine nodal, do you know how could I do?

Thank's a lot

Yves

Hi
You try to use a 2 years old fortran file which is no longer compatible with the version of Telemac you use...
Start with your version and adapt the modification done in the example.
Regards

Merci pour la réponse!

Effectivement, je suis censé pouvoir le faire directement avec
OPTION POUR LES CONDITIONS AUX LIMITES DE MAREE = 5

Mais j'ai 2 soucis,
Quand je mets OPTION POUR LES CONDITIONS AUX LIMITES DE MAREE = 5, j'ai cette erreur :
INITIALISATION BASEE SUR TPXO :
ONDES DISPONIBLES POUR MAREES SCHEMATIQUES :
C_ID( 1 ) = m2
C_ID( 2 ) = s2
C_ID( 11 ) = m4
DONNER UN NUMERO DE NOEUD DE FRONTIERE D UNE FRONTIERE LIQUIDEPOUR CALAGE DE PLEINE MER EN MAREES SCHEMATIQUES

Normalement, il faut choisir un noeud avec ICALW.le code détermine automatiquement un indice de nœud situé à la moitié de la numérotation des nœuds de frontière maritime, si il n'est pas déterminé. Pourquoi ça plante alors? Et ou puis-je le déterminer moi-même?


Deuxièmement, si je mets OPTION POUR LES CONDITIONS AUX LIMITES DE MAREE = 8
J'ai l'erreur :
MAUVAISE VALEUR POUR TIDALTYPE = 8
ELLE DOIT ETRE COMPRISE ENTRE 1 ET 6
AVEC LA BASE DE DONNEES DE MAREE TPXO


Alors que dans le code tpxo.f j'ai bien :
ELSEIF(TIDALTYPE.EQ.8) THEN
DO IC = 1,NCON
IF(C_ID(IC).EQ.'m2 ') THEN
I = I+1
WRITE(LU,*) 'C_ID(',IC,') = ',C_ID(IC)
ENDIF

Merci pour l'aide!

Hello, guys. Im sorry for renew this topic, but Im facing the same situation. I would like to simulate with TPXO database only for M2, M4, S2, O1 and K1. Or I would like to put my own tidal database, once I know phase, amplitude of the local tidal constituents... I know it is possible, but I did not see any examples like this.

Did you solve your problem? If yes, may I ask how?

Thanks in advance,