!------------------------------------------------------------------------------- ! Code_Saturne version 2.0.0-beta2 ! -------------------------- ! This file is part of the Code_Saturne Kernel, element of the ! Code_Saturne CFD tool. ! Copyright (C) 1998-2008 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 usproj & !================ ( idbia0 , idbra0 , & ndim , ncelet , ncel , nfac , nfabor , nfml , nprfml , & nnod , lndfac , lndfbr , ncelbr , & nvar , nscal , nphas , & nbpmax , nvp , nvep , nivep , ntersl , nvlsta , nvisbr , & nideve , nrdeve , nituse , nrtuse , & ifacel , ifabor , ifmfbr , ifmcel , iprfml , maxelt , lstelt , & ipnfac , nodfac , ipnfbr , nodfbr , itepa , & idevel , ituser , ia , & xyzcen , surfac , surfbo , cdgfac , cdgfbo , xyznod , volume , & dt , rtpa , rtp , propce , propfa , propfb , & coefa , coefb , & ettp , ettpa , tepa , statis , stativ , tslagr , parbor , & rdevel , rtuser , ra ) !=============================================================================== ! FONCTION : ! -------- ! MODIFICATION UTILISATEUR EN FIN DE PAS DE TEMPS ! TOUT EST POSSIBLE, ! ON DONNE ICI PLUSIEURS EXEMPLES : ! - CALCUL DE BILAN THERMIQUE ! (au besoin, voir "ADAPTATION A UN SCALAIRE QUELCONQUE") ! - CALCUL DES EFFORTS GLOBAUX SUR UN SOUS-ENSEMBLE DE FACES ! - MODIFICATION ARBITRAIRE D'UNE VARIABLE DE CALCUL ! - EXTRACTION D'UN PROFIL 1D ! - IMPRESSION D'UN MOMENT ! - EXEMPLES D'UTILISATION DES ROUTINES DE PARALLELISME ! CES EXEMPLES SONT DONNES EN SUPPOSANT UN CAS AVEC PERIODICITE ! (IPERIO .GT.0) ET PARALLELISME (IRANGP.GE.0). ! LE CALCUL DE BILAN THERMIQUE FOURNIT EN OUTRE UNE TRAME POUR ! PLUSIEURS CHOSES ! - CALCUL DE GRADIENT (AVEC LES PRECAUTIONS UTILES EN PARALLELE ET ! PERIODIQUE) ! - CALCUL DE GRANDEUR DEPENDANT DES VALEURS AUX CELLULES VOISINES ! D'UNE FACE (AVEC LES PRECAUTIONS A PRENDRE EN PARALLELE ET ! PERIODIQUE : VOIR L'ECHANGE DE DT ET DE CP) ! - CALCUL D'UNE SOMME SUR LES PROCESSEURS LORS D'UN CALCUL ! PARALLELE (PARSOM) ! IDENTIFICATION DES CELLULES/FACES DE BORD/FACES INTERNES ! ======================================================== ! Les commandes GETCEL, GETFBR et GETFAC permettent d'identifier ! respectivement les cellules, faces ou faces de bord en fonction ! de differents criteres. ! GETCEL(CHAINE,NLELT,LSTELT) : ! - CHAINE est une chaine de caractere fournie par l'utilisateur ! qui donne les criteres de selection ! - NLTELT est renvoye par la commande. C'est un entier qui ! correspond au nombre de cellules trouveees repondant au ! critere ! - LSTELT est renvoye par la commande. C'est un tableau d'entiers ! de taille NLTELT donnant la liste des cellules trouvees ! repondant au critere. ! CHAINE peut etre constitue de : ! - references de couleurs (ex. : 1, 8, 26, ... ! - references de groupes (ex. : entrees, groupe1, ...) ! - criteres geometriques (ex. X<0.1, Y>=0.25, ...) ! Ces criteres peuvent etre combines par des operateurs logiques ! (AND et OR) et des parentheses ! ex. : '1 AND (groupe2 OR groupe3) AND Y<1' permettra de recuperer ! les cellules de couleur 1, appartenant aux groupes 'groupe2' ! ou 'groupe3' et de coordonnee Y inferieure a 1. ! La syntaxe des commandes GETFBR et GETFAC est identique. !------------------------------------------------------------------------------- ! Arguments !__________________.____._____.________________________________________________. ! nom !type!mode ! role ! !__________________!____!_____!________________________________________________! ! idbia0 ! e ! <-- ! numero de la 1ere case libre dans ia ! ! idbra0 ! e ! <-- ! numero de la 1ere case libre dans ra ! ! ndim ! e ! <-- ! dimension de l'espace ! ! ncelet ! e ! <-- ! nombre d'elements halo compris ! ! ncel ! e ! <-- ! nombre d'elements actifs ! ! nfac ! e ! <-- ! nombre de faces internes ! ! nfabor ! e ! <-- ! nombre de faces de bord ! ! nfml ! e ! <-- ! nombre de familles d entites ! ! nprfml ! e ! <-- ! nombre de proprietese des familles ! ! nnod ! e ! <-- ! nombre de sommets ! ! lndfac ! e ! <-- ! longueur du tableau nodfac (optionnel ! ! lndfbr ! e ! <-- ! longueur du tableau nodfbr (optionnel ! ! ncelbr ! e ! <-- ! nombre d'elements ayant au moins une ! ! ! ! ! face de bord ! ! nvar ! e ! <-- ! nombre total de variables ! ! nscal ! e ! <-- ! nombre total de scalaires ! ! nphas ! e ! <-- ! nombre de phases ! ! nbpmax ! e ! <-- ! nombre max de particules autorise ! ! nvp ! e ! <-- ! nombre de variables particulaires ! ! nvep ! e ! <-- ! nombre info particulaires (reels) ! ! nivep ! e ! <-- ! nombre info particulaires (entiers) ! ! ntersl ! e ! <-- ! nbr termes sources de couplage retour ! ! nvlsta ! e ! <-- ! nombre de var statistiques lagrangien ! ! nvisbr ! e ! <-- ! nombre de statistiques aux frontieres ! ! nideve nrdeve ! e ! <-- ! longueur de idevel rdevel ! ! nituse nrtuse ! e ! <-- ! longueur de ituser rtuser ! ! ifacel ! te ! <-- ! elements voisins d'une face interne ! ! (2, nfac) ! ! ! ! ! ifabor ! te ! <-- ! element voisin d'une face de bord ! ! (nfabor) ! ! ! ! ! ifmfbr ! te ! <-- ! numero de famille d'une face de bord ! ! (nfabor) ! ! ! ! ! ifmcel ! te ! <-- ! numero de famille d'une cellule ! ! (ncelet) ! ! ! ! ! iprfml ! te ! <-- ! proprietes d'une famille ! ! nfml ,nprfml ! ! ! ! ! maxelt ! e ! <-- ! nb max d'elements (cell,fac,fbr) ! ! lstelt(maxelt) te ! --- ! tableau de travail ! ! ipnfac ! te ! <-- ! position du premier noeud de chaque ! ! (nfac+1) ! ! ! face interne dans nodfac (optionnel) ! ! nodfac ! te ! <-- ! connectivite faces internes/noeuds ! ! (lndfac) ! ! ! (optionnel) ! ! ipnfbr ! te ! <-- ! position du premier noeud de chaque ! ! (nfabor+1) ! ! ! face de bord dans nodfbr (optionnel) ! ! nodfbr ! te ! <-- ! connectivite faces de bord/noeuds ! ! (lndfbr ) ! ! ! (optionnel) ! ! itepa ! te ! <-- ! info particulaires (entiers) ! ! (nbpmax,nivep ! ! ! (cellule de la particule,...) ! ! idevel(nideve ! te ! <-- ! tab entier complementaire developemt ! ! ituser(nituse ! te ! <-- ! tab entier complementaire utilisateur ! ! ia(*) ! tr ! --- ! macro tableau entier ! ! xyzcen ! tr ! <-- ! point associes aux volumes de control ! ! (ndim,ncelet ! ! ! ! ! surfac ! tr ! <-- ! vecteur surface des faces internes ! ! (ndim,nfac) ! ! ! ! ! surfbo ! tr ! <-- ! vecteur surface des faces de bord ! ! (ndim,nfabor) ! ! ! ! ! cdgfac ! tr ! <-- ! centre de gravite des faces internes ! ! (ndim,nfac) ! ! ! ! ! cdgfbo ! tr ! <-- ! centre de gravite des faces de bord ! ! (ndim,nfabor) ! ! ! ! ! xyznod ! tr ! <-- ! coordonnes des noeuds (optionnel) ! ! (ndim,nnod) ! ! ! ! ! volume ! tr ! <-- ! volume d'un des ncelet elements ! ! (ncelet ! ! ! ! ! dt(ncelet) ! tr ! <-- ! pas de temps ! ! rtp, rtpa ! tr ! <-- ! variables de calcul au centre des ! ! (ncelet,*) ! ! ! cellules (instant courant ou prec) ! ! propce ! tr ! <-- ! proprietes physiques au centre des ! ! (ncelet,*) ! ! ! cellules ! ! propfa ! tr ! <-- ! proprietes physiques au centre des ! ! (nfac,*) ! ! ! faces internes ! ! propfb ! tr ! <-- ! proprietes physiques au centre des ! ! (nfabor,*) ! ! ! faces de bord ! ! coefa, coefb ! tr ! <-- ! conditions aux limites aux ! ! (nfabor,*) ! ! ! faces de bord ! ! ettp ! tr ! <-- ! tableaux des variables liees ! ! (nbpmax,nvp) ! ! ! aux particules etape courante ! ! ettpa ! tr ! <-- ! tableaux des variables liees ! ! (nbpmax,nvp) ! ! ! aux particules etape precedente ! ! tepa ! tr ! <-- ! info particulaires (reels) ! ! (nbpmax,nvep) ! ! ! (poids statistiques,...) ! ! statis ! tr ! <-- ! moyennes statistiques ! !(ncelet,nvlsta ! ! ! ! ! stativ ! tr ! <-- ! cumul pour les variances des ! !(ncelet, ! ! ! statistiques volumiques ! ! nvlsta-1) ! ! ! ! ! tslagr ! tr ! <-- ! terme de couplage retour du ! !(ncelet,ntersl ! ! ! lagrangien sur la phase porteuse ! ! parbor ! tr ! <-- ! infos sur interaction des particules ! !(nfabor,nvisbr ! ! ! aux faces de bord ! ! rdevel(nrdeve ! tr ! <-- ! tab reel complementaire developemt ! ! rtuser(nrtuse ! tr ! <-- ! tab reel complementaire utilisateur ! ! ra(*) ! tr ! --- ! macro tableau reel ! !__________________!____!_____!________________________________________________! ! TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU) ! L (LOGIQUE) .. ET TYPES COMPOSES (EX : TR TABLEAU REEL) ! MODE : <-- donnee, --> resultat, <-> Donnee modifiee ! --- tableau de travail !=============================================================================== implicit none !=============================================================================== ! DONNEES EN COMMON !=============================================================================== include "dimfbr.h" include "paramx.h" include "pointe.h" include "numvar.h" include "optcal.h" include "cstphy.h" include "cstnum.h" include "entsor.h" include "lagpar.h" include "lagran.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 nbpmax , nvp , nvep , nivep integer ntersl , nvlsta , nvisbr integer nideve , nrdeve , nituse , nrtuse integer ifacel(2,nfac) , ifabor(nfabor) integer ifmfbr(nfabor) , ifmcel(ncelet) integer iprfml(nfml,nprfml) integer maxelt, lstelt(maxelt) integer ipnfac(nfac+1), nodfac(lndfac) integer ipnfbr(nfabor+1), nodfbr(lndfbr) integer itepa(nbpmax,nivep) integer idevel(nideve), ituser(nituse) integer 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,*), rtpa(ncelet,*) double precision propce(ncelet,*) double precision propfa(nfac,*), propfb(ndimfb,*) double precision coefa(ndimfb,*), coefb(ndimfb,*) double precision ettp(nbpmax,nvp) , ettpa(nbpmax,nvp) double precision tepa(nbpmax,nvep) double precision statis(ncelet,nvlsta), stativ(ncelet,nvlsta-1) double precision tslagr(ncelet,ntersl) double precision parbor(nfabor,nvisbr) double precision rdevel(nrdeve), rtuser(nrtuse), ra(*) ! VARIABLES LOCALES integer idebia, idebra integer iel , ielg , ifac , ifacg , ivar integer iel1 , iel2 , ieltsm integer iortho , impout integer ifinia , ifinra integer igradx , igrady , igradz integer itravx , itravy , itravz , itreco integer inc , iccocg integer nswrgp , imligp , iphydp , iwarnp integer iutile , iphas , iclvar , iii integer ipcrom , ipcvst , iflmas , iflmab , ipccp, ipcvsl integer idimte , itenso , iscal integer ii , nbr , irangv , irang1 , npoint integer imom , ipcmom , idtcm integer itab(3), iun integer ncesmp , icesmp , ismacp , itpsmp integer ilelt , nlelt double precision xrtpa , xrtp double precision xbilan , xbilvl , xbilpa , xbilpt double precision xbilsy , xbilen , xbilso , xbildv double precision xbilmi , xbilma double precision epsrgp , climgp , extrap double precision xfluxf , xgamma double precision diipbx, diipby, diipbz, surfbn, distbr double precision visct, flumab , xcp , xvsl, cp0iph, rrr double precision xfor(3), xyz(3), xabs, xu, xv, xw, xk, xeps, XCOF(3) integer, parameter :: out_unit=20 !=============================================================================== !=============================================================================== ! 1. INITIALISATION !=============================================================================== ! ---> Gestion memoire idebia = idbia0 idebra = idbra0 !=============================================================================== ! 3. EXEMPLE : CALCUL DES EFFORTS GLOBAUX SUR UN SOUS-ENSEMBLE DE FACES ! A FAIRE AVEC PRECAUTIONS ... ! L'UTILISATEUR PREND SES RESPONSABILITES. !=============================================================================== ! ---------------------------------------------- ! Il est assez courant que l'on oublie d'eliminer cet exemple ! de ce sous-programme. ! On a donc prevu le test suivant pour eviter les mauvaises surprises iutile = 0 if(iutile.eq.0) return ! ---------------------------------------------- ! Si les efforts ont bien ete calcules : if (ineedf.eq.1) then do ii = 1, ndim xfor(ii) = 0.d0 XCOF(ii) = 0.d0 enddo ! open the file OPEN(FILE='drag_lift.dat',UNIT=IMPUSR(1),ACCESS='APPEND') CALL GETFBR('Perfil',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) do ii = 1, ndim xfor(ii) = xfor(ii) + ra(iforbr + (ifac-1)*ndim + ii-1) enddo enddo if (irangp.ge.0) then call parrsm(ndim,xfor) endif endif ! calculate x,y,zz coefficents DO II=1,NDIM XCOF(II) = 2.D0*xfor(II)/RO0(1)/UREF(1)/UREF(1) ENDDO ! print our coefficients for current timestep WRITE(IMPUSR(1),'(I12,1X,3(E12.5,1X))') NTCABS,XCOF(1),XCOF(2),XCOF(3) ! ! close CLOSE(UNIT=IMPUSR(1)) return end