prismdrv.F90

Go to the documentation of this file.

!------------------------------------------------------------------------
! Copyright 2006-2010, CERFACS, Toulouse, France.
! All rights reserved. Use is subject to OASIS4 license terms.
!------------------------------------------------------------------------
!BOP
!
!
! !MODULE: PRISMDrv
MODULE PRISMDrv

! !PUBLIC TYPES
  USE PSMILe_common
  USE PRISMDrv_Constants

  IMPLICIT NONE

! !DESCRIPTION
! This module describes the general variables used by the driver and the 
! driver's routines
!
! !REVISED HISTORY
!   Date      Programmer   Description
! ----------  ----------   -----------
! 03/10/2002  D. Declat    Creation
! 06/03/2008  J. Charles   Modifications added for the use of bundle fields
!
! EOP
!----------------------------------------------------------------------
! $Id: prismdrv.F90 2963 2011-02-17 14:52:53Z coquart $
! $Author: coquart $
!----------------------------------------------------------------------
!
! Structure for Process Management
! ================================
!
!  appli_name     : name of the application
!
!  comp_name      : name of the component
!
!  global_comp_id : global id of the component
!
!  global_rk      : global rank of the process
!
  TYPE Drv_Proc_manage

    Character(len=max_name)          :: appli_name    
    Character(len=max_name), Pointer :: comp_name(:)
    Integer, Pointer                 :: global_rank(:)
    Integer                          :: comps_per_rank
    Integer, Pointer                 :: global_comp_id(:)

  END TYPE Drv_Proc_manage
!
!----------------------------------------------------------------------
!
! Structure for Grid
! ==================
!
! The grids table
!
!  status        : Status of the grid information (see PSMILe statuses)
!
!  grid_name     : name of the grid 
!
!  grid_id       : Id of the grid for the component comp_id
!  comp_id       : Id of the component the grid belongs
! 
!  ***_units     : units of the ***
!  ***_min       : minimum of the ***
!  ***_max       : maximum of the *** 
!  ***_pointer   : *** array
!
!  vrt_stand_name 
!  vrt_formula_terms 
!
!  ig_grid_type  : type of structure for the grid
!  ig_volume_type : type of volume
!
!  pole_covered  : indicates if the pole are covered or not
!  extent        : extent of the grid in the dimension
!  size          : size of the grid (nb of points)
!  overlap       : number of overlaping points in the dimension 
!  periodic      : periodicity in the dimension 
!                  (PSMILe_true or PSMILe_false)
!  timedep       : time dependency of the dimension 
!                  (at max 3 dimensions i,j,k)
!
  TYPE Drv_Grid

    Integer                     :: grid_id
    Integer                     :: status

    Integer                     :: comp_id

    Integer                     :: lon_units
    Integer                     :: lat_units
    Integer                     :: z_units

    Integer                     :: grid_type

    Integer                     :: pole_covered     
    Integer                     :: size
    Integer, Dimension(3)       :: overlap    
    Integer, Dimension(3)       :: periodic         
    Integer, Dimension(3)       :: extent         
    Integer, Dimension(3)       :: timedep 

    Character(len=max_name)     :: grid_name

    Character(len=max_name)     :: vrt_stand_name
    Character(len=max_name)     :: vrt_formula_terms 

    Double Precision, Dimension(:), Pointer :: lon_pointer
    Double Precision, Dimension(:), Pointer :: lat_pointer
    Double Precision, Dimension(:), Pointer :: z_pointer


  END TYPE Drv_Grid

!
!======================================================================
!
! Structure for EPIO
! ==================
!
! The EPIO is the ensemble of points used in the interpolation operations
!
! The following structure defines the storage of the EPIOS and EPIOT for
! on one association source component/ target component
!
  TYPE Drv_EPIO

    Integer           :: epio_id     ! id set by the transformer (T)
    Integer           :: trans_rank  ! rank of the T process treating this epio
    Integer           :: status      ! indicates if this index has already been
                                     ! used (needed for the allocation process)
    Integer           :: src_comp_id ! source component global id
    Integer           :: tgt_comp_id ! target component global id
    Integer           :: src_process ! local rank of the source process
    Integer           :: tgt_process ! local rank of the target process
    Integer           :: src_status  ! status to inform if the src EPIO has 
                                     ! already been set or not
    Integer           :: tgt_status  ! status to inform if the tgt EPIO has 
                                     ! already been set or not
    Integer           :: src_size    ! number of source points in the EPIO
    Integer           :: src_lonlatz_size !size of source lat/lon/z compact 
                                     ! arrays (points or corners)
    Integer           :: src_nbr_corner ! number of corner for the source grid
    Integer           :: tgt_size    ! number of target points in the EPIO
    Integer           :: tgt_nbr_corner ! number of corner for the target grid
    Integer           :: weights_status ! indicates if the weights were already
                                        ! calculated or not
    Integer           :: src_coord_type ! PRISM_Real or PRISM_Double_Precision
    Integer           :: tgt_coord_type ! PRISM_Real or PRISM_Double_Precision
    INTEGER           :: max_links_map1 ! current size of link arrays
    INTEGER           :: num_links_map1 ! actual number of links for remapping
    INTEGER           :: resize_increment ! default to increase array size
    INTEGER           :: src_grid_type    ! source grid type
    INTEGER           :: gaussred_stride  ! stride in lat/lon corner array
             ! between cell minimum and maximum values.
    Integer, Pointer  :: src_mask_pointer(:) ! mask for source epio
    Integer, Pointer  :: tgt_mask_pointer(:) ! mask for source epio
    Integer, Pointer  :: index_array(:,:) ! for each EPIO target point, indices
             ! to access the compact arrays of corresponding source 
             ! point/corner latitudes, longitudes, and z
    Integer, Pointer  :: same_lat(:) ! array of integer to know if the bascule must
             ! be done or not for each target point for gaussian reduced source grid
    Integer, Pointer  :: nbsrccells_pertgtpt(:) ! for each EPIO target point, 
             ! the number of corresponding source cells (used for  conservative
             ! remapping only as this number varies for each tgt pt)
    INTEGER, POINTER  :: srcepio_add(:) ! array containing for each element of
             ! the lon/lat/z compact arrays the corresponding rank in the total
             ! number of epio source points
    INTEGER, POINTER  :: grid1_add_map1(:) ! for each 2d conservative remapping
             ! link: rank of the corresponding source point in the total
             ! number of epio source points
    INTEGER, POINTER  :: grid2_add_map1(:) ! for each 2d conservative remapping
             ! link: rank of the corresponding target point in the total
             ! number of epio target points
    REAL, POINTER     :: src_lat_pointer_real(:) ! source epio point/corner 
             ! latitude (compact array)
    Real, Pointer     :: src_lon_pointer_real(:) ! source epio point/corner 
             ! longitude (compact array)
    Real, Pointer     :: src_z_pointer_real(:)   ! source epio point/corner 
             ! position in the vertical dimension (compact array)
    Real, Pointer     :: tgt_lat_pointer_real(:) ! target epio point/corner 
             ! latitude (not compact array)
    Real, Pointer     :: tgt_lon_pointer_real(:) ! target epio point/corner 
             ! longitude (not compact array)
    Real, Pointer     :: tgt_z_pointer_real(:)   ! target epio point/corner 
             ! position in the vertical dimension (compact array)
    Double Precision, Pointer     :: src_lat_pointer_dble(:) ! see above
    Double Precision, Pointer     :: src_lon_pointer_dble(:) ! see above
    Double Precision, Pointer     :: src_z_pointer_dble(:)   ! see above
    Double Precision, Pointer     :: tgt_lat_pointer_dble(:) ! see above
    Double Precision, Pointer     :: tgt_lon_pointer_dble(:) ! see above
    Double Precision, Pointer     :: tgt_z_pointer_dble(:)   ! see above 
    Double Precision, Pointer     :: weights(:,:) ! interpolation weights
    DOUBLE PRECISION, POINTER     :: wts_map1(:,:)! 2d cons remapping weights
    DOUBLE PRECISION, POINTER     :: grid1_area(:)! area of source epio cell
              ! as provided by the component 
    DOUBLE PRECISION, POINTER     :: grid2_area(:)! area of target epio cell
              ! as calculated by the 2d conservative remappin    
  END TYPE Drv_EPIO

