Télécharger la plaquette Cast3M


fr en

Télécharger prija1.eso

Retour à la liste

Numérotation des lignes : 

prija1

  1. C PRIJA1    SOURCE    CB215821  20/11/25    13:36:49     10792          
  2.       SUBROUTINE PRIJA1(MELEMC,MLMCON,MLMPRI,IRN,IVN,IPN,IGAMN,IJACO)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  PRIJA1
  9. C
  10. C DESCRIPTION       :  Voir PRIJAC
  11. C                      Cas deux dimensions, gaz "calorically perfect"
  12. C
  13. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  14. C
  15. C AUTEUR            :  A. BECCANTINI, DM2S/SFME/LTMF
  16. C
  17. C************************************************************************
  18. C
  19. C ENTREES
  20. C
  21. C     MELEMC   : SPG des CHPOINTs
  22. C
  23. C     MLMPRI : liste des inconnues primitives (P,v,T)
  24. C
  25. C     MLMCON : liste des inconnues conservatives (rho, rhov, rhoet)
  26. C
  27. C     IRN    : CHPOINT CENTRE  contenant la densité
  28. C
  29. C     IVN    : CHPOINT CENTRE  contenant la vitesse
  30. C
  31. C     IPN    : CHPOINT CENTRE  contenant la pression
  32. C
  33. C     IGAM   : CHPOINT CENTRE  contenant le gamma
  34. C
  35. C SORTIES
  36. C
  37. C     IJACO  : pointeur de la MATRIK jacobienne
  38. C
  39. C************************************************************************
  40. C
  41. C HISTORIQUE (Anomalies et modifications éventuelles)
  42. C
  43. C HISTORIQUE : crée le 29.05.02
  44. C
  45. C************************************************************************
  46. C
  47. C
  48. C N.B.: On suppose qu'on a déjà controllé RO, P > 0
  49. C                                         GAMMA \in (1,3)
  50. C       Si non il faut le faire!!!
  51. C
  52. C************************************************************************
  53. C
  54. C
  55. C**** Variables de COOPTIO
  56. C
  57. C      INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  58. C     &        ,IOTER,   IOLEC,  IOIMP,     IOCAR,  IOACQ
  59. C     &        ,IOPER,   IOSGB,  IOGRA,     IOSAU,  IORES
  60. C     &        ,IECHO,   IIMPI,  IOSPI
  61. C     &        ,IDIM
  62. C     &        ,MCOORD
  63. C     &        ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  64. C     &        ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  65. C     &        ,NORINC,NORVAL,NORIND,NORVAD
  66. C     &        ,NUCROU, IPSAUV, IFICLE, IPREFI
  67. C
  68.       IMPLICIT INTEGER(I-N)
  69.       INTEGER  IRN, IVN, IPN, IGAMN, IJACO, IGEOM, ICEN
  70.      &       , MP, NP,  NBEL, NBME, NBSOUS
  71.      &       , NRIGE, NMATRI, NKID, NKMT
  72. C
  73.       REAL*8 RHO,UX, UY, GAMMA
  74.       CHARACTER*8 TYPE
  75. C
  76. C**** LES INCLUDES
  77. C
  78.  
  79. -INC PPARAM
  80. -INC CCOPTIO
  81. -INC SMCHPOI
  83. -INC SMELEME
  84. -INC SMLMOTS
  85. -INC SMLENTI
  86.       POINTEUR MPVN.MPOVAL, MPRN.MPOVAL, MPGAM.MPOVAL
  87.       POINTEUR MELEMC.MELEME
  88.       POINTEUR RR.IZAFM, RUX.IZAFM, RUY.IZAFM, RP.IZAFM,
  89.      &         RUXR.IZAFM, RUXUX.IZAFM, RUXUY.IZAFM, RUXP.IZAFM,
  90.      &         RUYR.IZAFM, RUYUX.IZAFM, RUYUY.IZAFM, RUYP.IZAFM,
  91.      &         RETR.IZAFM, RETUX.IZAFM, RETUY.IZAFM, RETP.IZAFM
  92.       POINTEUR MLMPRI.MLMOTS, MLMCON.MLMOTS
  93. C
  94.       CALL LICHT(IVN,MPVN,TYPE,IGEOM)
  95.       CALL LICHT(IRN,MPRN,TYPE,IGEOM)
  96.       CALL LICHT(IGAMN,MPGAM,TYPE,IGEOM)
  97. C
  98. C     SEGACT MPVN*MOD
  99. C     SEGACT MPPN*MOD
  100. C     SEGACT MPRN*MOD
  101. C
  102. C**** Maillage des inconnues primales = Maillage des inconnues primales
  103. C     = MELEMC
  104. C
  105. C
  106.       NRIGE = 7
  107.       NMATRI = 1
  108.       NKID = 9
  109.       NKMT = 7
  110. C
  111.       SEGINI MATRIK
  112.       IJACO = MATRIK
  113.       MATRIK.IRIGEL(1,1) = MELEMC
  114.       MATRIK.IRIGEL(2,1) = MELEMC
  115. C
  116. C**** Matrice non symetrique
  117. C
  118.       MATRIK.IRIGEL(7,1) = 2
  119. C
  120.       NBME = 16
  121.       NBSOUS = 1
  122.       SEGINI IMATRI
  123.       MATRIK.IRIGEL(4,1) = IMATRI
  124. C
  125.       SEGACT MLMPRI
  126.       IMATRI.LISPRI(1) = MLMPRI.MOTS(1)
  127.       IMATRI.LISPRI(2) = MLMPRI.MOTS(2)
  128.       IMATRI.LISPRI(3) = MLMPRI.MOTS(3)
  129.       IMATRI.LISPRI(4) = MLMPRI.MOTS(4)
  130.       IMATRI.LISPRI(5) = MLMPRI.MOTS(1)
  131.       IMATRI.LISPRI(6) = MLMPRI.MOTS(2)
  132.       IMATRI.LISPRI(7) = MLMPRI.MOTS(3)
  133.       IMATRI.LISPRI(8) = MLMPRI.MOTS(4)
  134.       IMATRI.LISPRI(9) = MLMPRI.MOTS(1)
  135.       IMATRI.LISPRI(10) = MLMPRI.MOTS(2)
  136.       IMATRI.LISPRI(11) = MLMPRI.MOTS(3)
  137.       IMATRI.LISPRI(12) = MLMPRI.MOTS(4)
  138.       IMATRI.LISPRI(13) = MLMPRI.MOTS(1)
  139.       IMATRI.LISPRI(14) = MLMPRI.MOTS(2)
  140.       IMATRI.LISPRI(15) = MLMPRI.MOTS(3)
  141.       IMATRI.LISPRI(16) = MLMPRI.MOTS(4)
  142.       SEGDES MLMPRI
  143. C
  144.       SEGACT MLMCON
  145.       IMATRI.LISDUA(1) = MLMCON.MOTS(1)
  146.       IMATRI.LISDUA(2) = MLMCON.MOTS(1)
  147.       IMATRI.LISDUA(3) = MLMCON.MOTS(1)
  148.       IMATRI.LISDUA(4) = MLMCON.MOTS(1)
  149.       IMATRI.LISDUA(5) = MLMCON.MOTS(2)
  150.       IMATRI.LISDUA(6) = MLMCON.MOTS(2)
  151.       IMATRI.LISDUA(7) = MLMCON.MOTS(2)
  152.       IMATRI.LISDUA(8) = MLMCON.MOTS(2)
  153.       IMATRI.LISDUA(9) = MLMCON.MOTS(3)
  154.       IMATRI.LISDUA(10) = MLMCON.MOTS(3)
  155.       IMATRI.LISDUA(11) = MLMCON.MOTS(3)
  156.       IMATRI.LISDUA(12) = MLMCON.MOTS(3)
  157.       IMATRI.LISDUA(13) = MLMCON.MOTS(4)
  158.       IMATRI.LISDUA(14) = MLMCON.MOTS(4)
  159.       IMATRI.LISDUA(15) = MLMCON.MOTS(4)
  160.       IMATRI.LISDUA(16) = MLMCON.MOTS(4)
  161.       SEGDES MLMCON
  162. C
  163.       SEGACT MELEMC
  164.       NBEL = MELEMC.NUM(/2)
  165.       SEGDES MELEMC
  166.       NBSOUS = 1
  167.       NP = 1
  168.       MP = 1
  169.       SEGINI   RR  , RUX , RUY , RP ,
  170.      &         RUXR , RUXUX , RUXUY , RUXP ,
  171.      &         RUYR , RUYUX , RUYUY , RUYP ,
  172.      &         RETR , RETUX , RETUY , RETP
  173. C
  174. C**** Duale = IMATRI.LISDUA(1) = 'RN'
  175. C     Primale = IMATRI.LISPRI(4) = 'P'
  176. C     -> IMATRI.LIZAFM(1,4) = RP
  177. C
  178.       IMATRI.LIZAFM(1,1) = RR
  179.       IMATRI.LIZAFM(1,2) = RUX
  180.       IMATRI.LIZAFM(1,3) = RUY
  181.       IMATRI.LIZAFM(1,4) = RP
  182.       IMATRI.LIZAFM(1,5) = RUXR
  183.       IMATRI.LIZAFM(1,6) = RUXUX
  184.       IMATRI.LIZAFM(1,7) = RUXUY
  185.       IMATRI.LIZAFM(1,8) = RUXP
  186.       IMATRI.LIZAFM(1,9) = RUYR
  187.       IMATRI.LIZAFM(1,10) = RUYUX
  188.       IMATRI.LIZAFM(1,11) = RUYUY
  189.       IMATRI.LIZAFM(1,12) = RUYP
  190.       IMATRI.LIZAFM(1,13) = RETR
  191.       IMATRI.LIZAFM(1,14) = RETUX
  192.       IMATRI.LIZAFM(1,15) = RETUY
  193.       IMATRI.LIZAFM(1,16) = RETP
  194. C
  195.       DO ICEN = 1, NBEL, 1
  196.          RHO=MPRN.VPOCHA(ICEN,1)
  197.          UX=MPVN.VPOCHA(ICEN,1)
  198.          UY=MPVN.VPOCHA(ICEN,2)
  199.          GAMMA=MPGAM.VPOCHA(ICEN,1)
  200. C
  201. C********AB.AM(IFAC,IPRIM,IDUAL)
  202. C        A = nom de l'inconnu duale (Ro,rUX,rUY,RET)
  203. C        B = nom de l'inconnu primale (Ro,UX,UY,P)
  204. C
  205. C******* Dual RN
  206. C
  207.          RR.AM(ICEN,1,1)  = 1.0D0
  208.          RUX.AM(ICEN,1,1) = 0.0D0
  209.          RUY.AM(ICEN,1,1) = 0.0D0
  210.          RP.AM(ICEN,1,1)  = 0.0D0
  211. C
  212. C******* Dual RUXN
  213. C
  214.          RUXR.AM(ICEN,1,1) = UX
  215.          RUXUX.AM(ICEN,1,1) = RHO
  216.          RUXUY.AM(ICEN,1,1) = 0.0D0
  217.          RUXP.AM(ICEN,1,1) = 0.0D0
  218. C
  219. C******* Dual RUYN
  220. C
  221.          RUYR.AM(ICEN,1,1) = UY
  222.          RUYUX.AM(ICEN,1,1) = 0.0D0
  223.          RUYUY.AM(ICEN,1,1) = RHO
  224.          RUYP.AM(ICEN,1,1) = 0.0D0
  225. C
  226. C********Dual RETN
  227. C
  228.          RETR.AM(ICEN,1,1) = 0.5D0*((UX*UX)+(UY*UY))
  229.          RETUX.AM(ICEN,1,1) = RHO*UX
  230.          RETUY.AM(ICEN,1,1) = RHO*UY
  231.          RETP.AM(ICEN,1,1) = 1.0D0/(GAMMA - 1.0D0)
  232. C
  233.       ENDDO
  234. C
  235.       SEGDES MPVN
  236.       SEGDES MPRN
  237.       SEGDES MPGAM
  238. C
  239.       SEGDES MATRIK
  240.       SEGDES IMATRI
  241. C
  242.       SEGDES   RR  , RUX , RUY , RP ,
  243.      &         RUXR , RUXUX , RUXUY , RUXP ,
  244.      &         RUYR , RUYUX , RUYUY , RUYP ,
  245.      &         RETR , RETUX , RETUY , RETP
  246. C
  247.       RETURN
  248.       END
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  

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