Télécharger la plaquette Cast3M


fr en

Télécharger saucon.eso

Retour à la liste

Numérotation des lignes : 

saucon

  1. C SAUCON    SOURCE    CHAT      05/01/13    03:11:11     5004
  2.       SUBROUTINE SAUCON(NDIM,ITY1,ITYG,NOEL1,NOUN
  3.      $     ,DIAM,UELEM,XDEP,XDEP2,DTREEL,XARI2,IZSH,UREEL,ITEST)
  4.  
  5. *************************************************************************
  6. *** SP 'SAUCON' : fait avancer la particule apres un pas de tps par
  7. *** convection explicite uniquement dans l'element reel considere
  8. ***
  9. *** APPELES 1 = aucun
  10. *** APPELES 2 = 'VITAPP'
  11. ***
  12. *** E=  'NDIM' dimension de l'espace
  13. ***     'ITY1' entier caracterisant le type de l'element
  14. ***     'ITYG' entier caracterisant la geometrie de l'element
  15. ***     'NOEL1' nombre de noeuds de l'element considere
  16. ***     'NOUN' nombre de flux consideres
  17. ***     'DIAM' "longueur caracteristique" de l'element
  18. ***     'UELEM' valeurs des flux aux faces
  19. ***     'XDEP' coordonnees reference du pt de depart
  20. ***     'XDEP2' coordonnees reelles du pt de depart
  21. ***     'DTREEL' pas de tps reel pour avancee de la particule
  22. ***
  23. *** S = 'XARI2' coordonnees reelles d'arrivee
  24. ***     'ITEST' vaut 0 si jacobien nul (ds approximation vitesse), 1 sinon
  25. ***     'UREEL' vitesse convective reelle
  26. ***
  27. *** E/S = 'IZSH' segment content pour l'elemt considere et au pt depart :
  28. ***             'XYZL' coordonnees reelles des noeuds (E)
  29. ***             'SHP' fonctions de forme et derivees ds elemt ref (S)
  30. ***             'SHY' fonctions de base et derivees ds elemt ref (S)
  31. ***
  32. ***          auteur Cyril Nou
  33. *************************************************************************
  34.  
  35.       IMPLICIT INTEGER(I-N)
  36.       IMPLICIT REAL*8 (A-H,O-Z)
  37. *-INC SMELEME
  38. *      POINTEUR IZCENT.MELEME,IELTFA.MELEME
  39.       SEGMENT IZSH
  40.           REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  41.       ENDSEGMENT
  42.       DIMENSION XDEP(3),XDEP2(3),UELEM(NOUN)
  43.       DIMENSION XARI2(3),UREF(3),UREEL(3)
  44.  
  45. ***   calcul de la vitesse reelle au pt de depart
  46.       CALL VITAPP(NDIM,ITY1,ITYG,NOEL1
  47.      $     ,NOUN,DIAM,UELEM,XDEP,UREF,UREEL,IZSH,ITEST)
  48. ***   calcul pt reel d'arrivee apres pas de tps 'DTREEL'
  49.       DO 10 I=1,NDIM
  50.          XARI2(I)=XDEP2(I)+UREEL(I)*DTREEL
  51.  10   CONTINUE
  52.       RETURN
  53.       END
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  

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