test_stencil_mod Module Reference

Functions/Subroutines
subroutine	test_create_stencil
subroutine	test_calculate_interpolated_cell_value

Function/Subroutine Documentation

test_calculate_interpolated_cell_value()

subroutine test_stencil_mod::test_calculate_interpolated_cell_value

Definition at line 49 of file test_stencil.F90.

    type(model_state_type), target :: current_state
    type(grid_stencil_type) :: star_stencil
    type(prognostic_field_ptr_type), dimension(3)  :: fields
    integer, dimension(3, 2) :: sizes
    type(prognostic_field_type), dimension(:), allocatable :: interpolated_fields
    integer :: i,num_fields,max_y_point
    logical, dimension(3) :: interpolate_in_dimension
    real(kind=default_precision) :: r_value
    real(kind=default_precision) :: e_value
 
    
    current_state%global_grid%size(z_index) = 2
    num_fields = 15
    current_state%local_grid%size(z_index)= 0
    current_state%local_grid%size(y_index) = 3
    current_state%local_grid%halo_size(y_index) = 4
 
    fields(num_fields)%ptr => current_state%u
    sizes(num_fields,:) = (/ 2, 2 /)
    
    allocate(interpolated_fields(0:num_fields))    
    allocate(interpolated_fields(num_fields)%data(current_state%global_grid%size(z_index),&
         -1:3, -1:3))
    max_y_point = (current_state%local_grid%size(y_index)+ &
         current_state%local_grid%halo_size(y_index) *2)-1
    
    interpolated_fields(num_fields)%active=.true.
    star_stencil = create_stencil(current_state%local_grid, fields, num_fields, 3, &
         sizes,.true., .false.)
    e_value = 10_default_precision
    star_stencil%interpolated_fields(1,1)%data(1,1,1) = e_value
   
    call calculate_interpolated_cell_value(current_state%local_grid, 1, 1, (/.true.,.true.,.true./),&
         star_stencil%fields(1), star_stencil%interpolated_fields(1,1)%data(:,1,1))  
    
    r_value =  star_stencil%interpolated_fields(1,1)%data(1,1,1)
    
    call assert_equals(e_value, r_value,"Test")

Here is the call graph for this function:

Here is the caller graph for this function:

test_create_stencil()

subroutine test_stencil_mod::test_create_stencil

Definition at line 16 of file test_stencil.F90.

    type(model_state_type), target :: current_state
    type(grid_stencil_type) :: star_stencil
    type(prognostic_field_ptr_type), dimension(3)  :: fields
    integer, dimension(3, 2) :: sizes
    type(prognostic_field_type), dimension(:), allocatable :: interpolated_fields
    integer :: i,num_fields,max_y_point
    
    call assert_equals(0,star_stencil%nfields,"Test a stencil has been created")
    current_state%global_grid%size(z_index) = 2
    num_fields = 15
    current_state%local_grid%size(y_index) = 3
    current_state%local_grid%halo_size(y_index) = 4
 
    fields(num_fields)%ptr => current_state%u
    sizes(num_fields,:) = (/ 2, 2 /)
    
    allocate(interpolated_fields(0:num_fields))    
    allocate(interpolated_fields(num_fields)%data(current_state%global_grid%size(z_index),&
         -1:3, -1:3))
    max_y_point = (current_state%local_grid%size(y_index)+ &
         current_state%local_grid%halo_size(y_index) *2)-1
    
    interpolated_fields(num_fields)%active=.true.
    star_stencil = create_stencil(current_state%local_grid, fields, num_fields, 3, &
         sizes,.true., .false.)
    call assert_equals(max_y_point,star_stencil%max_y_point,"Test max Y point within the stencil")
    call assert_equals(-1,star_stencil%max_x_point,"Test max X point within the stencil")
    call assert_equals(num_fields+1,size(interpolated_fields),"Test size of interpolated_fields")
    call assert_equals(num_fields,star_stencil%nfields,"Test a stencil has been created")

Here is the call graph for this function:

Here is the caller graph for this function: