Discussion :

[karimbazi]

Bonjour,
voila j'ai un tableaux qui contient dans le premier colonne différents character (sont des symboles des atomes) rangés d'une manière aléatoire , dans les trois autres colonnes contient les coordonnes x,y, z de ces atomes.
j'aimerais savoir comment extraire des sous tableaux qui contient chacun seulement un types d'atome et ces cordonnes . j'utilise fortran .

voici mon tableau.

C -5.14023279 -0.63308436 -0.00049947
C -3.73883287 -0.63318804 -0.00004484
C -3.03804286 0.58040785 -0.00062098
C -3.73865278 1.79410744 -0.00163128
C -5.14005270 1.79421111 -0.00209632
C -5.84084271 0.58061523 -0.00151219
H -5.67530152 -1.55969174 -0.00004605
H -3.20390165 -1.55987461 0.00072242
H -3.20358406 2.72071484 -0.00205950
H -5.67498391 2.72089768 -0.00287525
H -6.91084265 0.58069438 -0.00187164
O -1.60804293 0.58030204 -0.00017796
H -1.28776259 0.57968679 -0.90517527

donc par exemple je vous savoir un tableaux qui contient que les H en premier colonnes et dans le reste les cordonnes .

merci d'avoir rependu.

[bobbyboy]

Bonjour,

comme j'imagine que ton tableau est dans un fichier, une routine de ce style conviendrai :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
 
subroutine extraire_type(type,sous_tab)
   implicit none
 
   ! variables d'entrée/sortie
   ! -----------------------
   character*1, intent(in) :: type               ! type d'atome à extraire
   real*8, intent(out)     :: tableau(100,3)     ! tableau qui contiendra les résultats
 
   ! variables locales
   ! ----------------
   integer :: i,ios
   character*1 :: tmp_type
   real*8 :: tmp_coord(3)
 
   open(unit=10,file='fichier_contenant_la_tableau_total',status='old')
 
 1 format(A,3(E15.8))
 
   ios = 0
   i = 1
   sous_tab = 0.0D00
 
   do while (ios .eq. 0)
      read(11,1,iostat=ios) tmp_type , tmp_coord(:)
 
      if (tmp_type .eq. type) then
         sous_tab(i,:) = tmp_coord(:)
         i = i + 1
      endif
   enddo
 
   ! ici :
   ! la variable i contient le nombre d'atomes de type recherché
   ! la variable sous_tab contient les coordonnées de ces atomes
 
   return
end subroutine extraire_type