Local variables to be added
The following local variables need to be defined for the examples in this section:
integer iel, ii, jj
double precision x, pormin, pormax, hc, ll, dhc
double precision, dimension(:), pointer :: cpro_porosi
Test to remove for use of subroutine end
The following initialization block needs to be added for the following examples:
Retrieve porosity field
Example: fixe a linear by part porosity profile
if (x.le.(ll/2.d0)) then
hc = 1.d0 - 2.d0*dhc*x/ll
else
hc = 1.d0 - dhc + dhc*(2.d0*x-ll)/ll
endif
cpro_porosi(iel) = hc
pormin = min(pormin,cpro_porosi(iel))
pormax = max(pormax,cpro_porosi(iel))
enddo
Periodicity and parallelism treatment: