Discussion :

[saida29]

Bonjour
lors de compilation de mon code je trouve le message suivant dans le subroutine qui calcul la temperature. merci bien pour vos aide

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
subroutine col(u,v,g,geq,th,theta,rho,
     &omegat,t_in,tw,w,cx,cy,n,m)
	real g(0:8,0:n,0:m),geq(0:8,0:n,0:m),th(0:n,0:m)
	real w(0:8),cx(0:8),cy(0:8)
	real u(0:n,0:m),v(0:n,0:m),theta(0:n,0:m)
        real rho(0:n,0:m)
        real t_in(m/2),tw
        tw=0.0
        t_in(m/2)=1
         do i=0,n
	 do j=0,m
 
         t1=u(i,j)*u(i,j)+v(i,j)*v(i,j)
         g(0,i,j)= -(2.d0/3.d0)*th(i,j)*t1(0,i,j)*rho(i,j)
         do k=1,4
         t2=u(i,j)*cx(k)+v(i,j)*cy(k)
          g(k,i,j)= w(k)* th(i,j)*rho(i,j)*(3.d0/2.d0
     &+(3.d0/2.d0)*t1+(9.d0/2.d0)*t2*t2-(3.d0/2.d0)*t2)
         end do
         do k=5,8
         t2=u(i,j)*cx(k)+v(i,j)*cy(k)
          g(k,i,j)= w(k)* th(i,j)*rho(i,j)*(3+(2.d0/3.d0)*t1+
     &(9.d0/2.d0)*t2*t2-(3.d0/2.d0)*t2)
          end do
         end do
         end do
	return
	end

modthereco.for:278.9:

g(0,i,j)= -(2.d0/3.d0)*th(i,j)*t1(0,i,j)*rho(i,j)
1
Error: Unclassifiable statement at (1)
merci

[LittleOwl]

Bonjour,

Le type de la variable t1 n'a pas été déclaré. Le type real lui est donc affecté implicitement à la ligne 14 de votre code. Or à la ligne suivante, t1 apparait dans l'expression de g(0,i,j) comme une référence à un tableau. D'où le message d'erreur.

En supposant que t1 est une variable locale, de type real, qui sert de variable auxiliaire pour vos calculs, le code suivant se compile correctement :

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
 
      subroutine col(u,v,g,th,rho,t_in,tw,w,cx,cy,n,m)
 
      real g(0:8,0:n,0:m),th(0:n,0:m)
      real w(0:8),cx(0:8),cy(0:8)
      real u(0:n,0:m),v(0:n,0:m)
      real rho(0:n,0:m)
      real t_in(m/2),tw
      real t1
 
       tw=0.0
       t_in(m/2)=1
       do i=0,n
        do j=0,m
         t1=u(i,j)*u(i,j)+v(i,j)*v(i,j)         
    	 g(0,i,j) = -(2.d0/3.d0)*th(i,j)*t1*rho(i,j)
 
         do k=1,4
          t2=u(i,j)*cx(k)+v(i,j)*cy(k)
          g(k,i,j)= w(k)* th(i,j)*rho(i,j)*(3.d0/2.d0
     &+(3.d0/2.d0)*t1+(9.d0/2.d0)*t2*t2-(3.d0/2.d0)*t2)
         end do
 
         do k=5,8
          t2=u(i,j)*cx(k)+v(i,j)*cy(k)
          g(k,i,j)= w(k)* th(i,j)*rho(i,j)*(3+(2.d0/3.d0)*t1+
     &(9.d0/2.d0)*t2*t2-(3.d0/2.d0)*t2)
          end do
 
       end do
      end do
 
      end subroutine col

Par ailleurs, votre subroutine ne faisant pas appel aux paramètres geq, omegat et theta, je les ai retiré de la liste des paramètres et des déclarations.