LinearGetValues.interface.F90

Go to the documentation of this file.

! (C) Copyright 2020 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.
 
! --------------------------------------------------------------------------------------------------
 
module mpas_lineargetvalues_interface_mod
 
! intrinsic
use iso_c_binding
 
! fckit dependencies
use fckit_configuration_module, only: fckit_configuration
 
! oops dependencies
use datetime_mod
use kinds, only: kind_real
 
! ufo dependencies
use ufo_locations_mod
use ufo_geovals_mod
use ufo_geovals_mod_c, only: ufo_geovals_registry
 
! self dependency
use mpasjedi_lineargetvalues_mod, only: mpasjedi_lineargetvalues, mpas_lineargetvalues_registry
 
! mpas-jedi dependencies
use mpas_geom_mod, only: mpas_geom, mpas_geom_registry
use mpas_fields_mod, only: mpas_fields, mpas_fields_registry
 
implicit none
private
 
! --------------------------------------------------------------------------------------------------
 
contains
 
! -------------------------------------------------------------------------------------------------
 
subroutine mpas_lineargetvalues_create_c(c_key_self, c_key_geom, c_locs, c_conf) &
           bind(c,name='mpas_lineargetvalues_create_f90')
 
integer(c_int),     intent(inout) :: c_key_self    !< Key to self
integer(c_int),     intent(in)    :: c_key_geom    !< Key to geometry
type(c_ptr), value, intent(in)    :: c_locs        !< Observation locations
type(c_ptr), value, intent(in)    :: c_conf        !< Key to configuration
 
type(mpasjedi_lineargetvalues), pointer :: self
type(mpas_geom),                pointer :: geom
type(ufo_locations)                     :: locs
type(fckit_configuration)               :: f_conf
 
! Create object
call mpas_lineargetvalues_registry%init()
call mpas_lineargetvalues_registry%add(c_key_self)
call mpas_lineargetvalues_registry%get(c_key_self, self)
 
! Others
call mpas_geom_registry%get(c_key_geom, geom)
locs = ufo_locations(c_locs)
f_conf = fckit_configuration(c_conf)
 
! Call method
call self%create(geom, locs, f_conf)
 
end subroutine mpas_lineargetvalues_create_c
 
! --------------------------------------------------------------------------------------------------
 
subroutine mpas_lineargetvalues_delete_c(c_key_self) &
           bind(c,name='mpas_lineargetvalues_delete_f90')
 
integer(c_int), intent(inout) :: c_key_self !< Key to self
 
type(mpasjedi_lineargetvalues), pointer :: self
 
! Get object
call mpas_lineargetvalues_registry%get(c_key_self, self)
 
! Call method
call self%delete()
 
! Remove object
call mpas_lineargetvalues_registry%remove(c_key_self)
 
end subroutine mpas_lineargetvalues_delete_c
 
! --------------------------------------------------------------------------------------------------
 
subroutine mpas_lineargetvalues_set_trajectory_c(c_key_self, c_key_geom, c_key_state, c_t1, &
                                                 c_t2, c_locs, c_key_geovals) &
           bind(c,name='mpas_lineargetvalues_set_trajectory_f90')
 
integer(c_int),     intent(in) :: c_key_self
integer(c_int),     intent(in) :: c_key_geom
integer(c_int),     intent(in) :: c_key_state
type(c_ptr), value, intent(in) :: c_t1
type(c_ptr), value, intent(in) :: c_t2
type(c_ptr), value, intent(in) :: c_locs
integer(c_int),     intent(in) :: c_key_geovals
 
type(mpasjedi_lineargetvalues), pointer :: self
type(mpas_geom),                pointer :: geom
type(mpas_fields),              pointer :: fields
type(datetime)                          :: t1
type(datetime)                          :: t2
type(ufo_locations)                     :: locs
type(ufo_geovals),              pointer :: geovals
 
