Basilisk




Exemples simples avec Basilisk de résolution des équations de Saint Venant

Pour installer, voir la page http://basilisk.fr/src/INSTALL

ces exemples sont aussi dupliqués sur http://basilisk.fr/sandbox/M1EMN/README






Premier exemple simple 1D: déferlement sur une plage inclinée

Une pertubation initiale de déplace vers une plage inclinée, elle monte sur la plage et se retire mollement
Fichier pour cet  exemple.

















Deuxième exemple: déferlement sur une plage inclinée


Le même mais en 2D cette fois, pour varier, on met la plage de biais. Fichier pour cet  exemple.





Exemple:




Exemple:



Faire tourner le tsunami de 2004

attention, il faut la base de donnée ETOPO2:

Récupérer la topographie du monde en suivant pas à pas http://gfs.sourceforge.net/wiki/index.php/Xyz2kdt
suivre exactement la procédure du point "Example 1: building a global terrain topography database using the ETOPO2 dataset"
remplacer éventuellement
wget http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO2/ETOPO2v2-2006/ETOPO2v2c/raw_binary/ETOPO2v2c_i2_LSB.zip
par
curl  http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO2/ETOPO2v2-2006/ETOPO2v2c/raw_binary/ETOPO2v2c_i2_LSB.zip >> ETOPO2v2c_i2_LSB.zip


bien utiliser
 /Users/pyl/basilisk/src/kdt/xyz2kdt -v etopo2

Récupérez le code de tsunami de 2004, le modifier pour reproduire le tsunami qui a détruit la civilisation Minoenne

se déplacer en méditerrannée
 // the domain is 45 degrees squared
  L0 = 45.;
  // centered on 15 40 longitude,latitude
  X0 = 15 - L0/2.;
  Y0 = 40. - L0/2.;

dans
event init (i = 0)
changer le chemin de la position de votre fichier avec les données du monde récupéré plus haut
  terrain (zb, " .../WORD/etopo2", NULL);

modifier la position de la faille pour la mettre sur Santorin qui est en   36°25'48.00"N   et  25°24'E  soit 36,43 Nord  25,4 Est
  fault (x = 25.4, y = 36.43,
     depth = 1.e3,
     strike = 0, dip = 12, rake = 90,
     length = 15e3, width = 15e3,
     U = 18*8,
     iterate = adapt);
}


  modifier la fenêtre du zoom pour "eta-zoom.mpg"
          box = {{23.5,34.5},{27.5,37.5}});


lancez et obtenez ce type de résultat...
Fichier pour cet  exemple.