Subroutine to initialize a deterministic 'smooth' problem
subroutine init_fe_prob
implicit none
if (mod(nx, 2) /= 1) then
write(OPU,*) 'It is better to use an odd number of nodes in the x dir to have profile for comparison with analytical sol.'
endif
fe_f%m%nx = nx
fe_f%m%ny = ny
fe_f%m%lx = lx
fe_f%m%ly = ly
fe_f%m%zx = 0.
fe_f%m%zy = 0.
call create_rect_FE_film( data_f = data_f, &
num_p = num_pts, &
fe_f = fe_f)
write(OPU,*) 'film created'
return
endsubroutine init_fe_prob