Discussion :

[jphysic]

Je remercie ceux qui ont répondu à mon dernier message [LittleOwl et François L].

Cette fois-ci, je voudrai que mes conditions initiales [x(1,1),x(2,1),x(3,1),y(1,1),y(2,1),y(3,1)], soient prises sur un cercle, par exemple de rayon r=1. De sorte à avoir une figure de ce genre [voir la 1ère pièce jointe] .


Mais mes résultats tracés avec ORIGINLAB, me donnent ça[2ème pièce jointe].

Je le trouve que C’est une figure effrayante,
N’est-ce-pas ?!!!!
C’est fort possible c’est la manière dont j’ai déclaré mes conditions initiales. Si vous pouvez me guider sur la manière dont je dois les prendre sur un cercle.

Le voici mon code :

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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
PROGRAM VortexEuler
Implicit None
INTEGER :: N,i,j,k,M
REAL, Dimension(3,5):: x,y,som1, som2,r
REAL, Dimension(3):: w
DOUBLE PRECISION :: alph1,alph2,alph3
DOUBLE PRECISION :: dt
DATA w/1.d0,-1.d0,1.d-1/
!<<<<<<<<       Initialisation >>>>>>>>>
N=3
M=5
dt= 1.0e-2
w=1.d0
write(*,*)'alph1'
read(*,*) alph1
write(*,*)'alph2'
read(*,*) alph2
write(*,*)'alph3'
read(*,*) alph3
!<<<<<< conditions initiales >>>>>>
x(1,1)= sin(alph1)
x(2,1)= sin(alph2)
x(3,1)= sin(alph3)
y(1,1)= cos(alph1)
y(2,1)= cos(alph2)
y(3,1)= cos(alph3)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
som1(:,1)=0.d0
som2(:,1)=0.d0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
open(10,File='Vortex euler corrigé.dat')
do k=1,M 
do j=1,N
do i=1,N
if (.not.(i.ne.j)) cycle
r(i,k)=(((x(i,k)-x(j,k))**2)+ ((y(i,k)-y(j,k))**2))
r(i,k)=r(i,k)
!
som1(i,k)=som1(i,k)+ dt*((w(i)*(y(i,k)-y(j,k)))/r(i,k))
!
som1(i,k)=som1(i,k)
!
som2(i,k)=som2(i,k)+ dt*((w(i)*(x(j,k)-x(i,k)))/r(i,k))
!
som2(i,k)=som2(i,k)
!
x(j,k)= x(j,k) + dt* som1(i,k)
x(j,k)=x(j,k)
y(j,k+1)= y(j,k) + dt* som2(i,k)
y(j,k)=y(j,k)
	write(10,*) x(j,k),y(j,k)
  enddo
enddo
enddo
close(10)
END PROGRAM VortexEuler

Cordialement.

[Ehouarn]

Bonjour,

On ne le répètera jamais assez... quand on développe un code, on compile toujours avec toutes les options de débogage (vérification des initialisations de variables, de possibles dépassement de tableaux, etc. ; voir la doc du compilateur).

Dans ce cas tu verras immédiatement certains des défauts de ton programme, par ex. le dépassement de tableau à la ligne

Code :

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

1
2
 
y(j,k+1)= y(j,k) + dt* som2(i,k)

Bonne continuation.