!------------------------------------------------------------------------------- ! Code_Saturne version 3.0.1 ! -------------------------- subroutine usphyv & ( nvar , nscal , & ibrom , & dt , rtp , rtpa , & propce , propfa , propfb ) use paramx use pointe use numvar use optcal use cstphy use entsor use parall use period use field 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,*) ! Local variables integer ivart, iclvar, iel, ifac integer ipcrom, ipbrom, ipcvis, ipccp integer ipcvsl, ith, iscal, ii double precision vara, varb, varc, varam, varbm, varcm, vardm double precision varal, varbl, varcl, vardl double precision varac, varbc double precision xrtp double precision, dimension(:), pointer :: coefap, coefbp if (.true.) then ! Position of variables, coefficients ! ----------------------------------- ! --- Number of the thermal variable ! (and of its boundary conditions) ! To use user scalar 2 instead, write 'ivart = isca(2)' if (iscalt.gt.0) then ivart = isca(iscalt) else write(nfecra,9010) iscalt call csexit (1) endif ! --- Position of boundary conditions for variable 'ivart' iclvar = iclrtp(ivart,icoef) ! --- Rank of density ! in 'propce', physical properties at element centers: 'ipcrom' ! in 'propfb', physical properties at boundary face centers: 'ipbrom' ipcrom = ipproc(irom) ipbrom = ipprob(irom) ! --- Coefficients of laws chosen by the user ! Values given here are fictitious vara = -4.0668d-3 varb = -5.0754d-2 varc = 1000.9d0 ! Density at cell centers !------------------------ ! law rho = t * ( a * t + b) + c ! so propce(iel, ipcrom) = xrtp * (vara*xrtp+varb) + varc ! Volumic thermal expansion coefficient !-------------------------------------- ! law beta = -1/rho * (d rho / d T) ! so propce(iel, ipproc(ibeta)) = (-1.d0/propce(iel,ipcrom))*(2.d0*vara*xrtp+varb) do iel = 1, ncel xrtp = rtp(iel,ivart) propce(iel,ipcrom) = xrtp * (vara*xrtp+varb) + varc propce(iel,ipproc(ibeta))= (-1.d0/propce(iel,ipcrom))*(2.d0*vara*xrtp+varb) enddo endif ! --- Test on .false. !-------- ! Formats !-------- #if defined(_CS_LANG_FR) 1000 format( & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/,& '@ @@ ATTENTION : ARRET LORS DU CALCUL DES GRANDEURS PHYSIQUES',/,& '@ ========= ',/,& '@ DONNEES DE CALCUL INCOHERENTES ',/,& '@ ',/,& '@ usipph indique que la chaleur specifique est uniforme ',/,& '@ ICP = ',I10 ,' alors que ',/,& '@ usphyv impose une chaleur specifique variable. ',/,& '@ ',/,& '@ Le calcul ne sera pas execute. ',/,& '@ ',/,& '@ Modifier usipph ou usphyv. ',/,& '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/) 1010 format( & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@ ',/,& '@ @@ ATTENTION : ARRET LORS DU CALCUL DES GRANDEURS PHYSIQUES',/,& '@ ========= ',/,& '@ DONNEES DE CALCUL INCOHERENTES ',/,& '@ ',/,& '@ Pour le scalaire ',I10 ,/,& '@ usipsc indique que la diffusivite est uniforme ',/,& '@ IVISLS(',I10 ,') = ',I10 ,' alors que ',/,& '@ usphyv impose une diffusivite variable. ',/,& '@ ',/,& '@ Le calcul ne sera pas execute. ',/,& '@ ',/,& '@ Modifier usipsc 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 ',/,& '@ usipsu. 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',/, & '@',/, & '@ usipph specifies that the specific heat is uniform',/, & '@ icp = ',i10 ,' while',/, & '@ usphyv prescribes a variable specific heat.',/, & '@',/, & '@ The calculation will not be run.',/, & '@',/, & '@ Modify usipph or usphyv.',/, & '@ ',/,& '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/) 1010 format( & '@',/, & '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/, & '@ @@ WARNING: stop when computing physical quantities',/, & '@ =======',/, & '@ Inconsistent calculation data',/, & '@',/, & '@ For scalar', i10,/, & '@ usipsu specifies that the diffusivity is uniform',/, & '@ ivislc(',i10 ,') = ',i10 ,' while',/, & '@ usphyv prescribes a variable diffusivity.',/, & '@',/, & '@ The calculation will not be run.',/, & '@',/, & '@ Modify usipsu 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',/, & '@ usipsu. If a scalar should represent the,',/, & '@ temperature, check that iscalt has been defined',/, & '@',/, & '@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,& '@',/) #endif !---- ! End !---- return end subroutine usphyv