echange.f90

国外大名顶顶的“台风”并行计算流体力学CFD软件的早期版本的源代码

!------------------------------------------------------------------------------!
! Procedure : echange            	  Auteur : E. Radenac
!                                         Date   : Mai 2003
! Fonction                                Modif  : Juin 2003
!   Echange de donnees entre zones de calcul
!
! Defauts/Limitations/Divers :
!
!------------------------------------------------------------------------------!
subroutine echange(echdata1, echdata2, normale, vecinter, d1, d2, nfacelim, &
			typecalcul, typemethode, solvercoupling, boco1, &
                        boco2, connface2)

use TYPHMAKE
use OUTPUT
use GEO3D
use DEFFIELD
use VARCOM
use MENU_BOCO

implicit none

! -- Declaration des entrees --
type(st_genericfield)      :: echdata1, echdata2
integer                    :: nfacelim   ! nombre de faces limites sur l'interface
type(v3d), dimension(nfacelim) &
                           :: vecinter                ! vecteur unitaire "intercellules"                       
type(v3d), dimension(nfacelim) &
                           :: normale ! normales a l'interface

real(krp), dimension(nfacelim) &
		           :: d1, d2  ! distance entre les centres des cellules gauche,
		      		      ! droite et l'interface
integer                    :: typecalcul, typemethode
integer                    :: solvercoupling
integer, dimension(nfacelim) &
                           :: connface2

! -- Declaration des entrees/sorties --2
type(mnu_boco)             :: boco1, boco2

! -- Declaration des variables internes --

! -- Debut de la procedure --

select case(solvercoupling)

case(kdif_kdif)
call echange_kdif(echdata1, echdata2, normale, vecinter, d1, d2, nfacelim, &
 			typecalcul, typemethode, boco1%boco_kdif, &
                        boco2%boco_kdif, connface2)
case(kdif_ns)
call erreur("incoherence interne (echange)", "non implemente")

case(ns_ns)
call erreur("incoherence interne (echange)", "non implemente")

case default
call erreur("incoherence interne (echange)", "couplage de solvers inconnu")

endselect


! "Determination" du pas de temps d'echange (minimum, maximum, senseurs) avant l'echange suivant
! ??
!call calc_tps_echange()


endsubroutine echange