@@ -369,9 +369,9 @@ contains
369
369
do k = 0 , n
370
370
do j = 0 , m
371
371
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
375
375
end do
376
376
end do
377
377
end do
@@ -386,14 +386,14 @@ contains
386
386
do j = 0 , m
387
387
do q = 1 , nnode
388
388
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)
392
392
393
393
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)
397
397
end do
398
398
end do
399
399
end do
0 commit comments