parameter (m=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 des valeurs
  horizontales
  v(i,j) = cos(gjsize*(float(j)-1.0))	! et verticales des vecteurs u et  v
  enddo
  enddo
  idm = 0
  rdm = 0.0
  call stinit(u,m,v,m,rdm,idm,m,n,wrk,iwsize)
  call stream(u,v,rdm,idm,idm,wrk)	!Trace des lignes de courant