Télécharger pre4.eso

Retour à la liste

Numérotation des lignes :

pre4
  1. C PRE4 SOURCE OF166741 24/12/13 21:17:08 12097
  2. SUBROUTINE PRE4()
  3. C************************************************************************
  4. C
  5. C PROJET : CASTEM 2000
  6. C
  7. C NOM : PRE4
  8. C
  9. C DESCRIPTION : Transport de scalaires passifs (vitesse aux faces)
  10. C
  11. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  12. C
  13. C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
  14. C
  15. C************************************************************************
  16. C
  17. C
  18. C APPELES (Outils) : LIRENT, ERREUR
  19. C
  20. C APPELES (Calcul) :
  21. C
  22. C PRE4-------------- PRE41 (1er ordre en espace)
  23. C
  24. C************************************************************************
  25. C
  26. C PHRASE D'APPEL (GIBIANE) :
  27. C
  28. C a) 1er ordre en espace (1er ordre en temps)
  29. C
  30. C CHAM1S = 'PRET' 'CLAUDEIS' 'FACE' ENTI1 TAB1 CHPO1 ;
  31. C
  32. C
  33. C ENTREES :
  34. C
  35. C ENTI1 : ordre en espace (=1)
  36. C
  37. C MOD1 : OBJET MODELE de TYPE NAVIER_STOKES
  38. C
  39. C CHPO1 : CHPOINT "CENTRE" contenant les scalaires à tranporter
  40. C
  41. C SORTIES :
  42. C
  43. C CHAM1S : MCHAML "FACEL" contenant les s scalaires à tranporter
  44. C (a gauche et a droite de chaque face).
  45. C Autant de composantes que CHPO1
  46. C
  47. C
  48. C************************************************************************
  49. C
  50. C HISTORIQUE (Anomalies et modifications éventuelles)
  51. C
  52. C HISTORIQUE : Créée le 28.11.01
  53. C
  54. C************************************************************************
  55. C
  56. C**** Les variables
  57. C
  58. IMPLICIT INTEGER(I-N)
  59. INTEGER ICOND, IRETOU
  60. & ,ORDESP,NBOPT, NBOPT0,IP
  61. PARAMETER(NBOPT=1)
  62. CHARACTER*(8) TOPTIO(NBOPT),MOTOP
  63. C
  64. DATA TOPTIO/'FACE'/
  65. C
  66. C**** Les Includes
  67. C
  68.  
  69. -INC PPARAM
  70. -INC CCOPTIO
  71. C
  72. C**** Les mots clés
  73. C
  74. ICOND = 1
  75. CALL LIRCHA(MOTOP,ICOND,IRETOU)
  76. IF(IERR .NE. 0) GOTO 9999
  77. NBOPT0 = NBOPT
  78. CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT0)
  79. IF(IERR .NE. 0)GOTO 9999
  80. IF(IP.NE.1)THEN
  81. *
  82. **** Tentative d'utilisation d'une option non-implementée
  83. C
  84. CALL ERREUR(251)
  85. GOTO 9999
  86. ENDIF
  87. C
  88. C**** Initialisation des parametres d'erreur
  89. C
  90. MOTERR(1:40) = ' '
  91. C
  92. C**** Lecture ORDESP
  93. C
  94. ICOND = 1
  95. CALL LIRENT(ORDESP,ICOND,IRETOU)
  96. IF(IERR .NE. 0) GOTO 9999
  97. IF((ORDESP .LT. 1) .OR. (ORDESP .GT. 2))THEN
  98. C
  99. C******* Message d'erreur standard
  100. C -301 0 %m1:40
  101. C
  102. MOTERR(1:40) = 'ENTI1 = ??? '
  103. WRITE(IOIMP,*) MOTERR
  104. C
  105. C******* Message d'erreur standard
  106. C Entier valant: %i1
  107. C -2 0
  108. C
  109. INTERR(1) = ORDESP
  110. CALL ERREUR(-2)
  111. C
  112. C******* Message d'erreur standard
  113. C 21 2
  114. C Données incompatibles
  115. C
  116. CALL ERREUR(21)
  117. GOTO 9999
  118. ENDIF
  119. IF(ORDESP .EQ. 1)THEN
  120. C
  121. C******* ORDESP = 1
  122. C
  123. CALL PRE41()
  124. ELSE
  125. C
  126. C******* Message d'erreur standard
  127. C 5 3
  128. C Erreur anormale. Contactez votre support
  129. C
  130. CALL ERREUR(5)
  131. ENDIF
  132. C
  133. 9999 CONTINUE
  134. C
  135. RETURN
  136. END
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  

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