boites pour le FOF_6133 (le halo peut-être "elliptique").

Les repertoires avec les fichiers grafic correspondants sont sur la meso dans aubertdom/Simus/resim/FOF_6133. Chacun de ces repertoires contient un champ de densité de baryons (ic_deltab) et 3 trois champs de vitesses (ic_velcx,y,z) de la matière noire. Si les vitesses des baryons sont necessaires, il suffit de dupliquer les champs de vitesses. La structure des fichiers grafic est expliquée sur le site web, mais en cas de souci yaka demander. Sinon il y a aussi Grafic2Gadget.f90 dans la section conversion de format su site interne.

En pratique, j’ai procédé de la façon suivante:

1/ j’ai centré la 100Mpc/1024 sur le barycentre des nuages de points à z=40 et z=0 du FOF_6133. Ce point a pour coordonnées:

 [3021.77, 92677.7, 12936.8] kpc/h
 [0.0302177,0.926777,0.129368] en unités de taille de boîte
 [31, 949, 132] en pixels

2/ j’ai dégradé cette boite jusqu’à 323 (au sens de la moyenne sur 8 point) 3/ j’ai extrait des sous-boites pour chacune de ces 100 Mpc/h dégradées. Elles ont donc la même résolution, mais moins de particules. Toutes ces boites sont centrées sur ce barycentre dans les conditions initiales.

flowchart of zoom procedure

Flowchart

Tout ca est résumé dans l’organigrame ci dessus, qui donne aussi la structure des répertoires. Au passage, certains repertoires s’appellent 12Mpc mais ils correspondent à des boites de 12.5 Mpc/h.

zoom procedure in init

Zoom strategy

Ci dessus une image-exemple de CI imbriquées. Elle correspond à la diagonale des 1283 dans le diagramme ci dessus avec une résolution centrale correspondant à celle d’une 100Mpc/1024.

Au passage, on pourrait pousser jusqu’à une résolution equivalente de 20483 au centre dans une boite 6 Mpc/h puisque le halo rentre completement dedans dans les conditions initiales. Mais bon, on va peut-etre d’abord regarder ce qui se passe à plus basse résolution.