/DEBUGGER=1
/----------------------------------------------------------------------
/   TELEMAC 3D : ESTUAIRE GIRONDE
/                ***  9 plans  ***
/----------------------------------------------------------------------
/
/NUMERO DE VERSION : V8P2R1
/TITRE = 'TELEMAC 3D : ESTUAIRE GIRONDE'
/
/--------------------------------------------------
/     OPTIONS GENERALES
/--------------------------------------------------
FORTRAN FILE = 'user_fortranR45'
BOUNDARY CONDITIONS FILE = 'girgc93.cli'
GEOMETRY FILE = 'girxl_GC_MOR1808tm_93V5.slf'
HARMONIC CONSTANTS FILE = 'GRGCNEA.CL'

ASCII ATMOSPHERIC DATA FILE =
   'forcage_atmos_OM.dat'
LIQUID BOUNDARIES FILE =
   'forcage_FRLIQ_ST.dat'
FORMATTED DATA FILE 1 =
   'BATHY_ZET_nouvelle_bathy_L93.dat'

3D RESULT FILE  =
   'gir3dNH_ML1_RE45_t3d-waq-gaia_20250219_3j.slf'
2D RESULT FILE  =
   'gir2dNH_ML1_RE45_t3d-waq-gaia_20250219_3j.slf'

COMPUTATION CONTINUED = YES
INITIAL TIME SET TO ZERO = YES
PREVIOUS COMPUTATION FILE =
   'init_gironde.slf'

NUMBER OF TIME STEPS = 8640
TIME STEP = 30.
ORIGINAL DATE OF TIME = 2025;2;19

VARIABLES FOR 2D GRAPHIC PRINTOUTS = U,V,B,S,W,H,X,Y,P,TA*
VARIABLES FOR 3D GRAPHIC PRINTOUTS =  Z,U,V,TA* /,P1,P2
LISTING PRINTOUT PERIOD = 20
GRAPHIC PRINTOUT PERIOD = 120 /20

INFORMATION ABOUT MASS-BALANCE FOR EACH LISTING PRINTOUT = YES
MASS-BALANCE = YES
NUMBER OF PRIVATE ARRAYS = 2

/--------------------------------------------------
/     CONDITIONS AUX LIMITES / INITIALES
/--------------------------------------------------
INITIAL CONDITIONS = 'CONSTANT ELEVATION'
INITIAL ELEVATION = 3.5
MINIMAL VALUE FOR DEPTH = 0.1

PRESCRIBED FLOWRATES     				= 0.;1030.;500.
PRESCRIBED VELOCITIES					= 0;0;0
VELOCITY PROFILES 						= 1;1;1
VELOCITY VERTICAL PROFILES              = 2;2;2
OPTION FOR LIQUID BOUNDARIES 			= 2;1;1
OPTION FOR TIDAL BOUNDARY CONDITIONS	= 1;0;0
TIDAL DATA BASE							= 3
COEFFICIENT TO CALIBRATE SEA LEVEL		= 0.35 /0.469
COEFFICIENT TO CALIBRATE TIDAL RANGE	= 0.95 /0.90

TREATMENT OF FLUXES AT THE BOUNDARIES = 2;2;2

/--------------------------------------------------
/     DISCRETISATION VERTICALE
/--------------------------------------------------
NUMBER OF HORIZONTAL LEVELS = 9
MESH TRANSFORMATION = 2

/--------------------------------------------------
/     PARAMETRES PHYSIQUES
/--------------------------------------------------
/NUMBER OF BOTTOM SMOOTHINGS = 6
/TURBULENCE MODEL FOR THE BOTTOM = 3
FRICTION COEFFICIENT FOR THE BOTTOM = 60.
LAW OF BOTTOM FRICTION = 5
VERTICAL TURBULENCE MODEL = 2
/ LONG DE MELANGE AVEC JET (NOUVEAU CALCUL)
/ Ici, version Tsanis
MIXING LENGTH MODEL = 6 /1
DAMPING FUNCTION=3
HORIZONTAL TURBULENCE MODEL = 1
COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES = 1
/COEFFICIENT FOR HORIZONTAL DIFFUSION OF TRACERS = 5
/PRANDTL NUMBER=8

WIND = YES
COEFFICIENT OF WIND INFLUENCE = 1.0
OPTION FOR WIND = 2
COEFFICIENT OF WIND INFLUENCE VARYING WITH WIND SPEED = YES
THRESHOLD DEPTH FOR WIND = 1.

