Algorithme 3 | ateliersiencesvilge

ALGORITHME 3

Analyse

Le but du programme est d’estimer le lien entre l’arrivée de pluies importantes et une baisse significative de la pression atmosphérique les précédant.

Comme base d’information sur laquelle s’appuyer, nous avons un quadrillage de la zone étudiée qui nous fournit pour chaque case du quadrillage, un relevé horaire des précipitations et de la pression atmosphérique.

Considérant qu’une constante de référence définissant une forte pluie nous a été fournie, nous commençons par recenser l’ensemble des cases avec l’heure associée où le relevé de précipitations dépassera la constante de référence désignée sous le nom de "seuil d’averse".

Le recensement des cases se fera dans un tableau d'entités "averse" qui contiendra toutes les informations caractérisant chacune d'elles. Ce tableau les rangera de manière chronologique.

Après quoi, en parcourant le tableau chronologique, on peut identifier les cases contiguës qui ont reçu une averse avec un "écart de temps" spécifié. La contiguïté est vérifiée par le max des écarts entre les deux abscisses et les deux ordonnées qui ne doit pas dépasser 1, mais doit aussi valoir 1 pour éliminer la superposition des cases.

Ce décalage nous assure une relation de causalité entre la première case arrosée et la seconde. Il faut alors associer les deux averses. Nous obtenons alors une liste chronologique des relations de causalités possibles. On pourra établir d'autres listes de relations de causalité en faisant varier l'écart de temps vers 2 ou 3 heures.

Dans un premier temps, on va étudier la case sensée contenir Bordeaux pour identifier les relations de causalité qu'elle peuple.

Structures de Données :

Deux tableaux 3 dimensions :
- V(x,y,h) : volume de précipitations dans la case (x,y) à l’heure h
- P(x,y,h) : pression atmosphérique dans la case (x,y) à l’heure h
Tableau 1 dimension d’objets structurés Averse
- Averse comprend les champs :
  - X appartient aux entiers positifs
  - Y appartient aux entiers positifs
  - H appartient aux réels positifs
  - V appartient aux réels positifs
Tableau 1 dimension d’objets structurés Relation Causalité
- Relation Causalité comprend les objets :
  - Averse 1
  - Averse 2
  - Pression
Des constantes :
- Seuil d’Averse appartient aux réels positifs
- Seuil de Pression appartient aux réels positifs
- Ecart de temps appartient aux réels positifs

Procédures proposées

·Recherche des Averses

Compteur_averse = 0
Pour H allant de 0 à 23
- Pour X allant de 1 à 5:
  - Pour Y allant de 1 à 5:

//On teste pour chaques cases s’il y a eu une averse, si oui on incremente le compteur de 1 et on met les coordonées et l’heure de l’averse dans Averse

Si V(x,y,h) >= Seuil d’Averse alors
- Compteur_averse = Compteur_averse +1
- Averse (compteur_averse) reçoit x,y,h,V(x,y,h)
Fin si
Fin pour
Fin pour
Fin pour
Afficher Averse(1) ; Averse(2) ;… ;Averse(Compteur_averse)

·Recherche des Causalités

Compteur_causalité = 0
Ecart_de_Temps = 1 (heure)

//Si il y a causalité entre une averse dans une zone qui se deplace dans une zone proche, on increment le compteur et on rentre les coordonnées du deplacement et la pression dans relation Causalité

Pour i allant de 0 à Compteur_averse – 1:
- Pour j allant de i +1 à Compteur_averse :
  - Si averse(j).h = averse(i).h + ecart_de_temps et
    maximum ( val_abs ( averse(j).x - averse(i).x) , val_abs (averse(j).y – averse(i).y ) ) = 1 alors : //
  - Relation_Causalité ( Compteur_Causalité ).averse_1 = averse(i) ,
  - Relation_Causalité ( Compteur_Causalité ).averse_2 = averse(j) ,
  - Relation_Causalité ( Compteur_Causalité ).pression
    = P(averse(j).x,averse(j).y,averse(j).h – ecart_de_temps) ,
  - Compteur_Causalité = Compteur_Causalité + 1;
  - Fin si
- Fin pour
Fin de boucle i

· Etude à Bordeaux

X_Bordeaux = 1,
Y_Bordeaux = 3;
//Si il y a causalité a Bordeaux, on en imprime les parametres de pression et d’averse
Boucle de i=0 jusqu’à Compteur_causalité – 1
- Si Relation_Causalité(i).averse_1.x = X_Bordeaux et
  Relation_Causalité(i).averse_1.y = Y_Bordeaux alors
  - Impression des valeurs des averses et de pression de la relation de causalité x,y,
- Fin si
Fin de boucle i

Pistes possibles :

Faire varier l’écart de temps
Comptabiliser les couples i,j et étudier la valeur de la pression , puis faire une analyse statistique des pressions obtenues
Faire varier les seuils d’averse et de pression pour identifier les causalités statistiquement établies

Retour