!  
!======================================================================
!
! Structure for files
! ===================
!
  TYPE Drv_file_struct

    CHARACTER(len=max_name)    :: file_name
    INTEGER                    :: suffix
    INTEGER                    :: file_format  
    INTEGER                    :: file_set        
    INTEGER                    :: file_threading  
    INTEGER                    :: file_pack    
    DOUBLE PRECISION           :: file_scal  
    DOUBLE PRECISION           :: file_add
    DOUBLE PRECISION           :: fill_val   
  END TYPE Drv_file_struct
!  
!======================================================================
!
! Structure for interpolations
! ============================
!
! See prism/src/lib/common_oa4/src/psmile_smioc.F90 for comments
!
  TYPE Drv_Interp

    Integer                 :: interp_id

    Integer                 :: interp_type  
    Integer, Dimension(3)   :: interp_method
    Integer                 :: nb_neighbors 

    Integer, Dimension(3)   :: arg1
    Integer, Dimension(3)   :: arg2
    Integer, Dimension(3)   :: arg3
    Integer, Dimension(3)   :: arg4
    Integer, Dimension(3)   :: arg5
    Integer, Dimension(3)   :: arg6
    Integer, Dimension(3)   :: arg7
    Type (Drv_file_struct)  :: arg10

    Character(len=max_name) :: arg9

    Double Precision        :: arg8

  END TYPE Drv_Interp

!
!======================================================================
! 
! Structure for local transformations
! ===================================
!
! transf_id     : Global id of the transformation
!
! scat_gath     : indicates if a scattering or a gathering is needed 
!                    (PSMILe_scat or PSMILe_gather or PSMILe_undef)
!
! reduc_type    : type of reduction for each dimension (PSMILe_min,
!                    PSMILe_max, PSMILe_ave, PSMILe_undef)
!
! bndl_combi_type : type of bndl_combine 
!   ->PSMIle_var_mask: in the comb, the bundle is masked with its mask
!   ->PSMILe_ext_mask: in the combination, the bundle is masked with an 
!                           ext mask coming from a file 
!   ->PSMILe_undef: no bundle combination
! bndl_combi_ext_mask_name: meaningful if PSMILe_ext_mask
! bndl_combi_ext_mask_file: meaningful if PSMILe_ext_mask 
! 
! mask_type     : type of masking 
!   ->PSMIle_var_mask: the transient is masked with its mask
!   ->PSMILe_ext_mask: the transient is masked with an 
!                           ext mask coming from a file 
!   ->PSMILe_undef: no masking 
!
! ext_mask_name : meaningful if PSMILe_ext_mask
! ext_mask_file : meaningful if PSMILe_ext_mask
! mask_value    : mask value 
!
! extrap_nnbr   : nb of neighbours in case of extrapolation
!
! extrap_wfile  : name of the file the weights for the
!                               extrapolation are stored
!
! add_scalar    : scalar to add to the field
! mult_scalar   : scalar to multiply the field
!
! src_size      : size of the source partition
! src_field     : source field partition
! 
! trans_size    : size of the transformed partition 
! trans_field   : transformed field partition
!
  TYPE Drv_Transf

    Integer                   :: transf_id

    Integer                   :: status    
    Integer                   :: request    

    Integer                   :: scat_gather

    Integer                   :: bndl_combi_type
    Integer                   :: mask_type

    Integer                   :: extrap_nnbr

    Integer, Dimension(4)     :: reduc_type

    Character(len=max_name)   :: bndl_combi_ext_mask_name
    Type (Drv_file_struct)    :: bndl_combi_ext_mask_file

    Character(len=max_name)   :: ext_mask_name
    Type (Drv_file_struct)    :: ext_mask_file

    Type (Drv_file_struct)    :: extrap_wfile

    Double Precision          :: mask_value

    Double Precision          :: add_scalar
    Double Precision          :: mult_scalar

  END TYPE Drv_Transf
  
!
!======================================================================
!
! Structure for the Fields exchange
!==================================
!
! trans_out_id    : Global id for the transient out variable
! trans_in_id     : Global id for the transient in variable
!
! epio_id         : Global id of the corresponding epio
!
! interp_id       : Global id of the corresponding interpolation 
!                   (refers to the interpolation table)
! interp_state    : state of the interpolation (UNDO, DONE)
!
! transf_id       : Global id of the corresponding transformations
!                   (refers to the transformation table)
! transf_state    : state of the interpolation (UNDO, DONE)
!
! trans_in_field_size : size of the target field
! trans_in_field      : target field
! trans_in_nbr_allocated_fields
!                     : number of fields that would fit into the
!                       memory currently allocated for
!                       trans_in_field_[int,real,dble]
! trans_in_nbr_allocated_sums
!                     : number of sums that would fit into the
!                       memory currently allocated for
!                       trans_in_field_[int,dble]
! 
  TYPE Drv_Exchange

    Integer                    :: trans_out_id
    Integer                    :: trans_in_id
    Integer                    :: epio_id

    Integer                    :: conservation

    Integer                    :: interp_id
    Integer                    :: interp_status

    Integer                    :: transf_id
    Integer                    :: transf_status

    Integer                    :: trans_in_field_size
    Integer                    :: trans_in_status
    Integer                    :: trans_in_request

    Integer                    :: trans_in_nbr_allocated_fields
    Integer                    :: trans_in_nbr_allocated_sums

    Integer                    :: trans_in_field_type  

    Integer, Pointer           :: trans_in_field_int(:)
    Real, Pointer              :: trans_in_field_real(:)
    Double Precision, Pointer  :: trans_in_field_dble(:) 

    Integer, Pointer           :: global_sum_int(:)
    Double Complex, Pointer    :: global_sum_dble(:)

  END TYPE Drv_Exchange

