@@ -143,7 +143,7 @@ contains
143143 case (BC_PERIODIC)
144144 call s_periodic(q_prim_vf, 1 , 1 , k, l, pb_in, mv_in, q_T_sf)
145145 case (BC_SLIP_WALL)
146- call s_slip_wall(q_prim_vf, 1 , 1 , k, l,q_T_sf)
146+ call s_slip_wall(q_prim_vf, 1 , 1 , k, l, q_T_sf)
147147 case (BC_NO_SLIP_WALL)
148148 call s_no_slip_wall(q_prim_vf, 1 , 1 , k, l, q_T_sf)
149149 case (BC_DIRICHLET)
@@ -908,8 +908,8 @@ contains
908908 type(scalar_field), dimension (sys_size), intent (inout ) :: q_prim_vf
909909 integer , intent (in ) :: bc_dir, bc_loc
910910 integer , intent (in ) :: k, l
911- type(scalar_field), intent (inout ) :: q_T_sf
912- real (wp), dimension (num_species) :: Ys_in
911+ type(scalar_field), intent (inout ) :: q_T_sf
912+ real (wp), dimension (num_species) :: Ys_in
913913 real (wp) :: mix_mol_weight
914914
915915 integer :: j, i
@@ -927,7 +927,7 @@ contains
927927 end if
928928 end do
929929 end do
930- if (chemistry) then
930+ if (chemistry) then
931931 do j = 1 , buff_size
932932 do i = chemxb, chemxe
933933 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(0 , k, l)
@@ -945,10 +945,10 @@ contains
945945 else
946946 q_prim_vf(i)%sf(m + j, k, l) = &
947947 q_prim_vf(i)%sf(m, k, l)
948- end if
948+ end if
949949 end do
950950 end do
951- if (chemistry) then
951+ if (chemistry) then
952952 do j = 1 , buff_size
953953 do i = chemxb, chemxe
954954 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(m, k, l)
@@ -971,7 +971,7 @@ contains
971971 end if
972972 end do
973973 end do
974- if (chemistry) then
974+ if (chemistry) then
975975 do j = 1 , buff_size
976976 do i = chemxb, chemxe
977977 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, 0 , l)
@@ -992,15 +992,15 @@ contains
992992 end if
993993 end do
994994 end do
995- if (chemistry) then
995+ if (chemistry) then
996996 do j = 1 , buff_size
997997 do i = chemxb, chemxe
998998 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, n, l)
999999 end do
10001000 call get_mixture_molecular_weight(Ys_in, mix_mol_weight)
1001- q_T_sf%sf(k, n+ j, l) = q_prim_vf(E_idx)%sf(k, n, l)* mix_mol_weight/ (gas_constant* q_prim_vf(contxb)%sf(k, n, l))
1001+ q_T_sf%sf(k, n + j, l) = q_prim_vf(E_idx)%sf(k, n, l)* mix_mol_weight/ (gas_constant* q_prim_vf(contxb)%sf(k, n, l))
10021002 end do
1003- end if
1003+ end if
10041004 end if
10051005 elseif (bc_dir == 3 ) then !< z- direction
10061006 if (bc_loc == - 1 ) then !< bc_z%beg
@@ -1015,15 +1015,15 @@ contains
10151015 end if
10161016 end do
10171017 end do
1018- if (chemistry) then
1018+ if (chemistry) then
10191019 do j = 1 , buff_size
10201020 do i = chemxb, chemxe
10211021 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, l, 0 )
10221022 end do
10231023 call get_mixture_molecular_weight(Ys_in, mix_mol_weight)
10241024 q_T_sf%sf(k, l, - j) = q_prim_vf(E_idx)%sf(k, l, 0 )* mix_mol_weight/ (gas_constant* q_prim_vf(contxb)%sf(k, l, 0 ))
10251025 end do
1026- end if
1026+ end if
10271027 else !< bc_z%end
10281028 do i = 1 , sys_size
10291029 do j = 1 , buff_size
@@ -1036,15 +1036,15 @@ contains
10361036 end if
10371037 end do
10381038 end do
1039- if (chemistry) then
1039+ if (chemistry) then
10401040 do j = 1 , buff_size
10411041 do i = chemxb, chemxe
10421042 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, l, p)
10431043 end do
10441044 call get_mixture_molecular_weight(Ys_in, mix_mol_weight)
10451045 q_T_sf%sf(k, l, p + j) = q_prim_vf(E_idx)%sf(k, l, p)* mix_mol_weight/ (gas_constant* q_prim_vf(contxb)%sf(k, l, p))
10461046 end do
1047- end if
1047+ end if
10481048 end if
10491049 end if
10501050
@@ -1057,8 +1057,8 @@ contains
10571057 type(scalar_field), dimension (sys_size), intent (inout ) :: q_prim_vf
10581058 integer , intent (in ) :: bc_dir, bc_loc
10591059 integer , intent (in ) :: k, l
1060- type(scalar_field), intent (inout ) :: q_T_sf
1061- real (wp), dimension (num_species) :: Ys_in
1060+ type(scalar_field), intent (inout ) :: q_T_sf
1061+ real (wp), dimension (num_species) :: Ys_in
10621062 real (wp) :: mix_mol_weight
10631063
10641064 integer :: j, i
@@ -1082,7 +1082,7 @@ contains
10821082 end if
10831083 end do
10841084 end do
1085- if (chemistry) then
1085+ if (chemistry) then
10861086 do j = 1 , buff_size
10871087 do i = chemxb, chemxe
10881088 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(0 , k, l)
@@ -1109,7 +1109,7 @@ contains
11091109 end if
11101110 end do
11111111 end do
1112- if (chemistry) then
1112+ if (chemistry) then
11131113 do j = 1 , buff_size
11141114 do i = chemxb, chemxe
11151115 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(m, k, l)
@@ -1138,7 +1138,7 @@ contains
11381138 end if
11391139 end do
11401140 end do
1141- if (chemistry) then
1141+ if (chemistry) then
11421142 do j = 1 , buff_size
11431143 do i = chemxb, chemxe
11441144 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, 0 , l)
@@ -1165,7 +1165,7 @@ contains
11651165 end if
11661166 end do
11671167 end do
1168- if (chemistry) then
1168+ if (chemistry) then
11691169 do j = 1 , buff_size
11701170 do i = chemxb, chemxe
11711171 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, n, l)
@@ -1194,7 +1194,7 @@ contains
11941194 end if
11951195 end do
11961196 end do
1197- if (chemistry) then
1197+ if (chemistry) then
11981198 do j = 1 , buff_size
11991199 do i = chemxb, chemxe
12001200 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, l, 0 )
@@ -1222,7 +1222,7 @@ contains
12221222 end do
12231223 end do
12241224
1225- if (chemistry) then
1225+ if (chemistry) then
12261226 do j = 1 , buff_size
12271227 do i = chemxb, chemxe
12281228 Ys_in(i - chemxb + 1 ) = q_prim_vf(i)%sf(k, l, p)
0 commit comments