Télécharger hyva3d.eso

Retour à la liste

Numérotation des lignes :

hyva3d
  1. C HYVA3D SOURCE PV090527 23/01/27 21:15:38 11574
  2. c***********************************************************************
  3. subroutine hydrvari3d(vari0,vari1,hydra0,hydra1,hydras)
  4. c effets de l hydratation sur les varibles internes
  5. c declaration externes
  6. real*8 vari0,vari1,hydra0,hydra1,hydras
  7. integer erreur
  8. c varibles locales
  9. real*8 yy,yymin
  10. parameter (yymin=1.d-3)
  11. c print*,'hdratation actuelle:',hydra1
  12. c print*,'hydratation precedente:',hydra0
  13. c print*,'hydratation seuil:',hydras
  14.  
  15. c initialisation controle d erreur
  16. if (hydra1.lt. yymin) then
  17. c le materiau est non coherent les varibles internes
  18. c sont mises a zero a chaque debut de pas et chaque fin de pas
  19. c pour l affichage
  20. vari1=0.d0
  21. else
  22. hydra1=min(hydra1,1.d0)
  23. if ((hydra1.ge.hydra0).and.(hydra0.gt.hydras)) then
  24. c apres le seuil on conserve les vari ou on les attenue
  25. if(hydra1.eq.hydra0) then
  26. vari1=vari0
  27. else
  28. vari1=vari0*(hydra0-hydras)/(hydra1-hydras)
  29. end if
  30. else
  31. if (hydra1.le.hydras) then
  32. c avant le seuil on met les vari a zero
  33. vari1=0.d0
  34. else
  35. if((hydra1.gt.hydras).and.(hydra0.le.hydras)) then
  36. c la fin de pas a passé le seuil mais pas le debut de pas
  37. vari1=0.d0
  38. else if (hydra1.lt.hydra0) then
  39. c l hydratation diminue on conserve la vari
  40. vari1=vari0
  41. else
  42. print*,'Cas imprevu dans hyva3d de fldo3d'
  43. vari1=vari0
  44. end if
  45. end if
  46. end if
  47. end if
  48. return
  49. end
  50. c***********************************************************************
  51.  
  52.  
  53.  

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