Shoreline S Donnes dentre Structure gnrale Shoreline S

  • Slides: 16
Download presentation
Shoreline. S Données d’entrée

Shoreline. S Données d’entrée

Structure générale Shoreline. S • Modèle execute par un script principal par modèle •

Structure générale Shoreline. S • Modèle execute par un script principal par modèle • Données d’entrée dans une structure Matlab nommé S %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Mettre les fonctions Shoreline. S disponible dans Matlab addpath(genpath('D: OSSShoreline. S_trunkfunctions' )); clear S; warning off; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% MODEL INPUT PARAMETERS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% S=struct; S. LDBcoastline='. . dataxy 1967_St. N. txt' ; % LDB with initial coastline shape ([Nx 2] ASCII FILE WITHOUT HEADER) S. struct=1; % switch for using hard structures Paramètres du modèle S. LDBstructures='. . dataijmuiden_port. txt' ; % LDB with hard structures ([Nx 2] ASCII FILE WITHOUT HEADER) S. LDBplot = {'. . dataxy 1967. txt', '1967', 'k--'; . . . '. . dataxy 2007. txt', '2007', 'r--'}; % Filenames, dates and line codes of observed coastlines S. xlimits=[90000 110000]; % x-limits of plot area S. ylimits=[490000 505000]; % y-limits of plot area. . . %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% RUN SHORELINES MODEL %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Exécution du [S, O]=Shoreline. S(S); save('. . outputResults. mat', 'S', 'O'); modèle + enrégistrement de résultats

Structure générale Shoreline. S • Tous les paramètres disponibles sont précisés dans le fichier

