anchor.diff:    integer, intent(OUT) :: my_point_index, target_point_index
     real(kind=8) :: xa1, ya1, xa2, ya2
     type(coef_type), pointer :: c_ptr
  -  integer :: target_polygon_index
  -  integer :: i1, i2
  +  !integer :: target_polygon_index
  +  integer :: i1
  +  !integer :: i2
   
   
     c_ptr => cb
  
positive.code:   character(len=*), intent(IN) :: coef_cal_type !<
    integer, intent(IN) :: side_div_num           !<
    type(polygon_type), pointer :: target_polygon
    real(kind=8) :: my_area
    real(kind=8) :: overlap_area, overlap_sum
    real(kind=8) :: point_x(MAX_POINT), point_y(MAX_POINT)
    real(kind=8) :: temp_coef
  
negative.code:   integer, intent(IN) :: div_num
    real(kind=8) :: point_pos_x, point_pos_y
    integer :: i, j, p
    type(polygon_type), pointer :: rect_polygon, target_polygon
    real(kind=8), dimension(div_num+1) :: x12, x23, x43, x14, y12, y23, y43, y14
    real(kind=8) :: x_pos(4), y_pos(4), x1, y1, x2, y2, x3, y3
    real(kind=8) :: px1, px2, px3, px4, py1, py2, py3, py4
  
    ! set rectangular polygon
    if (pa%num_of_point == 4) then ! rectangular polygon
  