!
!======================================================================
!
! Structure for persistent variables 
! ==================================
!
!  local_name    : local name of the persistent variable
!
!  stand_name    : standard name of the persistent variable
!
!  comp_name     : name of the component the persistent belongs
!
!  comp_id       : component id the persistent belongs (set by the drv)
!
!  persist_type  : variable type (PSMILe_local or PSMILe_global)
!
!  units         : unit of the persistent variable
!  persist_min   : minimum value for the variable
!  persist_max   : maximum value for the variable
!
!  datatype      : numerical type of the variable 
!
!  persist       : value of the persistent variable
!
   TYPE Drv_persistent

     Integer                 :: comp_id
     Integer                 :: persist_type

     Integer                 :: units             
     Integer                 :: datatype 

     Character(len=max_name) :: local_name  
     Character(len=max_name) :: stand_name
 
     Character(len=max_name) :: comp_name   
              
     Double Precision        :: persist_min              
     Double Precision        :: persist_max
     Double Precision        :: persist 
   
   END TYPE Drv_persistent
!
!======================================================================
!----------------------------------------------------------------------
!
! Global information
!
!----------------------------------------------------------------------
!
! Name of the different applications of the simulation
  Character(LEN=max_name), Dimension(:), ALLOCATABLE :: cga_appli_name

! Name of the different application executables of the simulation
  Character(LEN=max_name), Dimension(:), ALLOCATABLE :: cga_appli_exe_name

! Argument of the different application executables of the simulation
  Character(LEN=max_name), Dimension(:), ALLOCATABLE :: cga_appli_args

! Name of the different hosts for all applications
  Character(LEN=max_name), Dimension(:), ALLOCATABLE :: cga_appli_hostname

! Name of the different components for each application
  Character(LEN=max_name), Dimension(:), ALLOCATABLE :: cga_appli_compname

! Number of hosts for each application
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_nb_hosts

! Number of processes for each application
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_nb_pes

! Number of components for each application
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_nb_comps

! Number of arguments for each application
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_nb_args

! Number of processes for each hosts of an application
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_hostnbprocs 

! Total number of rank sets for each application in the experiment 
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_nbtot_ranksets

! Total number of rank sets in the experiment 
  Integer :: ig_nbtot_ranksets

! Number of rank sets (min-max-inc) for all components for all applications
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_compnbranksets

! Array of rank sets (min-max-inc) for all components for all applications
  Integer, Dimension(:,:), ALLOCATABLE :: iga_appli_compranks

! Stdout redirected or not for all applications
  Integer, Dimension(:), ALLOCATABLE :: iga_appli_redirect

! Dates of the simulations
  TYPE(PRISM_Time_Struct) :: sga_experiment_start_date
  TYPE(PRISM_Time_Struct) :: sga_experiment_end_date
  TYPE(PRISM_Time_Struct) :: sga_run_start_date
  TYPE(PRISM_Time_Struct) :: sga_run_end_date
!
! * Dimensionning numbers
!
  Integer :: ig_driver_nb_pes   ! nb of pes for the driver
!
  Integer :: ig_nb_appl         ! nb of applications
  Integer :: ig_nb_tot_pes      ! nb of processes
  Integer :: ig_nb_tot_hosts    ! total number of hosts 
  Integer :: ig_nb_tot_comps    ! total number of comps 
  Integer :: ig_nb_tot_args     ! total number of args 

! Number of Fortran unit sets, grids, transients  and  persistents
  INTEGER :: ig_nb_tot_unitsets ! total number of Forrean unit sets
  INTEGER :: ig_nb_tot_grids    ! total number of grids
  INTEGER :: ig_nb_tot_transi   ! total number of transients
  INTEGER :: ig_nb_tot_persis   ! total number of presistents

! MPI implementation
  Integer :: ig_MPI

!
! * Communicators
!
  Integer :: comm_drv_global     ! Global communicator
  Integer :: comm_drv_local      ! Driver communicator
  Integer :: comm_drv_psmile     ! communicator that gathers psmile appli.

  Integer :: driver_rank

  Integer :: comm_drv_trans  ! Global communicator dedicated to the trans.

!
! ** comm_coupling        : Array - Coupler internal communicator for
!                           communication with individual applications
  Integer, Dimension (:), ALLOCATABLE :: comm_coupling

!
! Management of the applications name, components name and id and process rank
!
  TYPE (Drv_Proc_manage) , Dimension(:), POINTER   :: Drv_Procs

! Management of the grids information
!
! Grids table
  TYPE (Drv_Grid), Dimension(:), POINTER   :: Drv_Grids
! Dimensioning variables
  Integer :: Number_of_Grids_drv       ! Number of grids

! Dimensioning variables
  Integer :: Number_of_Epios_allocated  ! Number of allocated Epios
! EPIOS, EPIOT arrays
  TYPE (Drv_EPIO), Dimension(:), POINTER    :: Drv_Epios

!
! Management of the transients information
!
! Interpolations table
  TYPE (Drv_Interp), Dimension(:), POINTER  :: Drv_Interps
! Transformations table
  TYPE (Drv_Transf), Dimension(:), POINTER  :: Drv_Transfs
! Exchange table
  TYPE (Drv_Exchange),Dimension(:), POINTER :: Drv_Exchanges

! Dimensioning variables
  Integer :: Number_of_comms         ! Number of exchanges of Transients
  Integer :: Number_of_Interps       ! Number of interpolations 
  Integer :: Number_of_Transfs       ! Number of transformations
  Integer :: Number_of_Exchanges     ! Number of exchanges transiting
