@@ -594,11 +594,11 @@ contains
594594 eta(3 ) = q_prim_vf(eqn_idx%mom%beg + 1 )%sf(cell(1 ), cell(2 ) + 1 , cell(3 ))
595595 else if (i == 3 ) then
596596 xi(1 ) = z_cc(cell(3 ) - 1 )
597- eta(1 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ) - 1 )
597+ eta(1 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ) - 1 )
598598 xi(2 ) = z_cc(cell(3 ))
599- eta(2 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ))
599+ eta(2 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ))
600600 xi(3 ) = z_cc(cell(3 ) + 1 )
601- eta(3 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ) + 1 )
601+ eta(3 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ) + 1 )
602602 end if
603603
604604 L(1 ) = ((pos - xi(2 ))* (pos - xi(3 )))/ ((xi(1 ) - xi(2 ))* (xi(1 ) - xi(3 )))
@@ -631,15 +631,15 @@ contains
631631 eta(5 ) = q_prim_vf(eqn_idx%mom%beg + 1 )%sf(cell(1 ), cell(2 ) + 2 , cell(3 ))
632632 else if (i == 3 ) then
633633 xi(1 ) = z_cc(cell(3 ) - 2 )
634- eta(1 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ) - 2 )
634+ eta(1 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ) - 2 )
635635 xi(2 ) = z_cc(cell(3 ) - 1 )
636- eta(2 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ) - 1 )
636+ eta(2 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ) - 1 )
637637 xi(3 ) = z_cc(cell(3 ))
638- eta(3 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ))
638+ eta(3 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ))
639639 xi(4 ) = z_cc(cell(3 ) + 1 )
640- eta(4 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ) + 1 )
640+ eta(4 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ) + 1 )
641641 xi(5 ) = z_cc(cell(3 ) + 2 )
642- eta(5 ) = q_prim_vf(eqn_idx%mom%beg )%sf(cell(1 ), cell(2 ), cell(3 ) + 2 )
642+ eta(5 ) = q_prim_vf(eqn_idx%mom%end )%sf(cell(1 ), cell(2 ), cell(3 ) + 2 )
643643 end if
644644
645645 L(1 ) = ((pos - xi(2 ))* (pos - xi(3 ))* (pos - xi(4 ))* (pos - xi(5 )))/ ((xi(1 ) - xi(2 ))* (xi(1 ) - xi(3 ))* (xi(1 ) - xi(4 )) &
0 commit comments