module farray_alloc
!
! Dynamical allocation of f and df.
!
  real, dimension(:,:,:,:), allocatable :: f,df

  public :: f,df
  public :: initialize, finalize

  contains
!******************************************************************************
  subroutine initialize

  use Cdata
  use Messages, only: fatal_error

  integer :: stat

    !allocate( f(mx,my,mz,nvar+naux+nscratch+nglobal),STAT=stat)
    !if (stat>0) call fatal_error('farray_alloc','Could not allocate memory for f')
    !allocate(df(mx,my,mz,nvar),STAT=stat)
    !if (stat>0) call fatal_error('farray_alloc','Could not allocate memory for df')
    !!print*, 'stat,mfarry,nfarray=',stat,mvar,maux,mscratch,mglobal,nvar,naux,nscratch,nglobal

    !mvar=nvar; maux=naux; maux_com=naux_com; mscratch=nscratch; mglobal=nglobal

    allocate(f(mx,my,mz,mfarray),STAT=stat)
    if (stat>0) call fatal_error('farray_alloc','Could not allocate memory for f')
    if (nt>0.and..not.lgpu) then
      allocate(df(mx,my,mz,mvar),STAT=stat)
      if (stat>0) call fatal_error('farray_alloc','Could not allocate memory for df')
    endif

  endsubroutine initialize
!******************************************************************************
  subroutine finalize

    deallocate(f) 
    if (allocated(df)) deallocate(df) 

  endsubroutine finalize
!******************************************************************************
  endmodule