init_kdif_ust.f90

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

!------------------------------------------------------------------------------!
! Procedure : init_kdif_ust               Auteur : J. Gressier
!                                         Date   : Mars 2003
! Fonction                                Modif  : juin 2003 (cf historique)
!   Traitement des parametres du fichier menu principal
!   Parametres principaux du projet
!
! Defauts/Limitations/Divers :
!   ATTENTION : initialisation des variables primitives
!
!------------------------------------------------------------------------------!
subroutine init_kdif_ust(kdif, champ, unif, mesh)

use TYPHMAKE
use DEFFIELD
use MENU_KDIF
use MENU_INIT

implicit none

! -- Declaration des entrees --
type(st_init_kdif) :: kdif
integer            :: unif ! uniformite de la condition initiale
type(st_mesh)      :: mesh

! -- Declaration des sorties --
type(st_field) :: champ

! -- Declaration des variables internes --
integer :: ip, ic

! -- Debut de la procedure --

if (unif == init_unif) then
  do ip = 1, champ%nscal
    champ%etatprim%tabscal(ip)%scal(:) = kdif%temp
  enddo
else !provisoire
  do ip = 1, champ%nscal
    do ic=1, champ%ncell
      champ%etatprim%tabscal(ip)%scal(ic)=kdif%coef(1)*mesh%centre(ic,1,1)%x+&
                                          kdif%coef(2)*mesh%centre(ic,1,1)%y+&
                                          kdif%coef(3)*mesh%centre(ic,1,1)%z+&
                                          kdif%coef(4)
    enddo
  enddo
endif

! pas de de variables vectorielles attendues (pas de test)

!!if (champ%allocgrad) champ%gradient(:,:,:,:,:) = 0._krp


endsubroutine init_kdif_ust

!------------------------------------------------------------------------------!
! Historique des modifications
!
! mars 2003 (v0.0.1b) : creation de la routine
! juin 2003           : maj pour variables conservatives et primitives
!------------------------------------------------------------------------------!