Télécharger la plaquette Cast3M


fr en

Télécharger modcal.eso

Retour à la liste

Numérotation des lignes : 

modcal

  1. C MODCAL    SOURCE    CHAT      05/01/13    01:47:31     5004
  2.       SUBROUTINE MODCAL(EPSILO,NDIM,ITYG,XINTER,XNORMA
  3.      $     ,XIREB,XNREB,IREBCO,IREBDI,XDEP2,TDEP,UREEL,U
  4.      $     ,IZSH,DTREEL,XCON2,XARI2,TARI,DTNEW,ITEST1,ITEST2,LTEST)
  5. ***************************************************************************
  6. *** SP 'MODCAL' : traite les cas ou calibrage particule via 'CALIDT' est
  7. *** inefficace (cas ou pt calibré invariant et hors maille). 'MODCAL'
  8. *** recalibre le saut convectif seul sans tenir compte du saut dspersif.
  9. ***
  10. *** APPELES 1 = aucun
  11. *** APPELES 2 = 'CALIDT', 'LIEUPT'
  12. ***
  13. *** E = 'EPSILO' marge relative acceptée position % element
  14. ***     'NDIM' dimension de l'espace
  15. ***     'ITYG' entier caracterisant la geometrie de l'element
  16. ***     'XINTER' coord reelles du pt intersection du calibrage
  17. ***     'XNORMA' vecteur unitaire normal du calibrage
  18. ***     'XIREB' pt d'impact sur la face impermeable
  19. ***     'XNREB' vecteur normal à la face impermeable
  20. ***     'IREBCO' vaut 1 si rebond lors de la convection, 0 sinon
  21. ***     'IREBDI' vaut 1 si rebond lors de la diffusion, 0 sinon
  22. ***     'XDEP2' coord reelles depart particule avt calibrage
  23. ***     'TDEP' tps reel depart particule avt calibrage
  24. ***     'UREEL' vecteur "vitesse convective" du calibrage
  25. ***     'U' vecteur "vitesse diffusive" du calibrage
  26. ***     'IZSH' segment content coord reelles noeuds, fcts forme et base
  27. ***     'DTREEL' pas de tps réel de calcul avancée particule
  28. ***
  29. *** E/S = 'XCON2' nouvelles coord reelles apres convection
  30. ***       'XARI2' nouvelles coord reelles arrivee particule
  31. ***       'TARI' nouveau tps reel d'arrivee particule
  32. ***       'DTNEW' pas de tps avancee particule au cours traitement
  33. ***       'ITEST1' vaut 1 si calib sur conv seule donne dt coherent, 0 sinon
  34. ***       'ITEST2' vaut 1 si nveau pt arrivee appartient à maille, 0 sinon
  35. ***       'LTEST' vaut 1 si 1er essai calibrage, 0 si 2eme essai
  36. ***
  37. *** ORIGINE = CYRIL NOU
  38. ***************************************************************************
  39.  
  40.  
  41.       IMPLICIT INTEGER(I-N)
  42.       IMPLICIT REAL*8 (A-H,O-Z)
  43.       SEGMENT IZSH
  44.           REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  45.       ENDSEGMENT
  46.       DIMENSION XINTER(3),XNORMA(3),XDEP2(3),XCON2(3),XARI2(3)
  47.       DIMENSION UZERO(3),UREEL(3),U(3),XIREB(3),XNREB(3),xcon(3)
  48. ***   on cree un vecteur "vitesse diffusive" nul
  49.       DO 10 I=1,NDIM
  50.          UZERO(I)=0.D0
  51.  10   CONTINUE
  52. ***   traitement consistant a recalibrer sur convec seule (m vecteur vitesse)
  53.       CALL CALIDT(EPSILO,NDIM,XINTER,XNORMA
  54.      $     ,XIREB,XNREB,IREBCO,IREBDI,XDEP2,TDEP,UREEL
  55.      $     ,UZERO,DTREEL,LTEST,DTNEW,XCON2,XARI2,TARI,ITEST1)
  56. ***   saut diffusif avec m vecteur vitesse diff si pas d'echec calibrage
  57.       DO 20 I=1,NDIM
  58.          XARI2(I)=XCON2(I)+U(I)*SQRT(DTNEW)
  59.  20   CONTINUE
  60. ***   localisation nouveau pt arrivee apres traitement à 0 pres
  61.       CALL LIEUPT(0.D0,NDIM,ITYG,XARI2,IZSH,ITEST2)
  62.  
  63.       RETURN
  64.       END
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  
  71.  

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