@@ -369,9 +369,9 @@ contains
369369 do k = 0 , n
370370 do j = 0 , m
371371 q_cons_ts(index)%vf(i)%sf(j, k, l) = &
372- (scaler1 * q_cons_ts(1 )%vf(i)%sf(j, k, l) &
373- + scaler2 * q_cons_ts(2 )%vf(i)%sf(j, k, l) &
374- + scaler3 * dt * rhs_vf(i)%sf(j, k, l)) / (scaler1 + scaler2) !! TODO :: scaler1 + scaler2 should be called a normalization constant
372+ (scaler1* q_cons_ts(1 )%vf(i)%sf(j, k, l) &
373+ + scaler2* q_cons_ts(2 )%vf(i)%sf(j, k, l) &
374+ + scaler3* dt * rhs_vf(i)%sf(j, k, l))/ (scaler1 + scaler2) !! TODO :: scaler1 + scaler2 should be called a normalization constant
375375 end do
376376 end do
377377 end do
@@ -386,14 +386,14 @@ contains
386386 do j = 0 , m
387387 do q = 1 , nnode
388388 pb_ts(index)%sf(j, k, l, q, i) = &
389- (scaler1 * pb_ts(1 )%sf(j, k, l, q, i) &
390- + scaler2 * pb_ts(2 )%sf(j, k, l, q, i) &
391- + scaler3 * dt * rhs_pb(j, k, l, q, i)) / (scaler1 + scaler2)
389+ (scaler1* pb_ts(1 )%sf(j, k, l, q, i) &
390+ + scaler2* pb_ts(2 )%sf(j, k, l, q, i) &
391+ + scaler3* dt * rhs_pb(j, k, l, q, i))/ (scaler1 + scaler2)
392392
393393 mv_ts(index)%sf(j, k, l, q, i) = &
394- (scaler1 * mv_ts(1 )%sf(j, k, l, q, i) &
395- + scaler2 * mv_ts(2 )%sf(j, k, l, q, i) &
396- + scaler3 * dt * rhs_mv(j, k, l, q, i)) / (scaler1 + scaler2)
394+ (scaler1* mv_ts(1 )%sf(j, k, l, q, i) &
395+ + scaler2* mv_ts(2 )%sf(j, k, l, q, i) &
396+ + scaler3* dt * rhs_mv(j, k, l, q, i))/ (scaler1 + scaler2)
397397 end do
398398 end do
399399 end do
0 commit comments