/ Ajout de la pression pour les surcotes
AIR PRESSURE = YES
VALUE OF ATMOSPHERIC PRESSURE = 101300.

TIDAL FLATS = YES
TREATMENT OF NEGATIVE DEPTHS = 2
OPTION FOR THE TREATMENT OF TIDAL FLATS = 1

CORIOLIS=YES
CORIOLIS COEFFICIENT : 0.00010281

/--------------------------------------------------
/     PARAMETERES NUMERIQUES
/--------------------------------------------------
SCHEME FOR ADVECTION OF VELOCITIES = 1
SCHEME FOR ADVECTION OF TRACERS  = 14
SCHEME OPTION FOR ADVECTION OF TRACERS = 1
/SCHEME FOR DIFFUSION OF TRACERS  = 1
/SCHEME FOR ADVECTION OF K-EPSILON = 14
PRECONDITIONING FOR DIFFUSION OF TRACERS = 34 /2
/PRECONDITIONING FOR DIFFUSION OF K-EPSILON : 34
SCHEME FOR ADVECTION OF DEPTH = 5
NON-HYDROSTATIC VERSION = NO
FREE SURFACE GRADIENT COMPATIBILITY = 0.9

ACCURACY FOR DIFFUSION OF VELOCITIES = 1.E-6
ACCURACY FOR DIFFUSION OF TRACERS = 1.E-6
PRECONDITIONING FOR DIFFUSION OF VELOCITIES = 2
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF VELOCITIES = 25000
MAXIMUM NUMBER OF ITERATIONS FOR PPE = 25000
MAXIMUM NUMBER OF ITERATIONS FOR PROPAGATION  = 25000
INITIAL GUESS FOR DEPTH  = 1
ACCURACY FOR PROPAGATION = 1.E-5
ACCURACY FOR PPE = 1.E-5
IMPLICITATION FOR DEPTH = 1.
IMPLICITATION FOR VELOCITIES = 1.
IMPLICITATION FOR DIFFUSION = 1.
MASS-LUMPING FOR DEPTH = 1
MASS-LUMPING FOR VELOCITIES = 1
MASS-LUMPING FOR DIFFUSION = 1
SOLVER FOR PROPAGATION = 1
PRECONDITIONING FOR PROPAGATION = 2
/SOLVER FOR VERTICAL VELOCITY = 7
MATRIX STORAGE = 3
MAXIMUM NUMBER OF ITERATIONS FOR ADVECTION SCHEMES = 2000
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF TRACERS = 240
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF VELOCITIES = 240
/HYDROSTATIC INCONSISTENCY FILTER : YES

/--------------------------------------------------
/     TRACEURS
/--------------------------------------------------
NAMES OF TRACERS =
   'SALINITE        G/L             ';
   'TEMPERATURE     DEGC            '
/ 3 traceurs sont ensuite implicites: O2, L, NH4
/ les 2 traceurs de GAIA COH1 et COH2 sont aussi implicites
/ mais gérés dans GAIA.cas
DENSITY LAW = 3

/ Valeurs initiales pour S,T,O2,L,NH4 >> cf. user_condi3d_trac
/ Les valeurs init. pour COH1 et COH2 sont données dans GAIA.cas
INITIAL VALUES OF TRACERS	= 35.;10.;8.;0.;0.

/ frliq: OCEAN, GARONNE, DORDOGNE
/ traceurs: S,T,O2,L,NH4 (cf. FRLIQ)
/ les forçages pour COH1 et COH2 sont faits dans GAIA.cas
/PRESCRIBED TRACERS VALUES	= 35.;12.;8.;0.;0.00;
/							  0. ;15.;8.;3.;0.054;
/							  0. ;15.;8.;3.;0.054
PRESCRIBED TRACERS VALUES	= 0;0;0;0;0;0;0;0;0;0;0;0;0;0;0

TRACERS VERTICAL PROFILES = 1;1;1; 1;1;1; 1;1;1; 1;1;1; 1;1;1

/ Couplage WAQTEL pour module THERMIC + O2
COUPLING WITH = WAQTEL,GAIA
WAQTEL STEERING FILE = 'waq_GXL_O2.cas'
/'waq_GXL_thermic-O2-ray.cas'
WATER QUALITY PROCESS = 22 /22 = 11 (THEMIC) x 2 (O2)
GAIA STEERING FILE = 'sis2gaiaTH.cas'