Esope source files

   1 : C NORMVI    SOURCE    CHAT      05/01/13    02:00:47     5004           
   2 :       SUBROUTINE NORMVI(NDIM,UE,XNORM,UXY)
   3 : ***************************************************************************
   4 : *** SP 'NORMVI' : calcule de maniere explicite la vitesse, sa norme et
   5 : *** 'UXY' (cf. ci-desssous) au pt de depart. Module valable dans element
   6 : *** reference ou element reel.
   7 : *** 
   8 : *** APPELES 1 = aucun
   9 : *** APPELES 2 = aucun
  10 : ***
  11 : *** E = 'NDIM' dimension de l'espace
  12 : ***     'UE' vitesse dans element considere
  13 : ***                       
  14 : *** S =  XNORM' norme de la vitesse
  15 : ***     'UXY' -> racine carre de UX*UX+UY*UY (CAS 3D)
  16 : ***
  17 : ***         Auteur Cyril Nou
  18 : ****************************************************************************          
  19 : 
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  22 :       DIMENSION UE(3)
  23 :       UXY=0.D0
  24 :       XNORM=0.D0
  25 : ***   cas 2D
  26 :       DO 10 I=1,NDIM
  27 :          XNORM=XNORM+UE(I)**2
  28 :  10   CONTINUE
  29 :       XNORM=SQRT(XNORM)
  30 : ***   cas 3D
  31 :       IF (NDIM.EQ.3) THEN
  32 :          DO 20 I=1,NDIM-1
  33 :             UXY=UXY+UE(I)**2
  34 :  20      CONTINUE
  35 :          UXY=SQRT(UXY)
  36 :       ENDIF
  37 :       RETURN                                                                    
  38 :       END                                                                       
  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 :