Discussion :

[kefi racha]

salut,
j'ai ecrit ce probleme en fortran mais ilya des fautes
je ne sais pas comment corriger le
svp aidez moi et merciiii ^_^

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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
 program  stationnaire
         implicit none I,J,Nmax,XL,XH,S,NI,NJ,XK,DX,DY
         dimension X(Nmax),Y(Nmax),TETA(NXmax,NYmax),QX(NXF,NYF),QY(NXF,NYF)
         parameter (NXmax=200,NYmax=200)
         XL=40.0
         XH=30.0
         NI=12
         NJ=12
         S=0
         XK=1.0E3
         DX=XL/float(NI-1)
         DY=XH/float(NJ-1)
         write(*,*)'longueur:',XL
         write(*,*)'largeur:',XH
         pause
         X(1)=0.0
         Y(1)=0.0
         do I=2,NI
         X(I)=X(I-1)+DX
         enddo
         do J=2,NJ
         Y(J)=Y(J-1)+DY
         enddo
         !conditions initiales
         do I=1,NI
         do J=1,NJ
         TETA(I,J)=0.0
         enddo
         enddo
         ! conditions aux limites
         J=1 !sud
         do I=1,NI
         TETA(I,J)=10.0
         enddo
         J=NJ !nord
         do I=1,NI
         TETA(I,J)=10.0
         enddo
         I=NI !EST
         do J=1,NJ
         TETA(I,J)=10.0+1.0E-1*sin(3.14159*Y(J)/XH)
         enddo
         I=1  !ouest
         do J=1,NJ
         TETA(I,J)=TETA(I+1,J)
         enddo
         !write(*,*)' les conditions initiales'
         !do J=NJ,1,-1
         !write(*,*) (TETA(I,J),I=1,NI)
         !enddo
         !pause
         ! CALCUL DU FLUX THERMIQUE
          I=1
        DO J=2,NJ-1
		QX(I,J)=-XK*(TETA(I+1,J)-TETA(I,J))/(DX)
		QY(I,J)=-XK*(TETA(I,J+1)-TETA(I,J-1))/(2.*DY)
       ENDDO
       I=NI
        DO J=2,NJ-1
		QX(I,J)=-XK*(TETA(I,J)-TETA(I-1,J))/(DX)
		QY(I,J)=-XK*(TETA(I,J+1)-TETA(I,J-1))/(2.*DY)
        ENDDO
        J=1
        DO I=2,NI-1
		QX(I,J)=-XK*(TETA(I+1,J)-TETA(I-1,J))/(2.*DX)
		QY(I,J)=-XK*(TETA(I,J+1)-TETA(I,J))/(DY)
       ENDDO
        J=NJ
       DO I=2,NI-1
		QX(I,J)=-XK*(TETA(I+1,J)-TETA(I-1,J))/(2.*DX)
		QY(I,J)=-XK*(TETA(I,J)-TETA(I,J-1))/(DY)
        ENDDO
		QX(1,1)=-XK*(TETA(2,1)-TETA(1,1))/(DX)
		QY(1,1)=-XK*(TETA(1,2)-TETA(1,1))/(DY)
		QX(NI,1)=-XK*(TETA(NI,1)-TETA(NI-1,1))/(DX)
		QY(NI,1)=-XK*(TETA(NI,1+1)-TETA(NI,1))/(DY)
		QX(NI,NJ)=-XK*(TETA(NI,NJ)-TETA(NI-1,NJ))/(DX)
		QY(NI,NJ)=-XK*(TETA(NI,NJ)-TETA(NI,NJ-1))/(DY)
		QX(1,NJ)=-XK*(TETA(1+1,NJ)-TETA(1,NJ))/(DX)
		QY(1,NJ)=-XK*(TETA(1,NJ)-TETA(1,NJ-1))/(DY)
	DO I=2,NI-1
		DO J=2,NJ-1
				QX(I,J)=-XK*(TETA(I+1,J)-TETA(I-1,J))/(2.*DX)
				QY(I,J)=-XK*(TETA(I,J+1)-TETA(I,J-1))/(2.*DY)
		ENDDO
	ENDDO
	end

ok les erreurs
en couleurs rouge
ok svp aidez moiiiiii

[FR119492]

Salut!

j'ai ecrit ce probleme en fortran mais ilya des fautes

Alors, deux questions:

1. Ce programme sert à quoi?
2. Quelles sont les fautes?

Jean-Marc Blanc

[kefi racha]

Salut!

Alors, deux questions:

1. Ce programme sert à quoi?
2. Quelles sont les fautes?

Jean-Marc Blanc

probleme de diffusion à 2D
les fautes sur les parametres !!!

[Marlan]

Bonjour,

Je te conseilels de lire attentivement cette partie de la FAQ, la plupart de tes erreurs s'y trouvent.

http://fortran.developpez.com/faq/?page=section2

Bonne journée,

Marlan