! VVINIT Initialisation d'un tracé de champ de vecteurs
  ! VVECTR Tracé de lignes courants
  ! VVSETx **Définition des paramètres vectors_params

  	parameter (mize=21, n=25, iwsize = 2*m*n)
  	real u(m,n), v(m,n), wrk(iwsize
  	gisize = 2.0*pi/float(m)
  	gjsize = 2.0*pi/float(n)
  	do j = 1,n
  	do  i = 1,m
  	  u(i,j) = cos(gisize*(float(i)-1.0))	! Calcul valeurs horizontales
  	  v(i,j) = cos(gjsize*(float(j)-1.0))	! et verticales vecteurs u et v
  	enddo
  	enddo
  	idm = 0
  	rdm = 0.0
  	call vvinit(u,m,v,m,rdm,idm,m,n,rdm,idm)
  	call vvectr(u,v,rdm,idm,idm,rdm)	! Trace des vecteurs