code clean, removing some comments and print statements

Author: sbryngelson

src/simulation_code/m_bubbles.f90

@@ -174,7 +174,7 @@ SUBROUTINE s_compute_bubble_source(idir, q_prim_vf, q_cons_vf, mydivu, &
 
                     bub_v_src(q,j,k,l) = nbub(j,k,l) * rddot
 
-                    IF (alf < 1.d-10) THEN
+                    IF (alf < 1.d-11) THEN
                         bub_adv_src(j,k,l) = 0d0
                         bub_r_src(q,j,k,l) = 0d0
                         bub_v_src(q,j,k,l) = 0d0

src/simulation_code/m_data_output.f90

@@ -2185,7 +2185,7 @@ SUBROUTINE s_write_probe_files(t_step,q_cons_vf,accel_mag) ! -----------
                         tmp = ptot
                         CALL s_mpi_allreduce_sum(tmp,ptot)
 
-                        IF (QBMM) THEN
+                        IF (qbmm) THEN
                             tmp = varR
                             CALL s_mpi_allreduce_sum(tmp,varR)
                             tmp = varV
@@ -2279,8 +2279,7 @@ SUBROUTINE s_write_probe_files(t_step,q_cons_vf,accel_mag) ! -----------
                         END IF
                     ELSEIF (p == 0) THEN
                         IF (bubbles) THEN
-                            WRITE(i+30,'(6X,F12.6,F24.8,F24.8,F24.8,F24.8,' // &
-                                           'F24.8,F24.8,F24.8,F24.8)') &
+                            WRITE(i+30,'(6X,10F24.8)') &
                                 nondim_time, &
                                 rho, &
                                 vel(1), &

src/simulation_code/m_global_parameters.f90

@@ -1154,13 +1154,13 @@ SUBROUTINE s_comp_n_from_cons( vftmp,nRtmp,ntmp )
                         ! nRtmp(i) = small_alf
                     ! END IF
                 ! END DO
-                ! STOP 'nR3 is negative'
+                STOP 'nR3 is negative'
                 nR3 = 1.d-12
             END IF
             IF (vftmp < 0d0) THEN
                 PRINT*, vftmp, nR3, nRtmp(:)
                 ! vftmp = small_alf
-                ! STOP 'vf negative'
+                STOP 'vf negative'
                 ntmp = DSQRT( (4.d0*pi/3.d0)*nR3/1.d-12 )
             ELSE
                 ntmp = DSQRT( (4.d0*pi/3.d0)*nR3/vftmp )

src/simulation_code/m_qbmm.f90

@@ -182,19 +182,6 @@ SUBROUTINE s_coeff( pres,rho,c,coeffs )
                         IF (  Web /=dflt_real) coeffs(6,i1,i2) = 2*i2/Web/rho
                     ELSE IF (bubble_model==2) THEN
                         ! KM with approximation of 1/(1-V/C) = 1+V/C
-                        ! coeffs(1,i1,i2)  = -3d0*i2/2d0
-                        ! coeffs(2,i1,i2)  = -i2/c
-                        ! coeffs(3,i1,i2)  = i2/(2d0*c*c)
-                        ! coeffs(4,i1,i2)  = -i2*pres/rho
-                        ! coeffs(5,i1,i2)  = -2d0*i2*pres/(c*rho)
-                        ! coeffs(6,i1,i2)  = -i2*pres/(c*c*rho)
-                        ! coeffs(7,i1,i2)  = i2/rho
-                        ! coeffs(8,i1,i2)  = 2d0*i2/(c*rho)
-                        ! coeffs(9,i1,i2)  = i2/(c*c*rho)
-                        ! coeffs(10,i1,i2) = -3d0*i2*gam/(c*rho)
-                        ! coeffs(11,i1,i2) = -3d0*i2*gam/(c*c*rho)
-                        ! coeffs(12,i1,i2) = i1
-
                         coeffs(1,i1,i2)  = -3d0*i2/2d0
                         coeffs(2,i1,i2)  = -i2/c
                         coeffs(3,i1,i2)  = i2/(2d0*c*c)
@@ -207,7 +194,6 @@ SUBROUTINE s_coeff( pres,rho,c,coeffs )
                         coeffs(10,i1,i2) = -3d0*i2*gam/(c*rho)
                         coeffs(11,i1,i2) = -3d0*i2*gam/(c*c*rho)
                         coeffs(12,i1,i2) = i1
-
                     END IF
                 END IF
             END DO; END DO

src/simulation_code/m_rhs.f90

@@ -3739,9 +3739,8 @@ SUBROUTINE s_get_monopole(idir,q_prim_vf,t_step,mymono) ! ----------------------
 
             IF (idir == ndirs) THEN
                 mytime = t_step*dt
-                print*, 'time', mytime, 'delay', mymono%delay
+                if (proc_rank == 0) print*, 'time', mytime, 'delay', mymono%delay, dflt_real
                 IF ( (mytime < mymono%delay) .AND. mymono%delay /= dflt_real ) RETURN
-                PRINT*, 'getting mono'
 
                 DO j = 0,m; DO k = 0,n; DO l=0,p
                     CALL s_convert_to_mixture_variables( q_prim_vf, myRho, n_tait, B_tait, Re, We, j, k, l )