mean_spectrum Subroutine

  1. mod_functions.f90
  2. mod_functions
  3. mean_spectrum

public subroutine mean_spectrum(data, spectrum, dim_v, dim_y, dim_x)

Compute the MEAN spectrum of a cube along the spatial axis

Arguments

Type IntentOptional AttributesName
real(kind=xp), intent(in), dimension(:,:,:), allocatable:: data

initial fits data
real(kind=xp), intent(inout), dimension(:), allocatable:: spectrum

std_spectrum of the observation
integer, intent(in) :: dim_v

dimension along v axis
integer, intent(in) :: dim_y

dimension along spatial axis y

integer, intent(in) :: dim_x

dimension along spatial axis x

Calls

proc~~mean_spectrum~~CallsGraph proc~mean_spectrum mean_spectrum proc~mean_2d mean_2D proc~mean_spectrum->proc~mean_2d proc~ravel_2d ravel_2D proc~mean_2d->proc~ravel_2d

Contents

Source Code

mean_spectrum

Source Code

  subroutine mean_spectrum(data, spectrum, dim_v, dim_y, dim_x)
    !! Compute the MEAN spectrum of a cube along the spatial axis
    implicit none
    
    real(xp), intent(in), dimension(:,:,:), allocatable :: data !! initial fits data
    integer, intent(in) :: dim_v !! dimension along v axis
    integer, intent(in) :: dim_y !! dimension along spatial axis y 
    integer, intent(in) :: dim_x !! dimension along spatial axis x

    real(xp), intent(inout), dimension(:), allocatable :: spectrum !! std_spectrum of the observation
    real(xp), dimension(:,:), allocatable :: map !! 2D array

    integer :: i !! loop index

    do i=1,dim_v
       allocate(map(dim_y,dim_x))
       map = data(i,:,:)
       spectrum(i) = mean_2D(map, dim_y, dim_x)
       deallocate(map)
    end do
    
  end subroutine mean_spectrum