!
!======================================================================
!----------------------------------------------------------------------
! 
! Interfaces of PRISM-Driver/Transformer routines
!
!----------------------------------------------------------------------
! 

  INTERFACE

    subroutine prismdrv_init (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_init

    subroutine prismdrv_init_appl (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_init_appl

    subroutine prismdrv_def_mpi_comm (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_def_mpi_comm

    subroutine prismdrv_finalize (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_finalize

!=======

    subroutine prismdrv_set_scc_info (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_set_scc_info

!=======
    subroutine prismdrv_spawn_child ( exec, args, n_args, application_number, &
   n_hosts, hostnames, npes, intracomm, intercomm, ierror )
      USE PSMILe_common
      CHARACTER(len=*),INTENT(In)         :: exec
      INTEGER,INTENT(In)                  :: n_args
      CHARACTER(len=*),INTENT(In)         :: args (n_args)
      INTEGER,INTENT(In)                  :: application_number
      INTEGER,INTENT(In)                  :: n_hosts
      CHARACTER(len=*),INTENT(In)         :: hostnames (n_hosts)
      INTEGER,INTENT(In)                  :: npes (n_hosts)
      INTEGER, INTENT(In)                 :: intracomm
      INTEGER,INTENT(Out)                 :: intercomm
      INTEGER,INTENT(Out)                 :: ierror
    end subroutine prismdrv_spawn_child

!=======
    subroutine prismdrv_set_smioc_info (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_set_smioc_info

    subroutine prismdrv_init_smioc_struct (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_init_smioc_struct

    subroutine prismdrv_finalize_smioc_struct (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_finalize_smioc_struct

    subroutine prismdrv_get_udef_transients ( id_err )
      Integer, Intent (Out) :: id_err
    end subroutine prismdrv_get_udef_transients


!======================================================================

    subroutine prismtrs_main (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismtrs_main

!=======
 
    subroutine prismtrs_bcast2trs (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismtrs_bcast2trs

!=======
    subroutine prismtrs_conserv2d_weight(                          &
       id_epiosrc_size, id_epiotgt_size,                           &
       id_src_nbr_corner, id_tgt_nbr_corner,                       & 
       id_src_lonlatz_size,                                        &
       id_index_size1,                                             &
       dda_src_lat, dda_src_lon, dda_src_z,                        &
       dda_tgt_lat, dda_tgt_lon, dda_tgt_z,                        &
       ida_nbsrccells_pertgtpt,                                    &
       ida_index_array,                                            &
       ida_srcepio_add,                                            &
       id_epio_id,                                                 &
       id_interp_id,                                               &
       id_idim,                                                    &
       id_num_wts,                                                 & 
       id_err) 
  !
      INTEGER :: id_epiosrc_size, id_epiotgt_size
      INTEGER :: id_src_nbr_corner, id_tgt_nbr_corner
      INTEGER :: id_src_lonlatz_size
      INTEGER :: id_index_size1
      DOUBLE PRECISION, DIMENSION(id_src_lonlatz_size) :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_lonlatz_size) :: dda_src_lon 
      DOUBLE PRECISION, DIMENSION(id_src_lonlatz_size) :: dda_src_z
      DOUBLE PRECISION, DIMENSION(id_epiotgt_size,id_tgt_nbr_corner) :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_epiotgt_size,id_tgt_nbr_corner) :: dda_tgt_lon 
      DOUBLE PRECISION, DIMENSION(id_epiotgt_size,id_tgt_nbr_corner) :: dda_tgt_z
      INTEGER, DIMENSION(id_epiotgt_size)                  :: ida_nbsrccells_pertgtpt
      INTEGER, DIMENSION(id_index_size1,id_src_nbr_corner) :: ida_index_array
      INTEGER, DIMENSION(id_index_size1)                   :: ida_srcepio_add
      INTEGER :: id_epio_id     ! EPIO Id
      INTEGER :: id_interp_id   ! Interpolation Id
      INTEGER :: id_idim        ! Dimension concerned with the current interpolation
      INTEGER :: id_num_wts     ! Number of weights for 2D conservative remapping
      INTEGER :: id_err
    end subroutine prismtrs_conserv2d_weight

!=======
    subroutine prismtrs_remap_conserv (     &
       grid1_size, grid2_size,              &
       grid1_corners, grid2_corners,        &
       id_src_lonlatz_size,                 &
       id_index_size1,                      &
       grid1_corner_lat, grid1_corner_lon,  &
       grid2_corner_lat, grid2_corner_lon,  &
       ida_nbsrccells_pertgtpt,             &
       ida_index_array,                     &
       ida_srcepio_add,                     &
       id_epio_id,                          &
       id_interp_id,                        &
       id_idim,                             &
       num_wts,                             &
       id_err)
      INTEGER :: grid1_size, grid2_size
      INTEGER :: grid1_corners, grid2_corners
      INTEGER :: id_src_lonlatz_size
      INTEGER :: id_index_size1
      DOUBLE PRECISION, DIMENSION (id_src_lonlatz_size)       :: grid1_corner_lat
      DOUBLE PRECISION, DIMENSION (id_src_lonlatz_size)       :: grid1_corner_lon
      DOUBLE PRECISION, DIMENSION (grid2_size, grid2_corners) :: grid2_corner_lat
      DOUBLE PRECISION, DIMENSION (grid2_size, grid2_corners) :: grid2_corner_lon
      INTEGER, DIMENSION(grid2_size)                          :: ida_nbsrccells_pertgtpt
      INTEGER, DIMENSION(id_index_size1,grid1_corners)        :: ida_index_array
      INTEGER, DIMENSION(id_index_size1)                      :: ida_srcepio_add
      INTEGER :: id_epio_id   ! EPIO Id
      INTEGER :: id_interp_id ! Interpolation Id
      INTEGER :: id_idim      ! Dimension concerned with the current interpolation
      INTEGER :: num_wts      ! Number of weights for 2D conservative remapping
      INTEGER :: id_err
    end subroutine prismtrs_remap_conserv

!=======
 
    subroutine prismtrs_deallocate (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismtrs_deallocate

!=======

    subroutine prismtrs_resize_remap_vars(increment, id_epio_id, num_wts)
      INTEGER, INTENT(in) :: increment, id_epio_id, num_wts
    end subroutine prismtrs_resize_remap_vars

!=======

    subroutine prismtrs_intersection(location,intrsct_lat,intrsct_lon,lcoinc, &
                             beglat, beglon, endlat, endlon, begseg, &
                             lbegin, lrevers, &
                             num_srch_cells, id_nbcorners, &
                             srch_corner_lat, srch_corner_lon, srch_add)
      INTEGER, INTENT(in) :: num_srch_cells, id_nbcorners
      DOUBLE PRECISION, DIMENSION(id_nbcorners, num_srch_cells), INTENT(in) ::
         srch_corner_lat, srch_corner_lon
      INTEGER, DIMENSION(num_srch_cells), INTENT(in) :: srch_add
      LOGICAL, INTENT(in) :: lbegin, lrevers 
      DOUBLE PRECISION, INTENT(in) :: beglat, beglon, endlat, endlon
      DOUBLE PRECISION, dimension(2), intent(inout) :: begseg 
      INTEGER, INTENT(out) :: location 
      LOGICAL, INTENT(out) :: lcoinc
      DOUBLE PRECISION, intent(out) ::  intrsct_lat, intrsct_lon
    end subroutine prismtrs_intersection

!=======

    subroutine prismtrs_line_integral(lcl_weights, & 
                             in_phi1, in_phi2, theta1, theta2, &
                             grid1_lon, grid2_lon, num_wts) 
      DOUBLE PRECISION, INTENT(in) ::  in_phi1, in_phi2, theta1, theta2,  
         grid1_lon, grid2_lon 
      INTEGER, INTENT(in) :: num_wts 
      DOUBLE PRECISION, DIMENSION(2*num_wts), INTENT(out) ::  lcl_weights 
    end subroutine prismtrs_line_integral

!=======

    subroutine prismtrs_loop (id_err)
      Integer, Intent (Out) :: id_err
    end subroutine prismtrs_loop

!=======

    subroutine prismtrs_pole_intersection(location, &
                      intrsct_lat,intrsct_lon,lcoinc,lthresh, &
                      beglat, beglon, endlat, endlon, begseg, lrevers, &
                      num_srch_cells, id_nbcorners, &
                      srch_corner_lat, srch_corner_lon, srch_add)

      INTEGER, INTENT(in) :: num_srch_cells, id_nbcorners
      DOUBLE PRECISION, DIMENSION(id_nbcorners, num_srch_cells), INTENT(in) ::
           srch_corner_lat, srch_corner_lon
      INTEGER, DIMENSION(num_srch_cells), INTENT(in) :: srch_add
      DOUBLE PRECISION, INTENT(in) ::  beglat, beglon, endlat, endlon 
      DOUBLE PRECISION, DIMENSION(2), INTENT(inout) ::  begseg 
      LOGICAL, INTENT(in) :: lrevers   
!
      INTEGER, INTENT(inout) :: location 
      LOGICAL, INTENT(out) :: lcoinc
      LOGICAL, INTENT(inout) :: lthresh
      DOUBLE PRECISION, intent(out) ::  intrsct_lat, intrsct_lon
    end subroutine prismtrs_pole_intersection

!=======

    subroutine prismtrs_set_src_epio_dble(ida_loop, id_err)
      USE PSMILe_common
      INTEGER, DIMENSION(PSMILe_trans_Header_length), INTENT (IN) :: ida_loop
      INTEGER, INTENT (Out)            :: id_err 
    end subroutine prismtrs_set_src_epio_dble

    subroutine prismtrs_set_src_epio_real(ida_loop, id_err)
      USE PSMILe_common
      INTEGER, DIMENSION(PSMILe_trans_Header_length), INTENT (IN) :: ida_loop
      INTEGER, INTENT (Out)            :: id_err 
    end subroutine prismtrs_set_src_epio_real

    subroutine prismtrs_set_tgt_epio_dble(ida_loop, id_err)
      USE PSMILe_common
      INTEGER, DIMENSION(PSMILe_trans_Header_length), INTENT (IN) :: ida_loop
      INTEGER, INTENT (Out)            :: id_err 
    end subroutine prismtrs_set_tgt_epio_dble

    subroutine prismtrs_set_tgt_epio_real(ida_loop, id_err)
      USE PSMILe_common
      INTEGER, DIMENSION(PSMILe_trans_Header_length), INTENT (IN) :: ida_loop
      INTEGER, INTENT (Out)            :: id_err 
    end subroutine prismtrs_set_tgt_epio_real

!=======

    subroutine prismtrs_set_epio_trans(id_process_global_rank, id_err)
      INTEGER, INTENT (IN)             :: id_process_global_rank
      INTEGER, INTENT (Out)            :: id_err 
    end subroutine prismtrs_set_epio_trans

!=======

    subroutine prismtrs_store_link_cnsrv(add1, add2, lcl_weights, id_action, id_epio_id, num_wts)
      INTEGER, INTENT(in) :: add1, add2, id_action, id_epio_id, num_wts
      DOUBLE PRECISION, dimension(2*num_wts), intent(in) :: lcl_weights
    end subroutine prismtrs_store_link_cnsrv

!=======

    subroutine prismtrs_get_trans_rank (id_trans_rank,             &
                                        id_err)
      INTEGER, INTENT (Out)           :: id_trans_rank
      INTEGER, INTENT (Out)           :: id_err
    end subroutine prismtrs_get_trans_rank

    subroutine prismtrs_get_epio_handle (id_epio_id,               &
                                         id_err)
      INTEGER, INTENT (Out)           :: id_epio_id
      INTEGER, INTENT (Out)           :: id_err
    end subroutine prismtrs_get_epio_handle

    subroutine prismtrs_set_triple_links(id_transient_out_id,     &
                                         id_transient_in_id,      &
                                         id_epio_id,              &
                                         id_err)
      INTEGER, INTENT (In)            :: id_transient_out_id
      INTEGER, INTENT (In)            :: id_transient_in_id
      INTEGER, INTENT (IN)            :: id_epio_id
      INTEGER, INTENT (Out)           :: id_err
    end subroutine prismtrs_set_triple_links

!=======

    subroutine prismtrs_set_neighbors3d(id_process_global_rank,   &
                                        id_epio_id,               &
                                        id_epio_tgt_size,       &
                                        id_nb_neighbors,          &
                                        id_source_size,      &
                                        id_source_grid_type,      &
                                        id_err)
      INTEGER, INTENT (IN)         :: id_process_global_rank
      INTEGER, INTENT (IN)         :: id_epio_id
      INTEGER, INTENT (IN)         :: id_epio_tgt_size
      INTEGER, INTENT (IN)         :: id_nb_neighbors
      INTEGER, INTENT (IN)         :: id_source_size
      INTEGER, INTENT (IN)         :: id_source_grid_type
      INTEGER, INTENT (Out)        :: id_err
    end subroutine prismtrs_set_neighbors3d

!=======

    subroutine prismtrs_set_neighbors2d1d(id_process_global_rank,   &
                                          id_epio_id,               &
                                          id_epio_field_size,       &
                                          id_nb_neighbors,          &
                                          id_err)
      INTEGER, INTENT (IN)         :: id_process_global_rank
      INTEGER, INTENT (IN)         :: id_epio_id
      INTEGER, INTENT (IN)         :: id_epio_field_size
      INTEGER, INTENT (IN)         :: id_nb_neighbors
      INTEGER, INTENT (Out)        :: id_err 
    end subroutine prismtrs_set_neighbors2d1d

!=======

    subroutine prismtrs_mind_dble(id_process_global_rank,  &
                                  id_transient_out_id,     &
                                  id_epio_id,              &
                                  il_trans_out_size,       &      
                                  nbr_fields,              &      
                                  id_err)
      INTEGER, INTENT (In)     :: id_process_global_rank
      INTEGER, INTENT (In)     :: id_transient_out_id
      INTEGER, INTENT (IN)     :: id_epio_id
      INTEGER, INTENT (IN)     :: il_trans_out_size
      INTEGER, INTENT (IN)     :: nbr_fields
      INTEGER, INTENT (Out)    :: id_err 
    end subroutine prismtrs_mind_dble

!=======

    subroutine prismtrs_sort_add(add1, add2, weights, num_links, num_wts)
      INTEGER, INTENT(in) :: num_links, num_wts
      INTEGER, INTENT(inout), DIMENSION(num_links) :: add1, add2
      DOUBLE PRECISION, INTENT(inout), DIMENSION(num_wts, num_links) :: weights
    end subroutine prismtrs_sort_add

!=======

    subroutine prismtrs_mind_real(id_process_global_rank,  &
                                  id_transient_out_id,     &
                                  id_epio_id,              &
                                  il_trans_out_size,       &      
                                  nbr_fields,              &      
                                  id_err)
      INTEGER, INTENT (In)     :: id_process_global_rank
      INTEGER, INTENT (In)     :: id_transient_out_id
      INTEGER, INTENT (IN)     :: id_epio_id
      INTEGER, INTENT (IN)     :: il_trans_out_size
      INTEGER, INTENT (IN)     :: nbr_fields
      INTEGER, INTENT (Out)    :: id_err 
    end subroutine prismtrs_mind_real

!=======

    subroutine prismtrs_mind_int(id_process_global_rank,  &
                                 id_transient_out_id,     &
                                 id_epio_id,              &
                                 il_trans_out_size,       &      
                                 nbr_fields,              &      
                                 id_err)
      INTEGER, INTENT (In)    :: id_process_global_rank
      INTEGER, INTENT (In)    :: id_transient_out_id
      INTEGER, INTENT (IN)    :: id_epio_id
      INTEGER, INTENT (IN)    :: il_trans_out_size
      INTEGER, INTENT (IN)    :: nbr_fields
      INTEGER, INTENT (Out)   :: id_err 
    end subroutine prismtrs_mind_int

!=======

    subroutine prismtrs_interp(id_exchange_id,       &
                               id_epio_id,           &
                               id_trans_out_size,    &
                               dda_trans_out,        &
                               id_trans_in_size,     &
                               dda_trans_in,         &
                               nbr_fields,           &
                               id_err)
      INTEGER, INTENT (In)       :: id_exchange_id
      INTEGER, INTENT (In)       :: id_epio_id
      INTEGER, INTENT (In)       :: id_trans_out_size
      DOUBLE PRECISION,DIMENSION(id_trans_out_size),INTENT(IN) :: dda_trans_out
      INTEGER, INTENT (In)       :: id_trans_in_size
      DOUBLE PRECISION,DIMENSION(id_trans_in_size),INTENT(Out) :: dda_trans_in
      INTEGER, INTENT (In)       :: nbr_fields
      INTEGER, INTENT (Out)      :: id_err
    end subroutine prismtrs_interp

!=======

    subroutine prismtrs_apply_weights(il_src_size,       &
                                      dda_trans_out,     &
                      il_tgt_size,       &
                      dda_trans_in,      &
                      ida_mask,          &
                      id_nb_neighbors,   &  
                      ida_neighbors,     &
                      dda_weights,       &
                      nbr_fields,        &
                      id_err)
      INTEGER, INTENT (IN)         :: il_src_size
      INTEGER, INTENT (IN)         :: nbr_fields 
      DOUBLE PRECISION, DIMENSION(il_src_size*nbr_fields),INTENT(IN)  :: dda_trans_out
      INTEGER, INTENT (IN)         :: il_tgt_size
      DOUBLE PRECISION, DIMENSION(il_tgt_size*nbr_fields),INTENT(Out) :: dda_trans_in
      INTEGER, DIMENSION(il_tgt_size), INTENT(IN) :: ida_mask
      INTEGER, INTENT (IN)         :: id_nb_neighbors
      INTEGER, DIMENSION(il_tgt_size,id_nb_neighbors),INTENT(IN):: 
       ida_neighbors
      DOUBLE PRECISION, DIMENSION(il_tgt_size,id_nb_neighbors), INTENT(IN) :: dda_weights
      INTEGER, INTENT (Out)        :: id_err 
    end subroutine prismtrs_apply_weights

!=======

    subroutine prismtrs_apply_weights_2dcons(il_src_size,    &
                                             dda_trans_out,  &
                                             il_tgt_size,    &
                                             dda_trans_in,   &
                                             ida_mask,       &
                                             id_nb_links,    &  
                                             ida_grid1_add,  &
                                             ida_grid2_add,  &
                                             id_num_wts,     &
                                             dda_wts_map1,   &
                                             nbr_fields,     &
                                             id_err)
      INTEGER, INTENT (IN)                              :: il_src_size
      INTEGER, INTENT (IN)                              :: nbr_fields 
      DOUBLE PRECISION, DIMENSION(il_src_size*nbr_fields), INTENT(IN) :: dda_trans_out
      INTEGER, INTENT (IN)                              :: il_tgt_size
      INTEGER, DIMENSION(il_tgt_size), INTENT(IN)       :: ida_mask
      INTEGER, INTENT (IN)                              :: id_nb_links
      INTEGER, DIMENSION(id_nb_links), INTENT (IN)      :: ida_grid1_add
      INTEGER, DIMENSION(id_nb_links), INTENT (IN)      :: ida_grid2_add
      INTEGER, INTENT(in)                               :: id_num_wts
      DOUBLE PRECISION, DIMENSION(id_num_wts,id_nb_links), INTENT(IN) :: dda_wts_map1

      DOUBLE PRECISION, DIMENSION(il_tgt_size*nbr_fields), INTENT (Out) :: dda_trans_in
      INTEGER, INTENT (Out)               :: id_err   ! error value
    end subroutine prismtrs_apply_weights_2Dcons

!=======   

    subroutine prismtrs_apply_grads ( il_src_size,      &
                                      dda_trans_out,    &
                      ida_src_mask,     &
                      il_tgt_size,      &
                      dda_trans_in,     &
                      ida_tgt_mask,     &
                      id_nb_neighbors,  &  
                      ida_neighbors,    &
                      dda_weights,      &
                                      nbr_fields,       &
                      id_err)
      INTEGER, INTENT (IN)         :: il_src_size
      INTEGER, INTENT (IN)         :: nbr_fields 
      DOUBLE PRECISION, DIMENSION(il_src_size*nbr_fields),INTENT(IN) :: dda_trans_out
      INTEGER, DIMENSION(il_src_size), INTENT(IN) :: ida_src_mask
      INTEGER, INTENT (IN)         :: il_tgt_size
      DOUBLE PRECISION, DIMENSION(il_tgt_size*nbr_fields),INTENT(Out) :: dda_trans_in
      INTEGER, DIMENSION(il_tgt_size), INTENT(IN) :: ida_tgt_mask
      INTEGER, INTENT (IN)         :: id_nb_neighbors
      INTEGER, DIMENSION(il_tgt_size,id_nb_neighbors),INTENT(IN):: 
       ida_neighbors
      DOUBLE PRECISION, DIMENSION(il_tgt_size,id_nb_neighbors), INTENT(IN) :: dda_weights
      INTEGER, INTENT (Out)        :: id_err 
    end subroutine prismtrs_apply_grads
    
!=======

    subroutine prismtrs_target_real(id_process_global_rank,  &
                                    id_transient_in_id,      &
                                    id_epio_id,              &
                                    nbr_fields,              &
                                    id_err)
      INTEGER, INTENT (In)       :: id_process_global_rank
      INTEGER, INTENT (In)       :: id_transient_in_id
      INTEGER, INTENT (IN)       :: id_epio_id
      INTEGER, INTENT (IN)       :: nbr_fields      
      INTEGER, INTENT (Out)      :: id_err
    end subroutine prismtrs_target_real

    subroutine prismtrs_target_dble(id_process_global_rank,  &
                                    id_transient_in_id,      &
                                    id_epio_id,              &
                                    nbr_fields,              &
                                    id_err)
      INTEGER, INTENT (In)       :: id_process_global_rank
      INTEGER, INTENT (In)       :: id_transient_in_id
      INTEGER, INTENT (IN)       :: id_epio_id
      INTEGER, INTENT (IN)       :: nbr_fields      
      INTEGER, INTENT (Out)      :: id_err
    end subroutine prismtrs_target_dble

    subroutine prismtrs_target_int(id_process_global_rank,  &
                                   id_transient_in_id,      &
                                   id_epio_id,              &
                                   nbr_fields,              &
                                   id_err)
      INTEGER, INTENT (In)      :: id_process_global_rank
      INTEGER, INTENT (In)      :: id_transient_in_id
      INTEGER, INTENT (IN)      :: id_epio_id
      INTEGER, INTENT (IN)      :: nbr_fields      
      INTEGER, INTENT (Out)     :: id_err
    end subroutine prismtrs_target_int

!======================================================================

    subroutine prismtrs_distwght_weight_2d(id_src_size,         &
                                           dda_src_lat,         &
                                           dda_src_lon,         &
                                           ida_src_mask,        &
                                           id_tgt_size,         &
                                           dda_tgt_lat,         &
                                           dda_tgt_lon,         &
                                           ida_tgt_mask,        &
                                           id_nb_neighbors,     &
                                           ida_neighbor_index,  & 
                                           dda_weights,         &
                                           id_err)
      INTEGER, INTENT (IN)              :: id_src_size
      INTEGER, INTENT (IN)              :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lon
      INTEGER, DIMENSION(id_src_size), INTENT (IN) :: ida_src_mask
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lon
      INTEGER, DIMENSION(id_tgt_size), INTENT (IN) :: ida_tgt_mask
      INTEGER, INTENT (IN)              :: id_nb_neighbors
      INTEGER, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (IN) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (Out) :: dda_weights
      INTEGER, INTENT (Out)             :: id_err
    end subroutine prismtrs_distwght_weight_2d

    subroutine prismtrs_gauswght_weight_2d(id_src_size,         &
                                           dda_src_lat,         &
                                           dda_src_lon,         &
                                           ida_src_mask,        &
                                           id_tgt_size,         &
                                           dda_tgt_lat,         &
                                           dda_tgt_lon,         &
                                           ida_tgt_mask,        &
                                           dd_gaus_var,         &
                                           id_nb_neighbors,     &
                                           ida_neighbor_index,  & 
                                           dda_weights,         &
                                           id_err)
      INTEGER, INTENT (IN)              :: id_src_size
      INTEGER, INTENT (IN)              :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size)   :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size)   :: dda_src_lon
      INTEGER, DIMENSION(id_src_size), INTENT (IN) :: ida_src_mask
      DOUBLE PRECISION, DIMENSION(id_tgt_size)   :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size)   :: dda_tgt_lon
      INTEGER, DIMENSION(id_tgt_size), INTENT (IN) :: ida_tgt_mask
      DOUBLE PRECISION, INTENT (IN)     :: dd_gaus_var
      INTEGER, INTENT (IN)              :: id_nb_neighbors
      INTEGER, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (INOUT) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (Out) :: dda_weights
      INTEGER, INTENT (Out)             :: id_err
    end subroutine prismtrs_gauswght_weight_2d

    subroutine prismtrs_bilinear_weight_2d(id_src_size,         &
                                           dda_src_lat,         &
                       dda_src_lon,         &
                       id_tgt_size,         &
                       dda_tgt_lat,         &
                       dda_tgt_lon,         &
                       ida_tgt_mask,        &
                       ida_neighbor_index,  & 
                       dda_weights,         &
                       id_err)
      INTEGER, INTENT (IN)            :: id_src_size
      INTEGER, INTENT (IN)            :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lon
      INTEGER, DIMENSION(id_tgt_size), INTENT(IN) :: ida_tgt_mask
      INTEGER, DIMENSION(id_tgt_size,4), INTENT(IN) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,4), INTENT(Out) :: dda_weights
      INTEGER, INTENT (Out)           :: id_err
    end subroutine prismtrs_bilinear_weight_2d


    subroutine prismtrs_bicubic_weight_2d(id_src_size,         &
                                           dda_src_lat,         &
                       dda_src_lon,         &
                       id_tgt_size,         &
                       dda_tgt_lat,         &
                       dda_tgt_lon,         &
                       ida_tgt_mask,        &
                       ida_neighbor_index,  &
                                           ida_same_lat,        &
                       dda_weights,         &
                       id_err)
      INTEGER, INTENT (IN)            :: id_src_size
      INTEGER, INTENT (IN)            :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lon
      INTEGER, DIMENSION(id_tgt_size), INTENT(IN) :: ida_tgt_mask
      INTEGER, DIMENSION(id_tgt_size,16), INTENT(IN) :: ida_neighbor_index
      INTEGER, DIMENSION(id_tgt_size), INTENT(IN) :: ida_same_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size,16), INTENT(Out) :: dda_weights
      INTEGER, INTENT (Out)           :: id_err
    end subroutine prismtrs_bicubic_weight_2d


    subroutine prismtrs_bicubic_grad_2d (  id_src_size,         &
                                           dda_src_lat,         &
                       dda_src_lon,         &
                       ida_src_mask,        &
                       id_tgt_size,         &
                       dda_tgt_lat,         &
                       dda_tgt_lon,         &
                       ida_tgt_mask,        &
                       ida_neighbor_index,  & 
                       dda_weights,         &
                       id_err)
      INTEGER, INTENT (IN)            :: id_src_size
      INTEGER, INTENT (IN)            :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lon
      INTEGER, DIMENSION(id_src_size), INTENT(IN) :: ida_src_mask
      INTEGER, DIMENSION(id_tgt_size), INTENT(IN) :: ida_tgt_mask
      INTEGER, DIMENSION(id_tgt_size,16), INTENT(IN) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,16), INTENT(Out) :: dda_weights
      INTEGER, INTENT (Out)           :: id_err
    end subroutine prismtrs_bicubic_grad_2d


    subroutine prismtrs_linear_weight_for_2d1d(id_src_size,         &
                                               dda_src_z,           &
                           id_tgt_size,         &
                           dda_tgt_z,           &
                           ida_tgt_mask,        &
                           id_nb_neighbors,     & 
                           ida_neighbor_index,  &
                           dda_weights,         &
                           id_err)
      INTEGER, INTENT (IN)                     :: id_src_size
      INTEGER, INTENT (IN)                     :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_z
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_z
      INTEGER, DIMENSION(id_tgt_size), INTENT (IN) :: ida_tgt_mask
      INTEGER, INTENT (IN)                     :: id_nb_neighbors
      INTEGER, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT(INOUT) :: 
     ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size, id_nb_neighbors),INTENT (InOut)::
     dda_weights
      INTEGER, INTENT (Out)                    :: id_err
    end subroutine prismtrs_linear_weight_for_2d1d

    subroutine prismtrs_distwght_weight_2d1d(id_src_size,         &
                                             dda_src_lat,         &
                                             dda_src_lon,         &
                                             dda_src_z,           &
                                             ida_src_mask,        &
                                             id_tgt_size,         &
                                             dda_tgt_lat,         &
                                             dda_tgt_lon,         &
                                             dda_tgt_z,           &
                                             ida_tgt_mask,        &
                                             id_nb_neighbors,     &
                                             ida_neighbor_index,  & 
                                             dda_weights,         &
                                             id_err)
      INTEGER, INTENT (IN)                :: id_src_size
      INTEGER, INTENT (IN)                :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_z
      INTEGER, DIMENSION(id_src_size), INTENT (IN) :: ida_src_mask
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_z
      INTEGER, DIMENSION(id_tgt_size), INTENT (IN) :: ida_tgt_mask
      INTEGER, INTENT (IN) :: id_nb_neighbors
      INTEGER, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (InOut) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (Out) :: dda_weights
      INTEGER, INTENT (Out)               :: id_err
    end subroutine prismtrs_distwght_weight_2d1d

    subroutine prismtrs_bilinear_weight_2d1d(id_src_size,         &
                                             dda_src_lat,         &
                                             dda_src_lon,         &
                                             dda_src_z,           &
                                             id_tgt_size,         &
                                             dda_tgt_lat,         &
                                             dda_tgt_lon,         &
                                             dda_tgt_z,           &
                                             ida_neighbor_index,  & 
                                             dda_weights,         &
                                             id_err)
      INTEGER, INTENT (IN)                     :: id_src_size
      INTEGER, INTENT (IN)                     :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_z
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_z
      INTEGER, DIMENSION(id_tgt_size,10), INTENT(INOUT) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,8), INTENT (Out) :: dda_weights
      INTEGER, INTENT (Out)               :: id_err
    end subroutine prismtrs_bilinear_weight_2d1d


    subroutine prismtrs_distwght_weight_3d(id_src_size,         &
                                           dda_src_lat,         &
                                           dda_src_lon,         &
                                           dda_src_z,           &
                                           ida_src_mask,        &
                                           id_tgt_size,         &
                                           dda_tgt_lat,         &
                                           dda_tgt_lon,         &
                                           dda_tgt_z,           &
                                           ida_tgt_mask,        &
                                           id_nb_neighbors,     &
                                           ida_neighbor_index,  & 
                                           dda_weights,         &
                                           id_err)
      INTEGER, INTENT (IN)              :: id_src_size
      INTEGER, INTENT (IN)              :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_z
      INTEGER, DIMENSION(id_src_size), INTENT (IN) :: ida_src_mask
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_z
      INTEGER, DIMENSION(id_tgt_size), INTENT (IN) :: ida_tgt_mask
      INTEGER, INTENT (IN)              :: id_nb_neighbors
      INTEGER, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (IN) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (Out) :: dda_weights
      INTEGER, INTENT (Out)             :: id_err
    end subroutine prismtrs_distwght_weight_3d

    subroutine prismtrs_gauswght_weight_3d(id_src_size,         &
                                           dda_src_lat,         &
                                           dda_src_lon,         &
                                           dda_src_z,           &
                                           ida_src_mask,        &
                                           id_tgt_size,         &
                                           dda_tgt_lat,         &
                                           dda_tgt_lon,         &
                                           dda_tgt_z,           &
                                           ida_tgt_mask,        &
                                           dd_gaus_var,         &
                                           id_nb_neighbors,     &
                                           ida_neighbor_index,  & 
                                           dda_weights,         &
                                           id_err)
      INTEGER, INTENT (IN)              :: id_src_size
      INTEGER, INTENT (IN)              :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_src_size)  :: dda_src_z
      INTEGER, DIMENSION(id_src_size), INTENT (IN) :: ida_src_mask
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size)  :: dda_tgt_z
      INTEGER, DIMENSION(id_tgt_size), INTENT (IN) :: ida_tgt_mask
      DOUBLE PRECISION, INTENT (IN)     :: dd_gaus_var
      INTEGER, INTENT (IN)              :: id_nb_neighbors
      INTEGER, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (IN) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,id_nb_neighbors), INTENT (Out) :: dda_weights
      INTEGER, INTENT (Out)             :: id_err
    end subroutine prismtrs_gauswght_weight_3d

    subroutine prismtrs_trilinear_weight(id_src_size,         &
                                         dda_src_lat,         &
                                         dda_src_lon,         &
                                         dda_src_z,           &
                                         id_tgt_size,         &
                                         dda_tgt_lat,         &
                                         dda_tgt_lon,         &
                                         dda_tgt_z,           &
                                         ida_tgt_mask,        &
                                         ida_neighbor_index,  & 
                                         dda_weights,         &
                                         id_err)
      INTEGER, INTENT (IN)                     :: id_src_size
      INTEGER, INTENT (IN)                     :: id_tgt_size
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lat
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_lon
      DOUBLE PRECISION, DIMENSION(id_src_size) :: dda_src_z
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lat
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_lon
      DOUBLE PRECISION, DIMENSION(id_tgt_size) :: dda_tgt_z
      INTEGER, DIMENSION(id_tgt_size), INTENT(IN)   :: ida_tgt_mask
      INTEGER, DIMENSION(id_tgt_size,8), INTENT(IN) :: ida_neighbor_index
      DOUBLE PRECISION, DIMENSION(id_tgt_size,8), INTENT(Out) :: dda_weights
      INTEGER, INTENT (Out)           :: id_err
    end subroutine prismtrs_trilinear_weight

