cooling_towers example
Local variables to be added
The following local variables need to be defined for the examples in this section:
integer iel, iutile
integer ilelt, nlelt
double precision d2s3
double precision, dimension(:,:), pointer :: vel
integer, allocatable, dimension(:) :: lstelt
double precision, dimension(:), pointer :: cvar_temp4, cvar_humid
Allocation
Before user initialization, work arrays lstelt must be allocated, like in basic example.
Initialization
The following initialization block needs to be added for the following examples:
d2s3 = 2.d0/3.d0
cvar_temp4(iel) = 11.d0
cvar_humid(iel) = 0.0063d0
enddo
call
getcel(
'6', nlelt, lstelt)
do ilelt = 1, nlelt
iel = lstelt(ilelt)
vel(1,iel) = -0.5d0
cvar_temp4(iel) = 20.d0
cvar_humid(iel) = 0.012d0
enddo
endif
Finalization
At the end of the subroutine, it is recommended to deallocate the work array lstelt, like in basic example.