#define IF_NOTOK_RETURN(action) if (status/=0) then; write (gol,'("in ",a)') rname; call goErr; action; return; end if
#define IF_ERROR_RETURN(action) if (status> 0) then; write (gol,'("in ",a)') rname; call goErr; action; return; end if
!
#include "tm5.inc"
!
!-----------------------------------------------------------------------------
!                    TM5                                                     !
!-----------------------------------------------------------------------------
!BOP
!
! !MODULE:     M7_DATA
!
! !DESCRIPTION: holds m7 data and methods to allocate (INIT_M7_DATA), and
!               deallocate them (FREE_M7_DATA).
!\\
!\\
! !INTERFACE: 
!
MODULE M7_DATA
  !
  ! !USES:
  !
  USE GO,           ONLY : gol, goErr, goPr
  USE mo_aero_m7,   ONLY : nmod, nsol
  USE tm5_distgrid, ONLY : Get_DistGrid, dgrid
  USE chem_param,   ONLY : d3_data
  USE dims,         ONLY : im,jm,lm,nregions

  IMPLICIT NONE
  !
  ! !PUBLIC DATA MEMBERS:
  !
  INTEGER, PARAMETER :: nm7procs = 88

  TYPE(d3_data), DIMENSION(nregions, nmod), TARGET :: rw_mode
  TYPE(d3_data), DIMENSION(nregions, nmod), TARGET :: dens_mode
  TYPE(d3_data), DIMENSION(nregions, nsol), TARGET :: h2o_mode, rwd_mode
  !
  ! !REVISION HISTORY: 
  !   20 Jun 2012 - P. Le Sager - updated to Narcissa/Achim version : i.e.
  !                               added rwd_mode
  !                             - protex; all public by default
  !                             - adapted for lon-lat MPI domain decomposition
  !
  ! !REMARKS:
  !
  !EOP
  !------------------------------------------------------------------------

CONTAINS

  !--------------------------------------------------------------------------
  !                    TM5                                                  !
  !--------------------------------------------------------------------------
  !BOP
  !
  ! !IROUTINE:   INIT_M7_DATA
  !
  ! !DESCRIPTION: 
  !\\
  !\\
  ! !INTERFACE:
  !
  SUBROUTINE INIT_M7_DATA( status )
    !
    ! !OUTPUT PARAMETERS:
    !
    INTEGER, INTENT(out)   :: status
    !
    ! !REVISION HISTORY: 
    !   20 Jun 2012 - P. Le Sager - adapted for lon-lat MPI domain decomposition
    !
    !EOP
    !------------------------------------------------------------------------
    !BOC

    INTEGER   :: region, imode, i1, i2, j1, j2, lmr
    INTEGER   :: l_halo

    l_halo = 1 

    DO region=1,nregions

       call Get_DistGrid( dgrid(region), I_STRT=i1, I_STOP=i2, J_STRT=j1, J_STOP=j2 )
       lmr = lm(region)

       DO imode=1,nmod

          ALLOCATE(rw_mode  (region,imode)%d3(i1:i2, j1:j2, lmr))
          ALLOCATE(dens_mode(region,imode)%d3(i1:i2, j1:j2, lmr))

          rw_mode  (region,imode)%d3 = 0.1e-6  ! m
          dens_mode(region,imode)%d3 = 1800.0  ! kg/m3

          rw_mode  (region,imode)%halo = 0
          dens_mode(region,imode)%halo = 0
          
       ENDDO

       DO imode=1,nsol

          ALLOCATE(rwd_mode(region,imode)%d3(i1:i2, j1:j2, lmr))
          ALLOCATE(h2o_mode(region,imode)%d3(i1-l_halo:i2+l_halo, j1-l_halo:j2+l_halo, lmr))
          
          rwd_mode(region,imode)%d3 = 0.1e-6
          h2o_mode(region,imode)%d3 = 0.0
          
          rwd_mode(region,imode)%halo = 0
          h2o_mode(region,imode)%halo = l_halo
          
       ENDDO
       
    ENDDO
    ! ok
    status = 0
  END SUBROUTINE INIT_M7_DATA
  !EOC

  !--------------------------------------------------------------------------
  !                    TM5                                                  !
  !--------------------------------------------------------------------------
  !BOP
  !
  ! !IROUTINE:    FREE_M7_DATA
  !
  ! !DESCRIPTION: 
  !\\
  !\\
  ! !INTERFACE:
  !
  SUBROUTINE FREE_M7_DATA
    !
    ! !REVISION HISTORY: 
    !   20 Jun 2012 - P. Le Sager - no reference to MPI anymore
    !
    !EOP
    !------------------------------------------------------------------------
    !BOC

    INTEGER   :: region, imode

    DO region=1,nregions
       
       DO imode=1,nmod
          DEALLOCATE(rw_mode(region,imode)%d3)
          DEALLOCATE(dens_mode(region,imode)%d3)
       ENDDO

       DO imode=1,nsol
          DEALLOCATE(h2o_mode(region,imode)%d3)
          DEALLOCATE(rwd_mode(region,imode)%d3)
       ENDDO
       
    ENDDO
    
  END SUBROUTINE FREE_M7_DATA
  !EOC

END MODULE M7_DATA