comin_descrdata_types.F90

Go to the documentation of this file.

!> @file comin_descrdata_types.F90
!! @brief ComIn descriptive data structures.
!
!  @authors 10/2021 :: ICON Community Interface  <comin@icon-model.org>
!
!  SPDX-License-Identifier: BSD-3-Clause
!
!  See LICENSES for license information.
!  Where software is supplied by third parties, it is indicated in the
!  headers of the routines.
!
MODULE comin_descrdata_types
  USE iso_c_binding,         ONLY : c_signed_char, c_int, c_bool
  USE comin_setup_constants, ONLY : wp
  IMPLICIT NONE
 
  PRIVATE
 
  PUBLIC :: t_comin_descrdata_global_data
  PUBLIC :: t_comin_descrdata_domain_data
  PUBLIC :: t_comin_descrdata_domain_cells_data
  PUBLIC :: t_comin_descrdata_domain_edges_data
  PUBLIC :: t_comin_descrdata_domain_verts_data
  PUBLIC :: t_comin_descrdata_simulation_interval
  PUBLIC :: comin_glb2loc_index_lookup_fct
 
#include "comin_global.inc"
 
  ! ------------------------------------
  ! data types for descriptive data structures
  ! ------------------------------------
 
  ! Descriptive data from p_patch apart from few exceptions (e.g. lon and lat)
  ! handed to ComIn as pointer, all other descriptive data are copies from
  ! host model, which need to be updated during the model simulation if changed
 
  !> Global data is invariant wrt the computational grid and never changed
  !! or updated.
  !! @ingroup fortran_interface
  TYPE :: t_comin_descrdata_global_data
    !> @brief Number of logical domains
    !!
    !! Fortran indexing. Starting with 1.
    INTEGER(C_INT)        :: n_dom
    !> @brief Maximum number of model domains
    INTEGER(C_INT)        :: max_dom
    !> @brief Loop chunk length
    !!
    !! Often used as the first dimension of variables
    INTEGER(C_INT)        :: nproma
    !> @brief Working precision
    !!
    !! Fortan KIND VALUE (REAL)
    INTEGER(C_INT)        :: wp
    !> @brief Cell row ordering (block index): halo cells have a negative numbering ranging between `min_rlcell` and `min_rlcell_int-1`.
    INTEGER(C_INT)        :: min_rlcell_int
    !> @brief Cell row ordering (block index): halo cells have a negative numbering ranging between `min_rlcell` and `min_rlcell_int-1`.
    INTEGER(C_INT)        :: min_rlcell
    !> @brief Cell row ordering (block index): cell rows in the lateral domain boundary have positive `refin_c_ctrl` numbers.
    !!
    !! In total, `max_rlcell = 5` cell rows are part of this section.
    INTEGER(C_INT)        :: max_rlcell
    !> @brief Vertex row ordering (block index). See the corresponding documentation for cells.
    INTEGER(C_INT)        :: min_rlvert_int
    !> @brief Vertex row ordering (block index). See the corresponding documentation for cells.
    INTEGER(C_INT)        :: min_rlvert
    !> @brief Vertex row ordering (block index). See the corresponding documentation for cells.
    INTEGER(C_INT)        :: max_rlvert
    !> @brief Vertex row ordering (block index). See the corresponding documentation for cells.
    INTEGER(C_INT)        :: min_rledge_int
    !> @brief Edge row ordering (block index). See the corresponding documentation for cells.
    INTEGER(C_INT)        :: min_rledge
    !> @brief Edge row ordering (block index). See the corresponding documentation for cells.
    INTEGER(C_INT)        :: max_rledge
    !> @brief Cell row ordering (block index): width of lateral boundary region, not including the nudging zone
    INTEGER(C_INT)        :: grf_bdywidth_c
    !> @brief Edge row ordering (block index): See the corresponding documentation for cells.
    INTEGER(C_INT)        :: grf_bdywidth_e
    !> @brief Whether this is a restarted run
    LOGICAL(C_BOOL)       :: lrestartrun
    !> @brief Parameter A of the vertical coordinate (without influence of topography)
    !!
    !! index1=1,nlev+1
    REAL(wp), ALLOCATABLE :: vct_a(:)
    !> @brief Limited area setup where forcing comes from the sides
    LOGICAL(C_BOOL) :: l_limited_area
    !> @brief The yac instance id used by the host model
    !!
    !! See comments in the technical documentation and the examples
    INTEGER(C_INT)        :: yac_instance_id
    !> @brief Host model version information - Git remote URL
    CHARACTER(LEN=:), ALLOCATABLE :: host_git_remote_url
    !> @brief Host model version information - Git branch name
    CHARACTER(LEN=:), ALLOCATABLE :: host_git_branch
    !> @brief Host model version information - Git tag name
    CHARACTER(LEN=:), ALLOCATABLE :: host_git_tag
    !> @brief Host model version information - revision version number
    CHARACTER(LEN=:), ALLOCATABLE :: host_revision
    !> @brief Information about the computation device (accelerator)
    !!
    !! Whether the host model uses an acceleraor (e.g. GPU). In this case the
    !! memory can be accessed directly on the device. See technical
    !! documentation and examples for further details.
    LOGICAL(C_BOOL)               :: has_device
    !> @brief Name of the accelerator device
    CHARACTER(LEN=:), ALLOCATABLE :: device_name
    !> @brief Vendor of the accelerator device
    CHARACTER(LEN=:), ALLOCATABLE :: device_vendor
    !> @brief Driver of the accelerator device
    CHARACTER(LEN=:), ALLOCATABLE :: device_driver
  END TYPE t_comin_descrdata_global_data
 
  !> Cell information for grid data structures
  !! @ingroup fortran_interface
  TYPE :: t_comin_descrdata_domain_cells_data
    !> @brief Number of local cells
    INTEGER(C_INT) :: ncells
    !> @brief Number of global cells
    INTEGER(C_INT) :: ncells_global
    !> @brief Number of blocks for cells
    !! @see nproma
    INTEGER(C_INT) :: nblks
    !> @brief Maximum number of connected neighbors (usually 3)
    INTEGER(C_INT)          :: max_connectivity
    !> @brief Number of edges connected to cell
    !!
    !! index1=1,nproma, index2=1,nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: num_edges(:,:) => null()
    !> @brief Lateral boundary distance indices
    !!
    !! index1=1,nproma, index2=1,nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS             :: refin_ctrl(:,:) => null()
    !> @brief List of start indices for each refin_ctrl level
    !!
    !! index1=min_rlcell,max_rlcell (defined in mo_impl_constants)
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: start_index(:) => null()
    !> @brief List of end indices for each refin_ctrl level
    !!
    !! index1=min_rlcell,max_rlcell
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: end_index(:) => null()
    !> @brief List of start block for each refin_ctrl level
    !!
    !! index1=min_rlcell,max_rlcell
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: start_block(:) => null()
    !> @brief List of end block for each refin_ctrl level
    !!
    !! index1=min_rlcell,max_rlcell
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: end_block(:) => null()
    !> @brief Domain ID of child triangles
    !!
    !! index1=1,nproma, index2=1,nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: child_id(:,:) => null()
    !> @brief Line index of parent triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: parent_glb_idx(:,:) => null()
    !> @brief Block index of parent triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: parent_glb_blk(:,:) => null()
    !> @brief Line indices of verts of triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c, index3=1,3
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: vertex_idx(:,:,:) => null()
    !> @brief Block indices of verts of triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c, index3=1,3
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: vertex_blk(:,:,:) => null()
    !> @brief Line indices of neighbors of triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c, index3=1,3
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: neighbor_idx(:,:,:) => null()
    !> @brief Block indices of neighbors of triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c, index3=1,3
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: neighbor_blk(:,:,:) => null()
    !> @brief Line indices of edges of triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c, index3=1,3
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: edge_idx(:,:,:) => null()
    !> @brief Block indices of edges of triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c, index3=1,3
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: edge_blk(:,:,:) => null()
    !> @brief Longitude of centers of triangular cells
    !!
    !! index1=1,nproma, index2=1,nblks_c
    REAL(wp), ALLOCATABLE     :: clon(:,:)
    !> @brief Latitude of centers of triangular cells
    !!
    !! index1=1,nproma, index2=1,nblks_c
    REAL(wp), ALLOCATABLE     :: clat(:,:)
    !> @brief Area of triangle
    !!
    !! index1=1,nproma, index2=1,nblks_c
    REAL(wp), POINTER, CONTIGUOUS         :: area(:,:) => null()
    !> @brief Geometrical height of half levels at cell centre
    !!
    !! index1=1,nproma, index2=1,nlev+1, index3=1,nblks_c
    REAL(wp), ALLOCATABLE                 :: hhl(:,:,:)
    !> @brief Coefficient for bilinear interpolation from edges to cells for scalar quantities
    !! index1=nproma, index2=3, index3=nblks_c
    REAL(wp), POINTER, CONTIGUOUS :: e_bln_c_s(:,:,:)
    !> @brief Factor for divergence
    !! index1=nproma , index2=cell_type, index3=nblks_c
    REAL(wp), POINTER, CONTIGUOUS :: geofac_div(:,:,:)
    !> @brief Factor for nabla2-scalar
    !! index1=nproma, index2=cell_type+1, index3=nblks_c
    REAL(wp), POINTER, CONTIGUOUS :: geofac_n2s(:,:,:)
    !> @brief Factor for Green-Gauss gradient
    !! index1=nproma, index2=4, index3=nblks_c, index4=2
    REAL(wp), POINTER, CONTIGUOUS :: geofac_grg(:,:,:,:)
    !> @brief Indices of the stencil of surrounding edges for vector rbf interpolation at each
    !! cell center
    !! index1=rbf_vec_dim_c(=9), index2=nproma, inedx3=nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS :: rbf_vec_idx(:,:,:)
    !> @brief Block-Indices of the stencil of surrounding edges for vector rbf
    !> interpolation at each
    !! cell center
    !! index1=rbf_vec_dim_c(=9), index2=nproma, inedx3=nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS :: rbf_vec_blk(:,:,:)
    !> @brief Number of surrounding edges in the stencil for vector rbf
    !! interpolation at each cell center
    !! index1=nproma, index2=nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS :: rbf_vec_stencil(:,:)
    !> @brief coefficients used for vector rbf interpolation at
    !! each cell center
    !! index1=rbf_vec_dim_c(=9), index2=2, INdex3=nproma, index4=nblks_c
    REAL(wp), POINTER, CONTIGUOUS :: rbf_vec_coeff(:,:,:,:)
    !> @brief Global cell indices
    !!
    !! index1=1,ncells
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: glb_index(:) => null()
    !> @brief Domain decomposition flag
    !!
    !! * decomp_domain==0: inner domain
    !! * decomp_domain>0: boundary
    !! * decomp_domain<0: undefined
    !!
    !! For cells:
    !! * 0=owned
    !! * 1=shared edge with owned cell
    !! * 2=shared vertex with owned cell
    !!
    !! index1=1,nproma, index2=1,nblks_c
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: decomp_domain(:,:) => null()
  END TYPE t_comin_descrdata_domain_cells_data
 
  !> Vertex information for grid data structures
  !! @ingroup fortran_interface
  TYPE :: t_comin_descrdata_domain_verts_data
    !> @brief Number of local verts
    INTEGER(C_INT) :: nverts
    !> @brief Number of global verts
    INTEGER(C_INT) :: nverts_global
    !> @brief Number of blocks for verts
    !!
    !! @see nproma
    INTEGER(C_INT) :: nblks
    !> @brief Number of edges connected to vertex
    !!
    !! index1=1,nproma, index2=1,nblks_v
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: num_edges(:,:) => null()
    !> @brief Lateral boundary distance indices
    !!
    !! index1=1,nproma, index2=1,nblks_v
    INTEGER(C_INT), POINTER, CONTIGUOUS             :: refin_ctrl(:,:) => null()
    !> @brief List of start indices for each refin_ctrl level
    !!
    !! index1=min_rlvert,max_rlvert
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: start_index(:) => null()
    !> @brief List of end indices for each refin_ctrl level
    !!
    !! index1=min_rlvert,max_rlvert
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: end_index(:) => null()
    !> @brief List of start block for each refin_ctrl level
    !!
    !! index1=min_rlvert,max_rlvert
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: start_block(:) => null()
    !> @brief List of end block for each refin_ctrl level
    !!
    !! index1=min_rlvert,max_rlvert
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: end_block(:) => null()
    !> @brief Block indices of neighbor vertices
    !!
    !! index1=1,nproma, index2=1,nblks_v, index3=1,6
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: neighbor_blk(:,:,:) => null()
    !> @brief Line indices of neighbor vertices
    !!
    !! index1=1,nproma, index2=1,nblks_v, index3=1,6
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: neighbor_idx(:,:,:) => null()
    !> @brief Line indices of cells around each vertex
    !!
    !! index1=1,nproma, index2=1,nblks_v, index3=1,6
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: cell_idx(:,:,:) => null()
    !> @brief Block indices of cells around each vertex
    !!
    !! index1=1,nproma, index2=1,nblks_v, index3=1,6
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: cell_blk(:,:,:) => null()
    !> @brief Line indices of edges around a vertex
    !!
    !! index1=1,nproma, index2=1,nblks_v, index3=1,6
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: edge_idx(:,:,:) => null()
    !> @brief Block indices of edges around a vertex
    !!
    !! index1=1,nproma, index2=1,nblks_v, index3=1,6
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: edge_blk(:,:,:) => null()
    !> @brief Longitude of vertex
    !!
    !! index1=1,nproma, index2=1,nblks_v
    REAL(wp), ALLOCATABLE     :: vlon(:,:)
    !> @brief Latitude of vertex
    !!
    !! index1=1,nproma, index2=1,nblks_v
    REAL(wp), ALLOCATABLE     :: vlat(:,:)
    !> Indices of the stencil of surrounding edges for vector rbf interpolation
    !! at each triangle vertex
    !! index1=rbf_vec_dim_v(=6), index2=nproma, index3=nblks_v
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: rbf_vec_idx(:,:,:)
    !> Block-indices of the stencil of surrounding edges for vector rbf
    !! interpolation at each triangle vertex
    !! index1=rbf_vec_dim_v(=6), index2=nproma, index3=nblks_v
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: rbf_vec_blk(:,:,:)
    !> @brief Number of surrounding edges in the stencil for vector rbf
    !! interpolation at each triangle vertex
    !! index1=nproma, index2=nblks_v
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: rbf_vec_stencil(:,:)
    !> @brief Coefficients used for vector rbf interpolation at each tringle vertex.
    !! input is normal component
    !! index1=rbf_vec_dim_v(=6), index2=2, index3=nproma, index4=nblks_v
    REAL(wp), POINTER, CONTIGUOUS :: rbf_vec_coeff(:,:,:,:)
    !> @brief Global vertex indices
    !!
    !! index1=1,nverts
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: glb_index(:) => null()
    !> @brief Domain decomposition flag for vertices
    !!
    !! * decomp_domain==0: inner domain
    !! * decomp_domain>0: boundary
    !! * decomp_domain<0: undefined
    !!
    !! For vertices:
    !! * 0=owned
    !! * 1=on owned cell
    !! * 2=on level 1 cell
    !!
    !! index1=nproma, index2=1,nblks_v
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: decomp_domain(:,:) => null()
  END TYPE t_comin_descrdata_domain_verts_data
 
  !> Edge information for grid data structures
  !! @ingroup fortran_interface
  TYPE :: t_comin_descrdata_domain_edges_data
    !> @brief Number of local edges
    INTEGER(C_INT) :: nedges
    !> @brief Number of global edges
    INTEGER(C_INT) :: nedges_global
    !> @brief Number of blocks for edges
    INTEGER(C_INT) :: nblks
    !> @brief Lateral boundary distance indices
    !!
    !! index1=1,nproma, index2=1,nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS             :: refin_ctrl(:,:) => null()
    !> @brief List of start indices for each refin_ctrl level
    !!
    !! index1=min_rledge,max_rledge
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: start_index(:) => null()
    !> @brief List of end indices for each refin_ctrl level
    !!
    !! index1=min_rledge,max_rledge
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: end_index(:) => null()
    !> @brief List of start block for each refin_ctrl level
    !!
    !! index1=min_rledge,max_rledge
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: start_block(:) => null()
    !> @brief List of end block for each refin_ctrl level
    !!
    !! index1=min_rledge,max_rledge
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: end_block(:) => null()
    !> @brief Domain ID of child edges
    !!
    !! index1=1,nproma, index2=1,nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: child_id(:,:) => null()
    !> @brief Line index of parent edges
    !!
    !! index1=1,nproma, index2=1,nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: parent_glb_idx(:,:) => null()
    !> @brief block index of parent edges
    !!
    !! index1=1,nproma, index2=1,nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: parent_glb_blk(:,:) => null()
    !> @brief Line indices of adjacent cells
    !!
    !! index1=1,nproma, index2=1,nblks_e, index3=1,2
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: cell_idx(:,:,:) => null()
    !> @brief Block indices of adjacent cells
    !!
    !! index1=1,nproma, index2=1,nblks_e, index3=1,2
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: cell_blk(:,:,:) => null()
    !> @brief Line indices of edge vertices
    !!
    !! vertex indices 3 and 4 are the non-edge-aligned vertices of cells 1 and 2
    !!
    !! index1=1,nproma, index2=1,nblks_e, index3=1,4
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: vertex_idx(:,:,:) => null()
    !> @brief Block indices of edge vertices
    !!
    !! index1=1,nproma, index2=1,nblks_e, index3=1,4
    INTEGER(C_INT), POINTER, CONTIGUOUS          :: vertex_blk(:,:,:) => null()
    !> @brief Longitude of edge midpoint
    !!
    !! index=1,nproma, index2=1,nblks_e
    REAL(wp), ALLOCATABLE     :: elon(:,:)
    !> @brief Latitude of edge midpoint
    !!
    !! index=1,nproma, index2=1,nblks_e
    REAL(wp), ALLOCATABLE     :: elat(:,:)
    !> @brief Orientation of the edge
    !! =1 if vector product of vector from vertex1 to vertex 2 (v2-v1) by vector
    !! from cell c1 to cell c2 (c2-c1) goes outside the sphere
    !! =-1 if vector product ...       goes inside  the sphere
    !! index=1,nproma, index2=1,nblks_e
    REAL(wp), POINTER, CONTIGUOUS                :: tangent_orientation(:,:) => null()
    !> @brief Nudging coefficient for cells
    !! index1=nproma, index2=nblks_e
    REAL(wp), POINTER, CONTIGUOUS :: nudgecoeff_e(:,:)
    !> @brief Indices for stencil of surrounding edges for vector rbf
    !! interpolation at each triangle edge
    !! index1=rbf_vec_dim_e(=4), index2=nproma, index3=nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS :: rbf_vec_idx(:,:,:)
    !> @brief Block-indices for stencil of surrounding edges for vector rbf
    !! interpolation at each triangle edge
    !! index1=rbf_vec_dim_e(=4), index2=nproma, index3=nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS :: rbf_vec_blk(:,:,:)
    !> @brief Number of surrounding edges in the stencil for vector rbf
    !! interpolation at each triangle edge
    !! index1=nproma, index2=nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS :: rbf_vec_stencil(:,:)
    !> @brief Coefficients used for rbf inter- polation of the tangential velocity component
    !! (from the surrounding normals) at each triangle edge
    !! index1=rbf_vec_dim_e(=4), index2=nproma, index3=nblks_e
    REAL(wp), POINTER, CONTIGUOUS :: rbf_vec_coeff(:,:,:)
    ! @brief Coefficient for interpolation from adjacent cells onto edge
    ! index1=nproma, index2=2, index3=nblks_e
    REAL(wp), POINTER, CONTIGUOUS :: c_lin_e(:,:,:)
    !> @brief Global edge indices
    !!
    !! index1=1,nedges
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: glb_index(:) => null()
    !> @brief Domain decomposition flag for edges
    !!
    !! * decomp_domain==0: inner domain
    !! * decomp_domain>0: boundary
    !! * decomp_domain<0: undefined
    !!
    !! For edges:
    !! * 0=owned
    !! * 1=on owned cell
    !! * 2=exactly one shared vertex with owned cell
    !! index1=nproma, index2=1,nblks_e
    INTEGER(C_INT), POINTER, CONTIGUOUS  :: decomp_domain(:,:) => null()
  END TYPE t_comin_descrdata_domain_edges_data
 
  !> Patch grid data structure, gathering information on grids
  !! @ingroup fortran_interface
  TYPE :: t_comin_descrdata_domain_data
    !> @brief Horizontal grid filename
    CHARACTER(LEN=:), POINTER    :: grid_filename => null()
    !> @brief Alphanumerical hash of grid
    INTEGER(c_signed_char), POINTER, CONTIGUOUS     :: grid_uuid(:) => null()
    !> @brief Number of grid used (GRIB2 key)
    INTEGER(C_INT)                      :: number_of_grid_used
    !> @brief Domain ID of current domain
    INTEGER(C_INT)             :: id
    !> @brief Domain id of parent
    INTEGER(C_INT)             :: parent_id
    !> @brief IDs for all child domains
    !!
    !! index=1,max_dom (usually 10) but only has meaningful data in 1,n_dom
    INTEGER(C_INT), POINTER             :: child_id(:) => null()
    !> @brief Actual number of child domains
    INTEGER(C_INT)             :: n_childdom
    !> @brief Time at which execution of a (nested) model domain starts
    REAL(wp)                     :: dom_start
    !> @brief Time at which execution of a (nested) model domain terminates
    REAL(wp)                     :: dom_end
    !> @brief Number of vertical model levels
    INTEGER(C_INT)             :: nlev
    !> @brief Half level of parent domain that coincides with upper margin of current domain
    INTEGER(C_INT)             :: nshift
    !> @brief Total shift of model top w.r.t. global domain
    INTEGER(C_INT)             :: nshift_total
 
    TYPE(t_comin_descrdata_domain_cells_data)    :: cells
    TYPE(t_comin_descrdata_domain_verts_data)    :: verts
    TYPE(t_comin_descrdata_domain_edges_data)    :: edges
  END TYPE t_comin_descrdata_domain_data
 
  !> Simulation status information, sim_current contains current time step
  !! @ingroup fortran_interface
  TYPE :: t_comin_descrdata_simulation_interval
    !> @brief Start datetime of experiment
    CHARACTER(LEN=COMIN_MAX_DATETIME_STR_LEN) :: exp_start
    !> @brief Stop datetime of experiment
    CHARACTER(LEN=COMIN_MAX_DATETIME_STR_LEN) :: exp_stop
    !> @brief Start datetime of current run
    CHARACTER(LEN=COMIN_MAX_DATETIME_STR_LEN) :: run_start
    !> @brief Stop datetime of current run
    CHARACTER(LEN=COMIN_MAX_DATETIME_STR_LEN) :: run_stop
  END TYPE t_comin_descrdata_simulation_interval
 
  INTERFACE
    FUNCTION comin_glb2loc_index_lookup_fct(jg, glb) RESULT(loc) BIND(C)
      IMPORT c_int
      INTEGER(C_INT), VALUE :: jg
      INTEGER(C_INT), VALUE :: glb
      INTEGER(C_INT) :: loc
    END FUNCTION comin_glb2loc_index_lookup_fct
  END INTERFACE
 
END MODULE comin_descrdata_types