yac_module.F90

Go to the documentation of this file.

! Copyright (c) 2024 The YAC Authors
!
! SPDX-License-Identifier: BSD-3-Clause
 
#ifdef HAVE_CONFIG_H
#include "config.h"
#endif
 
module yac
  use, intrinsic :: iso_c_binding, only : yac_int_fc_kind, &
                                          yac_mpi_fint_fc_kind
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  public
 
  integer, parameter :: yac_max_charlen = 132
  integer, parameter :: yac_mpi_fint_kind = yac_mpi_fint_fc_kind
  integer, parameter :: yac_int_kind = yac_int_fc_kind
  double precision, parameter :: yac_frac_mask_no_value = 133713371337.0d0
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  enum, bind(c)
   enumerator :: yac_location_cell   = 0
   enumerator :: yac_location_corner = 1
   enumerator :: yac_location_edge   = 2
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  enum, bind(c)
   enumerator :: yac_exchange_type_none   = 0
   enumerator :: yac_exchange_type_source = 1
   enumerator :: yac_exchange_type_target = 2
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  enum, bind(c)
   enumerator :: yac_action_none               = 0
   enumerator :: yac_action_reduction          = 1
   enumerator :: yac_action_coupling           = 2
   enumerator :: yac_action_get_for_restart    = 4
   enumerator :: yac_action_put_for_restart    = 5
   enumerator :: yac_action_out_of_bound       = 8
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
  enum, bind(c)
   enumerator :: yac_time_unit_millisecond = 0
   enumerator :: yac_time_unit_second = 1
   enumerator :: yac_time_unit_minute = 2
   enumerator :: yac_time_unit_hour = 3
   enumerator :: yac_time_unit_day = 4
   enumerator :: yac_time_unit_month = 5
   enumerator :: yac_time_unit_year = 6
   enumerator :: yac_time_unit_iso_format = 7
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  enum, bind(c)
   enumerator :: yac_reduction_time_none       = 0
   enumerator :: yac_reduction_time_accumulate = 1
   enumerator :: yac_reduction_time_average    = 2
   enumerator :: yac_reduction_time_minimum    = 3
   enumerator :: yac_reduction_time_maximum    = 4
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
  enum, bind(c)
   enumerator :: yac_calendar_not_set = 0
   enumerator :: yac_proleptic_gregorian = 1
   enumerator :: yac_year_of_365_days = 2
   enumerator :: yac_year_of_360_days = 3
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  enum, bind(c)
     enumerator :: yac_avg_arithmetic = 0
     enumerator :: yac_avg_dist       = 1
     enumerator :: yac_avg_bary       = 2
  end enum
 
  enum, bind(c)
     enumerator :: yac_ncc_avg  = 0
     enumerator :: yac_ncc_dist = 1
  end enum
 
  enum, bind(c)
     enumerator :: yac_nnn_avg    = 0
     enumerator :: yac_nnn_dist   = 1
     enumerator :: yac_nnn_gauss  = 2
     enumerator :: yac_nnn_rbf    = 3
  end enum
 
  enum, bind(c)
     enumerator :: yac_dnn_avg   = 0
     enumerator :: yac_dnn_dist  = 1
     enumerator :: yac_dnn_gauss = 2
     enumerator :: yac_dnn_rbf   = 3
  end enum
 
  enum, bind(c)
     enumerator :: yac_dnn_search_distance_fixed     = 0
     enumerator :: yac_dnn_search_distance_cell_area = 1
  end enum
 
  enum, bind(c)
     enumerator :: yac_conserv_destarea = 0
     enumerator :: yac_conserv_fracarea = 1
  end enum
 
  enum, bind(c)
     enumerator :: yac_spmap_avg  = 0
     enumerator :: yac_spmap_dist = 1
  end enum
 
  enum, bind(c)
     enumerator :: yac_spmap_none        = 0
     enumerator :: yac_spmap_srcarea     = 1
     enumerator :: yac_spmap_invtgtarea  = 2
     enumerator :: yac_spmap_fracarea    = 3
  end enum
 
  enum, bind(c)
    enumerator :: yac_file_missing_error = 0
    enumerator :: yac_file_missing_cont  = 1
  end enum
 
  enum, bind(c)
    enumerator :: yac_file_success_stop = 0
    enumerator :: yac_file_success_cont = 1
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  enum, bind(c)
     enumerator :: yac_config_output_format_yaml = 0
     enumerator :: yac_config_output_format_json = 1
  end enum
 
  enum, bind(c)
     enumerator :: yac_config_output_sync_loc_def_comp = 0
     enumerator :: yac_config_output_sync_loc_sync_def = 1
     enumerator :: yac_config_output_sync_loc_enddef   = 2
  end enum
 
  !----------------------------------
  ! Handling of existing weight files
  !----------------------------------
 
  enum, bind(c)
   enumerator :: yac_wgt_on_existing_error     = 0
   enumerator :: yac_wgt_on_existing_keep      = 1
   enumerator :: yac_wgt_on_existing_overwrite = 2
  end enum
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  integer :: yac_yaml_emitter_default_f
  integer :: yac_yaml_emitter_json_f
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  type yac_real_ptr
    real, pointer :: p(:)
  end type yac_real_ptr
 
  type yac_dble_ptr
    double precision, pointer :: p(:)
 end type yac_dble_ptr
 
 type :: yac_string
     character(len=:), allocatable :: string
  end type  yac_string
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fmpi_handshake
    subroutine yac_fmpi_handshake ( comm, group_names, group_comms )
      import :: yac_max_charlen
      integer, intent(in) :: comm
      character(len=YAC_MAX_CHARLEN), intent(in) :: &
        group_names(:)
      integer, intent(out) :: &
        group_comms(SIZE(group_names))
    end subroutine yac_fmpi_handshake
  end interface yac_fmpi_handshake
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_finit_comm
 
     subroutine yac_finit_comm ( comm )
 
       integer, intent(in) :: comm
 
     end subroutine yac_finit_comm
 
     subroutine yac_finit_comm_instance ( comm, yac_instance_id )
 
       integer, intent(in)  :: comm
       integer, intent(out) :: yac_instance_id
 
     end subroutine yac_finit_comm_instance
 
  end interface yac_finit_comm
 
  interface yac_finit
 
     subroutine yac_finit ( )
 
     end subroutine yac_finit
 
     subroutine yac_finit_instance ( yac_instance_id)
 
       integer, intent(out) :: yac_instance_id
 
     end subroutine yac_finit_instance
 
  end interface yac_finit
 
  interface yac_finit_dummy
     subroutine yac_finit_dummy ()
     end subroutine yac_finit_dummy
 
     subroutine yac_finit_comm_dummy ( comm )
 
       integer, intent(in) :: comm
 
     end subroutine yac_finit_comm_dummy
  end interface yac_finit_dummy
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
    function yac_fget_default_instance_id()
 
      integer :: yac_fget_default_instance_id
 
    end function yac_fget_default_instance_id
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
    function yac_fdefault_instance_defined()
      logical :: yac_fdefault_instance_defined
    end function yac_fdefault_instance_defined
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fread_config_yaml
     subroutine yac_fread_config_yaml (yaml_filename)
       character(len=*), intent(in) :: yaml_filename
     end subroutine yac_fread_config_yaml
 
     subroutine yac_fread_config_yaml_instance (yac_instance_id, yaml_filename)
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: yaml_filename
     end subroutine yac_fread_config_yaml_instance
  end interface yac_fread_config_yaml
 
  interface yac_fread_config_json
     subroutine yac_fread_config_json (json_filename)
       character(len=*), intent(in) :: json_filename
     end subroutine yac_fread_config_json
 
     subroutine yac_fread_config_json_instance (yac_instance_id, json_filename)
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: json_filename
     end subroutine yac_fread_config_json_instance
  end interface yac_fread_config_json
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fset_config_output_file
     subroutine yac_fset_config_output_file ( &
       filename, fileformat, sync_location, include_definitions)
       character(len=*), intent(in) :: filename
       integer, intent(in) :: fileformat
       integer, intent(in) :: sync_location
       logical, intent(in), optional :: include_definitions
     end subroutine yac_fset_config_output_file
 
     subroutine yac_fset_config_output_file_instance ( &
       yac_instance_id, filename, fileformat, sync_location, &
       include_definitions)
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: filename
       integer, intent(in) :: fileformat
       integer, intent(in) :: sync_location
       logical, intent(in), optional :: include_definitions
     end subroutine yac_fset_config_output_file_instance
  end interface yac_fset_config_output_file
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fset_grid_output_file
     subroutine yac_fset_grid_output_file ( &
       gridname, filename)
       character(len=*), intent(in) :: gridname
       character(len=*), intent(in) :: filename
     end subroutine yac_fset_grid_output_file
 
     subroutine yac_fset_grid_output_file_instance ( &
       yac_instance_id, gridname, filename)
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: gridname
       character(len=*), intent(in) :: filename
     end subroutine yac_fset_grid_output_file_instance
  end interface yac_fset_grid_output_file
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fset_grids_output_file
     subroutine yac_fset_grids_output_file ( &
       filename)
       character(len=*), intent(in) :: filename
     end subroutine yac_fset_grids_output_file
 
     subroutine yac_fset_grids_output_file_instance ( &
       yac_instance_id, filename)
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: filename
     end subroutine yac_fset_grids_output_file_instance
  end interface yac_fset_grids_output_file
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fset_coordinates_check
     subroutine yac_fset_coordinates_check ( apply_check )
 
       logical, intent(in) :: apply_check
     end subroutine yac_fset_coordinates_check
 
     subroutine yac_fset_coordinates_check_instance ( yac_instance_id, apply_check )
       integer, intent(in) :: yac_instance_id
       logical, intent(in) :: apply_check
     end subroutine yac_fset_coordinates_check_instance
  end interface yac_fset_coordinates_check
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fset_coordinates_mismatch_is_fatal
     subroutine yac_fset_coordinates_mismatch_is_fatal ( &
       mismatch_is_fatal )
 
       logical, intent(in) :: mismatch_is_fatal
     end subroutine yac_fset_coordinates_mismatch_is_fatal
 
     subroutine yac_fset_coordinates_mismatch_is_fatal_instance ( &
       yac_instance_id, mismatch_is_fatal )
 
       integer, intent(in) :: yac_instance_id
       logical, intent(in) :: mismatch_is_fatal
     end subroutine yac_fset_coordinates_mismatch_is_fatal_instance
  end interface yac_fset_coordinates_mismatch_is_fatal
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fcleanup
 
     subroutine yac_fcleanup ()
     end subroutine yac_fcleanup
 
     subroutine yac_fcleanup_instance ( yac_instance_id )
 
      integer, intent(in) :: yac_instance_id
 
     end subroutine yac_fcleanup_instance
 
  end interface yac_fcleanup
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_ffinalize
 
     subroutine yac_ffinalize ()
     end subroutine yac_ffinalize
 
     subroutine yac_ffinalize_instance ( yac_instance_id )
 
      integer, intent(in) :: yac_instance_id
 
     end subroutine
 
  end interface yac_ffinalize
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
     function yac_fget_version () result (version_string)
     character(len=:), ALLOCATABLE :: version_string
     end function yac_fget_version
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
    function yac_fget_mpi_handshake_group_name () &
      result(mpi_handshake_group_name)
      character(len=:), ALLOCATABLE :: mpi_handshake_group_name
    end function yac_fget_mpi_handshake_group_name
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  INTERFACE yac_fpredef_comp
 
     SUBROUTINE yac_fpredef_comp ( comp_name, comp_id )
 
       character(len=*), intent(in) :: comp_name
       integer, intent(out)         :: comp_id
 
     END SUBROUTINE yac_fpredef_comp
 
     SUBROUTINE yac_fpredef_comp_instance ( yac_instance_id, comp_name, comp_id )
 
       integer, intent(in)          :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       integer, intent(out)         :: comp_id
 
     END SUBROUTINE yac_fpredef_comp_instance
 
  END INTERFACE yac_fpredef_comp
 
  interface yac_fdef_comp
 
     subroutine yac_fdef_comp ( comp_name, comp_id )
 
       character(len=*), intent(in) :: comp_name
       integer, intent(out)         :: comp_id
 
     end subroutine yac_fdef_comp
 
     subroutine yac_fdef_comp_instance ( yac_instance_id, comp_name, comp_id )
 
       integer, intent(in)          :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       integer, intent(out)         :: comp_id
 
     end subroutine yac_fdef_comp_instance
 
  end interface yac_fdef_comp
 
  interface yac_fdef_comps
 
     subroutine yac_fdef_comps ( comp_names, num_comps, comp_ids )
 
        use, intrinsic :: iso_c_binding, only : c_char
        import :: yac_max_charlen
 
        integer, intent(in)   :: num_comps
        character(kind=c_char, len=YAC_MAX_CHARLEN), intent(in) :: &
                                 comp_names(num_comps)
        integer, intent(out)  :: comp_ids(num_comps)
 
     end subroutine yac_fdef_comps
 
     subroutine yac_fdef_comps_instance ( yac_instance_id, &
                                          comp_names,      &
                                          num_comps,       &
                                          comp_ids )
 
        use, intrinsic :: iso_c_binding, only : c_char
        import :: yac_max_charlen
 
        integer, intent(in)   :: yac_instance_id
        integer, intent(in)   :: num_comps
        character(kind=c_char, len=YAC_MAX_CHARLEN), intent(in) :: &
                                 comp_names(num_comps)
        integer, intent(out)  :: comp_ids(num_comps)
 
     end subroutine yac_fdef_comps_instance
 
  end interface yac_fdef_comps
 
  interface yac_fdef_comp_dummy
 
    subroutine yac_fdef_comp_dummy ( )
    end subroutine yac_fdef_comp_dummy
 
    subroutine yac_fdef_comp_dummy_instance ( yac_instance_id )
 
      integer, intent(in) :: yac_instance_id
 
    end subroutine yac_fdef_comp_dummy_instance
 
  end interface yac_fdef_comp_dummy
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fdef_datetime
 
     subroutine yac_fdef_datetime ( start_datetime, end_datetime )
 
       character(len=*), intent(in), optional :: start_datetime
       character(len=*), intent(in), optional :: end_datetime
 
     end subroutine yac_fdef_datetime
 
     subroutine yac_fdef_datetime_instance ( yac_instance_id, &
                                             start_datetime,  &
                                             end_datetime )
 
       integer, intent(in)                    :: yac_instance_id
       character(len=*), intent(in), optional :: start_datetime
       character(len=*), intent(in), optional :: end_datetime
 
     end subroutine yac_fdef_datetime_instance
 
  end interface yac_fdef_datetime
 
  interface yac_fdef_calendar
 
     subroutine yac_fdef_calendar( calendar )
       integer, intent(in)     :: calendar
     end subroutine yac_fdef_calendar
 
  end interface yac_fdef_calendar
 
  interface yac_fget_calendar
 
    subroutine yac_fget_calendar( calendar )
      integer, intent(out) :: calendar
    end subroutine yac_fget_calendar
 
  end interface yac_fget_calendar
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fget_start_datetime
 
     function yac_fget_start_datetime () result (start_datetime_string)
     character(len=:), ALLOCATABLE :: start_datetime_string
     end function yac_fget_start_datetime
 
     function yac_fget_start_datetime_instance (yac_instance_id) &
       result(start_datetime_string)
     integer, intent(in)            :: yac_instance_id
     character(len=:), ALLOCATABLE :: start_datetime_string
     end function yac_fget_start_datetime_instance
 
  end interface yac_fget_start_datetime
 
  interface yac_fget_end_datetime
 
     function yac_fget_end_datetime () result (end_datetime_string)
     character(len=:), ALLOCATABLE :: end_datetime_string
     end function yac_fget_end_datetime
 
     function yac_fget_end_datetime_instance (yac_instance_id) &
       result(end_datetime_string)
     integer, intent(in)            :: yac_instance_id
     character(len=:), ALLOCATABLE :: end_datetime_string
     end function yac_fget_end_datetime_instance
 
  end interface yac_fget_end_datetime
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fdef_grid
 
     subroutine yac_fdef_grid_nonuniform_real ( grid_name,             &
                                                nbr_vertices,          &
                                                nbr_cells,             &
                                                nbr_connections,       &
                                                nbr_vertices_per_cell, &
                                                x_vertices,            &
                                                y_vertices,            &
                                                cell_to_vertex,        &
                                                grid_id,               &
                                                use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_connections
       integer, intent(in)           :: nbr_vertices_per_cell(nbr_cells)
       real, intent(in)              :: x_vertices(nbr_vertices)
       real, intent(in)              :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_vertex(nbr_connections)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_nonuniform_real
 
     subroutine yac_fdef_grid_nonuniform_dble ( grid_name,             &
                                                nbr_vertices,          &
                                                nbr_cells,             &
                                                nbr_connections,       &
                                                nbr_vertices_per_cell, &
                                                x_vertices,            &
                                                y_vertices,            &
                                                cell_to_vertex,        &
                                                grid_id,               &
                                                use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_connections
       integer, intent(in)           :: nbr_vertices_per_cell(nbr_cells)
       double precision, intent(in)  :: x_vertices(nbr_vertices)
       double precision, intent(in)  :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_vertex(nbr_connections)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_nonuniform_dble
 
     subroutine yac_fdef_grid_unstruct_real ( grid_name,             &
                                              nbr_vertices,          &
                                              nbr_cells,             &
                                              nbr_vertices_per_cell, &
                                              x_vertices,            &
                                              y_vertices,            &
                                              cell_to_vertex,        &
                                              grid_id,               &
                                              use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_vertices_per_cell
       real, intent(in)              :: x_vertices(nbr_vertices)
       real, intent(in)              :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_vertex( &
                                          nbr_vertices_per_cell,nbr_cells)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_unstruct_real
 
     subroutine yac_fdef_grid_unstruct_dble ( grid_name,             &
                                              nbr_vertices,          &
                                              nbr_cells,             &
                                              nbr_vertices_per_cell, &
                                              x_vertices,            &
                                              y_vertices,            &
                                              cell_to_vertex,        &
                                              grid_id,               &
                                              use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_vertices_per_cell
       double precision, intent(in)  :: x_vertices(nbr_vertices)
       double precision, intent(in)  :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_vertex( &
                                          nbr_vertices_per_cell,nbr_cells)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_unstruct_dble
 
     subroutine yac_fdef_grid_nonuniform_edge_real ( grid_name,          &
                                                     nbr_vertices,       &
                                                     nbr_cells,          &
                                                     nbr_edges,          &
                                                     nbr_connections,    &
                                                     nbr_edges_per_cell, &
                                                     x_vertices,         &
                                                     y_vertices,         &
                                                     cell_to_edge,       &
                                                     edge_to_vertex,     &
                                                     grid_id,            &
                                                     use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_edges
       integer, intent(in)           :: nbr_connections
       integer, intent(in)           :: nbr_edges_per_cell(nbr_cells)
       real, intent(in)              :: x_vertices(nbr_vertices)
       real, intent(in)              :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_edge(nbr_connections)
       integer, intent(in)           :: edge_to_vertex(2,nbr_edges)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_nonuniform_edge_real
 
     subroutine yac_fdef_grid_nonuniform_edge_dble ( grid_name,          &
                                                     nbr_vertices,       &
                                                     nbr_cells,          &
                                                     nbr_edges,          &
                                                     nbr_connections,    &
                                                     nbr_edges_per_cell, &
                                                     x_vertices,         &
                                                     y_vertices,         &
                                                     cell_to_edge,       &
                                                     edge_to_vertex,     &
                                                     grid_id,            &
                                                     use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_edges
       integer, intent(in)           :: nbr_connections
       integer, intent(in)           :: nbr_edges_per_cell(nbr_cells)
       double precision, intent(in)  :: x_vertices(nbr_vertices)
       double precision, intent(in)  :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_edge(nbr_connections)
       integer, intent(in)           :: edge_to_vertex(2,nbr_edges)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_nonuniform_edge_dble
 
     subroutine yac_fdef_grid_unstruct_edge_real ( grid_name,          &
                                                   nbr_vertices,       &
                                                   nbr_cells,          &
                                                   nbr_edges,          &
                                                   nbr_edges_per_cell, &
                                                   x_vertices,         &
                                                   y_vertices,         &
                                                   cell_to_edge,       &
                                                   edge_to_vertex,     &
                                                   grid_id,            &
                                                   use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_edges
       integer, intent(in)           :: nbr_edges_per_cell
       real, intent(in)              :: x_vertices(nbr_vertices)
       real, intent(in)              :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_edge( &
                                          nbr_edges_per_cell,nbr_cells)
       integer, intent(in)           :: edge_to_vertex(2,nbr_edges)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_unstruct_edge_real
 
     subroutine yac_fdef_grid_unstruct_edge_dble ( grid_name,          &
                                                   nbr_vertices,       &
                                                   nbr_cells,          &
                                                   nbr_edges,          &
                                                   nbr_edges_per_cell, &
                                                   x_vertices,         &
                                                   y_vertices,         &
                                                   cell_to_edge,       &
                                                   edge_to_vertex,     &
                                                   grid_id,            &
                                                   use_ll_edges)
 
       character(len=*), intent(in)  :: grid_name
       integer, intent(in)           :: nbr_vertices
       integer, intent(in)           :: nbr_cells
       integer, intent(in)           :: nbr_edges
       integer, intent(in)           :: nbr_edges_per_cell
       double precision, intent(in)  :: x_vertices(nbr_vertices)
       double precision, intent(in)  :: y_vertices(nbr_vertices)
       integer, intent(in)           :: cell_to_edge( &
                                          nbr_edges_per_cell,nbr_cells)
       integer, intent(in)           :: edge_to_vertex(2,nbr_edges)
       integer, intent(out)          :: grid_id
       logical, optional, intent(in) :: use_ll_edges
 
     end subroutine yac_fdef_grid_unstruct_edge_dble
 
     subroutine yac_fdef_grid_curve2d_real ( grid_name,             &
                                             nbr_vertices,          &
                                             cyclic,                &
                                             x_vertices,            &
                                             y_vertices,            &
                                             grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_vertices(2)
       integer, intent(in)          :: cyclic(2)
       real, intent(in)             :: &
         x_vertices(nbr_vertices(1),nbr_vertices(2))
       real, intent(in)             :: &
         y_vertices(nbr_vertices(1),nbr_vertices(2))
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_curve2d_real
 
     subroutine yac_fdef_grid_curve2d_dble ( grid_name,             &
                                             nbr_vertices,          &
                                             cyclic,                &
                                             x_vertices,            &
                                             y_vertices,            &
                                             grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_vertices(2)
       integer, intent(in)          :: cyclic(2)
       double precision, intent(in) :: &
         x_vertices(nbr_vertices(1),nbr_vertices(2))
       double precision, intent(in) :: &
         y_vertices(nbr_vertices(1),nbr_vertices(2))
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_curve2d_dble
 
     subroutine yac_fdef_grid_reg2d_real ( grid_name,             &
                                           nbr_vertices,          &
                                           cyclic,                &
                                           x_vertices,            &
                                           y_vertices,            &
                                           grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_vertices(2)
       integer, intent(in)          :: cyclic(2)
       real, intent(in)             :: x_vertices(nbr_vertices(1))
       real, intent(in)             :: y_vertices(nbr_vertices(2))
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_reg2d_real
 
     subroutine yac_fdef_grid_reg2d_dble ( grid_name,             &
                                           nbr_vertices,          &
                                           cyclic,                &
                                           x_vertices,            &
                                           y_vertices,            &
                                           grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_vertices(2)
       integer, intent(in)          :: cyclic(2)
       double precision, intent(in) :: x_vertices(nbr_vertices(1))
       double precision, intent(in) :: y_vertices(nbr_vertices(2))
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_reg2d_dble
 
     subroutine yac_fdef_grid_cloud_real ( grid_name,  &
                                           nbr_points, &
                                           x_points,   &
                                           y_points,   &
                                           grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_points
       real, intent(in)             :: x_points(nbr_points)
       real, intent(in)             :: y_points(nbr_points)
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_cloud_real
 
     subroutine yac_fdef_grid_cloud_dble ( grid_name,  &
                                           nbr_points, &
                                           x_points,   &
                                           y_points,   &
                                           grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_points
       double precision, intent(in) :: x_points(nbr_points)
       double precision, intent(in) :: y_points(nbr_points)
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_cloud_dble
 
     subroutine yac_fdef_grid_reg2d_rot_real ( grid_name,    &
                                               nbr_vertices, &
                                               cyclic,       &
                                               x_vertices,   &
                                               y_vertices,   &
                                               x_north_pole, &
                                               y_north_pole, &
                                               grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_vertices(2)
       integer, intent(in)          :: cyclic(2)
       real, intent(in)             :: x_vertices(nbr_vertices(1))
       real, intent(in)             :: y_vertices(nbr_vertices(2))
       real, intent(in)             :: x_north_pole
       real, intent(in)             :: y_north_pole
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_reg2d_rot_real
 
     subroutine yac_fdef_grid_reg2d_rot_dble ( grid_name,    &
                                               nbr_vertices, &
                                               cyclic,       &
                                               x_vertices,   &
                                               y_vertices,   &
                                               x_north_pole, &
                                               y_north_pole, &
                                               grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer, intent(in)          :: nbr_vertices(2)
       integer, intent(in)          :: cyclic(2)
       double precision, intent(in) :: x_vertices(nbr_vertices(1))
       double precision, intent(in) :: y_vertices(nbr_vertices(2))
       double precision, intent(in) :: x_north_pole
       double precision, intent(in) :: y_north_pole
       integer, intent(out)         :: grid_id
 
     end subroutine yac_fdef_grid_reg2d_rot_dble
 
  end interface yac_fdef_grid
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fdef_points
 
     !--------------------------------------------------------------
     ! Interfaces for unnamed point definitions
     !--------------------------------------------------------------
 
     subroutine yac_fdef_points_reg2d_real ( grid_id,       &
                                             nbr_points,    &
                                             location,      &
                                             x_points_real, &
                                             y_points_real, &
                                             point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       real, intent(in)     :: x_points_real(nbr_points(1))
       real, intent(in)     :: y_points_real(nbr_points(2))
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_real
 
     subroutine yac_fdef_points_reg2d_dble ( grid_id,       &
                                             nbr_points,    &
                                             location,      &
                                             x_points,      &
                                             y_points,      &
                                             point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       double precision, intent(in) :: x_points(nbr_points(1))
       double precision, intent(in) :: y_points(nbr_points(2))
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_dble
 
     subroutine yac_fdef_points_curve2d_real ( grid_id,       &
                                               nbr_points,    &
                                               location,      &
                                               x_points_real, &
                                               y_points_real, &
                                               point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       real, intent(in)     :: &
         x_points_real(nbr_points(1),nbr_points(2))
       real, intent(in)     :: &
         y_points_real(nbr_points(1),nbr_points(2))
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_curve2d_real
 
     subroutine yac_fdef_points_curve2d_dble ( grid_id,       &
                                               nbr_points,    &
                                               location,      &
                                               x_points,      &
                                               y_points,      &
                                               point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       double precision, intent(in) :: &
         x_points(nbr_points(1),nbr_points(2))
       double precision, intent(in) :: &
         y_points(nbr_points(1),nbr_points(2))
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_curve2d_dble
 
     subroutine yac_fdef_points_unstruct_real ( grid_id,       &
                                                nbr_points,    &
                                                location,      &
                                                x_points_real, &
                                                y_points_real, &
                                                point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       real, intent(in)     :: x_points_real(nbr_points)
       real, intent(in)     :: y_points_real(nbr_points)
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_unstruct_real
 
     subroutine yac_fdef_points_unstruct_dble ( grid_id,       &
                                                nbr_points,    &
                                                location,      &
                                                x_points,      &
                                                y_points,      &
                                                point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       double precision, intent(in) :: x_points(nbr_points)
       double precision, intent(in) :: y_points(nbr_points)
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_unstruct_dble
 
     subroutine yac_fdef_points_reg2d_rot_real ( grid_id,           &
                                                 nbr_points,        &
                                                 location,          &
                                                 x_points_real,     &
                                                 y_points_real,     &
                                                 x_north_pole_real, &
                                                 y_north_pole_real, &
                                                 point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       real, intent(in)     :: x_points_real(nbr_points(1))
       real, intent(in)     :: y_points_real(nbr_points(2))
       real, intent(in)     :: x_north_pole_real
       real, intent(in)     :: y_north_pole_real
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_rot_real
 
     subroutine yac_fdef_points_reg2d_rot_dble ( grid_id,      &
                                                 nbr_points,   &
                                                 location,     &
                                                 x_points,     &
                                                 y_points,     &
                                                 x_north_pole, &
                                                 y_north_pole, &
                                                 point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       double precision, intent(in) :: x_points(nbr_points(1))
       double precision, intent(in) :: y_points(nbr_points(2))
       double precision, intent(in) :: x_north_pole
       double precision, intent(in) :: y_north_pole
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_rot_dble
 
     !--------------------------------------------------------------
     ! Interfaces for named point definitions
     !--------------------------------------------------------------
 
     subroutine yac_fdef_points_reg2d_named_real ( grid_id,       &
                                                   nbr_points,    &
                                                   location,      &
                                                   x_points_real, &
                                                   y_points_real, &
                                                   name,          &
                                                   point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       real, intent(in)     :: x_points_real(nbr_points(1))
       real, intent(in)     :: y_points_real(nbr_points(2))
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_named_real
 
     subroutine yac_fdef_points_reg2d_named_dble ( grid_id,    &
                                                   nbr_points, &
                                                   location,   &
                                                   x_points,   &
                                                   y_points,   &
                                                   name,       &
                                                   point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       double precision, intent(in) :: x_points(nbr_points(1))
       double precision, intent(in) :: y_points(nbr_points(2))
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_named_dble
 
     subroutine yac_fdef_points_curve2d_named_real ( grid_id,       &
                                                     nbr_points,    &
                                                     location,      &
                                                     x_points_real, &
                                                     y_points_real, &
                                                     name,          &
                                                     point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       real, intent(in)     :: &
         x_points_real(nbr_points(1),nbr_points(2))
       real, intent(in)     :: &
         y_points_real(nbr_points(1),nbr_points(2))
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_curve2d_named_real
 
     subroutine yac_fdef_points_curve2d_named_dble ( grid_id,    &
                                                     nbr_points, &
                                                     location,   &
                                                     x_points,   &
                                                     y_points,   &
                                                     name,       &
                                                     point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       double precision, intent(in) :: &
         x_points(nbr_points(1),nbr_points(2))
       double precision, intent(in) :: &
         y_points(nbr_points(1),nbr_points(2))
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_curve2d_named_dble
 
     subroutine yac_fdef_points_unstruct_named_real ( grid_id,       &
                                                      nbr_points,    &
                                                      location,      &
                                                      x_points_real, &
                                                      y_points_real, &
                                                      name,          &
                                                      point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       real, intent(in)     :: x_points_real(nbr_points)
       real, intent(in)     :: y_points_real(nbr_points)
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_unstruct_named_real
 
     subroutine yac_fdef_points_unstruct_named_dble ( grid_id,    &
                                                      nbr_points, &
                                                      location,   &
                                                      x_points,   &
                                                      y_points,   &
                                                      name,       &
                                                      point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       double precision, intent(in) :: x_points(nbr_points)
       double precision, intent(in) :: y_points(nbr_points)
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_unstruct_named_dble
 
     subroutine yac_fdef_points_reg2d_rot_named_real ( grid_id,           &
                                                       nbr_points,        &
                                                       location,          &
                                                       x_points_real,     &
                                                       y_points_real,     &
                                                       x_north_pole_real, &
                                                       y_north_pole_real, &
                                                       name,              &
                                                       point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       real, intent(in)     :: x_points_real(nbr_points(1))
       real, intent(in)     :: y_points_real(nbr_points(2))
       real, intent(in)     :: x_north_pole_real
       real, intent(in)     :: y_north_pole_real
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_rot_named_real
 
     subroutine yac_fdef_points_reg2d_rot_named_dble ( grid_id,      &
                                                       nbr_points,   &
                                                       location,     &
                                                       x_points,     &
                                                       y_points,     &
                                                       x_north_pole, &
                                                       y_north_pole, &
                                                       name,         &
                                                       point_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points(2)
       integer, intent(in)  :: location
       double precision, intent(in) :: x_points(nbr_points(1))
       double precision, intent(in) :: y_points(nbr_points(2))
       double precision, intent(in) :: x_north_pole
       double precision, intent(in) :: y_north_pole
       character(len=*), intent(in) :: name
       integer, intent(out) :: point_id
 
     end subroutine yac_fdef_points_reg2d_rot_named_dble
 
  end interface yac_fdef_points
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
 
     subroutine yac_fset_global_index ( global_index, &
                                        location,     &
                                        grid_id )
 
       integer, intent(in) :: global_index(*)
       integer, intent(in) :: location
       integer, intent(in) :: grid_id
 
     end subroutine yac_fset_global_index
 
  end interface
 
  interface
 
     subroutine yac_fset_global_index_yac_int ( global_index, &
                                                location,     &
                                                grid_id )
 
       import :: yac_int_kind
 
       integer(kind=YAC_INT_KIND), intent(in)  :: global_index(*)
       integer, intent(in)                     :: location
       integer, intent(in)                     :: grid_id
 
     end subroutine yac_fset_global_index_yac_int
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fset_core_mask
 
     subroutine yac_fset_core_imask ( is_core,  &
                                      location, &
                                      grid_id )
 
       integer, intent(in)  :: is_core(*)
       integer, intent(in)  :: location
       integer, intent(in)  :: grid_id
 
     end subroutine yac_fset_core_imask
 
     subroutine yac_fset_core_lmask ( is_core,  &
                                      location, &
                                      grid_id )
 
       logical, intent(in)  :: is_core(*)
       integer, intent(in)  :: location
       integer, intent(in)  :: grid_id
 
     end subroutine yac_fset_core_lmask
 
  end interface yac_fset_core_mask
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fset_mask
 
     subroutine yac_fset_imask ( is_valid, points_id )
 
       integer, intent(in)  :: is_valid(*)
       integer, intent(in)  :: points_id
 
     end subroutine yac_fset_imask
 
     subroutine yac_fset_lmask ( is_valid, points_id )
 
       logical, intent(in)  :: is_valid(*)
       integer, intent(in)  :: points_id
 
     end subroutine yac_fset_lmask
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fdef_mask
 
     subroutine yac_fdef_imask ( grid_id,    &
                                 nbr_points, &
                                 location,   &
                                 is_valid,   &
                                 mask_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       integer, intent(in)  :: is_valid(*)
       integer, intent(out) :: mask_id
 
     end subroutine yac_fdef_imask
 
     subroutine yac_fdef_lmask ( grid_id,    &
                                 nbr_points, &
                                 location,   &
                                 is_valid,   &
                                 mask_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       logical, intent(in)  :: is_valid(*)
       integer, intent(out) :: mask_id
 
     end subroutine yac_fdef_lmask
 
  end interface
 
  interface yac_fdef_mask_named
 
     subroutine yac_fdef_imask_named ( grid_id,    &
                                       nbr_points, &
                                       location,   &
                                       is_valid,   &
                                       name,       &
                                       mask_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       integer, intent(in)  :: is_valid(*)
       character(len=*), intent(in) :: name
       integer, intent(out) :: mask_id
 
     end subroutine yac_fdef_imask_named
 
     subroutine yac_fdef_lmask_named ( grid_id,    &
                                       nbr_points, &
                                       location,   &
                                       is_valid,   &
                                       name,       &
                                       mask_id )
 
       integer, intent(in)  :: grid_id
       integer, intent(in)  :: nbr_points
       integer, intent(in)  :: location
       logical, intent(in)  :: is_valid(*)
       character(len=*), intent(in) :: name
       integer, intent(out) :: mask_id
 
     end subroutine yac_fdef_lmask_named
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
 
     subroutine yac_fdef_field ( field_name,     &
                                 component_id,   &
                                 point_ids,      &
                                 num_pointsets,  &
                                 collection_size,&
                                 timestep,       &
                                 time_unit,      &
                                 field_id )
 
       character(len=*), intent (in) :: field_name
       integer, intent (in)          :: component_id
       integer, intent (in)          :: point_ids(*)
       integer, intent (in)          :: num_pointsets
       integer, intent (in)          :: collection_size
       character(len=*), intent (in) :: timestep
       integer, intent (in)          :: time_unit
       integer, intent (out)         :: field_id
 
     end subroutine yac_fdef_field
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
 
     subroutine yac_fdef_field_mask ( field_name,     &
                                      component_id,   &
                                      point_ids,      &
                                      mask_ids,       &
                                      num_pointsets,  &
                                      collection_size,&
                                      timestep,       &
                                      time_unit,      &
                                      field_id )
 
       character(len=*), intent (in) :: field_name
       integer, intent (in)          :: component_id
       integer, intent (in)          :: point_ids(*)
       integer, intent (in)          :: mask_ids(*)
       integer, intent (in)          :: num_pointsets
       integer, intent (in)          :: collection_size
       character(len=*), intent (in) :: timestep
       integer, intent (in)          :: time_unit
       integer, intent (out)         :: field_id
 
     end subroutine yac_fdef_field_mask
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
 
    subroutine yac_fcheck_field_dimensions( field_id,          &
                                            collection_size,   &
                                            num_interp_fields, &
                                            interp_field_sizes )
 
      integer, intent (in) :: field_id
      integer, intent (in) :: collection_size
      integer, intent (in) :: num_interp_fields
      integer, intent (in) :: interp_field_sizes(num_interp_fields)
 
    end subroutine yac_fcheck_field_dimensions
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
 
    subroutine yac_fcheck_src_field_buffer_size( field_id,        &
                                                 collection_size, &
                                                 src_field_buffer_size )
 
      integer, intent (in) :: field_id
      integer, intent (in) :: collection_size
      integer, intent (in) :: src_field_buffer_size
 
    end subroutine yac_fcheck_src_field_buffer_size
 
    subroutine yac_fcheck_src_field_buffer_sizes( field_id,        &
                                                  num_src_fields,  &
                                                  collection_size, &
                                                  src_field_buffer_sizes )
 
      integer, intent (in) :: field_id
      integer, intent (in) :: num_src_fields
      integer, intent (in) :: collection_size
      integer, intent (in) :: src_field_buffer_sizes(num_src_fields)
                                              
 
    end subroutine yac_fcheck_src_field_buffer_sizes
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
 
    subroutine yac_fget_raw_interp_weights_data ( &
      field_id,                                   &
      frac_mask_fallback_value,                   &
      scaling_factor,                             &
      scaling_summand,                            &
      num_fixed_values,                           &
      fixed_values,                               &
      num_tgt_per_fixed_value,                    &
      tgt_idx_fixed,                              &
      num_wgt_tgt,                                &
      wgt_tgt_idx,                                &
      num_src_per_tgt,                            &
      weights,                                    &
      src_field_idx,                              &
      src_idx,                                    &
      num_src_fields,                             &
      src_field_buffer_sizes )
 
      integer, intent (in) :: field_id
      double precision, intent(out) :: frac_mask_fallback_value
      double precision, intent(out) :: scaling_factor
      double precision, intent(out) :: scaling_summand
      integer, intent(out) :: num_fixed_values
      double precision, allocatable, intent(out) :: fixed_values(:)
      integer, allocatable, intent(out) :: num_tgt_per_fixed_value(:)
      integer, allocatable, intent(out) :: tgt_idx_fixed(:)
      integer, intent(out) :: num_wgt_tgt
      integer, allocatable, intent(out) :: wgt_tgt_idx(:)
      integer, allocatable, intent(out) :: num_src_per_tgt(:)
      double precision, allocatable, intent(out) :: weights(:)
      integer, allocatable, intent(out) :: src_field_idx(:)
      integer, allocatable, intent(out) :: src_idx(:)
      integer, intent(out) :: num_src_fields
      integer, allocatable, intent(out) :: src_field_buffer_sizes(:)
    end subroutine yac_fget_raw_interp_weights_data
 
    subroutine yac_fget_raw_interp_weights_data_csr ( &
      field_id,                                       &
      frac_mask_fallback_value,                       &
      scaling_factor,                                 &
      scaling_summand,                                &
      num_fixed_values,                               &
      fixed_values,                                   &
      num_tgt_per_fixed_value,                        &
      tgt_idx_fixed,                                  &
      src_indptr,                                     &
      weights,                                        &
      src_field_idx,                                  &
      src_idx,                                        &
      num_src_fields,                                 &
      src_field_buffer_sizes )
 
      integer, intent (in) :: field_id
      double precision, intent(out) :: frac_mask_fallback_value
      double precision, intent(out) :: scaling_factor
      double precision, intent(out) :: scaling_summand
      integer, intent(out) :: num_fixed_values
      double precision, allocatable, intent(out) :: fixed_values(:)
      integer, allocatable, intent(out) :: num_tgt_per_fixed_value(:)
      integer, allocatable, intent(out) :: tgt_idx_fixed(:)
      integer, allocatable, intent(out) :: src_indptr(:)
      double precision, allocatable, intent(out) :: weights(:)
      integer, allocatable, intent(out) :: src_field_idx(:)
      integer, allocatable, intent(out) :: src_idx(:)
      integer, intent(out) :: num_src_fields
      integer, allocatable, intent(out) :: src_field_buffer_sizes(:)
    end subroutine yac_fget_raw_interp_weights_data_csr
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fput
 
     subroutine yac_fput_real ( field_id,         &
                                nbr_hor_points,   &
                                nbr_pointsets,    &
                                collection_size,  &
                                send_field,       &
                                info,             &
                                ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: nbr_pointsets
       integer, intent (in)          :: collection_size
       real, intent (in)             :: send_field(nbr_hor_points, &
                                                   nbr_pointsets,  &
                                                   collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_real
 
     subroutine yac_fput_real_ptr ( field_id,         &
                                    nbr_pointsets,    &
                                    collection_size,  &
                                    send_field,       &
                                    info,             &
                                    ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)     :: field_id
       integer, intent (in)     :: nbr_pointsets
       integer, intent (in)     :: collection_size
       type(yac_real_ptr), intent (in) :: send_field(nbr_pointsets, collection_size)
                                                   
       integer, intent (out)    :: info
       integer, intent (out)    :: ierror
 
     end subroutine yac_fput_real_ptr
 
     subroutine yac_fput_single_pointset_real ( field_id,         &
                                                nbr_hor_points,   &
                                                collection_size,  &
                                                send_field,       &
                                                info,             &
                                                ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: collection_size
       real, intent (in)             :: send_field(nbr_hor_points, collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_single_pointset_real
 
     subroutine yac_fput_dble ( field_id,         &
                                nbr_hor_points,   &
                                nbr_pointsets,    &
                                collection_size,  &
                                send_field,       &
                                info,             &
                                ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: nbr_pointsets
       integer, intent (in)          :: collection_size
       double precision, intent (in) :: send_field(nbr_hor_points, nbr_pointsets, collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_dble
 
     subroutine yac_fput_dble_ptr ( field_id,         &
                                    nbr_pointsets,    &
                                    collection_size,  &
                                    send_field,       &
                                    info,             &
                                    ierror )
 
         import :: yac_dble_ptr
 
         integer, intent (in)     :: field_id
         integer, intent (in)     :: nbr_pointsets
         integer, intent (in)     :: collection_size
         type(yac_dble_ptr), intent (in) :: send_field(nbr_pointsets, collection_size)
                                                     
         integer, intent (out)    :: info
         integer, intent (out)    :: ierror
 
     end subroutine yac_fput_dble_ptr
 
     subroutine yac_fput_single_pointset_dble ( field_id,         &
                                                nbr_hor_points,   &
                                                collection_size,  &
                                                send_field,       &
                                                info,             &
                                                ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: collection_size
       double precision, intent (in) :: send_field(nbr_hor_points, collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_single_pointset_dble
 
     subroutine yac_fput_frac_real ( field_id,         &
                                     nbr_hor_points,   &
                                     nbr_pointsets,    &
                                     collection_size,  &
                                     send_field,       &
                                     send_frac_mask,   &
                                     info,             &
                                     ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: nbr_pointsets
       integer, intent (in)          :: collection_size
       real, intent (in)             :: send_field(nbr_hor_points, &
                                                   nbr_pointsets,  &
                                                   collection_size)
                                                        
       real, intent (in)             :: send_frac_mask(nbr_hor_points, &
                                                       nbr_pointsets,  &
                                                       collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_frac_real
 
     subroutine yac_fput_frac_real_ptr ( field_id,         &
                                         nbr_pointsets,    &
                                         collection_size,  &
                                         send_field,       &
                                         send_frac_mask,   &
                                         info,             &
                                         ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)     :: field_id
       integer, intent (in)     :: nbr_pointsets
       integer, intent (in)     :: collection_size
       type(yac_real_ptr), intent (in) :: send_field(nbr_pointsets, collection_size)
                                                   
       type(yac_real_ptr), intent (in) :: send_frac_mask(nbr_pointsets, collection_size)
                                                   
       integer, intent (out)    :: info
       integer, intent (out)    :: ierror
 
     end subroutine yac_fput_frac_real_ptr
 
     subroutine yac_fput_frac_single_pointset_real ( field_id,         &
                                                     nbr_hor_points,   &
                                                     collection_size,  &
                                                     send_field,       &
                                                     send_frac_mask,   &
                                                     info,             &
                                                     ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: collection_size
       real, intent (in)             :: send_field(nbr_hor_points, collection_size)
                                                        
       real, intent (in)             :: send_frac_mask(nbr_hor_points, collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_frac_single_pointset_real
 
     subroutine yac_fput_frac_dble ( field_id,         &
                                     nbr_hor_points,   &
                                     nbr_pointsets,    &
                                     collection_size,  &
                                     send_field,       &
                                     send_frac_mask,   &
                                     info,             &
                                     ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: nbr_pointsets
       integer, intent (in)          :: collection_size
       double precision, intent (in) :: send_field(nbr_hor_points, nbr_pointsets, collection_size)
                                                        
       double precision, intent (in) :: send_frac_mask(nbr_hor_points, nbr_pointsets, collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_frac_dble
 
     subroutine yac_fput_frac_dble_ptr ( field_id,         &
                                         nbr_pointsets,    &
                                         collection_size,  &
                                         send_field,       &
                                         send_frac_mask,   &
                                         info,             &
                                         ierror )
 
         import :: yac_dble_ptr
 
         integer, intent (in)     :: field_id
         integer, intent (in)     :: nbr_pointsets
         integer, intent (in)     :: collection_size
         type(yac_dble_ptr), intent (in) :: send_field(nbr_pointsets, collection_size)
                                                     
         type(yac_dble_ptr), intent (in) :: send_frac_mask(nbr_pointsets, collection_size)
                                                     
         integer, intent (out)    :: info
         integer, intent (out)    :: ierror
 
     end subroutine yac_fput_frac_dble_ptr
 
     subroutine yac_fput_frac_single_pointset_dble ( field_id,         &
                                                     nbr_hor_points,   &
                                                     collection_size,  &
                                                     send_field,       &
                                                     send_frac_mask,   &
                                                     info,             &
                                                     ierror )
 
       integer, intent (in)          :: field_id
       integer, intent (in)          :: nbr_hor_points
       integer, intent (in)          :: collection_size
       double precision, intent (in) :: send_field(nbr_hor_points, collection_size)
                                                        
       double precision, intent (in) :: send_frac_mask(nbr_hor_points, collection_size)
                                                        
       integer, intent (out)         :: info
       integer, intent (out)         :: ierror
 
     end subroutine yac_fput_frac_single_pointset_dble
 
  end interface yac_fput
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fget
 
     subroutine yac_fget_real ( field_id,         &
                                nbr_hor_points,   &
                                collection_size,  &
                                recv_field,       &
                                info,             &
                                ierror )
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: nbr_hor_points
       integer, intent (in)  :: collection_size
       real, intent (inout)  :: recv_field(nbr_hor_points, collection_size)
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_real
 
     subroutine yac_fget_real_ptr ( field_id,         &
                                    collection_size,  &
                                    recv_field,       &
                                    info,             &
                                    ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: collection_size
       type(yac_real_ptr)      :: recv_field(collection_size)
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_real_ptr
 
     subroutine yac_fget_dble ( field_id,         &
                                nbr_hor_points,   &
                                collection_size,  &
                                recv_field,       &
                                info,             &
                                ierror )
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: nbr_hor_points
       integer, intent (in)  :: collection_size
       double precision, intent (inout) :: &
          recv_field(nbr_hor_points, collection_size)
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_dble
 
     subroutine yac_fget_dble_ptr ( field_id,         &
                                    collection_size,  &
                                    recv_field,       &
                                    info,             &
                                    ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr)    :: recv_field(collection_size)
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_dble_ptr
 
  end interface yac_fget
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fget_raw
 
     subroutine yac_fget_raw_real ( field_id,              &
                                    src_field_buffer_size, &
                                    collection_size,       &
                                    src_field_buffer,      &
                                    info,                  &
                                    ierror )
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: src_field_buffer_size
       integer, intent (in)  :: collection_size
       real, intent (out)    :: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                      
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_real
 
     subroutine yac_fget_raw_real_ptr ( field_id,         &
                                        num_src_fields,   &
                                        collection_size,  &
                                        src_field_buffer, &
                                        info,             &
                                        ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: num_src_fields
       integer, intent (in)  :: collection_size
       type(yac_real_ptr), intent(inout) :: &
         src_field_buffer(num_src_fields, collection_size)
                                                
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_real_ptr
 
     subroutine yac_fget_raw_dble ( field_id,              &
                                    src_field_buffer_size, &
                                    collection_size,       &
                                    src_field_buffer,      &
                                    info,                  &
                                    ierror )
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: src_field_buffer_size
       integer, intent (in)  :: collection_size
       double precision, intent (out) :: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                      
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_dble
 
     subroutine yac_fget_raw_dble_ptr ( field_id,         &
                                        num_src_fields,   &
                                        collection_size,  &
                                        src_field_buffer, &
                                        info,             &
                                        ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: num_src_fields
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr), intent(inout) :: &
         src_field_buffer(num_src_fields, collection_size)
                                                
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_dble_ptr
 
     subroutine yac_fget_raw_frac_real ( field_id,              &
                                         src_field_buffer_size, &
                                         collection_size,       &
                                         src_field_buffer,      &
                                         src_frac_mask_buffer,  &
                                         info,                  &
                                         ierror )
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: src_field_buffer_size
       integer, intent (in)  :: collection_size
       real, intent (out)    :: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                      
       real, intent (out)    :: &
         src_frac_mask_buffer(src_field_buffer_size, collection_size)
                                                      
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_frac_real
 
     subroutine yac_fget_raw_frac_real_ptr ( field_id,             &
                                             num_src_fields,       &
                                             collection_size,      &
                                             src_field_buffer,     &
                                             src_frac_mask_buffer, &
                                             info,                 &
                                             ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: num_src_fields
       integer, intent (in)  :: collection_size
       type(yac_real_ptr), intent(inout) :: &
         src_field_buffer(num_src_fields, collection_size)
                                                
       type(yac_real_ptr), intent(inout) :: &
         src_frac_mask_buffer(num_src_fields, collection_size)
                                                
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_frac_real_ptr
 
     subroutine yac_fget_raw_frac_dble ( field_id,              &
                                         src_field_buffer_size, &
                                         collection_size,       &
                                         src_field_buffer,      &
                                         src_frac_mask_buffer,  &
                                         info,                  &
                                         ierror )
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: src_field_buffer_size
       integer, intent (in)  :: collection_size
       double precision, intent (out) :: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                      
       double precision, intent (out) :: &
         src_frac_mask_buffer(src_field_buffer_size, collection_size)
                                                      
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_frac_dble
 
     subroutine yac_fget_raw_frac_dble_ptr ( field_id,             &
                                             num_src_fields,       &
                                             collection_size,      &
                                             src_field_buffer,     &
                                             src_frac_mask_buffer, &
                                             info,                 &
                                             ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: num_src_fields
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr), intent(inout) :: &
         src_field_buffer(num_src_fields, collection_size)
                                                
       type(yac_dble_ptr), intent(inout) :: &
         src_frac_mask_buffer(num_src_fields, collection_size)
                                                
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_frac_dble_ptr
 
  end interface yac_fget_raw
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fget_async
 
     subroutine yac_fget_async_dble_ptr ( field_id,         &
                                          collection_size,  &
                                          recv_field,       &
                                          info,             &
                                          ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr)    :: recv_field(collection_size)
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_async_dble_ptr
 
  end interface yac_fget_async
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fget_raw_async
 
     subroutine yac_fget_raw_async_dble_ptr ( field_id,         &
                                              num_src_fields,   &
                                              collection_size,  &
                                              src_field_buffer, &
                                              info,             &
                                              ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: num_src_fields
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr), intent(inout) :: &
         src_field_buffer(num_src_fields, collection_size)
                                                
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_async_dble_ptr
 
     subroutine yac_fget_raw_frac_async_dble_ptr ( field_id,             &
                                                   num_src_fields,       &
                                                   collection_size,      &
                                                   src_field_buffer,     &
                                                   src_frac_mask_buffer, &
                                                   info,                 &
                                                   ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: field_id
       integer, intent (in)  :: num_src_fields
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr), intent(inout) :: &
         src_field_buffer(num_src_fields, collection_size)
                                                
       type(yac_dble_ptr), intent(inout) :: &
         src_frac_mask_buffer(num_src_fields, collection_size)
                                                
       integer, intent (out) :: info
       integer, intent (out) :: ierror
 
     end subroutine yac_fget_raw_frac_async_dble_ptr
 
  end interface yac_fget_raw_async
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fexchange
 
     subroutine yac_fexchange_real ( send_field_id,       &
                                     recv_field_id,       &
                                     send_nbr_hor_points, &
                                     send_nbr_pointsets,  &
                                     recv_nbr_hor_points, &
                                     collection_size,     &
                                     send_field,          &
                                     recv_field,          &
                                     send_info,           &
                                     recv_info,           &
                                     ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       real, intent (in)     :: send_field(send_nbr_hor_points, &
                                           send_nbr_pointsets,  &
                                           collection_size)
                                                    
       real, intent (inout)  :: recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_real
 
     subroutine yac_fexchange_real_ptr ( send_field_id,      &
                                         recv_field_id,      &
                                         send_nbr_pointsets, &
                                         collection_size,    &
                                         send_field,         &
                                         recv_field,         &
                                         send_info,          &
                                         recv_info,          &
                                         ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: collection_size
       type(yac_real_ptr), intent (in) ::                        &
                                send_field(send_nbr_pointsets, &
                                           collection_size)
                                                    
       type(yac_real_ptr)      :: recv_field(collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_real_ptr
 
     subroutine yac_fexchange_single_pointset_real ( send_field_id,       &
                                                     recv_field_id,       &
                                                     send_nbr_hor_points, &
                                                     recv_nbr_hor_points, &
                                                     collection_size,     &
                                                     send_field,          &
                                                     recv_field,          &
                                                     send_info,           &
                                                     recv_info,           &
                                                     ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       real, intent (in)     :: send_field(send_nbr_hor_points, &
                                           collection_size)
                                                     
       real, intent (inout)  :: recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                     
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_single_pointset_real
 
     subroutine yac_fexchange_dble ( send_field_id,       &
                                     recv_field_id,       &
                                     send_nbr_hor_points, &
                                     send_nbr_pointsets,  &
                                     recv_nbr_hor_points, &
                                     collection_size,     &
                                     send_field,          &
                                     recv_field,          &
                                     send_info,           &
                                     recv_info,           &
                                     ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       double precision, intent (in) :: &
                                send_field(send_nbr_hor_points, &
                                           send_nbr_pointsets,  &
                                           collection_size)
                                                    
       double precision, intent (inout):: &
                                recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_dble
 
     subroutine yac_fexchange_dble_ptr ( send_field_id,      &
                                         recv_field_id,      &
                                         send_nbr_pointsets, &
                                         collection_size,    &
                                         send_field,         &
                                         recv_field,         &
                                         send_info,          &
                                         recv_info,          &
                                         ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr), intent (in) ::                        &
                                send_field(send_nbr_pointsets, &
                                           collection_size)
                                                   
       type(yac_dble_ptr)      :: recv_field(collection_size)
                                                   
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_dble_ptr
 
     subroutine yac_fexchange_single_pointset_dble ( send_field_id,       &
                                                     recv_field_id,       &
                                                     send_nbr_hor_points, &
                                                     recv_nbr_hor_points, &
                                                     collection_size,     &
                                                     send_field,          &
                                                     recv_field,          &
                                                     send_info,           &
                                                     recv_info,           &
                                                     ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       double precision, intent (in) ::                         &
                                send_field(send_nbr_hor_points, &
                                           collection_size)
                                                    
       double precision, intent (inout)::                       &
                                recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_single_pointset_dble
 
     subroutine yac_fexchange_frac_real ( send_field_id,       &
                                          recv_field_id,       &
                                          send_nbr_hor_points, &
                                          send_nbr_pointsets,  &
                                          recv_nbr_hor_points, &
                                          collection_size,     &
                                          send_field,          &
                                          send_frac_mask,      &
                                          recv_field,          &
                                          send_info,           &
                                          recv_info,           &
                                          ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       real, intent (in)     :: send_field(send_nbr_hor_points, &
                                           send_nbr_pointsets,  &
                                           collection_size)
                                                    
       real, intent (in)     :: send_frac_mask(send_nbr_hor_points, &
                                               send_nbr_pointsets,  &
                                               collection_size)
                                                    
       real, intent (inout)  :: recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_frac_real
 
     subroutine yac_fexchange_frac_real_ptr ( send_field_id,      &
                                              recv_field_id,      &
                                              send_nbr_pointsets, &
                                              collection_size,    &
                                              send_field,         &
                                              send_frac_mask,     &
                                              recv_field,         &
                                              send_info,          &
                                              recv_info,          &
                                              ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: collection_size
       type(yac_real_ptr), intent (in) ::                      &
                                send_field(send_nbr_pointsets, &
                                           collection_size)
                                                    
       type(yac_real_ptr), intent (in) ::                          &
                                send_frac_mask(send_nbr_pointsets, &
                                               collection_size)
                                                    
       type(yac_real_ptr)    :: recv_field(collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_frac_real_ptr
 
     subroutine yac_fexchange_frac_single_pointset_real ( send_field_id,       &
                                                          recv_field_id,       &
                                                          send_nbr_hor_points, &
                                                          recv_nbr_hor_points, &
                                                          collection_size,     &
                                                          send_field,          &
                                                          send_frac_mask,      &
                                                          recv_field,          &
                                                          send_info,           &
                                                          recv_info,           &
                                                          ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       real, intent (in)     :: send_field(send_nbr_hor_points, &
                                           collection_size)
                                                     
       real, intent (in)     :: send_frac_mask(send_nbr_hor_points, &
                                               collection_size)
                                                     
       real, intent (inout)  :: recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                     
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_frac_single_pointset_real
 
     subroutine yac_fexchange_frac_dble ( send_field_id,       &
                                          recv_field_id,       &
                                          send_nbr_hor_points, &
                                          send_nbr_pointsets,  &
                                          recv_nbr_hor_points, &
                                          collection_size,     &
                                          send_field,          &
                                          send_frac_mask,      &
                                          recv_field,          &
                                          send_info,           &
                                          recv_info,           &
                                          ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       double precision, intent (in) :: &
                                send_field(send_nbr_hor_points, &
                                           send_nbr_pointsets,  &
                                           collection_size)
                                                    
       double precision, intent (in) :: &
                                send_frac_mask(send_nbr_hor_points, &
                                               send_nbr_pointsets,  &
                                               collection_size)
                                                    
       double precision, intent (inout):: &
                                recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_frac_dble
 
     subroutine yac_fexchange_frac_dble_ptr ( send_field_id,      &
                                              recv_field_id,      &
                                              send_nbr_pointsets, &
                                              collection_size,    &
                                              send_field,         &
                                              send_frac_mask,     &
                                              recv_field,         &
                                              send_info,          &
                                              recv_info,          &
                                              ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: collection_size
       type(yac_dble_ptr), intent (in) ::                      &
                                send_field(send_nbr_pointsets, &
                                           collection_size)
                                                   
       type(yac_dble_ptr), intent (in) ::                          &
                                send_frac_mask(send_nbr_pointsets, &
                                               collection_size)
                                                   
       type(yac_dble_ptr)      :: recv_field(collection_size)
                                                   
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_frac_dble_ptr
 
     subroutine yac_fexchange_frac_single_pointset_dble ( send_field_id,       &
                                                          recv_field_id,       &
                                                          send_nbr_hor_points, &
                                                          recv_nbr_hor_points, &
                                                          collection_size,     &
                                                          send_field,          &
                                                          send_frac_mask,      &
                                                          recv_field,          &
                                                          send_info,           &
                                                          recv_info,           &
                                                          ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: recv_nbr_hor_points
       integer, intent (in)  :: collection_size
       double precision, intent (in) ::                         &
                                send_field(send_nbr_hor_points, &
                                           collection_size)
                                                    
       double precision, intent (in) ::                             &
                                send_frac_mask(send_nbr_hor_points, &
                                               collection_size)
                                                    
       double precision, intent (inout)::                       &
                                recv_field(recv_nbr_hor_points, &
                                           collection_size)
                                                    
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_frac_single_pointset_dble
 
  end interface yac_fexchange
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fexchange_raw
 
     subroutine yac_fexchange_raw_real ( send_field_id,         &
                                         recv_field_id,         &
                                         send_nbr_hor_points,   &
                                         send_nbr_pointsets,    &
                                         src_field_buffer_size, &
                                         collection_size,       &
                                         send_field,            &
                                         src_field_buffer,      &
                                         send_info,             &
                                         recv_info,             &
                                         ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: src_field_buffer_size
       integer, intent (in)  :: collection_size
       real, intent (in)     :: &
        send_field(send_nbr_hor_points, send_nbr_pointsets, collection_size)
                                                      
       real, intent (out)    :: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                      
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_raw_real
 
     subroutine yac_fexchange_raw_real_ptr ( send_field_id,      &
                                             recv_field_id,      &
                                             send_nbr_pointsets, &
                                             collection_size,    &
                                             send_field,         &
                                             src_field_buffer,   &
                                             send_info,          &
                                             recv_info,          &
                                             ierror )
 
       import :: yac_real_ptr
 
       integer, intent (in)               :: send_field_id
       integer, intent (in)               :: recv_field_id
       integer, intent (in)               :: send_nbr_pointsets
       integer, intent (in)               :: collection_size
       type(yac_real_ptr), intent (in)    :: &
         send_field(send_nbr_pointsets, collection_size)
       type(yac_real_ptr), intent (inout) :: &
         src_field_buffer(send_nbr_pointsets, collection_size)
       integer, intent (out)              :: send_info
       integer, intent (out)              :: recv_info
       integer, intent (out)              :: ierror
 
     end subroutine yac_fexchange_raw_real_ptr
 
     subroutine yac_fexchange_raw_single_pointset_real ( send_field_id,         &
                                                         recv_field_id,         &
                                                         send_nbr_hor_points,   &
                                                         src_field_buffer_size, &
                                                         collection_size,       &
                                                         send_field,            &
                                                         src_field_buffer,      &
                                                         send_info,             &
                                                         recv_info,             &
                                                         ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: src_field_buffer_size
       integer, intent (in)  :: collection_size
       real, intent (in)     :: send_field(send_nbr_hor_points, &
                                           collection_size)
                                                      
       real, intent (out) :: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                      
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_raw_single_pointset_real
 
     subroutine yac_fexchange_raw_dble ( send_field_id,         &
                                         recv_field_id,         &
                                         send_nbr_hor_points,   &
                                         send_nbr_pointsets,    &
                                         src_field_buffer_size, &
                                         collection_size,       &
                                         send_field,            &
                                         src_field_buffer,      &
                                         send_info,             &
                                         recv_info,             &
                                         ierror )
 
       integer, intent (in)  :: send_field_id
       integer, intent (in)  :: recv_field_id
       integer, intent (in)  :: send_nbr_hor_points
       integer, intent (in)  :: send_nbr_pointsets
       integer, intent (in)  :: src_field_buffer_size
       integer, intent (in)  :: collection_size
       double precision, intent (in)  ::                        &
                                send_field(send_nbr_hor_points, &
                                           send_nbr_pointsets,  &
                                           collection_size)
                                                      
       double precision, intent (out) :: &
        src_field_buffer(src_field_buffer_size, collection_size)
                                                      
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_raw_dble
 
     subroutine yac_fexchange_raw_dble_ptr ( send_field_id,      &
                                             recv_field_id,      &
                                             send_nbr_pointsets, &
                                             collection_size,    &
                                             send_field,         &
                                             src_field_buffer,   &
                                             send_info,          &
                                             recv_info,          &
                                             ierror )
 
       import :: yac_dble_ptr
 
       integer, intent (in)               :: send_field_id
       integer, intent (in)               :: recv_field_id
       integer, intent (in)               :: send_nbr_pointsets
       integer, intent (in)               :: collection_size
       type(yac_dble_ptr), intent (in)    :: &
         send_field(send_nbr_pointsets, collection_size)
       type(yac_dble_ptr), intent (inout) :: &
         src_field_buffer(send_nbr_pointsets, collection_size)
       integer, intent (out)              :: send_info
       integer, intent (out)              :: recv_info
       integer, intent (out)              :: ierror
 
     end subroutine yac_fexchange_raw_dble_ptr
 
     subroutine yac_fexchange_raw_single_pointset_dble ( send_field_id,         &
                                                         recv_field_id,         &
                                                         send_nbr_hor_points,   &
                                                         src_field_buffer_size, &
                                                         collection_size,       &
                                                         send_field,            &
                                                         src_field_buffer,      &
                                                         send_info,             &
                                                         recv_info,             &
                                                         ierror )
 
       integer, intent (in)          :: send_field_id
       integer, intent (in)          :: recv_field_id
       integer, intent (in)          :: send_nbr_hor_points
       integer, intent (in)          :: src_field_buffer_size
       integer, intent (in)          :: collection_size
       double precision, intent (in) :: &
         send_field(send_nbr_hor_points, collection_size)
       double precision, intent (out) :: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                              
       integer, intent (out) :: send_info
       integer, intent (out) :: recv_info
       integer, intent (out) :: ierror
 
     end subroutine yac_fexchange_raw_single_pointset_dble
 
    subroutine yac_fexchange_raw_frac_real ( send_field_id,         &
                                             recv_field_id,         &
                                             send_nbr_hor_points,   &
                                             send_nbr_pointsets,    &
                                             src_field_buffer_size, &
                                             collection_size,       &
                                             send_field,            &
                                             send_frac_mask,        &
                                             src_field_buffer,      &
                                             src_frac_mask_buffer,  &
                                             send_info,             &
                                             recv_info,             &
                                             ierror )
 
      integer, intent (in)  :: send_field_id
      integer, intent (in)  :: recv_field_id
      integer, intent (in)  :: send_nbr_hor_points
      integer, intent (in)  :: send_nbr_pointsets
      integer, intent (in)  :: src_field_buffer_size
      integer, intent (in)  :: collection_size
      real, intent (in)     :: &
        send_field(send_nbr_hor_points, send_nbr_pointsets, collection_size)
                                                     
      real, intent (in)     :: &
        send_frac_mask(send_nbr_hor_points, send_nbr_pointsets, collection_size)
                                                     
      real, intent (out)    :: &
        src_field_buffer(src_field_buffer_size, collection_size)
                                                     
      real, intent (out)    :: &
        src_frac_mask_buffer(src_field_buffer_size, collection_size)
                                                     
      integer, intent (out) :: send_info
      integer, intent (out) :: recv_info
      integer, intent (out) :: ierror
 
    end subroutine yac_fexchange_raw_frac_real
 
    subroutine yac_fexchange_raw_frac_real_ptr ( send_field_id,        &
                                                 recv_field_id,        &
                                                 send_nbr_pointsets,   &
                                                 collection_size,      &
                                                 send_field,           &
                                                 send_frac_mask,       &
                                                 src_field_buffer,     &
                                                 src_frac_mask_buffer, &
                                                 send_info,            &
                                                 recv_info,            &
                                                 ierror )
 
      import :: yac_real_ptr
 
      integer, intent (in)               :: send_field_id
      integer, intent (in)               :: recv_field_id
      integer, intent (in)               :: send_nbr_pointsets
      integer, intent (in)               :: collection_size
      type(yac_real_ptr), intent (in)    :: &
        send_field(send_nbr_pointsets, collection_size)
      type(yac_real_ptr), intent (in)    :: &
        send_frac_mask(send_nbr_pointsets, collection_size)
      type(yac_real_ptr), intent(inout) :: &
        src_field_buffer(send_nbr_pointsets, collection_size)
      type(yac_real_ptr), intent(inout) :: &
        src_frac_mask_buffer(send_nbr_pointsets, collection_size)
                                                               
      integer, intent (out)              :: send_info
      integer, intent (out)              :: recv_info
      integer, intent (out)              :: ierror
 
    end subroutine yac_fexchange_raw_frac_real_ptr
 
    subroutine yac_fexchange_raw_frac_single_pointset_real ( send_field_id,         &
                                                             recv_field_id,         &
                                                             send_nbr_hor_points,   &
                                                             src_field_buffer_size, &
                                                             collection_size,       &
                                                             send_field,            &
                                                             send_frac_mask,        &
                                                             src_field_buffer,      &
                                                             src_frac_mask_buffer,  &
                                                             send_info,             &
                                                             recv_info,             &
                                                             ierror )
 
      integer, intent (in)  :: send_field_id
      integer, intent (in)  :: recv_field_id
      integer, intent (in)  :: send_nbr_hor_points
      integer, intent (in)  :: src_field_buffer_size
      integer, intent (in)  :: collection_size
      real, intent (in)     :: &
        send_field(send_nbr_hor_points, collection_size)
                                                     
      real, intent (in)     :: &
        send_frac_mask(send_nbr_hor_points, collection_size)
                                                     
      real, intent (out)    :: &
        src_field_buffer(src_field_buffer_size, collection_size)
                                                     
      real, intent (out)    :: &
        src_frac_mask_buffer(src_field_buffer_size, collection_size)
                                                     
      integer, intent (out) :: send_info
      integer, intent (out) :: recv_info
      integer, intent (out) :: ierror
 
    end subroutine yac_fexchange_raw_frac_single_pointset_real
 
    subroutine yac_fexchange_raw_frac_dble ( send_field_id,         &
                                             recv_field_id,         &
                                             send_nbr_hor_points,   &
                                             send_nbr_pointsets,    &
                                             src_field_buffer_size, &
                                             collection_size,       &
                                             send_field,            &
                                             send_frac_mask,        &
                                             src_field_buffer,      &
                                             src_frac_mask_buffer,  &
                                             send_info,             &
                                             recv_info,             &
                                             ierror )
 
      integer, intent (in)  :: send_field_id
      integer, intent (in)  :: recv_field_id
      integer, intent (in)  :: send_nbr_hor_points
      integer, intent (in)  :: send_nbr_pointsets
      integer, intent (in)  :: src_field_buffer_size
      integer, intent (in)  :: collection_size
      double precision, intent (in) :: &
        send_field(send_nbr_hor_points, send_nbr_pointsets, collection_size)
                                                     
      double precision, intent (in) :: &
        send_frac_mask(send_nbr_hor_points, send_nbr_pointsets, collection_size)
                                                     
      double precision, intent (out):: &
         src_field_buffer(src_field_buffer_size, collection_size)
                                                     
      double precision, intent (out):: &
         src_frac_mask_buffer(src_field_buffer_size, collection_size)
                                                     
      integer, intent (out) :: send_info
      integer, intent (out) :: recv_info
      integer, intent (out) :: ierror
 
    end subroutine yac_fexchange_raw_frac_dble
 
    subroutine yac_fexchange_raw_frac_dble_ptr ( send_field_id,        &
                                                 recv_field_id,        &
                                                 send_nbr_pointsets,   &
                                                 collection_size,      &
                                                 send_field,           &
                                                 send_frac_mask,       &
                                                 src_field_buffer,     &
                                                 src_frac_mask_buffer, &
                                                 send_info,            &
                                                 recv_info,            &
                                                 ierror )
 
      import :: yac_dble_ptr
 
      integer, intent (in)              :: send_field_id
      integer, intent (in)              :: recv_field_id
      integer, intent (in)              :: send_nbr_pointsets
      integer, intent (in)              :: collection_size
      type(yac_dble_ptr), intent (in)   :: &
        send_field(send_nbr_pointsets, collection_size)
      type(yac_dble_ptr), intent (in)   :: &
        send_frac_mask(send_nbr_pointsets, collection_size)
      type(yac_dble_ptr), intent(inout) :: &
        src_field_buffer(send_nbr_pointsets, collection_size)
      type(yac_dble_ptr), intent(inout) :: &
        src_frac_mask_buffer(send_nbr_pointsets, collection_size)
                                                              
      integer, intent (out)             :: send_info
      integer, intent (out)             :: recv_info
      integer, intent (out)             :: ierror
 
    end subroutine yac_fexchange_raw_frac_dble_ptr
 
    subroutine yac_fexchange_raw_frac_single_pointset_dble ( send_field_id,         &
                                                             recv_field_id,         &
                                                             send_nbr_hor_points,   &
                                                             src_field_buffer_size, &
                                                             collection_size,       &
                                                             send_field,            &
                                                             send_frac_mask,        &
                                                             src_field_buffer,      &
                                                             src_frac_mask_buffer,  &
                                                             send_info,             &
                                                             recv_info,             &
                                                             ierror )
 
      integer, intent (in)           :: send_field_id
      integer, intent (in)           :: recv_field_id
      integer, intent (in)           :: send_nbr_hor_points
      integer, intent (in)           :: src_field_buffer_size
      integer, intent (in)           :: collection_size
      double precision, intent (in)  :: &
        send_field(send_nbr_hor_points, collection_size)
      double precision, intent (in)  :: &
        send_frac_mask(send_nbr_hor_points, collection_size)
      double precision, intent (out) :: &
        src_field_buffer(src_field_buffer_size, collection_size)
                                                              
      double precision, intent (out) :: &
        src_frac_mask_buffer(src_field_buffer_size, collection_size)
                                                              
      integer, intent (out)          :: send_info
      integer, intent (out)          :: recv_info
      integer, intent (out)          :: ierror
 
    end subroutine yac_fexchange_raw_frac_single_pointset_dble
 
  end interface yac_fexchange_raw
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_ftest
 
     subroutine yac_ftest_i ( field_id, flag )
 
       integer, intent (in)  :: field_id
       integer, intent (out) :: flag
 
     end subroutine yac_ftest_i
 
     subroutine yac_ftest_l ( field_id, flag )
 
       integer, intent (in)  :: field_id
       logical, intent (out) :: flag
 
     end subroutine yac_ftest_l
 
  end interface yac_ftest
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fwait
 
     subroutine yac_fwait ( field_id )
 
       integer, intent (in)  :: field_id
 
     end subroutine yac_fwait
 
  end interface yac_fwait
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface
 
     subroutine yac_fget_comp_comm ( comp_id, comp_comm )
 
       integer, intent (in)  :: comp_id
       integer, intent (out) :: comp_comm
 
     end subroutine yac_fget_comp_comm
 
  end interface
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fget_comps_comm
 
     subroutine yac_fget_comps_comm ( comp_names, &
                                      num_comps,  &
                                      comps_comm)
 
       use, intrinsic :: iso_c_binding, only : c_char
       import :: yac_max_charlen
 
       integer, intent(in)   :: num_comps
       character(kind=c_char, len=YAC_MAX_CHARLEN), intent(in) :: &
                                comp_names(num_comps)
       integer, intent (out) :: comps_comm
 
     end subroutine yac_fget_comps_comm
 
     subroutine yac_fget_comps_comm_instance ( yac_instance_id, &
                                               comp_names,      &
                                               num_comps,       &
                                               comps_comm)
 
       use, intrinsic :: iso_c_binding, only : c_char
       import :: yac_max_charlen
 
       integer, intent(in)          :: yac_instance_id
       integer, intent(in)   :: num_comps
       character(kind=c_char, len=YAC_MAX_CHARLEN), intent(in) :: &
                                comp_names(num_comps)
       integer, intent (out) :: comps_comm
 
     end subroutine yac_fget_comps_comm_instance
 
  end interface yac_fget_comps_comm
 
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fsync_def
 
     subroutine yac_fsync_def ( )
 
     end subroutine yac_fsync_def
 
     subroutine yac_fsync_def_instance ( yac_instance_id )
 
       integer, intent(in)  :: yac_instance_id
 
     end subroutine yac_fsync_def_instance
 
  end interface yac_fsync_def
 
  interface yac_fsync_def_comps
 
     subroutine yac_fsync_def_comps ( comp_names, num_comps )
 
       use, intrinsic :: iso_c_binding, only : c_char
       import :: yac_max_charlen
 
       integer, intent(in)   :: num_comps
       character(kind=c_char, len=YAC_MAX_CHARLEN), intent(in) :: &
                                comp_names(num_comps)
 
     end subroutine yac_fsync_def_comps
 
     subroutine yac_fsync_def_comps_instance ( yac_instance_id, comp_names, num_comps )
 
       use, intrinsic :: iso_c_binding, only : c_char
       import :: yac_max_charlen
 
       integer, intent(in) :: yac_instance_id
       integer, intent(in) :: num_comps
       character(kind=c_char, len=YAC_MAX_CHARLEN), intent(in) :: &
                              comp_names(num_comps)
 
     end subroutine yac_fsync_def_comps_instance
 
  end interface yac_fsync_def_comps
 
 !----------------------------------------------------------------------
 !----------------------------------------------------------------------
 
 interface
 
    subroutine yac_fget_interp_stack_config(interp_stack_config_id)
      integer, intent(out) :: interp_stack_config_id
    end subroutine yac_fget_interp_stack_config
 
    subroutine yac_ffree_interp_stack_config(interp_stack_config_id)
      integer, intent(in) :: interp_stack_config_id
    end subroutine yac_ffree_interp_stack_config
 
    subroutine yac_fadd_interp_stack_config_average(interp_stack_config_id, &
         reduction_type, partial_coverage)
      integer, intent(in)     :: interp_stack_config_id
      integer, intent(in)     :: reduction_type
      integer, intent(in)     :: partial_coverage
    end subroutine yac_fadd_interp_stack_config_average
 
    subroutine yac_fadd_interp_stack_config_ncc(interp_stack_config_id, &
         weight_type, partial_coverage)
      integer, intent(in)     :: interp_stack_config_id
      integer, intent(in)     :: weight_type
      integer, intent(in)     :: partial_coverage
    end subroutine yac_fadd_interp_stack_config_ncc
 
    subroutine yac_fadd_interp_stack_config_nnn(interp_stack_config_id, &
         type, n, max_search_distance, scale)
      integer, intent(in)          :: interp_stack_config_id
      integer, intent(in)          :: type
      integer, intent(in)          :: n
      double precision, intent(in) :: max_search_distance
      double precision, intent(in) :: scale
    end subroutine yac_fadd_interp_stack_config_nnn
 
    subroutine yac_fadd_interp_stack_config_dnn(interp_stack_config_id, &
         weight_type, search_distance_type, search_distance, scale)
      integer, intent(in)          :: interp_stack_config_id
      integer, intent(in)          :: weight_type
      integer, intent(in)          :: search_distance_type
      double precision, intent(in) :: search_distance
      double precision, intent(in) :: scale
    end subroutine yac_fadd_interp_stack_config_dnn
 
    subroutine yac_fadd_interp_stack_config_rbf(interp_stack_config_id, &
         n, max_search_distance, scale)
      integer, intent(in)                     :: interp_stack_config_id
      integer, optional, intent(in)          :: n
      double precision, optional, intent(in) :: max_search_distance
      double precision, optional, intent(in) :: scale
    end subroutine yac_fadd_interp_stack_config_rbf
 
    subroutine yac_fadd_interp_stack_config_conservative(interp_stack_config_id, &
         order, enforced_conserv, partial_coverage, normalization)
      integer, intent(in)     :: interp_stack_config_id
      integer, intent(in)     :: order
      integer, intent(in)     :: enforced_conserv
      integer, intent(in)     :: partial_coverage
      integer, intent(in)     :: normalization
    end subroutine yac_fadd_interp_stack_config_conservative
 
    subroutine yac_fget_spmap_overwrite_config_bnd_circle ( &
      center_lon, center_lat, inc_angle, spmap_overwrite_config_id, &
      spread_distance, max_search_distance, weight_type )
      double precision, intent(in)           :: center_lon
      double precision, intent(in)           :: center_lat
      double precision, intent(in)           :: inc_angle
      integer, intent(out)                   :: spmap_overwrite_config_id
      double precision, optional, intent(in) :: spread_distance
      double precision, optional, intent(in) :: max_search_distance
      integer, optional, intent(in)          :: weight_type
    end subroutine yac_fget_spmap_overwrite_config_bnd_circle
 
    subroutine yac_ffree_spmap_overwrite_config ( &
      spmap_overwrite_config_id )
      integer, intent(in) :: spmap_overwrite_config_id
    end subroutine yac_ffree_spmap_overwrite_config
 
    subroutine yac_fadd_interp_stack_config_spmap(interp_stack_config_id, &
         search_distance, max_search_distance, weight_type, scale_type, &
         src_sphere_radius, src_filename, src_varname, src_min_global_id, &
         tgt_sphere_radius, tgt_filename, tgt_varname, tgt_min_global_id, &
         spmap_overwrite_config_ids)
      integer, intent(in)           :: interp_stack_config_id
      double precision, intent(in)  :: search_distance
      double precision, intent(in)  :: max_search_distance
      integer, intent(in)           :: weight_type
      integer, intent(in)           :: scale_type
      double precision, intent(in)  :: src_sphere_radius
      character (len=*), intent(in) :: src_filename
      character (len=*), intent(in) :: src_varname
      integer, intent(in)           :: src_min_global_id
      double precision, intent(in)  :: tgt_sphere_radius
      character (len=*), intent(in) :: tgt_filename
      character (len=*), intent(in) :: tgt_varname
      integer, intent(in)           :: tgt_min_global_id
      integer, optional, intent(in) :: spmap_overwrite_config_ids(:)
    end subroutine yac_fadd_interp_stack_config_spmap
 
    subroutine yac_fadd_interp_stack_config_hcsbb(interp_stack_config_id)
      integer, intent(in)     :: interp_stack_config_id
    end subroutine yac_fadd_interp_stack_config_hcsbb
 
    subroutine yac_fadd_interp_stack_config_user_file( &
         interp_stack_config_id, filename, on_missing_file, on_success)
      integer, intent(in)     :: interp_stack_config_id
      character (len=*), intent(in) :: filename
      integer, optional, intent(in) :: on_missing_file
      integer, optional, intent(in) :: on_success
    end subroutine yac_fadd_interp_stack_config_user_file
 
    subroutine yac_fadd_interp_stack_config_fixed(interp_stack_config_id, &
         val)
      integer, intent(in)          :: interp_stack_config_id
      double precision, intent(in) :: val
    end subroutine yac_fadd_interp_stack_config_fixed
 
    subroutine yac_fadd_interp_stack_config_creep(interp_stack_config_id, &
         creep_distance)
      integer, intent(in) :: interp_stack_config_id
      integer, intent(in) :: creep_distance
    end subroutine yac_fadd_interp_stack_config_creep
 
    subroutine yac_fget_interp_stack_config_from_string_yaml(&
        interp_stack_config, interp_stack_config_id)
      character ( len=* ), intent(in) :: interp_stack_config
      integer, intent(out) :: interp_stack_config_id
    end subroutine yac_fget_interp_stack_config_from_string_yaml
 
    subroutine yac_fget_interp_stack_config_from_string_json(&
        interp_stack_config, interp_stack_config_id)
      character ( len=* ), intent(in) :: interp_stack_config
      integer, intent(out) :: interp_stack_config_id
    end subroutine yac_fget_interp_stack_config_from_string_json
 
 end interface
 
 !----------------------------------------------------------------------
 !----------------------------------------------------------------------
 
 interface yac_fdef_couple
 
    subroutine yac_fdef_couple (                     &
       src_comp_name, src_grid_name, src_field_name, &
       tgt_comp_name, tgt_grid_name, tgt_field_name, &
       coupling_timestep, time_unit, time_reduction, &
       interp_stack_config_id, src_lag, tgt_lag,     &
       weight_file, weight_file_on_existing,         &
       mapping_side, scale_factor, scale_summand,    &
       src_mask_names, tgt_mask_name,                &
       yaxt_exchanger_name, collection_indices,      &
       use_raw_exchange )
 
      import :: yac_string
      character ( len=* ), intent(in)           :: src_comp_name
      character ( len=* ), intent(in)           :: src_grid_name
      character ( len=* ), intent(in)           :: src_field_name
      character ( len=* ), intent(in)           :: tgt_comp_name
      character ( len=* ), intent(in)           :: tgt_grid_name
      character ( len=* ), intent(in)           :: tgt_field_name
      character ( len=* ), intent(in)           :: coupling_timestep
      integer, intent(in)                       :: time_unit
      integer, intent(in)                       :: time_reduction
      integer, intent(in)                       :: interp_stack_config_id
      integer, intent(in), optional             :: src_lag
      integer, intent(in), optional             :: tgt_lag
      character ( len=* ), intent(in), optional :: weight_file
      integer, intent(in), optional             :: weight_file_on_existing
      integer, intent(in), optional             :: mapping_side
      double precision, intent(in), optional    :: scale_factor
      double precision, intent(in), optional    :: scale_summand
      type(yac_string), intent(in), optional    :: src_mask_names(:)
      character ( len=* ), intent(in), optional :: tgt_mask_name
      character ( len=* ), intent(in), optional :: yaxt_exchanger_name
      integer, intent(in), optional             :: collection_indices(:)
      logical, intent(in), optional             :: use_raw_exchange
    end subroutine yac_fdef_couple
 
    subroutine yac_fdef_couple_instance ( instance_id, &
         src_comp_name, src_grid_name, src_field_name, &
         tgt_comp_name, tgt_grid_name, tgt_field_name, &
         coupling_timestep, time_unit, time_reduction, &
         interp_stack_config_id, src_lag, tgt_lag,     &
         weight_file, weight_file_on_existing,         &
         mapping_side, scale_factor, scale_summand,    &
         src_mask_names, tgt_mask_name,                &
         yaxt_exchanger_name, collection_indices,      &
         use_raw_exchange )
 
      import :: yac_string
      integer, intent(in)                       :: instance_id
      character ( len=* ), intent(in)           :: src_comp_name
      character ( len=* ), intent(in)           :: src_grid_name
      character ( len=* ), intent(in)           :: src_field_name
      character ( len=* ), intent(in)           :: tgt_comp_name
      character ( len=* ), intent(in)           :: tgt_grid_name
      character ( len=* ), intent(in)           :: tgt_field_name
      character ( len=* ), intent(in)           :: coupling_timestep
      integer, intent(in)                       :: time_unit
      integer, intent(in)                       :: time_reduction
      integer, intent(in)                       :: interp_stack_config_id
      integer, intent(in), optional             :: src_lag
      integer, intent(in), optional             :: tgt_lag
      character ( len=* ), intent(in), optional :: weight_file
      integer, intent(in), optional             :: weight_file_on_existing
      integer, intent(in), optional             :: mapping_side
      double precision, intent(in), optional    :: scale_factor
      double precision, intent(in), optional    :: scale_summand
      type(yac_string), intent(in), optional    :: src_mask_names(:)
      character ( len=* ), intent(in), optional :: tgt_mask_name
      character ( len=* ), intent(in), optional :: yaxt_exchanger_name
      integer, intent(in), optional             :: collection_indices(:)
      logical, intent(in), optional             :: use_raw_exchange
    end subroutine yac_fdef_couple_instance
 
 end interface yac_fdef_couple
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fenddef
 
    subroutine yac_fenddef ( )
 
    end subroutine yac_fenddef
 
    subroutine yac_fenddef_instance ( yac_instance_id )
 
      integer, intent(in)  :: yac_instance_id
 
    end subroutine yac_fenddef_instance
 
    subroutine yac_fenddef_and_emit_config(emit_flags, config)
 
      integer, intent (in)           :: emit_flags
      character (len=:), ALLOCATABLE :: config
 
    end subroutine yac_fenddef_and_emit_config
 
    subroutine yac_fenddef_and_emit_config_instance( &
      yac_instance_id, emit_flags, config)
 
      integer, intent (in)           :: yac_instance_id
      integer, intent (in)           :: emit_flags
      character (len=:), ALLOCATABLE :: config
 
    end subroutine yac_fenddef_and_emit_config_instance
 
  end interface yac_fenddef
 
  !----------------------------------------------------------------------
  !----------------------------------------------------------------------
 
  interface yac_fget_grid_size
 
    function yac_fget_grid_size( location, grid_id ) result (grid_size)
 
      integer, intent(in) :: location
      integer, intent(in) :: grid_id
      integer :: grid_size
 
    end function yac_fget_grid_size
 
  end interface yac_fget_grid_size
 
  ! ---------------------------------------------------------------------
 
  interface yac_fcompute_grid_cell_areas
 
    subroutine yac_fcompute_grid_cell_areas_real( &
      grid_id, nbr_cells, cell_areas )
 
      integer, intent(in) :: grid_id
      integer, intent(in) :: nbr_cells
      real, intent(out)   :: cell_areas(nbr_cells)
    end subroutine yac_fcompute_grid_cell_areas_real
 
    subroutine yac_fcompute_grid_cell_areas_dble( &
      grid_id, nbr_cells, cell_areas )
 
      integer, intent(in)           :: grid_id
      integer, intent(in)           :: nbr_cells
      double precision, intent(out) :: cell_areas(nbr_cells)
    end subroutine yac_fcompute_grid_cell_areas_dble
 
  end interface
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_points_size
 
    function yac_fget_points_size( point_id ) result (points_size)
 
      integer, intent(in) :: point_id
      integer :: points_size
 
    end function yac_fget_points_size
 
  end interface yac_fget_points_size
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_comp_id
 
     function yac_fget_comp_id(comp_name) result(comp_id)
 
       character(len=*), intent(in) :: comp_name
       integer :: comp_id
 
     end function yac_fget_comp_id
 
     function yac_fget_comp_id_instance(yac_id, comp_name) result(comp_id)
 
       integer, intent(in) :: yac_id
       character(len=*), intent(in) :: comp_name
       integer :: comp_id
 
     end function yac_fget_comp_id_instance
  end interface yac_fget_comp_id
 
  interface yac_fget_grid_id
 
     function yac_fget_grid_id(grid_name) result(grid_id)
 
       character(len=*), intent(in) :: grid_name
       integer :: grid_id
 
     end function yac_fget_grid_id
 
  end interface yac_fget_grid_id
 
  interface yac_fget_points_id
 
     function yac_fget_points_id(grid_id, location, point_name) &
       result(point_id)
 
       integer, intent(in) :: grid_id
       integer, intent(in) :: location
       character(len=*), intent(in) :: point_name
       integer :: point_id
 
     end function yac_fget_points_id
 
  end interface yac_fget_points_id
 
  interface yac_fget_mask_id
 
     function yac_fget_mask_id(grid_id, location, mask_name) result(mask_id)
 
       integer, intent(in) :: grid_id
       integer, intent(in) :: location
       character(len=*), intent(in) :: mask_name
       integer :: mask_id
 
     end function yac_fget_mask_id
 
  end interface yac_fget_mask_id
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_is_defined
 
     function yac_fget_field_is_defined ( comp_name, grid_name, field_name ) &
       result(field_is_defined)
 
       character(len=*), intent(in)  :: comp_name
       character(len=*), intent(in)  :: grid_name
       character(len=*), intent(in)  :: field_name
       logical :: field_is_defined
 
     end function yac_fget_field_is_defined
 
     function yac_fget_field_is_defined_instance ( yac_id,&
                                           comp_name,     &
                                           grid_name,     &
                                           field_name )   &
       result(field_is_defined)
 
       integer, intent(in) :: yac_id
       character(len=*), intent(in)  :: comp_name
       character(len=*), intent(in)  :: grid_name
       character(len=*), intent(in)  :: field_name
       logical :: field_is_defined
 
     end function yac_fget_field_is_defined_instance
  end interface yac_fget_field_is_defined
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_id
 
     function yac_fget_field_id ( comp_name, grid_name, field_name ) result(field_id)
 
       character(len=*), intent(in)  :: comp_name
       character(len=*), intent(in)  :: grid_name
       character(len=*), intent(in)  :: field_name
       integer :: field_id
 
     end function yac_fget_field_id
 
     function yac_fget_field_id_instance ( yac_id,    &
                                           comp_name, &
                                           grid_name, &
                                           field_name ) result(field_id)
 
       integer, intent(in) :: yac_id
       character(len=*), intent(in)  :: comp_name
       character(len=*), intent(in)  :: grid_name
       character(len=*), intent(in)  :: field_name
       integer :: field_id
 
     end function yac_fget_field_id_instance
  end interface yac_fget_field_id
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_component_name
 
     function yac_fget_component_name_from_field_id ( field_id ) &
          result( comp_name )
 
       integer, intent (in)           :: field_id
       character(len=:), ALLOCATABLE :: comp_name
 
     end function yac_fget_component_name_from_field_id
 
  end interface yac_fget_component_name
 
  interface yac_fget_grid_name
 
     function yac_fget_grid_name_from_field_id ( field_id ) &
          result( grid_name )
 
       integer, intent (in)           :: field_id
       character(len=:), ALLOCATABLE :: grid_name
 
     end function yac_fget_grid_name_from_field_id
 
  end interface yac_fget_grid_name
 
  interface yac_fget_field_name
 
     function yac_fget_field_name_from_field_id ( field_id ) &
          result( field_name )
 
       integer, intent (in)           :: field_id
       character(len=:), ALLOCATABLE :: field_name
 
     end function yac_fget_field_name_from_field_id
 
  end interface yac_fget_field_name
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_comp_names
 
     function yac_fget_comp_names( ) result(comp_names)
       import :: yac_string
       type(yac_string), allocatable :: comp_names(:)
     end function  yac_fget_comp_names
 
     function yac_fget_comp_names_instance( yac_instance_id) &
         result( comp_names )
       import :: yac_string
       integer, intent(in) :: yac_instance_id
       type(yac_string), allocatable :: comp_names(:)
     end function yac_fget_comp_names_instance
 
  end interface yac_fget_comp_names
 
  interface yac_fget_grid_names
 
     function yac_fget_grid_names( ) result( grid_names )
       import :: yac_string
       type(yac_string), allocatable :: grid_names(:)
     end function yac_fget_grid_names
 
     function yac_fget_grid_names_instance( yac_instance_id ) &
         result( grid_names )
       import :: yac_string
       integer, intent(in) :: yac_instance_id
       type(yac_string), allocatable :: grid_names(:)
     end function yac_fget_grid_names_instance
 
  end interface yac_fget_grid_names
 
  interface yac_fget_comp_grid_names
 
     function yac_fget_comp_grid_names( comp_name ) result( grid_names )
       import :: yac_string
       character(len=*), intent(in)  :: comp_name
       type(yac_string), allocatable :: grid_names(:)
     end function yac_fget_comp_grid_names
 
     function yac_fget_comp_grid_names_instance( yac_instance_id, comp_name ) &
         result( grid_names )
       import :: yac_string
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in)  :: comp_name
       type(yac_string), allocatable :: grid_names(:)
     end function yac_fget_comp_grid_names_instance
 
  end interface yac_fget_comp_grid_names
 
  interface yac_fget_field_names
 
     function yac_fget_field_names( comp_name, grid_name ) result ( field_names )
       import :: yac_string
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       type(yac_string), allocatable :: field_names(:)
     end function yac_fget_field_names
 
     function yac_fget_field_names_instance( yac_instance_id, &
                                             comp_name,       &
                                             grid_name )      &
          result( field_names )
       import :: yac_string
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       type(yac_string), allocatable :: field_names(:)
     end function yac_fget_field_names_instance
 
  end interface yac_fget_field_names
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_role
 
     function yac_fget_role_from_field_id ( field_id )
 
       integer, intent (in)           :: field_id
       integer                        :: yac_fget_role_from_field_id
 
     end function yac_fget_role_from_field_id
 
     function yac_fget_field_role ( comp_name, grid_name, field_name )
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       integer :: yac_fget_field_role
     end function yac_fget_field_role
 
     function yac_fget_field_role_instance ( yac_instance_id, &
                                             comp_name,       &
                                             grid_name,       &
                                             field_name )
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       integer :: yac_fget_field_role_instance
     end function yac_fget_field_role_instance
 
  end interface yac_fget_field_role
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_source
 
     subroutine yac_fget_field_source ( tgt_comp_name,   &
                                        tgt_grid_name,   &
                                        tgt_field_name,  &
                                        src_comp_name,   &
                                        src_grid_name,   &
                                        src_field_name )
       character(len=*), intent(in) :: tgt_comp_name
       character(len=*), intent(in) :: tgt_grid_name
       character(len=*), intent(in) :: tgt_field_name
       character(len=:), ALLOCATABLE :: src_comp_name
       character(len=:), ALLOCATABLE :: src_grid_name
       character(len=:), ALLOCATABLE :: src_field_name
     end subroutine yac_fget_field_source
 
     subroutine yac_fget_field_source_instance ( yac_instance_id, &
                                                 tgt_comp_name,   &
                                                 tgt_grid_name,   &
                                                 tgt_field_name,  &
                                                 src_comp_name,   &
                                                 src_grid_name,   &
                                                 src_field_name )
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: tgt_comp_name
       character(len=*), intent(in) :: tgt_grid_name
       character(len=*), intent(in) :: tgt_field_name
       character(len=:), ALLOCATABLE :: src_comp_name
       character(len=:), ALLOCATABLE :: src_grid_name
       character(len=:), ALLOCATABLE :: src_field_name
     end subroutine yac_fget_field_source_instance
 
  end interface yac_fget_field_source
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_timestep
 
 
     function yac_fget_timestep_from_field_id ( field_id ) result( timestep )
       integer, intent (in)           :: field_id
       character(len=:), ALLOCATABLE :: timestep
 
     end function yac_fget_timestep_from_field_id
 
     function yac_fget_field_timestep ( comp_name, grid_name, field_name ) &
          result( timestep )
 
       character(len=*), intent(in)   :: comp_name
       character(len=*), intent(in)   :: grid_name
       character(len=*), intent(in)   :: field_name
       character(len=:), ALLOCATABLE :: timestep
 
     end function  yac_fget_field_timestep
 
     function yac_fget_field_timestep_instance ( yac_instance_id, &
                                                 comp_name,       &
                                                 grid_name,       &
                                                field_name )      &
          result( timestep )
 
       integer, intent(in)            :: yac_instance_id
       character(len=*), intent(in)   :: comp_name
       character(len=*), intent(in)   :: grid_name
       character(len=*), intent(in)   :: field_name
       character(len=:), ALLOCATABLE :: timestep
 
     end function  yac_fget_field_timestep_instance
 
  end interface yac_fget_field_timestep
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_datetime
 
     function yac_fget_field_datetime ( field_id ) &
          result( datetime )
       integer, intent(in) :: field_id
       character(len=:), allocatable :: datetime
     end function yac_fget_field_datetime
 
  end interface yac_fget_field_datetime
 
  ! ---------------------------------------------------------------------
 
  interface yac_fenable_field_frac_mask
 
     subroutine yac_fenable_field_frac_mask ( &
       comp_name, grid_name, field_name, frac_mask_fallback_value )
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       double precision, intent(in) :: frac_mask_fallback_value
     end subroutine yac_fenable_field_frac_mask
 
     subroutine yac_fenable_field_frac_mask_instance ( &
       yac_instance_id, comp_name, grid_name, field_name, &
       frac_mask_fallback_value )
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       double precision, intent(in) :: frac_mask_fallback_value
     end subroutine yac_fenable_field_frac_mask_instance
 
  end interface yac_fenable_field_frac_mask
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_frac_mask_fallback_value
 
     function yac_fget_field_frac_mask_fallback_value ( &
       comp_name, grid_name, field_name )
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       double precision :: yac_fget_field_frac_mask_fallback_value
     end function yac_fget_field_frac_mask_fallback_value
 
     function yac_fget_field_frac_mask_fallback_value_instance ( &
       yac_instance_id, comp_name, grid_name, field_name )
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       double precision :: yac_fget_field_frac_mask_fallback_value_instance
     end function yac_fget_field_frac_mask_fallback_value_instance
 
  end interface yac_fget_field_frac_mask_fallback_value
 
  ! ---------------------------------------------------------------------
 
  interface yac_fget_field_collection_size
 
     function yac_fget_collection_size_from_field_id ( field_id ) &
          result(collection_size)
 
       integer, intent (in)  :: field_id
       integer :: collection_size
 
     end function  yac_fget_collection_size_from_field_id
 
     function yac_fget_field_collection_size ( comp_name, grid_name, field_name )
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       integer :: yac_fget_field_collection_size
     end function yac_fget_field_collection_size
 
     function yac_fget_field_collection_size_instance ( yac_instance_id, &
                                                        comp_name,       &
                                                        grid_name,       &
                                                        field_name )
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       integer :: yac_fget_field_collection_size_instance
     end function yac_fget_field_collection_size_instance
 
  end interface yac_fget_field_collection_size
 
  ! ---------------------------------------------------------------------
 
  interface yac_fdef_component_metadata
     subroutine yac_fdef_component_metadata(comp_name, metadata)
       character(len=*), intent(in)    :: comp_name
       character(len=*), intent(in)    :: metadata
     end subroutine yac_fdef_component_metadata
 
     subroutine yac_fdef_component_metadata_instance(yac_instance_id, comp_name, &
          metadata)
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in)    :: comp_name
       character(len=*), intent(in)    :: metadata
     end subroutine yac_fdef_component_metadata_instance
  end interface yac_fdef_component_metadata
 
  interface yac_fcomponent_has_metadata
    function yac_fcomponent_has_metadata(comp_name) result( has_metadata )
      character(len=*), intent(in) :: comp_name
      logical :: has_metadata
    end function yac_fcomponent_has_metadata
 
    function yac_fcomponent_has_metadata_instance(yac_instance_id, comp_name) &
      result( has_metadata )
      integer, intent(in) :: yac_instance_id
      character(len=*), intent(in) :: comp_name
      logical :: has_metadata
    end function yac_fcomponent_has_metadata_instance
  end interface yac_fcomponent_has_metadata
 
  interface yac_fget_component_metadata
     function yac_fget_component_metadata(comp_name) result( metadata )
       character(len=*), intent(in)  :: comp_name
       character(len=:), allocatable :: metadata
     end function yac_fget_component_metadata
 
     function yac_fget_component_metadata_instance(yac_instance_id, comp_name) &
          result( metadata )
       integer, intent(in)           :: yac_instance_id
       character(len=*), intent(in)  :: comp_name
       character(len=:), allocatable :: metadata
     end function yac_fget_component_metadata_instance
  end interface yac_fget_component_metadata
 
  interface yac_fdef_grid_metadata
     subroutine yac_fdef_grid_metadata(grid_name, metadata)
       character(len=*), intent(in)    :: grid_name
       character(len=*), intent(in)    :: metadata
     end subroutine yac_fdef_grid_metadata
 
     subroutine yac_fdef_grid_metadata_instance(yac_instance_id, grid_name, metadata)
       integer, intent(in)             :: yac_instance_id
       character(len=*), intent(in)    :: grid_name
       character(len=*), intent(in)    :: metadata
     end subroutine yac_fdef_grid_metadata_instance
  end interface yac_fdef_grid_metadata
 
  interface yac_fgrid_has_metadata
    function yac_fgrid_has_metadata(grid_name) result( has_metadata )
      character(len=*), intent(in) :: grid_name
      logical :: has_metadata
    end function yac_fgrid_has_metadata
 
    function yac_fgrid_has_metadata_instance(yac_instance_id, grid_name) &
      result( has_metadata )
      integer, intent(in) :: yac_instance_id
      character(len=*), intent(in) :: grid_name
      logical :: has_metadata
    end function yac_fgrid_has_metadata_instance
  end interface yac_fgrid_has_metadata
 
  interface yac_fget_grid_metadata
     function yac_fget_grid_metadata(grid_name) result( metadata )
       character(len=*), intent(in) :: grid_name
       character(len=:), allocatable :: metadata
     end function yac_fget_grid_metadata
 
     function yac_fget_grid_metadata_instance(yac_instance_id, grid_name) &
          result( metadata )
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: grid_name
       character(len=:), allocatable :: metadata
     end function yac_fget_grid_metadata_instance
  end interface yac_fget_grid_metadata
 
  interface yac_fdef_field_metadata
     subroutine yac_fdef_field_metadata(comp_name, grid_name, field_name, metadata)
       character(len=*), intent(in)    :: comp_name
       character(len=*), intent(in)    :: grid_name
       character(len=*), intent(in)    :: field_name
       character(len=*), intent(in)    :: metadata
     end subroutine yac_fdef_field_metadata
 
     subroutine yac_fdef_field_metadata_instance(yac_instance_id, comp_name, &
          grid_name, field_name, metadata)
       integer, intent(in)             :: yac_instance_id
       character(len=*), intent(in)    :: comp_name
       character(len=*), intent(in)    :: grid_name
       character(len=*), intent(in)    :: field_name
       character(len=*), intent(in)    :: metadata
     end subroutine yac_fdef_field_metadata_instance
  end interface yac_fdef_field_metadata
 
  interface yac_ffield_has_metadata
    function yac_ffield_has_metadata(comp_name, grid_name, field_name) &
        result( has_metadata )
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       logical :: has_metadata
    end function yac_ffield_has_metadata
 
    function yac_ffield_has_metadata_instance( &
        yac_instance_id, comp_name, grid_name, field_name) &
      result( has_metadata )
      integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       logical :: has_metadata
    end function yac_ffield_has_metadata_instance
  end interface yac_ffield_has_metadata
 
  interface yac_fget_field_metadata
     function yac_fget_field_metadata(comp_name, grid_name, field_name) &
          result( metadata )
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       character(len=:), allocatable :: metadata
     end function yac_fget_field_metadata
 
     function yac_fget_field_metadata_instance(yac_instance_id, comp_name, &
          grid_name, field_name) &
          result( metadata )
       integer, intent(in) :: yac_instance_id
       character(len=*), intent(in) :: comp_name
       character(len=*), intent(in) :: grid_name
       character(len=*), intent(in) :: field_name
       character(len=:), allocatable :: metadata
     end function yac_fget_field_metadata_instance
  end interface yac_fget_field_metadata
 
  ! ---------------------------------------------------------------------
 
  interface
 
    subroutine yac_fget_action ( field_id, action )
 
      integer, intent (in)  :: field_id
      integer, intent (out) :: action
 
    end subroutine yac_fget_action
 
  end interface
 
  ! ---------------------------------------------------------------------
 
  interface
 
    subroutine yac_fupdate ( field_id )
 
      integer, intent (in)  :: field_id
 
    end subroutine yac_fupdate
 
 end interface
 
! --- ISO_C interface -------------------------------------------------
 
  interface
 
     subroutine yac_abort_message ( text, file, line ) &
           bind( c, name='yac_abort_message' )
 
       use, intrinsic :: iso_c_binding, only : c_int, c_char
 
       character ( kind=c_char ), dimension(*) :: text
       character ( kind=c_char ), dimension(*) :: file
       integer ( kind=c_int ), value           :: line
 
     end subroutine yac_abort_message
 
  end interface
 
end module yac