mod_atmos_mesh.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!> module Atmosphere / Mesh
!!
!! @par Description
!!          Base module for mesh with atmospheric model
!!
!! @author Yuta kawai, Team SCALE
!!
!<
!-------------------------------------------------------------------------------
#include "scaleFElib.h"
module mod_atmos_mesh
  !-----------------------------------------------------------------------------
  !
  !++ Used modules
  !
  use scale_precision
  use scale_io
  use scale_prc
 
  use scale_meshfield_base, only: meshfield3d
  use scale_mesh_base2d, only: meshbase2d
  use scale_mesh_base3d, only: meshbase3d
  use scale_element_base, only: elementbase3d
  use scale_element_line, only: lineelement
  use scale_element_hexahedral, only: hexahedralelement
  use scale_localmesh_2d, only: localmesh2d
  use scale_localmesh_3d, only: localmesh3d
  use scale_sparsemat, only: sparsemat
  use scale_meshfield_base, only: meshfield2d
  use scale_mesh_topography, only: meshtopography
 
  use scale_file_restart_meshfield, only: file_restart_meshfield_component
  use scale_model_var_manager, only: modelvarmanager
  use scale_model_mesh_manager, only: modelmesh3d
  
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public type & procedures
  !
 
  !> Derived type to manage a computational mesh (base class)
  !!
  type, abstract, extends(modelmesh3d), public :: atmosmesh
    type(hexahedralelement) :: element !< Object to manage 3D reference element
    type(lineelement) :: element_v1d   !< Object to manage 1D reference element for the vertical direction
 
    type(meshtopography) :: topography  !< Object to manage topography
    integer :: vcoord_type_id           !< ID of vertical coordinate type
 
    logical :: comm_use_mpi_pc              !< Flag whether persistent communication in MPI is used
    logical :: comm_use_mpi_pc_fujitsu_ext  !< Flag whether Fujitsu extension in MPI persistent communication is used
  contains
    procedure :: atmosmesh_init
    procedure :: atmosmesh_final
    procedure(atmosmesh_setup_restartfile1), public, deferred :: setup_restartfile1
    procedure(atmosmesh_setup_restartfile2), public, deferred :: setup_restartfile2
    procedure(atmosmesh_calc_uvmet), public, deferred :: calc_uvmet
    generic :: setup_restartfile => setup_restartfile1, setup_restartfile2
    procedure(atmosmesh_setup_vcoord), public, deferred :: setup_vcoordinate
  end type atmosmesh
 
  interface
    subroutine atmosmesh_setup_restartfile1( this, restart_file, var_num )
      import atmosmesh
      import file_restart_meshfield_component
      class(atmosmesh), target, intent(inout) :: this
      class(file_restart_meshfield_component), intent(inout) :: restart_file
      integer, intent(in) :: var_num  
    end subroutine atmosmesh_setup_restartfile1
  end interface
  interface
    subroutine atmosmesh_setup_restartfile2( this, restart_file, &
      in_basename, in_postfix_timelabel,                         &
      out_basename, out_postfix_timelabel,                       &
      out_dtype, out_title, var_num                              )
      import atmosmesh
      import file_restart_meshfield_component
      class(atmosmesh), target, intent(inout) :: this
      class(file_restart_meshfield_component), intent(inout) :: restart_file
      character(*), intent(in) :: in_basename
      logical, intent(in) :: in_postfix_timelabel
      character(*), intent(in) :: out_basename
      logical, intent(in) :: out_postfix_timelabel
      character(*), intent(in) :: out_title
      character(*), intent(in) :: out_dtype  
      integer, intent(in) :: var_num  
    end subroutine atmosmesh_setup_restartfile2
  end interface
  interface
    subroutine atmosmesh_calc_uvmet( this, U, V, &
        Umet, Vmet )
        import atmosmesh
        import meshfield3d
        class(atmosmesh), target, intent(in) :: this
        type(meshfield3d), intent(in) :: U
        type(meshfield3d), intent(in) :: V
        type(meshfield3d), intent(inout) :: Umet
        type(meshfield3d), intent(inout) :: Vmet
    end subroutine atmosmesh_calc_uvmet
  end interface
  interface 
    subroutine atmosmesh_setup_vcoord( this )
      import atmosmesh
      class(atmosmesh), target, intent(inout) :: this
    end subroutine atmosmesh_setup_vcoord
  end interface
  integer, parameter, public :: atm_mesh_max_commnuicator_num = 10
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private procedures
  !
  !-------------------
 
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
 
contains
 
!> Setup an object to manage a computational mesh
  subroutine atmosmesh_init( this, mesh )
 
    use scale_file_monitor_meshfield, only: &
      file_monitor_meshfield_set_dim
    
    implicit none
    class(atmosmesh), target, intent(inout) :: this
    class(meshbase3d), intent(in) :: mesh            !< Object to manage 3D computational mesh
 
    class(meshbase2d), pointer :: mesh2d
    !-------------------------------------------
 
    call this%ModelMesh3D_Init( mesh )
 
    !-
    call file_monitor_meshfield_set_dim( mesh, 'ATM3D' )
    
    !-
    call mesh%GetMesh2D( mesh2d )
    call this%topography%Init( "topo", mesh2d )
 
    return
  end subroutine atmosmesh_init
 
!> Finalize an object to manage a computational mesh
  subroutine atmosmesh_final(this)
    implicit none
 
    class(atmosmesh), intent(inout) :: this
    !-------------------------------------------
 
    call this%topography%Final()
    call this%ModelMesh3D_Final()
 
    call this%element%Final()
    call this%element_v1D%Final()
 
    return
  end subroutine atmosmesh_final
 
end module mod_atmos_mesh