init_fe_prob Subroutine

private subroutine init_fe_prob()


Arguments

None

Calls

proc~~init_fe_prob~~CallsGraph proc~init_fe_prob init_fe_prob proc~create_rect_fe_film create_rect_FE_film proc~init_fe_prob->proc~create_rect_fe_film proc~create_rect_x_ymesh create_rect_x_ymesh proc~create_rect_fe_film->proc~create_rect_x_ymesh

Called by

proc~~init_fe_prob~~CalledByGraph proc~init_fe_prob init_fe_prob proc~test_rough_fe test_rough_fe proc~test_rough_fe->proc~init_fe_prob proc~test_bearing_x_fe test_bearing_x_fe proc~test_bearing_x_fe->proc~init_fe_prob proc~test_slider_fe test_slider_fe proc~test_slider_fe->proc~init_fe_prob proc~test_pocket_fe test_pocket_fe proc~test_pocket_fe->proc~init_fe_prob proc~test_bearing_y_fe test_bearing_y_fe proc~test_bearing_y_fe->proc~init_fe_prob proc~run_test run_test proc~run_test->proc~test_rough_fe proc~run_test->proc~test_bearing_x_fe proc~run_test->proc~test_slider_fe proc~run_test->proc~test_pocket_fe proc~run_test->proc~test_bearing_y_fe program~main main program~main->proc~run_test

Contents

Source Code


Source Code

   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