New subroutine for storing past element material number on cloud

src/store_ematnump_on_cloud.f90

+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+!                                                                              |
+!              ||===\\                                                         | 
+!              ||    \\                                                        |
+!              ||     ||   //==\\   ||  ||   //==||  ||/==\\                   |
+!              ||     ||  ||    ||  ||  ||  ||   ||  ||    ||                  |
+!              ||    //   ||    ||  ||  ||  ||   ||  ||                        |
+!              ||===//     \\==//    \\==\\  \\==\\  ||                        |
+!                                                                              |
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+!                                                                              |
+!             STORE EMATNUMP ON CLOUD    OCT. 2012                         |
+!                                                                              |
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+
+subroutine store_ematnump_on_cloud(cl,osolve)
+
+use definitions
+
+!------------------------------------------------------------------------------|
+!(((((((((((((((( Purpose of the routine  ))))))))))))))))))))))))))))))))))))))
+!------------------------------------------------------------------------------|
+! This routine loops over all of the cloud particles and stores the past       |
+! material number for each element on the cloud                                |
+!------------------------------------------------------------------------------|
+!((((((((((((((((  declaration of the subroutine arguments  ))))))))))))))))))))
+!------------------------------------------------------------------------------|
+
+implicit none
+
+include 'mpif.h'
+
+type (octreesolve) osolve
+type (cloud) cl
+
+!------------------------------------------------------------------------------|
+!(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
+!------------------------------------------------------------------------------|
+
+double precision :: x0,y0,z0,dxyz
+integer :: i,iproc,nproc,ierr,leaf,level,locc
+
+!-------------------------------------------------------------------------------
+!-------------------------------------------------------------------------------
+
+call mpi_comm_size (mpi_comm_world,nproc,ierr)
+call mpi_comm_rank (mpi_comm_world,iproc,ierr)
+
+! Loop over all cloud particles, find associated element and store past element mat number
+do i = 1,cl%np
+  call octree_find_leaf (osolve%octree,osolve%noctree,cl%x(i),cl%y(i),cl%z(i),  &
+                        leaf,level,locc,x0,y0,z0,dxyz)
+  cl%ematnum(i)=osolve%matnum(leaf)
+enddo
+
+end subroutine store_ematnump_on_cloud
+
+!-------------------------------------------------------------------------------
+!-------------------------------------------------------------------------------