PROGRAM conv_mesh

  implicit none

  include 'mpif.h'

  integer, save :: nprocs, my_rank, ierr, world, lastproc

  logical, save :: master

  INTEGER, PARAMETER :: infl=666

  INTEGER :: fldidx, &
       o_nG1, o_nG2dnG3nL1, &
       n_nG1, n_nG2dnG3nL1, &
       reclength, i

  REAL, DIMENSION(:), ALLOCATABLE, TARGET :: field

  REAL :: o_dx, o_AX

  CHARACTER (len=6), PARAMETER :: datarep='native'

  INTEGER :: fh, amode, etype, filetype, disp
  INTEGER, DIMENSION(MPI_STATUS_SIZE) :: status
  INTEGER(kind=MPI_OFFSET_KIND) :: offset

  INTEGER :: count

  character (len=128) :: flname

  integer :: nG1, nG2, nG3, nG2d

  !> START CODE ----------------------------------------------------------
  !

  call MPI_INIT(ierr)
  world=MPI_COMM_WORLD

  call MPI_COMM_RANK(world, my_rank, ierr)
  call MPI_COMM_SIZE(world, nprocs, ierr)

  ! NO TIME FOR CODING FOR PARALLEL

  IF(nprocs.GT.1) STOP

  ! NO TIME FOR CODING FOR PARALLEL

  if(my_rank.eq.0) master = .true.
  lastproc = nprocs-1

  ! NO TIME FOR CODING IN PARALLEL

  o_nG1 = 1250
  nG2 = 1250
  nG3 = 120
  nG2d = nG2 +2
  
  o_nG2dnG3nL1 = nG2d*nG3*o_nG1
  n_nG2dnG3nL1 = nG2d*nG3*n_nG1

  if(master) print '("Allocs")'

  ALLOCATE(field(o_nG2dnG3nL1), stat=ierr)
  if(ierr.gt.0) PRINT *, "bad alloc"

  do i=1,o_nG2dnG3nL1
     field(i) = 69.
  enddo

  ! WRITE STUFF

  flname = 'x_mpi.neu'

  inquire (iolength=reclength) field
  print '("reclength ",I0)', reclength
  print '("size(field) ",I0)', size(field)

  call system('du -sh *neu ')

  if(master)then
     call MPI_File_delete(trim(flname), MPI_INFO_NULL, ierr)
     print *, "ierr ", ierr
  endif

  call system('du -sh *neu ')

  etype = MPI_REAL
  filetype = MPI_REAL

  amode = ior(MPI_MODE_CREATE,MPI_MODE_WRONLY)

  CALL MPI_file_OPEN(world, flname, &
       amode, MPI_INFO_NULL, fh, ierr)

  call system('du -sh *neu ')
     
  disp = 0
  offset = 0
  
  CALL MPI_File_set_view(fh, disp, etype, filetype, &
       trim(datarep),  MPI_INFO_NULL, ierr)

  print *, "ierr ", ierr
  print *, "o_nG2dnG3nL1 ", o_nG2dnG3nL1
  print *, "dim() ", size(field)

  CALL MPI_file_write_at_all(fh, field, o_nG2dnG3nL1, &
       etype, status, ierr)
  
  print '("Wrote ")'
  print *, "ierr ", ierr

  call MPI_Get_Count(status, etype, count, ierr)
  PRINT *, "<COUNT ", count
  
  call MPI_File_get_size(fh, count, ierr);
  PRINT *, "<SIZE ", count
     
  CALL MPI_Barrier(world, ierr)
  CALL MPI_File_sync(fh, ierr)
  
  CALL MPI_file_CLOSE(fh, ierr)

  call system('du -sh *neu ')
  
  CALL MPI_Barrier(world, ierr)



  flname='x_normal.neu'
  open(infl,file=flname,access='stream')
  write(infl) field
  close(infl)


  DEALLOCATE(field)

  !     --- SYNC and out ----------------------------------------------- /
  call MPI_Barrier(world, ierr)
  call MPI_Finalize(ierr)


END PROGRAM conv_mesh