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_lmdz
    USE lmdz_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_lmdz
    USE lmdz_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_lmdz
    USE lmdz_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