!------------------------------------------------------------------------------- ! Code_Saturne version 2.0.0-rc1 ! -------------------------- ! This file is part of the Code_Saturne Kernel, element of the ! Code_Saturne CFD tool. ! Copyright (C) 1998-2009 EDF S.A., France ! contact: saturne-support@edf.fr ! The Code_Saturne Kernel is free software; you can redistribute it ! and/or modify it under the terms of the GNU General Public License ! as published by the Free Software Foundation; either version 2 of ! the License, or (at your option) any later version. ! The Code_Saturne Kernel is distributed in the hope that it will be ! useful, but WITHOUT ANY WARRANTY; without even the implied warranty ! of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the ! GNU General Public License for more details. ! You should have received a copy of the GNU General Public License ! along with the Code_Saturne Kernel; if not, write to the ! Free Software Foundation, Inc., ! 51 Franklin St, Fifth Floor, ! Boston, MA 02110-1301 USA !------------------------------------------------------------------------------- subroutine uscfxi & !================ ( idbia0 , idbra0 , & ndim , ncelet , ncel , nfac , nfabor , nfml , nprfml , & nnod , lndfac , lndfbr , ncelbr , & nvar , nscal , nphas , & nideve , nrdeve , nituse , nrtuse , & ifacel , ifabor , ifmfbr , ifmcel , iprfml , maxelt , lstelt , & ipnfac , nodfac , ipnfbr , nodfbr , & idevel , ituser , ia , & xyzcen , surfac , surfbo , cdgfac , cdgfbo , xyznod , volume , & dt , rtp , propce , propfa , propfb , coefa , coefb , & w1 , w2 , w3 , w4 , & rdevel , rtuser , ra ) !=============================================================================== ! Purpose: ! ------- ! User subroutine. ! Initialize the unknown variables for the compressible flow scheme. ! Description ! =========== ! This subroutine is similar to the user subroutine 'usiniv', but ! is dedicated to the compressible flow scheme. ! It is called at the beginning of the computation (only if it is ! not a restart), just before the time marching loop starts. ! It allows to initialize all the unknown variables. ! The standard initialization has been reproduced here as an example. ! More examples can be found in 'usiniv'. ! Physical properties ! =================== ! The physical properties (viscosity, specific heat, thermal ! conductivity, Schmidt number) that are stored in the arrays propce, ! propfa and propfb must not be modified here: if it is necessary to ! do so, it must be done in the dedicated user programme 'uscfpv'. ! Cells identification ! ==================== ! Cells may be identified using the 'getcel' subroutine. ! The syntax of this subroutine is described in the 'usclim' subroutine, ! but a more thorough description can be found in the user guide. ! Arguments !__________________.____._____.________________________________________________. ! nom !type!mode ! role ! !__________________!____!_____!________________________________________________! ! idbia0 ! i ! <-- ! number of first free position in ia ! ! idbra0 ! i ! <-- ! number of first free position in ra ! ! ndim ! i ! <-- ! spatial dimension ! ! ncelet ! i ! <-- ! number of extended (real + ghost) cells ! ! ncel ! i ! <-- ! number of cells ! ! nfac ! i ! <-- ! number of interior faces ! ! nfabor ! i ! <-- ! number of boundary faces ! ! nfml ! i ! <-- ! number of families (group classes) ! ! nprfml ! i ! <-- ! number of properties per family (group class) ! ! nnod ! i ! <-- ! number of vertices ! ! lndfac ! i ! <-- ! size of nodfac indexed array ! ! lndfbr ! i ! <-- ! size of nodfbr indexed array ! ! ncelbr ! i ! <-- ! number of cells with faces on boundary ! ! nvar ! i ! <-- ! total number of variables ! ! nscal ! i ! <-- ! total number of scalars ! ! nphas ! i ! <-- ! number of phases ! ! nideve, nrdeve ! i ! <-- ! sizes of idevel and rdevel arrays ! ! nituse, nrtuse ! i ! <-- ! sizes of ituser and rtuser arrays ! ! ifacel(2, nfac) ! ia ! <-- ! interior faces -> cells connectivity ! ! ifabor(nfabor) ! ia ! <-- ! boundary faces -> cells connectivity ! ! ifmfbr(nfabor) ! ia ! <-- ! boundary face family numbers ! ! ifmcel(ncelet) ! ia ! <-- ! cell family numbers ! ! iprfml ! ia ! <-- ! property numbers per family ! ! (nfml, nprfml) ! ! ! ! ! maxelt ! i ! <-- ! max number of cells and faces (int/boundary) ! ! lstelt(maxelt) ! ia ! --- ! work array ! ! ipnfac(nfac+1) ! ia ! <-- ! interior faces -> vertices index (optional) ! ! nodfac(lndfac) ! ia ! <-- ! interior faces -> vertices list (optional) ! ! ipnfbr(nfabor+1) ! ia ! <-- ! boundary faces -> vertices index (optional) ! ! nodfac(lndfbr) ! ia ! <-- ! boundary faces -> vertices list (optional) ! ! idevel(nideve) ! ia ! <-> ! integer work array for temporary developpement ! ! ituser(nituse ! ia ! <-> ! user-reserved integer work array ! ! ia(*) ! ia ! --- ! main integer work array ! ! xyzcen ! ra ! <-- ! cell centers ! ! (ndim, ncelet) ! ! ! ! ! surfac ! ra ! <-- ! interior faces surface vectors ! ! (ndim, nfac) ! ! ! ! ! surfbo ! ra ! <-- ! boundary faces surface vectors ! ! (ndim, nfavor) ! ! ! ! ! cdgfac ! ra ! <-- ! interior faces centers of gravity ! ! (ndim, nfac) ! ! ! ! ! cdgfbo ! ra ! <-- ! boundary faces centers of gravity ! ! (ndim, nfabor) ! ! ! ! ! xyznod ! ra ! <-- ! vertex coordinates (optional) ! ! (ndim, nnod) ! ! ! ! ! volume(ncelet) ! ra ! <-- ! cell volumes ! ! dt(ncelet) ! ra ! <-- ! time step (per cell) ! ! rtp(ncelet,*) ! ra ! <-> ! calculated variables at cell centers ! ! propce(ncelet, *)! ra ! <-- ! physical properties at cell centers ! ! propfa(nfac, *) ! ra ! <-- ! physical properties at interior face centers ! ! propfb(nfabor, *)! ra ! <-- ! physical properties at boundary face centers ! ! coefa, coefb ! ra ! <-- ! boundary conditions ! ! (nfabor, *) ! ! ! ! ! w1..4(ncelet) ! tr ! --- ! work arrays ! ! rdevel(nrdeve) ! ra ! <-> ! real work array for temporary developpement ! ! rtuser(nituse ! ra ! <-> ! user-reserved real work array ! ! ra(*) ! ra ! --- ! main real work array ! !__________________!____!_____!________________________________________________! ! Type: i (integer), r (real), s (string), a (array), l (logical), ! and composite types (ex: ra real array) ! mode: <-- input, --> output, <-> modifies data, --- work array !=============================================================================== implicit none !=============================================================================== ! Common blocks !=============================================================================== include "paramx.h" include "pointe.h" include "numvar.h" include "optcal.h" include "cstphy.h" include "cstnum.h" include "entsor.h" include "parall.h" include "period.h" include "ppppar.h" include "ppthch.h" include "ppincl.h" !=============================================================================== ! Arguments integer idbia0 , idbra0 integer ndim , ncelet , ncel , nfac , nfabor integer nfml , nprfml integer nnod , lndfac , lndfbr , ncelbr integer nvar , nscal , nphas integer nideve , nrdeve , nituse , nrtuse integer ifacel(2,nfac) , ifabor(nfabor) integer ifmfbr(nfabor) , ifmcel(ncelet) integer iprfml(nfml,nprfml), maxelt, lstelt(maxelt) integer ipnfac(nfac+1), nodfac(lndfac) integer ipnfbr(nfabor+1), nodfbr(lndfbr) integer idevel(nideve), ituser(nituse), ia(*) double precision xyzcen(ndim,ncelet) double precision surfac(ndim,nfac), surfbo(ndim,nfabor) double precision cdgfac(ndim,nfac), cdgfbo(ndim,nfabor) double precision xyznod(ndim,nnod), volume(ncelet) double precision dt(ncelet), rtp(ncelet,*), propce(ncelet,*) double precision propfa(nfac,*), propfb(nfabor,*) double precision coefa(nfabor,*), coefb(nfabor,*) double precision w1(ncelet), w2(ncelet), w3(ncelet), w4(ncelet) double precision rdevel(nrdeve), rtuser(nrtuse), ra(*) ! Local variables integer idebia, idebra integer iel, iphas integer iccfth, iscal, imodif, iutile !=============================================================================== !=============================================================================== ! 1. Control print !=============================================================================== write(nfecra,9001) !=============================================================================== ! 2. Initialization of local variables !=============================================================================== idebia = idbia0 idebra = idbra0 imodif = 1 !=============================================================================== ! 3. Unknown variable initialization ! for initial calculations (not in case of restart) !=============================================================================== if ( isuite.eq.0 ) then iphas = 1 ! --- Velocity components do iel = 1, ncel rtp(iel,iu(iphas)) = 67.d0 rtp(iel,iv(iphas)) = 0.d0 rtp(iel,iw(iphas)) = 0.d0 enddo ! --- User defined scalars ! If there are user defined scalars if(nscaus.gt.0) then ! For each scalar do iscal = 1, nscaus ! If the scalar is associated to the considered phase iphas if(iphsca(iscal).eq.iphas) then ! Initialize each cell value do iel = 1, ncel rtp(iel,isca(iscal)) = 0.d0 enddo endif enddo endif ! --- Pressure, Density, Temperature, Total Energy ! Only 2 out of these 4 variables are independent: one may choose to ! initialize any pair of variables picked out of these 4, except ! (Temperature-Energy). The remaining 2 variables will be deduced ! automatically. ! Initialize 2 and only 2 variables ! To do so, set iutile=1 for each of the 2 selected variables ! and iutile=0 for each of the 2 others ! In the example provided below, Pressure and Temperature are ! initialized. ! iccfth indicates which variables have been set: ! it is completed automatically for each variable and ! it must not be modified. iccfth = 10000 ! 1. Pressure (Pa) iutile = 1 if(iutile.eq.1) then iccfth = iccfth*2 do iel = 1, ncel rtp(iel,ipr (iphas) ) = 101325.d0 enddo endif ! 2. Density (kg/m3) iutile = 0 if(iutile.eq.1) then iccfth = iccfth*3 do iel = 1, ncel rtp(iel,isca(irho (iphas))) = ro0(iphas) enddo endif ! 3. Temperature (K -- Warning: Kelvin) iutile = 1 if(iutile.eq.1) then iccfth = iccfth*5 do iel = 1, ncel rtp(iel,isca(itempk(iphas))) = 293.d0 enddo endif ! 4. Total Energy (J/kg) iutile = 0 if(iutile.eq.1) then iccfth = iccfth*7 do iel = 1, ncel rtp(iel,isca(ienerg(iphas))) = cv0(iphas)*t0(iphas) enddo endif ! ** The following subroutine returns automatically the values for the ! two remaining variables that need to be computed, using the ! indicator iccfth. call uscfth & !========== ( idebia , idebra , & ndim , ncelet , ncel , nfac , nfabor , nfml , nprfml , & nnod , lndfac , lndfbr , ncelbr , & nvar , nscal , nphas , & iccfth , imodif , iphas , & nideve , nrdeve , nituse , nrtuse , & ifacel , ifabor , ifmfbr , ifmcel , iprfml , & ipnfac , nodfac , ipnfbr , nodfbr , & idevel , ituser , ia , & xyzcen , surfac , surfbo , cdgfac , cdgfbo , xyznod , volume , & dt , rtp , rtp , propce , propfa , propfb , & coefa , coefb , & w1 , w2 , w3 , w4 , & rdevel , rtuser , ra ) endif !---- ! Formats !---- 9001 format( & /, & ' uscfxi: User defined initialization of the variables.',/, & /) !---- ! End !---- return end subroutine