!======================================================================

    subroutine prismdrv_get_smioc_file_name(id_comp_id,       &
                                            cda_file_name,    &
                                            id_file_name_len, &
                                            id_err)
      INTEGER, INTENT (In)                :: id_comp_id
      CHARACTER(LEN=256), INTENT (Out)    :: cda_file_name
      INTEGER, INTENT (Out)               :: id_file_name_len
      INTEGER, INTENT (Out)               :: id_err
    end subroutine prismdrv_get_smioc_file_name

!======================================================================

subroutine prismtrs_enqueue_in_field_dble(field, field_size, exchange_id, ierror)
   integer, intent(in)                   :: field_size
   double precision, intent(in)          :: field(field_size)
   integer, intent(in)                   :: exchange_id
   integer, intent(out)                  :: ierror
end subroutine prismtrs_enqueue_in_field_dble

subroutine prismtrs_enqueue_in_field_real(field, field_size, exchange_id, ierror)
   integer, intent(in)                   :: field_size
   real, intent(in)                      :: field(field_size)
   integer, intent(in)                   :: exchange_id
   integer, intent(out)                  :: ierror
end subroutine prismtrs_enqueue_in_field_real

subroutine prismtrs_enqueue_in_field_int(field, field_size, exchange_id, ierror)
   integer, intent(in)                   :: field_size
   integer, intent(in)                   :: field(field_size)
   integer, intent(in)                   :: exchange_id
   integer, intent(out)                  :: ierror
