allocate_cut_cell_arrays.f

Go to the documentation of this file.

      SUBROUTINE allocate_cut_cell_arrays

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
!
!  Module name: ALLOCATE_ARRAYS
!  Purpose: allocate arrays
!                                                                      C
!  Author: Jeff Dietiker                              Date: 21-Feb-08  C
!  Reviewer:
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      Use indices

      USE cutcell
      USE stl
      USE discretelement

      IMPLICIT NONE

      dimension_max_cut_cell = int(fac_dim_max_cut_cell*dimension_3g)

      Allocate(  interior_cell_at(dimension_3) )

      Allocate( xg_e(0:dimension_i) )
      Allocate( yg_n(0:dimension_j) )
      Allocate( zg_t(0:dimension_k) )

      Allocate(  x_u(dimension_3) )
      Allocate(  y_u(dimension_3) )
      Allocate(  z_u(dimension_3) )

      Allocate(  x_v(dimension_3) )
      Allocate(  y_v(dimension_3) )
      Allocate(  z_v(dimension_3) )

      Allocate(  x_w(dimension_3) )
      Allocate(  y_w(dimension_3) )
      Allocate(  z_w(dimension_3) )

      Allocate(  intersect_x(dimension_3) )
      Allocate(  intersect_y(dimension_3) )
      Allocate(  intersect_z(dimension_3) )

      Allocate(  x_int(dimension_3) )
      Allocate(  y_int(dimension_3) )
      Allocate(  z_int(dimension_3) )

      Allocate(  x_new_point(dimension_max_cut_cell) )
      Allocate(  y_new_point(dimension_max_cut_cell) )
      Allocate(  z_new_point(dimension_max_cut_cell) )

      Allocate(  x_new_u_point(dimension_max_cut_cell) )
      Allocate(  y_new_u_point(dimension_max_cut_cell) )
      Allocate(  z_new_u_point(dimension_max_cut_cell) )

      Allocate(  x_new_v_point(dimension_max_cut_cell) )
      Allocate(  y_new_v_point(dimension_max_cut_cell) )
      Allocate(  z_new_v_point(dimension_max_cut_cell) )

      Allocate(  x_new_w_point(dimension_max_cut_cell) )
      Allocate(  y_new_w_point(dimension_max_cut_cell) )
      Allocate(  z_new_w_point(dimension_max_cut_cell) )

      Allocate(  number_of_nodes(dimension_3) )
      Allocate(  number_of_u_nodes(dimension_3) )
      Allocate(  number_of_v_nodes(dimension_3) )
      Allocate(  number_of_w_nodes(dimension_3) )

     number_of_nodes   = 0
     number_of_u_nodes = 0
     number_of_v_nodes = 0
     number_of_w_nodes = 0

      Allocate(  connectivity(dimension_3,15) )
      Allocate(  connectivity_u(dimension_3,15) )
      Allocate(  connectivity_v(dimension_3,15) )
      Allocate(  connectivity_w(dimension_3,15) )

      Allocate(  partition(dimension_3) )

      Allocate(  wall_u_at(dimension_3) )
      Allocate(  wall_v_at(dimension_3) )
      Allocate(  wall_w_at(dimension_3) )

      wall_u_at = .false.
      wall_v_at = .false.
      wall_w_at = .false.

      Allocate( area_cut(dimension_3) )
      Allocate( area_u_cut(dimension_3) )
      Allocate( area_v_cut(dimension_3) )
      Allocate( area_w_cut(dimension_3) )


      Allocate( delx_ue(dimension_3) )
      Allocate( delx_uw(dimension_3) )
      Allocate( dely_un(dimension_3) )
      Allocate( dely_us(dimension_3) )
      Allocate( delz_ut(dimension_3) )
      Allocate( delz_ub(dimension_3) )

      Allocate( delx_ve(dimension_3) )
      Allocate( delx_vw(dimension_3) )
      Allocate( dely_vn(dimension_3) )
      Allocate( dely_vs(dimension_3) )
      Allocate( delz_vt(dimension_3) )
      Allocate( delz_vb(dimension_3) )

      Allocate( delx_we(dimension_3) )
      Allocate( delx_ww(dimension_3) )
      Allocate( dely_wn(dimension_3) )
      Allocate( dely_ws(dimension_3) )
      Allocate( delz_wt(dimension_3) )
      Allocate( delz_wb(dimension_3) )

      Allocate( x_u_ec(dimension_3) )
      Allocate( y_u_ec(dimension_3) )
      Allocate( z_u_ec(dimension_3) )
      Allocate( x_u_nc(dimension_3) )
      Allocate( y_u_nc(dimension_3) )
      Allocate( z_u_nc(dimension_3) )
      Allocate( x_u_tc(dimension_3) )
      Allocate( y_u_tc(dimension_3) )
      Allocate( z_u_tc(dimension_3) )

      Allocate( x_v_ec(dimension_3) )
      Allocate( y_v_ec(dimension_3) )
      Allocate( z_v_ec(dimension_3) )
      Allocate( x_v_nc(dimension_3) )
      Allocate( y_v_nc(dimension_3) )
      Allocate( z_v_nc(dimension_3) )
      Allocate( x_v_tc(dimension_3) )
      Allocate( y_v_tc(dimension_3) )
      Allocate( z_v_tc(dimension_3) )

      Allocate( x_w_ec(dimension_3) )
      Allocate( y_w_ec(dimension_3) )
      Allocate( z_w_ec(dimension_3) )
      Allocate( x_w_nc(dimension_3) )
      Allocate( y_w_nc(dimension_3) )
      Allocate( z_w_nc(dimension_3) )
      Allocate( x_w_tc(dimension_3) )
      Allocate( y_w_tc(dimension_3) )
      Allocate( z_w_tc(dimension_3) )

      Allocate( delh_scalar(dimension_3) )

      Allocate( delh_u(dimension_3) )
      Allocate( theta_ue(dimension_3) )
      Allocate( theta_ue_bar(dimension_3) )
      Allocate( theta_u_ne(dimension_3) )
      Allocate( theta_u_nw(dimension_3) )
      Allocate( theta_u_te(dimension_3) )
      Allocate( theta_u_tw(dimension_3) )
      Allocate( alpha_ue_c(dimension_3) )
      Allocate( noc_u_e(dimension_3) )
      Allocate( theta_un(dimension_3) )
      Allocate( theta_un_bar(dimension_3) )
      Allocate( alpha_un_c(dimension_3) )
      Allocate( noc_u_n(dimension_3) )
      Allocate( theta_ut(dimension_3) )
      Allocate( theta_ut_bar(dimension_3) )
      Allocate( alpha_ut_c(dimension_3) )
      Allocate( noc_u_t(dimension_3) )
      Allocate( a_upg_e(dimension_3) )
      Allocate( a_upg_w(dimension_3) )

      Allocate( delh_v(dimension_3) )
      Allocate( theta_v_ne(dimension_3) )
      Allocate( theta_v_se(dimension_3) )
      Allocate( theta_vn(dimension_3) )
      Allocate( theta_vn_bar(dimension_3) )
      Allocate( theta_v_nt(dimension_3) )
      Allocate( theta_v_st(dimension_3) )
      Allocate( theta_ve(dimension_3) )
      Allocate( theta_ve_bar(dimension_3) )
      Allocate( alpha_ve_c(dimension_3) )
      Allocate( noc_v_e(dimension_3) )
      Allocate( alpha_vn_c(dimension_3) )
      Allocate( noc_v_n(dimension_3) )
      Allocate( theta_vt(dimension_3) )
      Allocate( theta_vt_bar(dimension_3) )
      Allocate( alpha_vt_c(dimension_3) )
      Allocate( noc_v_t(dimension_3) )
      Allocate( a_vpg_n(dimension_3) )
      Allocate( a_vpg_s(dimension_3) )

      Allocate( delh_w(dimension_3) )
      Allocate( theta_w_te(dimension_3) )
      Allocate( theta_w_be(dimension_3) )
      Allocate( theta_w_tn(dimension_3) )
      Allocate( theta_w_bn(dimension_3) )
      Allocate( theta_wt(dimension_3) )
      Allocate( theta_wt_bar(dimension_3) )
      Allocate( theta_we(dimension_3) )
      Allocate( theta_we_bar(dimension_3) )
      Allocate( alpha_we_c(dimension_3) )
      Allocate( noc_w_e(dimension_3) )
      Allocate( theta_wn(dimension_3) )
      Allocate( theta_wn_bar(dimension_3) )
      Allocate( alpha_wn_c(dimension_3) )
      Allocate( noc_w_n(dimension_3) )
      Allocate( alpha_wt_c(dimension_3) )
      Allocate( noc_w_t(dimension_3) )
      Allocate( a_wpg_t(dimension_3) )
      Allocate( a_wpg_b(dimension_3) )


      Allocate( normal_s(dimension_3,3) )
      Allocate( normal_u(dimension_3,3) )
      Allocate( normal_v(dimension_3,3) )
      Allocate( normal_w(dimension_3,3) )

      Allocate( refp_s(dimension_3,3) )
      Allocate( refp_u(dimension_3,3) )
      Allocate( refp_v(dimension_3,3) )
      Allocate( refp_w(dimension_3,3) )

      Allocate(  oneodx_e_u(dimension_3) )
      Allocate(  oneody_n_u(dimension_3) )
      Allocate(  oneodz_t_u(dimension_3) )

      Allocate(  oneodx_e_v(dimension_3) )
      Allocate(  oneody_n_v(dimension_3) )
      Allocate(  oneodz_t_v(dimension_3) )

      Allocate(  oneodx_e_w(dimension_3) )
      Allocate(  oneody_n_w(dimension_3) )
      Allocate(  oneodz_t_w(dimension_3) )

      Allocate(  xn_int(dimension_3) )
      Allocate(  xn_u_int(dimension_3) )
      Allocate(  xn_v_int(dimension_3) )
      Allocate(  xn_w_int(dimension_3) )

      Allocate(  ye_int(dimension_3) )
      Allocate(  ye_u_int(dimension_3) )
      Allocate(  ye_v_int(dimension_3) )
      Allocate(  ye_w_int(dimension_3) )

      Allocate(  zt_int(dimension_3) )
      Allocate(  zt_u_int(dimension_3) )
      Allocate(  zt_v_int(dimension_3) )
      Allocate(  zt_w_int(dimension_3) )

      Allocate(  snap(dimension_3) )

      snap = .false.

      Allocate(  cut_treatment_at(dimension_3) )
      Allocate(  cut_u_treatment_at(dimension_3) )
      Allocate(  cut_v_treatment_at(dimension_3) )
      Allocate(  cut_w_treatment_at(dimension_3) )


      cut_treatment_at = .false.
      cut_u_treatment_at = .false.
      cut_v_treatment_at = .false.
      cut_w_treatment_at = .false.

      Allocate(  cut_cell_at(dimension_3) )
      Allocate(  cut_u_cell_at(dimension_3) )
      Allocate(  cut_v_cell_at(dimension_3) )
      Allocate(  cut_w_cell_at(dimension_3) )

      cut_cell_at = .false.
      cut_u_cell_at = .false.
      cut_v_cell_at = .false.
      cut_w_cell_at = .false.

      Allocate( small_cell_at(dimension_3) )
      small_cell_at = .false.

      Allocate( small_cell_flag(dimension_3) )
      small_cell_flag = 0

      Allocate(  blocked_cell_at(dimension_3) )
      Allocate(  blocked_u_cell_at(dimension_3) )
      Allocate(  blocked_v_cell_at(dimension_3) )
      Allocate(  blocked_w_cell_at(dimension_3) )

      blocked_cell_at   = .false.
      blocked_u_cell_at = .false.
      blocked_v_cell_at = .false.
      blocked_w_cell_at = .false.

      Allocate(  standard_cell_at(dimension_3) )
      Allocate(  standard_u_cell_at(dimension_3) )
      Allocate(  standard_v_cell_at(dimension_3) )
      Allocate(  standard_w_cell_at(dimension_3) )


      Allocate(  vorticity(dimension_3) )
      Allocate(  lambda2(dimension_3) )

      Allocate(  trd_g_out(dimension_3) )
      Allocate(  pp_g_out(dimension_3) )
      Allocate(  epp_out(dimension_3) )

      Allocate(  dudx_out(dimension_3) )
      Allocate(  dvdy_out(dimension_3) )
      Allocate(  delv_out(dimension_3) )

      Allocate(  u_master_of(dimension_3) )
      Allocate(  v_master_of(dimension_3) )
      Allocate(  w_master_of(dimension_3) )

      Allocate(  bc_id(dimension_3) )
      Allocate(  bc_u_id(dimension_3) )
      Allocate(  bc_v_id(dimension_3) )
      Allocate(  bc_w_id(dimension_3) )

      bc_id   = 0
      bc_u_id = 0
      bc_v_id = 0
      bc_w_id = 0

      Allocate(  debug_cg(dimension_3,15) )

      Allocate(  u_g_cc(dimension_3) )
      Allocate(  v_g_cc(dimension_3) )
      Allocate(  w_g_cc(dimension_3) )

      Allocate(  u_s_cc(dimension_3, dimension_m) )
      Allocate(  v_s_cc(dimension_3, dimension_m) )
      Allocate(  w_s_cc(dimension_3, dimension_m) )

      ALLOCATE(n_facet_at(dimension_3))
      n_facet_at = 0

      ALLOCATE(list_facet_at(dimension_3,dim_facets_per_cell))

      ALLOCATE(potential_cut_cell_at(dimension_3))


      Allocate(  f_at(dimension_3) )

      Allocate(  dwall(dimension_3) )


      RETURN
      END SUBROUTINE allocate_cut_cell_arrays