Esope source files

   1 : C VSUR      SOURCE    PV        09/03/12    21:37:01     6325           
   2 :       SUBROUTINE VSUR
   3 : C=====================================================================
   4 : C    OPERATEUR POUR CALCULER LES VECTEURS SURFACIQUES AUX POINTS 
   5 : C               DE GAUSS DES COQUES
   6 : C
   7 : C    CHAM = VSUR  MODL  (MOT1)
   8 : C              MODL : MODELE de calcul , type MMODEL
   9 : C              CHAM : CHAMELEM de sous type VECTEUR contenant
  10 : C                       les JACOBIENS aux POINTS DE GAUSS
  11 : C              MOT1 : MOT 'NORM' : les vecteurs sont normes
  12 : C
  13 : C=====================================================================
  14 : C
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8(A-H,O-Z)
  17 :       PARAMETER(XZER=0.D0,UNDEMI=.5D0,UN=1.D0)
  18 : -INC CCOPTIO
  19 : -INC SMCHAML
  20 : C
  21 : C
  22 : C
  23 : C
  24 : C           IPMODL          MODELE MMODEL
  25 : C           IPCHE           CHAMELEM DE SOUS TYPE SCALAIRE
  26 : C                           CONTENANT LES VECTEURS SURFACIQUES
  27 : C
  28 : C
  29 :       CHARACTER*4 MCLE(1)  
  30 :       DATA NCLE /1/
  31 :       DATA MCLE /'NORM'/
  32 :       IPMODL=0
  33 :       IPCHE=0
  34 :       IRET=0
  35 : C
  36 : C     LECTURE DU MODELE
  37 : C
  37 :       CALL LIROBJ('MMODEL',IPMODL,1,IRTM)
  38 :       IF (IERR.NE.0) RETURN
  39 : C
  40 : C     LECTURE DU MOT
  41 : C
  41 :       CALL LIRMOT(MCLE,NCLE,ICLE,0)
  42 :       IF (IERR.NE.0) RETURN
  43 :       IF (ICLE.NE.1) THEN
  44 :          INORM = 0
  45 :       ELSE
  46 :          INORM = 1
  47 :       ENDIF
  48 : C
  49 : C     CALCUL DU JACOBIEN
  50 : C
  50 :       CALL JACONO(IPMODL,INORM,IPCHE,IRET)
  51 : C
  51 :       IF(IRET.NE.0) CALL ECROBJ('MCHAML',IPCHE)
  52 :       RETURN
  53 :       END
  54 :  
  55 :  
  56 :  
  57 :