/usr/share/code_saturne/user_examples/cs_user_extra_operations-global_efforts.f90 is in code-saturne-data 4.3.3+repack-1build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 | !-------------------------------------------------------------------------------
! Code_Saturne version 4.3.3
! --------------------------
! This file is part of Code_Saturne, a general-purpose CFD tool.
!
! Copyright (C) 1998-2016 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.
!-------------------------------------------------------------------------------
!===============================================================================
! Purpose:
! -------
!> \file cs_user_extra_operations-global_efforts.f90
!>
!> \brief This function is called at the end of each time step, and has a very
!> general purpose
!> (i.e. anything that does not have another dedicated user subroutine)
!>
!> See \subpage cs_user_extra_operations_examples and
!> \subpage cs_user_extra_operations-nusselt_calculation for examples.
!>
!> This is an example of cs_user_extra_operations.f90 which
!> performs global efforts
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
! Arguments
!______________________________________________________________________________.
! mode name role !
!______________________________________________________________________________!
!> \param[in] nvar total number of variables
!> \param[in] nscal total number of scalars
!> \param[in] dt time step (per cell)
!_______________________________________________________________________________
subroutine cs_f_user_extra_operations &
( nvar , nscal , &
dt )
!===============================================================================
!===============================================================================
! Module files
!===============================================================================
use paramx
use dimens, only: ndimfb
use pointe
use numvar
use optcal
use cstphy
use cstnum
use entsor
use lagran
use parall
use period
use ppppar
use ppthch
use ppincl
use mesh
use field
use field_operator
use turbomachinery
!===============================================================================
implicit none
! Arguments
integer nvar , nscal
double precision dt(ncelet)
! Local variables
!< [loc_var_dec]
integer ifac
integer ii
integer ilelt , nlelt
double precision xfor(3)
double precision, dimension(:,:), pointer :: bfprp_for
integer, allocatable, dimension(:) :: lstelt
!< [loc_var_dec]
!===============================================================================
!===============================================================================
! Initialization
!===============================================================================
if (iforbr.ge.0) call field_get_val_v(iforbr, bfprp_for)
! Allocate a temporary array for cells or interior/boundary faces selection
allocate(lstelt(max(ncel,nfac,nfabor)))
!===============================================================================
! Example: compute global efforts on a subset of faces
!===============================================================================
! If efforts have been calculated correctly:
!< [example_1]
if (iforbr.ge.0) then
do ii = 1, ndim
xfor(ii) = 0.d0
enddo
call getfbr('2 or 3', nlelt, lstelt)
!==========
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
do ii = 1, ndim
xfor(ii) = xfor(ii) + bfprp_for(ii, ifac)
enddo
enddo
if (irangp.ge.0) then
call parrsm(ndim,xfor)
endif
endif
!< [example_1]
! Deallocate the temporary array
deallocate(lstelt)
return
end subroutine cs_f_user_extra_operations
|