Télécharger la plaquette Cast3M


fr en

Télécharger fgrf3d.eso

Retour à la liste

Numérotation des lignes : 

fgrf3d

  1. C FGRF3D    SOURCE    FD218221  24/02/07    21:15:11     11834          
  2.         subroutine fgrf3d(NDIMG,NDIMPL,NTMAX,DPFDS,icr,JSP,PF,
  3.      #  PF_GRAD_PF,affiche)
  4.  
  5. c       calcul du gradient d une fonction scalaire de l espace generalise de l ecoulement (36)    
  6.  
  7.         implicit real*8 (a-h,o-z)
  8.         implicit integer (i-n)
  9.  
  10. c       i numero du critere en cours 
  11.         integer NDIMG,NTMAX,icr
  12.         integer lnumcrt(NTMAX)
  13.         real*8 DPFDS(NTMAX,NDIMG)
  14.         real*8 jsp(NDIMG,NDIMPL)
  15.         real*8 PF(NTMAX),PF_GRAD_PF(NDIMPL)
  16.         logical affiche, affiche_local
  17.  
  18.  
  19.         integer i,j,k
  20.  
  21.         affiche_local=affiche
  22.  
  23.         if(affiche_local) then
  24.             print*,'On est ds fgrf3d pour le critere:',icr
  25.         end if
  26.  
  27.  
  28.         do i=1,NDIMPL
  29. c           grad f
  30.             PF_GRAD_PF(i)=0.d0            
  31.             do j=1,NDIMG
  32.                 PF_GRAD_PF(i)=PF_GRAD_PF(i)+DPFDS(icr,j)*JSP(j,i)
  33. C                 print*,'DPFDS(',icr,j,')=',DPFDS(icr,j),
  34. C      #                 'JSP(',j,i,')=',JSP(j,i)
  35.             end do
  36. c           f grad f
  37.             PF_GRAD_PF(i)=PF(icr)*PF_GRAD_PF(i)
  38.         end do
  39.  
  40.         if (affiche_local) then
  41.              do i=1,NDIMG
  42.                 write(*,'(A11,2(1x,I2),1X,A1,E10.3)') 'PF.Grad(PF)',
  43.      #          icr,i,'=',PF_GRAD_PF(i)
  44.              end do
  45.              read*
  46.         end if
  47.         return
  48.         end
  49.  
  50.  
  51.  
  52.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales