clean

Author: sbryngelson

src/simulation/inline_riemann.fpp

@@ -46,74 +46,3 @@ if (avg_state == 2) then
 end if
 
 #:enddef compute_average_state
-
-! #:def compute_wave_speeds()
-
-! if (wave_speeds == 1) then
-!     s_L = min(vel_L(dir_idx(1)) - c_L, vel_R(dir_idx(1)) - c_R)
-!     s_R = max(vel_R(dir_idx(1)) + c_R, vel_L(dir_idx(1)) + c_L)
-
-!     s_S = (pres_R - pres_L + rho_L*vel_L(dir_idx(1))* &
-!            (s_L - vel_L(dir_idx(1))) - &
-!            rho_R*vel_R(dir_idx(1))* &
-!            (s_R - vel_R(dir_idx(1)))) &
-!           /(rho_L*(s_L - vel_L(dir_idx(1))) - &
-!             rho_R*(s_R - vel_R(dir_idx(1))))
-! elseif (wave_speeds == 2) then
-!     pres_SL = 5d-1*(pres_L + pres_R + rho_avg*c_avg* &
-!                     (vel_L(dir_idx(1)) - &
-!                      vel_R(dir_idx(1))))
-
-!     pres_SR = pres_SL
-
-!     Ms_L = max(1d0, sqrt(1d0 + ((5d-1 + gamma_L)/(1d0 + gamma_L))* &
-!                          (pres_SL/pres_L - 1d0)*pres_L/ &
-!                          ((pres_L + pi_inf_L/(1d0 + gamma_L)))))
-!     Ms_R = max(1d0, sqrt(1d0 + ((5d-1 + gamma_R)/(1d0 + gamma_R))* &
-!                          (pres_SR/pres_R - 1d0)*pres_R/ &
-!                          ((pres_R + pi_inf_R/(1d0 + gamma_R)))))
-
-!     s_L = vel_L(dir_idx(1)) - c_L*Ms_L
-!     s_R = vel_R(dir_idx(1)) + c_R*Ms_R
-
-!     s_S = 5d-1*((vel_L(dir_idx(1)) + vel_R(dir_idx(1))) + &
-!                 (pres_L - pres_R)/ &
-!                 (rho_avg*c_avg))
-! end if
-
-
-!         if (wave_speeds == 1) then
-!             s_L = min(vel_L(dir_idx(1)) - c_L, vel_R(dir_idx(1)) - c_R)
-!             s_R = max(vel_R(dir_idx(1)) + c_R, vel_L(dir_idx(1)) + c_L)
-
-!             s_S = (pres_R - pres_L + rho_L*vel_L(dir_idx(1))* &
-!                    (s_L - vel_L(dir_idx(1))) - &
-!                    rho_R*vel_R(dir_idx(1))* &
-!                    (s_R - vel_R(dir_idx(1)))) &
-!                   /(rho_L*(s_L - vel_L(dir_idx(1))) - &
-!                     rho_R*(s_R - vel_R(dir_idx(1))))
-!         elseif (wave_speeds == 2) then
-!             pres_SL = 5d-1*(pres_L + pres_R + rho_avg*c_avg* &
-!                             (vel_L(dir_idx(1)) - &
-!                              vel_R(dir_idx(1))))
-
-!             pres_SR = pres_SL
-
-!             Ms_L = max(1d0, sqrt(1d0 + ((5d-1 + gamma_L)/(1d0 + gamma_L))* &
-!                                  (pres_SL/pres_L - 1d0)*pres_L/ &
-!                                  ((pres_L + pi_inf_L/(1d0 + gamma_L)))))
-!             Ms_R = max(1d0, sqrt(1d0 + ((5d-1 + gamma_R)/(1d0 + gamma_R))* &
-!                                  (pres_SR/pres_R - 1d0)*pres_R/ &
-!                                  ((pres_R + pi_inf_R/(1d0 + gamma_R)))))
-
-!             s_L = vel_L(dir_idx(1)) - c_L*Ms_L
-!             s_R = vel_R(dir_idx(1)) + c_R*Ms_R
-
-!             s_S = 5d-1*((vel_L(dir_idx(1)) + vel_R(dir_idx(1))) + &
-!                         (pres_L - pres_R)/ &
-!                         (rho_avg*c_avg))
-!         end if
-
-! #:enddef compute_wave_speeds
-
-