sensitivity.f90

Go to the documentation of this file.

 
 
module sensitivity
  use simulation_component, only: simulation_component_t
  use json_file_module, only: json_file
  use case, only: case_t
  use num_types, only: rp
  use field, only: field_t
  use field_registry, only: neko_field_registry
  use json_utils, only: json_get_or_default
  implicit none
  private
 
  type, public :: sensitivity_t
 
     type(field_t), pointer :: test
 
   contains
 
     procedure, pass(this) :: init => sensitivity_init
     procedure, pass(this) :: free => sensitivity_free
     procedure, pass(this) :: update => update_sensitivity
  end type sensitivity_t
 
contains
 
  subroutine sensitivity_init(this, json, neko_case)
    class(sensitivity_t), intent(inout) :: this
    class(json_file), intent(inout) :: json
    class(case_t), intent(inout) :: neko_case
 
    ! Allocate new fields if they don't exist yet.
    if (.not. neko_field_registry%field_exists("test")) then
       call neko_field_registry%add_field(neko_case%fluid%u%dof, "test")
    end if
 
    this%test => neko_field_registry%get_field_by_name("test")
 
    call neko_case%f_out%fluid%append(this%test)
 
  end subroutine sensitivity_init
 
  subroutine sensitivity_free(this)
    class(sensitivity_t), intent(inout) :: this
 
  end subroutine sensitivity_free
 
  subroutine update_sensitivity(this, t, tstep)
    class(sensitivity_t), intent(inout) :: this
    real(kind=rp), intent(in) :: t
    integer, intent(in) :: tstep
 
    call compute_sensitivity(this)
 
  end subroutine update_sensitivity
 
  ! ========================================================================== !
  ! Each part of the simulation component is implemented as a subroutine.
 
  subroutine compute_sensitivity(this)
    class(sensitivity_t), intent(inout) :: this
 
  end subroutine compute_sensitivity
 
  ! subroutine compute_adjoint(this)
  !   class(sensitivity_t), intent(inout) :: this
 
  ! end subroutine compute_adjoint
 
end module sensitivity