init_champ_ust.f90

来自「国外大名顶顶的“台风”并行计算流体力学CFD软件的早期版本的源代码」· F90 代码 · 共 57 行

!------------------------------------------------------------------------------!
! Procedure : init_champ_ust              Auteur : J. Gressier
!                                         Date   : Mars 2003
! Fonction                                Modif  : see history
!   Initialization of fields according to the solver
!
! Defauts/Limitations/Divers :
!
!------------------------------------------------------------------------------!
!subroutine init_champ_ust(defsolver, ust_mesh, champ, grid)
subroutine init_champ_ust(defsolver, ust_mesh, grid)

use TYPHMAKE
use VARCOM
use OUTPUT
use USTMESH
use DEFFIELD
use MENU_SOLVER

implicit none

! -- Declaration des entrees --
type(mnu_solver) :: defsolver            ! parametres du solveur
type(st_ustmesh) :: ust_mesh             ! maillage et connectivites
type(st_grid)    :: grid                 ! grille

! -- Declaration des sorties --
!type(st_field), pointer :: champ                ! champ d'etat et de gradients

! -- Declaration des variables internes --
integer :: i
type(st_field), pointer :: champ

! -- Debut de la procedure --

call print_info(8, ". initialisation et allocation des champs")

! allocation des champs

select case(defsolver%typ_solver)
case(solNS)
  champ=>newfield(grid, 2, 1, ust_mesh%ncell, ust_mesh%nface)
case(solKDIF)
  champ=>newfield(grid, 1, 0, ust_mesh%ncell, ust_mesh%nface)  ! A MODIFIER SELON GAZ
case(solVORTEX)
  champ=>newfield(grid, 1, 0, ust_mesh%ncell, ust_mesh%nface)
case default
  call erreur("Incoherence interne (init_champ_ust)","type de solveur inconnu")
endselect 

call alloc_prim(champ)

! Boucle sur les definitions de champ

do i = 1, defsolver%ninit

  write(str_w,'(a,i3)') "    initialisation n