! Get objects
call mpas_lineargetvalues_registry%get(c_key_self, self)
call mpas_geom_registry%get(c_key_geom, geom)
call mpas_fields_registry%get(c_key_state, fields)
call c_f_datetime(c_t1, t1)
call c_f_datetime(c_t2, t2)
locs = ufo_locations(c_locs)
call ufo_geovals_registry%get(c_key_geovals, geovals)
 
! Call method
call self%set_trajectory(geom, fields, t1, t2, locs, geovals)
 
end subroutine mpas_lineargetvalues_set_trajectory_c
 
 
! --------------------------------------------------------------------------------------------------
 
subroutine mpas_lineargetvalues_fill_geovals_tl_c(c_key_self, c_key_geom, c_key_inc, c_t1, &
                                                  c_t2, c_locs, c_key_geovals) &
           bind(c,name='mpas_lineargetvalues_fill_geovals_tl_f90')
 
integer(c_int),     intent(in) :: c_key_self
integer(c_int),     intent(in) :: c_key_geom
integer(c_int),     intent(in) :: c_key_inc
type(c_ptr), value, intent(in) :: c_t1
type(c_ptr), value, intent(in) :: c_t2
type(c_ptr), value, intent(in) :: c_locs
integer(c_int),     intent(in) :: c_key_geovals
 
type(mpasjedi_lineargetvalues), pointer :: self
type(mpas_geom),                pointer :: geom
type(mpas_fields),              pointer :: fields
type(datetime)                          :: t1
type(datetime)                          :: t2
type(ufo_locations)                     :: locs
type(ufo_geovals),              pointer :: geovals
 
! Get objects
call mpas_lineargetvalues_registry%get(c_key_self, self)
call mpas_geom_registry%get(c_key_geom, geom)
call mpas_fields_registry%get(c_key_inc, fields)
call c_f_datetime(c_t1, t1)
call c_f_datetime(c_t2, t2)
locs = ufo_locations(c_locs)
call ufo_geovals_registry%get(c_key_geovals, geovals)
 
! Call method
call self%fill_geovals_tl(geom, fields, t1, t2, locs, geovals)
 
end subroutine mpas_lineargetvalues_fill_geovals_tl_c
 
! --------------------------------------------------------------------------------------------------
 
subroutine mpas_lineargetvalues_fill_geovals_ad_c(c_key_self, c_key_geom, c_key_inc, c_t1, &
                                                  c_t2, c_locs, c_key_geovals) &
           bind(c,name='mpas_lineargetvalues_fill_geovals_ad_f90')
 
integer(c_int),     intent(in) :: c_key_self
integer(c_int),     intent(in) :: c_key_geom
integer(c_int),     intent(in) :: c_key_inc
type(c_ptr), value, intent(in) :: c_t1
type(c_ptr), value, intent(in) :: c_t2
type(c_ptr), value, intent(in) :: c_locs
integer(c_int),     intent(in) :: c_key_geovals
 
type(mpasjedi_lineargetvalues), pointer :: self
type(mpas_geom),                pointer :: geom
type(mpas_fields),              pointer :: fields
type(datetime)                          :: t1
type(datetime)                          :: t2
type(ufo_locations)                     :: locs
type(ufo_geovals),              pointer :: geovals
 
! Get objects
call mpas_lineargetvalues_registry%get(c_key_self, self)
call mpas_geom_registry%get(c_key_geom, geom)
call mpas_fields_registry%get(c_key_inc, fields)
call c_f_datetime(c_t1, t1)
call c_f_datetime(c_t2, t2)
locs = ufo_locations(c_locs)
call ufo_geovals_registry%get(c_key_geovals, geovals)
 
! Call method
call self%fill_geovals_ad(geom, fields, t1, t2, locs, geovals)
 
end subroutine mpas_lineargetvalues_fill_geovals_ad_c
 
! --------------------------------------------------------------------------------------------------
 
end module mpas_lineargetvalues_interface_mod