I need to know how can i affect my scalar temperature. My subroutine doesnt work.
How is the best practice ? what do i need to know ?
Any help is welcome.
Code Saturne 4.09
no specific module
scalar transported : temperature
Code: Select all
subroutine usphyv &
( nvar , nscal , &
mbrom , &
dt )
!===============================================================================
!===============================================================================
! Module files
!===============================================================================
use paramx
use pointe
use numvar
use optcal
use cstphy
use entsor
use parall
use period
use ppppar
use ppthch
use ppincl
use field
use mesh
!===============================================================================
implicit none
! Arguments
integer nvar , nscal
integer mbrom
double precision dt(ncelet)
! Local variables
integer ivart, iel, ifac
double precision, dimension(:), pointer :: cvar_scalt
if (iscalt.gt.0) then
ivart = isca(iscalt)
call field_get_val_s(ivarfl(ivart), cvar_scalt)
else
write(nfecra,9010) iscalt
call csexit (1)
endif
xvart= 0.5d0
do iel = 1, ncel
cvar_scalt(iel)= xvart + 0.5d0
enddo
return
end subroutine usphyv