Debug SamplingMonitor

Author: Himpu Marbona

Solver/src/MultiphaseSolver/SpatialDiscretization.f90

@@ -190,17 +190,18 @@ subroutine Initialize_SpaceAndTimeMethods(controlVariables, mesh)
                error stop 
 
             end select
-   
+			
+            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
-               write(STD_OUT,'(A)') "  Implementing artificial compressibility with a constant speed of sound in each phase"
+               write(STD_OUT,'(A)') "  Implementing artificial compressibility with a constant ACM factor in each phase"
             endif
-
-            call CHDiscretization % Construct(controlVariables, ELLIPTIC_CH)
-            call CHDiscretization % Describe
-
 
          end if
 

Solver/src/libs/monitors/Samplings.f90

@@ -267,7 +267,7 @@ subroutine Sampling_Destruct (self)
 
       end subroutine
 
-      elemental subroutine Sampling_Assign ( to, from )
+      impure elemental subroutine Sampling_Assign ( to, from )
          implicit none
          !-arguments--------------------------------------
          class(Sampling_t), intent(inout)  :: to
@@ -290,6 +290,7 @@ elemental subroutine Sampling_Assign ( to, from )
          safedeallocate ( to % surfaceSamplings )
 		 safedeallocate ( to % planeSamplings )
 		 safedeallocate ( to % spatialMeanNodes )
+		 
          allocate ( to % surfaceSamplings ( size(from % surfaceSamplings) ) )
 		 allocate ( to % planeSamplings   ( size(from % planeSamplings  ) ) )
 		 allocate ( to % spatialMeanNodes ( size(from % spatialMeanNodes  ) ) )