write_field_p.F90

Go to the documentation of this file.

module write_field_p
implicit none
  
  interface writefield_p
    module procedure write_field3d_p,write_field2d_p,write_field1d_p
  end interface writefield_p
  
  contains
  
  subroutine write_field1d_p(name,Field)
    USE parallel
    USE write_field
    implicit none
  
    integer, parameter :: maxdim=1
    character(len=*)   :: name
    real, dimension(:) :: field
    real, dimension(:),allocatable :: new_field
    integer, dimension(MaxDim) :: dim
    
    
    dim=shape(field)
    allocate(new_field(dim(1)))
    new_field(:)=field(:)
    call gather_field(new_field,dim(1),1,0)
    
    if (mpi_rank==0) call writefield(name,new_field)
    
    end subroutine write_field1d_p

  subroutine write_field2d_p(name,Field)
    USE parallel
    USE write_field
    implicit none
  
    integer, parameter :: maxdim=2
    character(len=*)   :: name
    real, dimension(:,:) :: field
    real, dimension(:,:),allocatable :: new_field
    integer, dimension(MaxDim) :: dim
    
    dim=shape(field)
    allocate(new_field(dim(1),dim(2)))
    new_field(:,:)=field(:,:)
    call gather_field(new_field(1,1),dim(1)*dim(2),1,0)
    
    if (mpi_rank==0) call writefield(name,new_field)
    
     
  end subroutine write_field2d_p
  
  subroutine write_field3d_p(name,Field)
    USE parallel
    USE write_field
    implicit none
  
    integer, parameter :: maxdim=3
    character(len=*)   :: name
    real, dimension(:,:,:) :: field
    real, dimension(:,:,:),allocatable :: new_field
    integer, dimension(MaxDim) :: dim
    
    dim=shape(field)
    allocate(new_field(dim(1),dim(2),dim(3)))
    new_field(:,:,:)=field(:,:,:)
    call gather_field(new_field(1,1,1),dim(1)*dim(2),dim(3),0)
    
   if (mpi_rank==0) call writefield(name,new_field)
    
  end subroutine write_field3d_p  

end module write_field_p