Télécharger vsur.eso

Retour à la liste

Numérotation des lignes :

  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
  38. CALL LIROBJ('MMODEL',IPMODL,1,IRTM)
  39. IF (IERR.NE.0) RETURN
  40. C
  41. C LECTURE DU MOT
  42. C
  43. CALL LIRMOT(MCLE,NCLE,ICLE,0)
  44. IF (IERR.NE.0) RETURN
  45. IF (ICLE.NE.1) THEN
  46. INORM = 0
  47. ELSE
  48. INORM = 1
  49. ENDIF
  50. C
  51. C CALCUL DU JACOBIEN
  52. C
  53. CALL JACONO(IPMODL,INORM,IPCHE,IRET)
  54. C
  55. IF(IRET.NE.0) CALL ECROBJ('MCHAML',IPCHE)
  56. RETURN
  57. END
  58.  
  59.  
  60.  
  61.  
  62.  

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