!------------------------------------------------------------------------------- ! Code_Saturne version 2.1.0 ! -------------------------- ! This file is part of Code_Saturne, a general-purpose CFD tool. ! ! Copyright (C) 1998-2011 EDF S.A. ! ! This program 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. ! ! This program 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 ! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin ! Street, Fifth Floor, Boston, MA 02110-1301, USA. !------------------------------------------------------------------------------- subroutine usphyv & !================ ( nvar , nscal , & ibrom , & dt , rtp , rtpa , & propce , propfa , propfb , & coefa , coefb ) !=============================================================================== ! Module files !=============================================================================== use paramx use pointe use numvar use optcal use cstphy use entsor use parall use period use mesh !=============================================================================== implicit none ! Arguments integer nvar , nscal integer ibrom double precision dt(ncelet), rtp(ncelet,*), rtpa(ncelet,*) double precision propce(ncelet,*) double precision propfa(nfac,*), propfb(nfabor,*) double precision coefa(nfabor,*), coefb(nfabor,*) ! Local variables integer ivart, iclvar, iel integer ipcrom, ipbrom, ipcvis, ipccp integer ipcvsl, ith, iscal, ii integer iutile double precision vara, varb, varc, varam, varbm, varcm, vardm double precision varal, varbl, varcl, vardl double precision varac, varbc double precision xrtp iutile = 1 if(iutile.eq.1) then ipcrom = ipproc(irom) do iel = 1, ncel propce(iel,ipcrom) = 140.00 enddo endif ! --- Test on 'iutile' !=============================================================================== !=============================================================================== ! Formats !---- #if defined(_CS_LANG_FR) 1000 format( & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/,& '@ @@ ATTENTION : ARRET LORS DU CALCUL DES GRANDEURS PHYSIQUES',/,& '@ ========= ',/,& '@ DONNEES DE CALCUL INCOHERENTES ',/,& '@ ',/,& '@ usini1 indique que la chaleur specifique est uniforme ',/,& '@ ICP = ',I10 ,' alors que ',/,& '@ usphyv impose une chaleur specifique variable. ',/,& '@ ',/,& '@ Le calcul ne sera pas execute. ',/,& '@ ',/,& '@ Modifier usini1 ou usphyv. ',/,& '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/) 1010 format( & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/,& '@ @@ ATTENTION : ARRET LORS DU CALCUL DES GRANDEURS PHYSIQUES',/,& '@ ========= ',/,& '@ DONNEES DE CALCUL INCOHERENTES ',/,& '@ ',/,& '@ Pour le scalaire ',I10 ,/,& '@ usini1 indique que la diffusivite est uniforme ',/,& '@ IVISLS(',I10 ,') = ',I10 ,' alors que ',/,& '@ usphyv impose une diffusivite variable. ',/,& '@ ',/,& '@ Le calcul ne sera pas execute. ',/,& '@ ',/,& '@ Modifier usini1 ou usphyv. ',/,& '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/) 9010 format( & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/,& '@ @@ ATTENTION : ARRET LORS DU CALCUL DES GRANDEURS PHYSIQUES',/,& '@ ========= ',/,& '@ APPEL A csexit DANS LE SOUS PROGRAMME usphyv ',/,& '@ ',/,& '@ La variable dont dependent les proprietes physiques ne ',/,& '@ semble pas etre une variable de calcul. ',/,& '@ En effet, on cherche a utiliser la temperature alors que',/,& '@ ISCALT = ',I10 ,/,& '@ Le calcul ne sera pas execute. ',/,& '@ ',/,& '@ Verifier le codage de usphyv (et le test lors de la ',/,& '@ definition de IVART). ',/,& '@ Verifier la definition des variables de calcul dans ',/,& '@ usini1. Si un scalaire doit jouer le role de la ',/,& '@ temperature, verifier que ISCALT a ete renseigne. ',/,& '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/) #else 1000 format( & '@',/, & '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/, & '@ @@ WARNING: stop when computing physical quantities',/, & '@ =======',/, & '@ Inconsistent calculation data',/, & '@',/, & '@ usini1 specifies that the specific heat is uniform',/, & '@ icp = ',i10 ,' while',/, & '@ usphyv prescribes a variable specific heat.',/, & '@',/, & '@ The calculation will not be run.',/, & '@',/, & '@ Modify usini1 or usphyv.',/, & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/) 1010 format( & '@',/, & '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/, & '@ @@ WARNING: stop when computing physical quantities',/, & '@ =======',/, & '@ Inconsistent calculation data',/, & '@',/, & '@ For scalar', i10,/, & '@ usini1 specifies that the diffusivity is uniform',/, & '@ ivislc(',i10 ,') = ',i10 ,' while',/, & '@ usphyv prescribes a variable diffusivity.',/, & '@',/, & '@ The calculation will not be run.',/, & '@',/, & '@ Modify usini1 or usphyv.',/, & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/) 9010 format( & '@',/, & '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/, & '@ @@ WARNING: stop when computing physical quantities',/, & '@ =======',/, & '@',/, & '@ The variable on which physical properties depend does',/, & '@ seem to be a calculation variable.',/, & '@ Indeed, we are trying to use the temperature while',/, & '@ iscalt = ',i10 ,/,& '@',/, & '@ The calculation will not be run.',/, & '@',/, & '@ Check the programming in usphyv (and the test when',/, & '@ defining ivart).',/, & '@ Check the definition of calculation variables in',/, & '@ usini1. If a scalar should represent the,',/, & '@ temperature, check that iscalt has been defined',/, & '@',/, & '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/) #endif !---- ! End !---- return end subroutine