calc_tempinter_kdif.f90
上传用户:goto8899
上传日期:2014-01-20
资源大小:303k
文件大小:1k
源码类别:

并行计算

开发平台:

MultiPlatform

  1. !------------------------------------------------------------------------------!
  2. ! Procedure :  calc_tempinter_kdif        Auteur : E. Radenac
  3. !                                         Date   : Avril 2003
  4. ! Fonction                                Modif  : Juin 2003
  5. !   Calcul de la temperature de l'interface entre deux cellules appartenant a 
  6. !   deux zones differentes
  7. !
  8. ! Defauts/Limitations/Divers :
  9. !
  10. !------------------------------------------------------------------------------!
  11. subroutine calc_tempinter_kdif(temp1, temp2, conduct1, conduct2, d1, d2, &
  12. temp_inter)
  13. use TYPHMAKE
  14. use OUTPUT
  15. implicit none
  16. ! -- Declaration des entrees --
  17. integer               :: if
  18. real(krp)             :: temp1, temp2   ! temperatures echangees
  19. real(krp)             :: conduct1, conduct2  ! conductivites echangees
  20. real(krp)             :: d1, d2  ! distance entre les centres des cellules et l'interface          
  21.         
  22. ! -- Declaration des entrees/sorties --
  23. real(krp) :: temp_inter
  24. ! -- Declaration des variables internes --
  25. real(krp) :: a, b
  26. ! -- Debut de la procedure --
  27. a = conduct1/d1
  28. b = conduct2/d2
  29. temp_inter = (a*temp1+b*temp2)/(a+b)
  30. !print*, "temperatures CALC_TEMPINTER_KDIF :", temp1, " et ", temp2
  31. !print*, "temperature a l'interface :", temp_inter
  32. endsubroutine calc_tempinter_kdif