Structure générale Shoreline. S • Tous les paramètres disponibles sont précisés dans le fichier Shoreline. S. m %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% DEFAULT INPUT PARAMETERS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% S. Hso=1; % S. phiw 0=330*pi/180; % S. spread=90; % S. WVCfile=''; % S. Waveclimfile=''; % S. Wavematfile=''; S. Wavecorr=0; S. d=10; % S. ddeep=25; % S. dnearshore=8; % S. ob=0; % S. LDBcoastline=''; % S. phif=[]; % S. trform='CERC'; % S. b=1 e 6; % S. qscal=1; % S. tper=6; % S. d 50=2. 0 e-4; % S. porosity=0. 4; % S. tanbeta=0. 03; % S. rhos=2650; % S. rhow=1025; . . . wave deep wave height [m] water wave angle [? N] spreading [? ] (wave_dir from range: S. phiw 0 +/- 0. 5*S. spread) time-series <-leave empty to use wave parameters ('S. Hso', 'S. phiw 0' climate file active profile height [m] Waterdepth the location of wave climate, corresponding with S. Hso [m] Waterdepth at the 'dynamic boundary', corresponding with S. phif [m] Switch for activation of observation point at which wave height will be c LDB with initial coastline shape ([Nx 2] ASCII FILE WITHOUT HEADER) <- lea Orientation of the foreshore [? N] (in radians) <- only relevant for 'KAMP switch for transport formulation (e. g. S. trform='CERC', 'KAMP', 'MILH' or CERC : coeff in simple cerc formula Calibration factor of the transport (works for all transport formulas) KAMP & MILH : peak wave period [s] KAMP & MILH & VR 14 : median grain diameter [m] KAMP & MILH & VR 14 : S. porosity (typically 0. 4) [-] KAMP & MILH & VR 14 : mean bed slope [ratio 1/slope] KAMP & MILH & VR 14 : density of sand [kg/m 3]

Données d’entrée • Géométrie • Trait de côte initiale • Trait de côte observés

Données d’entrée • Géométrie • Trait de côte initiale • Trait de côte observés • Structures • Conditions aux limites • Climat / série temporelle de houle • Paramètres d’exécution • Temps • Transport littoral • Numériques

Géométrie • Trait de côte initial • S. LDBcoastline=‘spit. txt’; • Fichier text avec

Géométrie • Trait de côte initial • S. LDBcoastline=‘spit. txt’; • Fichier text avec 2 colonnes X, Y; sans en-tête Spit. txt 0. 0. 0. 5000. 10000.

Géométrie • Traits de côte observés • S. LDBplot = {'xy 1967. txt', '1967',

Géométrie • Traits de côte observés • S. LDBplot = {'xy 1967. txt', '1967', 'k--'; 'xy 2007. txt', '2007', 'r--'}; • {<nom fichier>, <entrée légende>, <style de courbe>} • Permet de visualiser le résultat du modèle par rapport aux mesures • Fichier text avec 2 colonnes X, Y; sans en-tête xy 1967. txt 0. 0. 0. 5000. 10000.

Géométrie • Structures: épis, brise-lames; côtes rocheuses • S. struct=1 • S. LDBstructures=‘briselame. txt’

Géométrie • Structures: épis, brise-lames; côtes rocheuses • S. struct=1 • S. LDBstructures=‘briselame. txt’ • Fichier text avec 2 colonnes X, Y; sans en-tête • Points dans le sens horaire briselames. txt 0. 0. 0. 200. 10. 0. Nan nan 200. 0. 200. 210. 0. Nan nan 400. 0. 400. 200. 410. 0. Épi 1 Épi 2 Épi 3

Géométrie Epi 3 Epi 2 Epi 1

Géométrie Epi 3 Epi 2 Epi 1

Conditions aux limites • Conditions constantes • S. Hso=1; • S. phiw 0=pi; •

Conditions aux limites • Conditions constantes • S. Hso=1; • S. phiw 0=pi; • S. spread=90; % hauteur de vague en m % direction de vague en degrees % angle diffusion houle en degrees

Conditions aux limites • Série temporelle • S. WVCfile=‘houle. txt’ % Hs 198401010000 198401010100

Conditions aux limites • Série temporelle • S. WVCfile=‘houle. txt’ % Hs 198401010000 198401010100 198401010200 198401010300 198401010400 198401010500 198401010600 198401010700 198401010800. . . Tp 0. 33 0. 315 0. 327 0. 362 0. 469 0. 602 0. 711 0. 942 1. 095 Dirm 3. 2 4. 8 4. 7 2. 6 3 3. 9 4. 2 3. 4 358. 5 342 315. 8 284. 3 279. 5 280. 6 280

Conditions aux limites • Climat de houle 2 5 8 10 • Séries temporelles

Conditions aux limites • Climat de houle 2 5 8 10 • Séries temporelles de la houle converties en valeurs de flux d’énergie • Subdivision en blocks avec contribution égale • Conditions moyennes de chaque block 7 4 1 3 6 • Probabilité d’occurrence 9 • Script EFShow. m

Conditions aux limites • Climat de houle • S. Waveclimfile=‘wavecon. txt’ % wavecon. txt:

Conditions aux limites • Climat de houle • S. Waveclimfile=‘wavecon. txt’ % wavecon. txt: <hauteur vague>, <periode>, <direction>, <probabilité> 8. 5367413 e-01 1. 0295535 e+00 1. 1296605 e+00 1. 2189080 e+00 1. 2891085 e+00 1. 3670497 e+00 1. 4380877 e+00 1. 4096041 e+00 1. 4837256 e+00 1. 4988507 e+00 1. 4126613 e+00 4. 7146347 e+00 5. 2167002 e+00 5. 2935060 e+00 5. 5028207 e+00 5. 5672058 e+00 5. 6226354 e+00 5. 8144669 e+00 5. 6367023 e+00 5. 7628238 e+00 5. 6689913 e+00 5. 4645597 e+00 2. 4666401 e+02 2. 5784366 e+02 2. 6148712 e+02 2. 6387059 e+02 2. 6590064 e+02 2. 6768376 e+02 2. 6940951 e+02 2. 7131450 e+02 2. 7361454 e+02 2. 7653712 e+02 2. 8119845 e+02 3. 0420079 e+01 1. 9330986 e+01 1. 5659597 e+01 1. 3336882 e+01 1. 1613577 e+01 1. 0339830 e+01 9. 0286195 e+00 9. 7778825 e+00 8. 4666722 e+00 8. 5415985 e+00 1. 0002661 e+01

Conditions aux limites Q_end Q_start • Type de conditions aux limites de transport •

Conditions aux limites Q_end Q_start • Type de conditions aux limites de transport • S. boundary_condition_start='FIXD 2'; (mur) ou FIXD 2 (Neumann) • S. boundary_condition_end='FIXD 2'; % CTAN (angle côte fixe), FUNC (Dirichlet), FIXD • Conditions fixes de transport (FUNC) • S. QS_start=50000; • S. QS_end=30000; % [m 3/année] • Conditions variables dans le temps (FUNC; fichier texte, 3 colonnes, pas d’entête) • S. BCfile=''; % Conditions en function de temps (<temps> , <QS_start>, <Qs_end>)

Paramétres d’exécution • Gestion de temps • • S. tc=1; % pas de temps

Paramétres d’exécution • Gestion de temps • • S. tc=1; % pas de temps adaptif S. reftime='2020 -01 -01'; % temps de reference S. endofsimulation='2040 -01 -01'; % point finale de simulation S. start=0; % debut en jours depuis S. refdate start reftime endofsimulation t

Paramétres d’exécution • Gestion de géométrie • S. ds 0=100; • S. d =

Paramétres d’exécution • Gestion de géométrie • S. ds 0=100; • S. d = 10; % resolution spatial initial [m] % hauteur actif du profil de côte transversal • Gestion de transport sédimentaire • S. trform='CERC'; • S. qscale=1. 0; • S. b=0. 5 e 6; % formule de transport (‘CERC', 'KAMP', 'MILH' ou 'VR 14') % facteur de calibration transports % B dans la formule CERC

Mise en page / enrégistrement des résultats S. plotvisible=1; S. outputdir='Output'; % output directory

Mise en page / enrégistrement des résultats S. plotvisible=1; S. outputdir='Output'; % output directory S. SLplot={'1970 -08 -01', '3 m', 'k--‘ % <date>, <entrée légende>, <style courbe> '1970 -11 -01', '6 m', 'b-' '1971 -04 -29', '1 yr', 'g-' ‘ 1972 -04 -29', '2 yr', 'K-' '1973 -04 -29', '3 yr', 'r-'}; S. usefill = 0; % switch that can be used to only plot lines instead of the fill S. fignryear=12; % number of saved png files per year S. plotinterval=1; % Time interval of plotting updated coastline on screen [day] S. storageinterval=50; % Time interval of storage of output file ('output. mat'; [day])