end subroutine prismtrs_enqueue_in_field_int

subroutine prismtrs_dequeue_in_field_dble(field_size, exchange_id, ierror)
   integer, intent(in)       :: field_size
   integer, intent(in)       :: exchange_id
   integer, intent(out)      :: ierror
end subroutine prismtrs_dequeue_in_field_dble

subroutine prismtrs_dequeue_in_field_real(field_size, exchange_id, ierror)
   integer, intent(in)       :: field_size
   integer, intent(in)       :: exchange_id
   integer, intent(out)      :: ierror
end subroutine prismtrs_dequeue_in_field_real

subroutine prismtrs_dequeue_in_field_int(field_size, exchange_id, ierror)
   integer, intent(in)       :: field_size
   integer, intent(in)       :: exchange_id
   integer, intent(out)      :: ierror
end subroutine prismtrs_dequeue_in_field_int

subroutine prismtrs_enqueue_glob_sum_dble(global_sum, nbr_fields, exchange_id, ierror)
   integer, intent(in)          :: nbr_fields
   double complex, intent(in)   :: global_sum(nbr_fields)
   integer, intent(in)          :: exchange_id
   integer, intent(out)         :: ierror
end subroutine prismtrs_enqueue_glob_sum_dble

subroutine prismtrs_enqueue_glob_sum_int(global_sum, nbr_fields, exchange_id, ierror)
   integer, intent(in)  :: nbr_fields
   integer, intent(in)  :: global_sum(nbr_fields) 
   integer, intent(in)  :: exchange_id
   integer, intent(out) :: ierror
end subroutine prismtrs_enqueue_glob_sum_int

subroutine prismtrs_dequeue_glob_sum_dble(nbr_fields, exchange_id, ierror)
   integer, intent(in)  :: nbr_fields
   integer, intent(in)  :: exchange_id
   integer, intent(out) :: ierror
end subroutine prismtrs_dequeue_glob_sum_dble

subroutine prismtrs_dequeue_glob_sum_int(nbr_fields, exchange_id, ierror)
   integer, intent(in)  :: nbr_fields
   integer, intent(in)  :: exchange_id
   integer, intent(out) :: ierror
end subroutine prismtrs_dequeue_glob_sum_int

subroutine prismdrv_init_drv_exchange(exchange_id)
   integer, intent(in) :: exchange_id
end subroutine prismdrv_init_drv_exchange

!======================================================================
  END INTERFACE

END MODULE PRISMDrv