ObsADT.interface.F90

Go to the documentation of this file.

! (C) Copyright 2017-2018 UCAR
!
! This software is licensed under the terms of the Apache Licence Version 2.0
! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
 
!> Fortran adt module for functions on the interface between C++ and Fortran
!  to handle observation operators
 
module ufo_adt_mod_c
 
  use fckit_configuration_module, only: fckit_configuration 
  use iso_c_binding
  use ufo_adt_mod 
  implicit none
  private
 
  ! ------------------------------------------------------------------------------
#define LISTED_TYPE ufo_adt
 
  !> Linked list interface - defines registry_t type
#include "oops/util/linkedList_i.f"
 
  !> Global registry
  type(registry_t) :: ufo_adt_registry
 
  ! ------------------------------------------------------------------------------
 
contains
 
  ! ------------------------------------------------------------------------------
  !> Linked list implementation
#include "oops/util/linkedList_c.f"
 
! ------------------------------------------------------------------------------
 
subroutine ufo_adt_setup_c(c_key_self, c_conf) bind(c,name='ufo_adt_setup_f90')
implicit none
integer(c_int), intent(inout)  :: c_key_self
type(c_ptr), value, intent(in) :: c_conf
 
type(ufo_adt), pointer :: self
type(fckit_configuration) :: f_conf
 
call ufo_adt_registry%setup(c_key_self, self)
f_conf = fckit_configuration(c_conf)
 
call self%setup(f_conf)
 
end subroutine ufo_adt_setup_c
 
! ------------------------------------------------------------------------------
 
subroutine ufo_adt_delete_c(c_key_self) bind(c,name='ufo_adt_delete_f90')
implicit none
integer(c_int), intent(inout) :: c_key_self
    
type(ufo_adt), pointer :: self
 
call ufo_adt_registry%get(c_key_self, self)
 
call self%delete()
 
call ufo_adt_registry%remove(c_key_self)
 
end subroutine ufo_adt_delete_c
 
! ------------------------------------------------------------------------------
 
subroutine ufo_adt_simobs_c(c_key_self, c_key_geovals, c_obsspace, c_nobs, c_hofx) bind(c,name='ufo_adt_simobs_f90')
 
implicit none
integer(c_int), intent(in) :: c_key_self
integer(c_int), intent(in) :: c_key_geovals
type(c_ptr), value, intent(in) :: c_obsspace
integer(c_int), intent(in) :: c_nobs
real(c_double), intent(inout) :: c_hofx(c_nobs)
 
type(ufo_adt), pointer :: self
 
call ufo_adt_registry%get(c_key_self, self)
call self%opr_simobs(c_key_geovals, c_obsspace, c_hofx)
 
end subroutine ufo_adt_simobs_c
 
! ------------------------------------------------------------------------------
 
end module ufo_adt_mod_c