Debug SamplingMonitor

Author: Himpu Marbona

Solver/src/MultiphaseSolver/SpatialDiscretization.f90

@@ -190,13 +190,14 @@ subroutine Initialize_SpaceAndTimeMethods(controlVariables, mesh)
                error stop 
 
             end select
-			
+   
+            use_non_constant_speed_of_sound = controlVariables % ContainsKey(FLUID1_COMPRESSIBILITY_KEY)
+
             call CHDiscretization % Construct(controlVariables, ELLIPTIC_CH)
             call CHDiscretization % Describe
 
 			if ( .not. MPI_Process % isRoot ) return
 
-            use_non_constant_speed_of_sound = controlVariables % ContainsKey(FLUID1_COMPRESSIBILITY_KEY)
             if(use_non_constant_speed_of_sound) then
                write(STD_OUT,'(A)') "  Implementing artificial compressibility with a non-constant speed of sound in each phase"
             else

Solver/src/libs/monitors/Samplings.f90

@@ -20,9 +20,9 @@ module Samplings
 !  
    type Sampling_t
       character(len=LINE_LENGTH)           :: solution_file
-      integer                              :: no_of_surfaceSamplings
-	  integer  							   :: no_of_planeSamplings
-	  integer  							   :: no_of_spatialMeanNodes
+      integer                              :: no_of_surfaceSamplings =0
+	  integer  							   :: no_of_planeSamplings   =0
+	  integer  							   :: no_of_spatialMeanNodes =0
       integer                              :: dt_restriction
       logical                              :: write_dt_restriction
       class(SurfaceSampling_t), allocatable :: surfaceSamplings(:)

Solver/src/libs/timeintegrator/pAdaptationClass.f90

@@ -363,7 +363,7 @@ subroutine OverEnriching_Initialize(this,oID)
 !        -----------------------------------------------         
          write(STD_OUT,'(/)')
          call SubSection_Header("Initialize Adaptation based on Overenriching Box")
-		    write(STD_OUT,'(30X,A,A27,A20)') "->" , "ID: " , this % ID
+		    write(STD_OUT,'(30X,A,A27,I4)') "->" , "ID: " , this % ID
             write(STD_OUT,'(30X,A,A27,A20)') "->" , "x span: " , x_span
 			write(STD_OUT,'(30X,A,A27,A20)') "->" , "y span: " , y_span
 			write(STD_OUT,'(30X,A,A27,A20)') "->" , "z span: " , z_span

Solver/src/libs/timeintegrator/pAdaptationClassRL.f90

@@ -341,12 +341,12 @@ subroutine pAdaptation_Construct(this, controlVariables, t0, mesh)
             case ("rho")
                this % acoustic_variable = 8
             case default
-               WRITE(STD_OUT,*) 'Not recognized acoustic variable. Using pressure by default.'
                this % acoustic_variable = 7
+			   if ( MPI_Process % isRoot ) WRITE(STD_OUT,*) 'Undefined acoustic variable. Using pressure by default.'
             end select
          else
-            WRITE(STD_OUT,*) 'Undefined acoustic variable. Using pressure by default.'
             this % acoustic_variable = 7
+			if ( MPI_Process % isRoot ) WRITE(STD_OUT,*) 'Undefined acoustic variable. Using pressure by default.'
          end if
 
          call mesh % DefineAcousticElements(this % observer, this % acoustic_sources, this % acoustic_distance, surfacesMesh % zones)