diffusion_mod Module Reference

Diffusion on the TH and Q fields. More...

Functions/Subroutines
type(component_descriptor_type) function, public	diffusion_get_descriptor ()
	Provides the descriptor back to the caller and is used in component registration. More...
subroutine	field_information_retrieval_callback (current_state, name, field_information)
	Field information retrieval callback, this returns information for a specific components published field. More...
subroutine	field_value_retrieval_callback (current_state, name, field_value)
	Field value retrieval callback, this returns the value of a specific published field. More...
subroutine	initialisation_callback (current_state)
	Sets up the stencil_mod (used in interpolation) and allocates data for the flux fields. More...
subroutine	finalisation_callback (current_state)
subroutine	timestep_callback (current_state)
	At each timestep will compute the diffusion source terms for TH and Q fields per column if these fields are active. More...
subroutine	perform_q_diffusion (current_state, local_y, local_x)
	Computes the diffusion source terms for each Q field. More...
subroutine	perform_th_diffusion (current_state, local_y, local_x)
	Computes the diffusion source terms for the theta field. More...
subroutine	general_diffusion (current_state, field, source_field, local_y, local_x, diagnostics)
	General diffusion computation for any field which is provided as arguments. Works in a column. More...
subroutine	perform_local_data_copy_for_diff (current_state, halo_depth, involve_corners, source_data)
	Does local data copying for diffusion coefficient variable halo swap. More...
subroutine	copy_halo_buffer_to_diff (current_state, neighbour_description, dim, target_index, neighbour_location, current_page, source_data)
	Copies the halo buffer to halo location for the diffusion coefficient field. More...
subroutine	copy_halo_buffer_to_diff_corners (current_state, neighbour_description, corner_loc, x_target_index, y_target_index, neighbour_location, current_page, source_data)
	Copies the corner halo buffer to the diffusion coefficient field corners. More...
subroutine	save_precomponent_tendencies (current_state, cxn, cyn, txn, tyn)
	Save the 3d tendencies coming into the component. More...
subroutine	compute_component_tendencies (current_state, cxn, cyn, txn, tyn)
	Computation of component tendencies. More...
subroutine	set_published_field_value (field_value, real_1d_field, real_2d_field, real_3d_field)
	Sets the published field value from the temporary diagnostic values held by this component. More...

Variables
real(kind=default_precision), dimension(:), allocatable	th_diffusion
real(kind=default_precision), dimension(:,:), allocatable	q_diffusion
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_th
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_qv
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_ql
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_qi
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_qr
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_qs
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_qg
real(kind=default_precision), dimension(:,:,:), allocatable	tend_3d_tabs
logical	l_tend_3d_th
logical	l_tend_3d_qv
logical	l_tend_3d_ql
logical	l_tend_3d_qi
logical	l_tend_3d_qr
logical	l_tend_3d_qs
logical	l_tend_3d_qg
logical	l_tend_3d_tabs
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_th
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_qv
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_ql
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_qi
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_qr
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_qs
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_qg
real(kind=default_precision), dimension(:), allocatable	tend_pr_tot_tabs
logical	l_tend_pr_tot_th
logical	l_tend_pr_tot_qv
logical	l_tend_pr_tot_ql
logical	l_tend_pr_tot_qi
logical	l_tend_pr_tot_qr
logical	l_tend_pr_tot_qs
logical	l_tend_pr_tot_qg
logical	l_tend_pr_tot_tabs
integer	iqv =0
integer	iql =0
integer	iqr =0
integer	iqi =0
integer	iqs =0
integer	iqg =0
integer	diagnostic_generation_frequency

Detailed Description

Diffusion on the TH and Q fields.

Function/Subroutine Documentation

compute_component_tendencies()

subroutine diffusion_mod::compute_component_tendencies ( type(model_state_type), intent(inout), target  current_state, integer, intent(in)  cxn, integer, intent(in)  cyn, integer, intent(in)  txn, integer, intent(in)  tyn  )

private

Computation of component tendencies.

Parameters

current_state	Current model state
cxn	The current slice, x, index
cyn	The current column, y, index.
txn	target_x_index
tyn	target_y_index

Definition at line 626 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
    integer, intent(in) ::  cxn, cyn, txn, tyn
 
    ! Calculate change in tendency due to component
    if (l_tend_3d_th) then
      tend_3d_th(:,tyn,txn)=current_state%sth%data(:,cyn,cxn)     - tend_3d_th(:,tyn,txn)
    endif
    if (l_tend_3d_qv) then
      tend_3d_qv(:,tyn,txn)=current_state%sq(iqv)%data(:,cyn,cxn) - tend_3d_qv(:,tyn,txn)
    endif
    if (l_tend_3d_ql) then
      tend_3d_ql(:,tyn,txn)=current_state%sq(iql)%data(:,cyn,cxn) - tend_3d_ql(:,tyn,txn)
    endif
    if (l_tend_3d_qi) then
      tend_3d_qi(:,tyn,txn)=current_state%sq(iqi)%data(:,cyn,cxn) - tend_3d_qi(:,tyn,txn)
    endif
    if (l_tend_3d_qr) then
      tend_3d_qr(:,tyn,txn)=current_state%sq(iqr)%data(:,cyn,cxn) - tend_3d_qr(:,tyn,txn)
    endif
    if (l_tend_3d_qs) then
      tend_3d_qs(:,tyn,txn)=current_state%sq(iqs)%data(:,cyn,cxn) - tend_3d_qs(:,tyn,txn)
    endif
    if (l_tend_3d_qg) then
      tend_3d_qg(:,tyn,txn)=current_state%sq(iqg)%data(:,cyn,cxn) - tend_3d_qg(:,tyn,txn)
    endif
    if (l_tend_3d_tabs) then
      tend_3d_tabs(:,tyn,txn)=   &
       current_state%sth%data(:,cyn,cxn) * current_state%global_grid%configuration%vertical%rprefrcp(:)   &
        - tend_3d_tabs(:,tyn,txn)
    endif
 
   ! Add local tendency fields to the profile total
    if (l_tend_pr_tot_th) then
      tend_pr_tot_th(:)=tend_pr_tot_th(:) + tend_3d_th(:,tyn,txn)
    endif
    if (l_tend_pr_tot_qv) then
      tend_pr_tot_qv(:)=tend_pr_tot_qv(:) + tend_3d_qv(:,tyn,txn)
    endif
    if (l_tend_pr_tot_ql) then
      tend_pr_tot_ql(:)=tend_pr_tot_ql(:) + tend_3d_ql(:,tyn,txn)
    endif
    if (l_tend_pr_tot_qi) then
      tend_pr_tot_qi(:)=tend_pr_tot_qi(:) + tend_3d_qi(:,tyn,txn)
    endif
    if (l_tend_pr_tot_qr) then
      tend_pr_tot_qr(:)=tend_pr_tot_qr(:) + tend_3d_qr(:,tyn,txn)
    endif
    if (l_tend_pr_tot_qs) then
      tend_pr_tot_qs(:)=tend_pr_tot_qs(:) + tend_3d_qs(:,tyn,txn)
    endif
    if (l_tend_pr_tot_qg) then
      tend_pr_tot_qg(:)=tend_pr_tot_qg(:) + tend_3d_qg(:,tyn,txn)
    endif
    if (l_tend_pr_tot_tabs) then
      tend_pr_tot_tabs(:)=tend_pr_tot_tabs(:) + tend_3d_tabs(:,tyn,txn)
    endif
 

Here is the caller graph for this function:

copy_halo_buffer_to_diff()

subroutine diffusion_mod::copy_halo_buffer_to_diff ( type(model_state_type), intent(inout)  current_state, type(neighbour_description_type), intent(inout)  neighbour_description, integer, intent(in)  dim, integer, intent(in)  target_index, integer, intent(in)  neighbour_location, integer, dimension(:), intent(inout)  current_page, type(field_data_wrapper_type), dimension(:), intent(in), optional  source_data  )

private

Copies the halo buffer to halo location for the diffusion coefficient field.

Parameters

current_state	The current model state
neighbour_description	The halo swapping description of the neighbour we are accessing the buffer of
dim	The dimension we receive for
target_index	The target index for the dimension we are receiving for
neighbour_location	The location in the local neighbour data stores of this neighbour
current_page	The current, next, halo swap page to read from (all previous have been read and copied already)
source_data	Optional source data which is written into

Definition at line 544 of file diffusion.F90.

    type(model_state_type), intent(inout) :: current_state
    integer, intent(in) :: dim, target_index, neighbour_location
    integer, intent(inout) :: current_page(:)
    type(neighbour_description_type), intent(inout) :: neighbour_description
    type(field_data_wrapper_type), dimension(:), intent(in), optional :: source_data
 
    call copy_buffer_to_field(current_state%local_grid, neighbour_description%recv_halo_buffer, &
         current_state%diff_coefficient%data, dim, target_index, current_page(neighbour_location))
 
    current_page(neighbour_location)=current_page(neighbour_location)+1

Here is the caller graph for this function:

copy_halo_buffer_to_diff_corners()

subroutine diffusion_mod::copy_halo_buffer_to_diff_corners ( type(model_state_type), intent(inout)  current_state, type(neighbour_description_type), intent(inout)  neighbour_description, integer, intent(in)  corner_loc, integer, intent(in)  x_target_index, integer, intent(in)  y_target_index, integer, intent(in)  neighbour_location, integer, dimension(:), intent(inout)  current_page, type(field_data_wrapper_type), dimension(:), intent(in), optional  source_data  )

private

Copies the corner halo buffer to the diffusion coefficient field corners.

Parameters

current_state	The current model state
neighbour_description	The halo swapping description of the neighbour we are accessing the buffer of
corner_loc	The corner location
x_target_index	The X target index for the dimension we are receiving for
y_target_index	The Y target index for the dimension we are receiving for
neighbour_location	The location in the local neighbour data stores of this neighbour
current_page	The current, next, halo swap page to read from (all previous have been read and copied already)
source_data	Optional source data which is written into

Definition at line 567 of file diffusion.F90.

    type(model_state_type), intent(inout) :: current_state
    integer, intent(in) :: corner_loc, x_target_index, y_target_index, neighbour_location
    integer, intent(inout) :: current_page(:)
    type(neighbour_description_type), intent(inout) :: neighbour_description
    type(field_data_wrapper_type), dimension(:), intent(in), optional :: source_data
 
    call copy_buffer_to_corner(current_state%local_grid, neighbour_description%recv_corner_buffer, &
         current_state%diff_coefficient%data, corner_loc, x_target_index, y_target_index, current_page(neighbour_location))
 
    current_page(neighbour_location)=current_page(neighbour_location)+1

Here is the caller graph for this function:

diffusion_get_descriptor()

type(component_descriptor_type) function, public diffusion_mod::diffusion_get_descriptor

Provides the descriptor back to the caller and is used in component registration.

Returns

The termination check component descriptor

Definition at line 48 of file diffusion.F90.

    diffusion_get_descriptor%name="diffusion"
    diffusion_get_descriptor%version=0.1
    diffusion_get_descriptor%initialisation=>initialisation_callback
    diffusion_get_descriptor%timestep=>timestep_callback
    diffusion_get_descriptor%finalisation=>finalisation_callback
 
    diffusion_get_descriptor%field_value_retrieval=>field_value_retrieval_callback
    diffusion_get_descriptor%field_information_retrieval=>field_information_retrieval_callback
    allocate(diffusion_get_descriptor%published_fields(2+8+8))
 
    diffusion_get_descriptor%published_fields(1)="th_diffusion"
    diffusion_get_descriptor%published_fields(2)="q_diffusion"
 
    diffusion_get_descriptor%published_fields(2+1)="tend_th_diffusion_3d_local"
    diffusion_get_descriptor%published_fields(2+2)="tend_qv_diffusion_3d_local"
    diffusion_get_descriptor%published_fields(2+3)="tend_ql_diffusion_3d_local"
    diffusion_get_descriptor%published_fields(2+4)="tend_qi_diffusion_3d_local"
    diffusion_get_descriptor%published_fields(2+5)="tend_qr_diffusion_3d_local"
    diffusion_get_descriptor%published_fields(2+6)="tend_qs_diffusion_3d_local"
    diffusion_get_descriptor%published_fields(2+7)="tend_qg_diffusion_3d_local"
    diffusion_get_descriptor%published_fields(2+8)="tend_tabs_diffusion_3d_local"
 
    diffusion_get_descriptor%published_fields(2+8+1)="tend_th_diffusion_profile_total_local"
    diffusion_get_descriptor%published_fields(2+8+2)="tend_qv_diffusion_profile_total_local"
    diffusion_get_descriptor%published_fields(2+8+3)="tend_ql_diffusion_profile_total_local"
    diffusion_get_descriptor%published_fields(2+8+4)="tend_qi_diffusion_profile_total_local"
    diffusion_get_descriptor%published_fields(2+8+5)="tend_qr_diffusion_profile_total_local"
    diffusion_get_descriptor%published_fields(2+8+6)="tend_qs_diffusion_profile_total_local"
    diffusion_get_descriptor%published_fields(2+8+7)="tend_qg_diffusion_profile_total_local"
    diffusion_get_descriptor%published_fields(2+8+8)="tend_tabs_diffusion_profile_total_local"
 

Here is the call graph for this function:

field_information_retrieval_callback()

subroutine diffusion_mod::field_information_retrieval_callback ( type(model_state_type), intent(inout), target  current_state, character(len=*), intent(in)  name, type(component_field_information_type), intent(out)  field_information  )

private

Field information retrieval callback, this returns information for a specific components published field.

Parameters

current_state	Current model state
name	The name of the field to retrieve information for
field_information	Populated with information about the field

Definition at line 87 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
    character(len=*), intent(in) :: name
    type(component_field_information_type), intent(out) :: field_information
    integer :: strcomp
 
    ! Field description is the same regardless of the specific field being retrieved
    field_information%field_type=component_array_field_type
    field_information%data_type=component_double_data_type
    if (name .eq. "q_diffusion") then
      field_information%number_dimensions=2
    else
      field_information%number_dimensions=1
    end if
    field_information%dimension_sizes(1)=current_state%local_grid%size(z_index)
    if (name .eq. "q_diffusion") field_information%dimension_sizes(2)=current_state%number_q_fields
    field_information%enabled=.true.
 
    ! Field information for 3d
    strcomp=index(name, "diffusion_3d_local")
    if (strcomp .ne. 0) then
      field_information%field_type=component_array_field_type
      field_information%number_dimensions=3
      field_information%dimension_sizes(1)=current_state%local_grid%size(z_index)
      field_information%dimension_sizes(2)=current_state%local_grid%size(y_index)
      field_information%dimension_sizes(3)=current_state%local_grid%size(x_index)
      field_information%data_type=component_double_data_type
 
      if      (name .eq. "tend_th_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_th
      else if (name .eq. "tend_qv_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_qv
      else if (name .eq. "tend_ql_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_ql
      else if (name .eq. "tend_qi_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_qi
      else if (name .eq. "tend_qr_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_qr
      else if (name .eq. "tend_qs_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_qs
      else if (name .eq. "tend_qg_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_qg
      else if (name .eq. "tend_tabs_diffusion_3d_local") then
        field_information%enabled=l_tend_3d_tabs
      else
        field_information%enabled=.true.
      end if
 
    end if !end 3d check
 
    ! Field information for profiles
    strcomp=index(name, "diffusion_profile_total_local")
    if (strcomp .ne. 0) then
      field_information%field_type=component_array_field_type
      field_information%number_dimensions=1
      field_information%dimension_sizes(1)=current_state%local_grid%size(z_index)
      field_information%data_type=component_double_data_type
 
      if      (name .eq. "tend_th_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_th
      else if (name .eq. "tend_qv_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qv
      else if (name .eq. "tend_ql_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_ql
      else if (name .eq. "tend_qi_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qi
      else if (name .eq. "tend_qr_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qr
      else if (name .eq. "tend_qs_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qs
      else if (name .eq. "tend_qg_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qg
      else if (name .eq. "tend_tabs_diffusion_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_tabs
      else
        field_information%enabled=.true.
      end if
 
    end if !end profile check
 

Here is the caller graph for this function:

field_value_retrieval_callback()

subroutine diffusion_mod::field_value_retrieval_callback ( type(model_state_type), intent(inout), target  current_state, character(len=*), intent(in)  name, type(component_field_value_type), intent(out)  field_value  )

private

Field value retrieval callback, this returns the value of a specific published field.

Parameters

current_state	Current model state
name	The name of the field to retrieve the value for
field_value	Populated with the value of the field

Definition at line 173 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
    character(len=*), intent(in) :: name
    type(component_field_value_type), intent(out) :: field_value
    
    if (name .eq. "th_diffusion") then
      call set_published_field_value(field_value, real_1d_field=th_diffusion)
    else if (name .eq. "q_diffusion") then
      call set_published_field_value(field_value, real_2d_field=q_diffusion)
 
    ! 3d Tendency Fields
    else if (name .eq. "tend_th_diffusion_3d_local" .and. allocated(tend_3d_th)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_th)
    else if (name .eq. "tend_qv_diffusion_3d_local" .and. allocated(tend_3d_qv)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_qv)
    else if (name .eq. "tend_ql_diffusion_3d_local" .and. allocated(tend_3d_ql)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_ql)
    else if (name .eq. "tend_qi_diffusion_3d_local" .and. allocated(tend_3d_qi)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_qi)
    else if (name .eq. "tend_qr_diffusion_3d_local" .and. allocated(tend_3d_qr)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_qr)
    else if (name .eq. "tend_qs_diffusion_3d_local" .and. allocated(tend_3d_qs)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_qs)
    else if (name .eq. "tend_qg_diffusion_3d_local" .and. allocated(tend_3d_qg)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_qg)
    else if (name .eq. "tend_tabs_diffusion_3d_local" .and. allocated(tend_3d_tabs)) then
      call set_published_field_value(field_value, real_3d_field=tend_3d_tabs)
 
    ! Profile Tendency Fields
    else if (name .eq. "tend_th_diffusion_profile_total_local" .and. allocated(tend_pr_tot_th)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_th)
    else if (name .eq. "tend_qv_diffusion_profile_total_local" .and. allocated(tend_pr_tot_qv)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_qv)
    else if (name .eq. "tend_ql_diffusion_profile_total_local" .and. allocated(tend_pr_tot_ql)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_ql)
    else if (name .eq. "tend_qi_diffusion_profile_total_local" .and. allocated(tend_pr_tot_qi)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_qi)
    else if (name .eq. "tend_qr_diffusion_profile_total_local" .and. allocated(tend_pr_tot_qr)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_qr)
    else if (name .eq. "tend_qs_diffusion_profile_total_local" .and. allocated(tend_pr_tot_qs)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_qs)
    else if (name .eq. "tend_qg_diffusion_profile_total_local" .and. allocated(tend_pr_tot_qg)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_qg)
    else if (name .eq. "tend_tabs_diffusion_profile_total_local" .and. allocated(tend_pr_tot_tabs)) then
      call set_published_field_value(field_value, real_1d_field=tend_pr_tot_tabs)
    end if
 

Here is the call graph for this function:

Here is the caller graph for this function:

finalisation_callback()

subroutine diffusion_mod::finalisation_callback ( type(model_state_type), intent(inout), target  current_state )

private

Definition at line 342 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
 
    if (allocated(th_diffusion)) deallocate(th_diffusion)
    if (allocated(q_diffusion)) deallocate(q_diffusion)
 
    if (allocated(tend_3d_th)) deallocate(tend_3d_th)
    if (allocated(tend_3d_qv)) deallocate(tend_3d_qv)
    if (allocated(tend_3d_ql)) deallocate(tend_3d_ql)
    if (allocated(tend_3d_qi)) deallocate(tend_3d_qi)
    if (allocated(tend_3d_qr)) deallocate(tend_3d_qr)
    if (allocated(tend_3d_qs)) deallocate(tend_3d_qs)
    if (allocated(tend_3d_qg)) deallocate(tend_3d_qg)
    if (allocated(tend_3d_tabs)) deallocate(tend_3d_tabs)
 
    if (allocated(tend_pr_tot_th)) deallocate(tend_pr_tot_th)
    if (allocated(tend_pr_tot_qv)) deallocate(tend_pr_tot_qv)
    if (allocated(tend_pr_tot_ql)) deallocate(tend_pr_tot_ql)
    if (allocated(tend_pr_tot_qi)) deallocate(tend_pr_tot_qi)
    if (allocated(tend_pr_tot_qr)) deallocate(tend_pr_tot_qr)
    if (allocated(tend_pr_tot_qs)) deallocate(tend_pr_tot_qs)
    if (allocated(tend_pr_tot_qg)) deallocate(tend_pr_tot_qg)
    if (allocated(tend_pr_tot_tabs)) deallocate(tend_pr_tot_tabs)
 

Here is the caller graph for this function:

general_diffusion()

subroutine diffusion_mod::general_diffusion ( type(model_state_type), intent(inout), target  current_state, type(prognostic_field_type), intent(inout)  field, type(prognostic_field_type), intent(inout)  source_field, integer, intent(in)  local_y, integer, intent(in)  local_x, real(kind=default_precision), dimension(:), intent(inout), optional  diagnostics  )

private

General diffusion computation for any field which is provided as arguments. Works in a column.

Parameters

current_state	The current model state
field	The field to take values from, typically zth or zq(n)
source_field	The source target field to update, typically sth or sq(n)
local_y	Local Y index
local_x	Local X index

Definition at line 477 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
    type(prognostic_field_type), intent(inout) :: field, source_field
    integer, intent(in) :: local_y, local_x
    real(kind=default_precision), dimension(:), intent(inout), optional :: diagnostics
 
    real(kind=default_precision) :: term
    integer :: k
 
    do k=2, current_state%local_grid%size(z_index)
      term=current_state%global_grid%configuration%horizontal%cx2*0.25_default_precision*&
           (((current_state%diff_coefficient%data(k, local_y, local_x)+&
           current_state%diff_coefficient%data(k, local_y, local_x-1))+&
           (current_state%diff_coefficient%data(k-1, local_y, local_x)+&
           current_state%diff_coefficient%data(k-1, local_y, local_x-1)))&
           *(field%data(k, local_y, local_x-1)-field%data(k, local_y, local_x)) -&
           ((current_state%diff_coefficient%data(k, local_y, local_x+1)+&
           current_state%diff_coefficient%data(k, local_y, local_x))+&
           (current_state%diff_coefficient%data(k-1, local_y, local_x+1)+&
           current_state%diff_coefficient%data(k-1, local_y, local_x)))&
           *(field%data(k, local_y, local_x)-field%data(k, local_y, local_x+1)) )&
           +current_state%global_grid%configuration%horizontal%cy2*0.25_default_precision*(&
           ((current_state%diff_coefficient%data(k, local_y, local_x)+&
           current_state%diff_coefficient%data(k, local_y-1, local_x))+&
           (current_state%diff_coefficient%data(k-1, local_y, local_x)+&
           current_state%diff_coefficient%data(k-1, local_y-1, local_x)))&
           *(field%data(k, local_y-1, local_x)-field%data(k, local_y, local_x)) -&
           ((current_state%diff_coefficient%data(k, local_y+1, local_x)+&
           current_state%diff_coefficient%data(k, local_y, local_x))+&
           (current_state%diff_coefficient%data(k-1, local_y+1, local_x)+&
           current_state%diff_coefficient%data(k-1, local_y, local_x)))&
           *(field%data(k, local_y, local_x)-field%data(k, local_y+1, local_x)) )&
           +( current_state%global_grid%configuration%vertical%czb(k)*&
           current_state%diff_coefficient%data(k-1, local_y, local_x)*&
           (field%data(k-1, local_y, local_x)-field%data(k, local_y, local_x)))
 
      if (k .lt. current_state%local_grid%size(z_index)) then
        term=term - current_state%global_grid%configuration%vertical%cza(k)*&
             current_state%diff_coefficient%data(k, local_y, local_x)*&
             (field%data(k, local_y, local_x)-field%data(k+1, local_y, local_x))
      end if
      source_field%data(k, local_y, local_x)=source_field%data(k, local_y, local_x)+term
      if (present(diagnostics)) diagnostics(k)=term
    end do

Here is the caller graph for this function:

initialisation_callback()

subroutine diffusion_mod::initialisation_callback ( type(model_state_type), intent(inout), target  current_state )

private

Sets up the stencil_mod (used in interpolation) and allocates data for the flux fields.

Parameters

current_state

The current model state_mod

Definition at line 224 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
 
    integer :: z_size, y_size, x_size
    logical :: l_qdiag
 
    z_size=current_state%local_grid%size(z_index) + current_state%local_grid%halo_size(z_index) * 2
    y_size=current_state%local_grid%size(y_index) + current_state%local_grid%halo_size(y_index) * 2
    x_size=current_state%local_grid%size(x_index) + current_state%local_grid%halo_size(x_index) * 2
    allocate(current_state%diff_coefficient%data(z_size, y_size, x_size))
 
    z_size=current_state%global_grid%size(z_index)
    allocate(th_diffusion(z_size))
    allocate(q_diffusion(z_size, current_state%number_q_fields))
 
    ! Set tendency diagnostic logicals based on availability
    ! Need to use 3d tendencies to compute the profiles, so they will be allocated
    !      in the case where profiles are available
    l_qdiag =  (.not. current_state%passive_q .and. current_state%number_q_fields .gt. 0)
 
    l_tend_pr_tot_th  = current_state%th%active 
    l_tend_pr_tot_qv  = l_qdiag .and. current_state%number_q_fields .ge. 1
    l_tend_pr_tot_ql  = l_qdiag .and. current_state%number_q_fields .ge. 2
    l_tend_pr_tot_qi  = l_qdiag .and. current_state%number_q_fields .ge. 11
    l_tend_pr_tot_qr  = l_qdiag .and. current_state%number_q_fields .ge. 11
    l_tend_pr_tot_qs  = l_qdiag .and. current_state%number_q_fields .ge. 11
    l_tend_pr_tot_qg  = l_qdiag .and. current_state%number_q_fields .ge. 11
    l_tend_pr_tot_tabs  = l_tend_pr_tot_th
 
    l_tend_3d_th  = current_state%th%active .or. l_tend_pr_tot_th
    l_tend_3d_qv  = (l_qdiag .and. current_state%number_q_fields .ge. 1) .or. l_tend_pr_tot_qv
    l_tend_3d_ql  = (l_qdiag .and. current_state%number_q_fields .ge. 2) .or. l_tend_pr_tot_ql
    l_tend_3d_qi  = (l_qdiag .and. current_state%number_q_fields .ge. 11) .or. l_tend_pr_tot_qi
    l_tend_3d_qr  = (l_qdiag .and. current_state%number_q_fields .ge. 11) .or. l_tend_pr_tot_qr
    l_tend_3d_qs  = (l_qdiag .and. current_state%number_q_fields .ge. 11) .or. l_tend_pr_tot_qs
    l_tend_3d_qg  = (l_qdiag .and. current_state%number_q_fields .ge. 11) .or. l_tend_pr_tot_qg
    l_tend_3d_tabs = l_tend_3d_th
 
    ! Allocate 3d tendency fields upon availability
    if (l_tend_3d_th) then
      allocate( tend_3d_th(current_state%local_grid%size(z_index),  &
                           current_state%local_grid%size(y_index),  &
                           current_state%local_grid%size(x_index)   )    )
    endif
    if (l_tend_3d_qv) then
      iqv=get_q_index(standard_q_names%VAPOUR, 'diffusion')
      allocate( tend_3d_qv(current_state%local_grid%size(z_index),  &
                           current_state%local_grid%size(y_index),  &
                           current_state%local_grid%size(x_index)   )    )
    endif
    if (l_tend_3d_ql) then
      iql=get_q_index(standard_q_names%CLOUD_LIQUID_MASS, 'diffusion')
      allocate( tend_3d_ql(current_state%local_grid%size(z_index),  &
                           current_state%local_grid%size(y_index),  &
                           current_state%local_grid%size(x_index)   )    )
    endif
    if (l_tend_3d_qi) then
      iqi=get_q_index(standard_q_names%ICE_MASS, 'diffusion')
      allocate( tend_3d_qi(current_state%local_grid%size(z_index),  &
                           current_state%local_grid%size(y_index),  &
                           current_state%local_grid%size(x_index)   )    )
    endif
    if (l_tend_3d_qr) then
      iqr=get_q_index(standard_q_names%RAIN_MASS, 'diffusion')
      allocate( tend_3d_qr(current_state%local_grid%size(z_index),  &
                           current_state%local_grid%size(y_index),  &
                           current_state%local_grid%size(x_index)   )    )
    endif
    if (l_tend_3d_qs) then
      iqs=get_q_index(standard_q_names%SNOW_MASS, 'diffusion')
      allocate( tend_3d_qs(current_state%local_grid%size(z_index),  &
                           current_state%local_grid%size(y_index),  &
                           current_state%local_grid%size(x_index)   )    )
    endif
    if (l_tend_3d_qg) then
      iqg=get_q_index(standard_q_names%GRAUPEL_MASS, 'diffusion')
      allocate( tend_3d_qg(current_state%local_grid%size(z_index),  &
                           current_state%local_grid%size(y_index),  &
                           current_state%local_grid%size(x_index)   )    )
    endif
    if (l_tend_3d_tabs) then
      allocate( tend_3d_tabs(current_state%local_grid%size(z_index),  &
                             current_state%local_grid%size(y_index),  &
                             current_state%local_grid%size(x_index)   )    )
    endif
 
    ! Allocate profile tendency fields upon availability
    if (l_tend_pr_tot_th) then
      allocate( tend_pr_tot_th(current_state%local_grid%size(z_index)) )
    endif
    if (l_tend_pr_tot_qv) then
      allocate( tend_pr_tot_qv(current_state%local_grid%size(z_index)) )
    endif
    if (l_tend_pr_tot_ql) then
      allocate( tend_pr_tot_ql(current_state%local_grid%size(z_index)) )
    endif
    if (l_tend_pr_tot_qi) then
      allocate( tend_pr_tot_qi(current_state%local_grid%size(z_index)) )
    endif
    if (l_tend_pr_tot_qr) then
      allocate( tend_pr_tot_qr(current_state%local_grid%size(z_index)) )
    endif
    if (l_tend_pr_tot_qs) then
      allocate( tend_pr_tot_qs(current_state%local_grid%size(z_index)) )
    endif
    if (l_tend_pr_tot_qg) then
      allocate( tend_pr_tot_qg(current_state%local_grid%size(z_index)) )
    endif
    if (l_tend_pr_tot_tabs) then
      allocate( tend_pr_tot_tabs(current_state%local_grid%size(z_index)) )
    endif
 
    ! Save the sampling_frequency to force diagnostic calculation on select time steps
    diagnostic_generation_frequency=options_get_integer(current_state%options_database, "sampling_frequency")
 

Here is the caller graph for this function:

perform_local_data_copy_for_diff()

subroutine diffusion_mod::perform_local_data_copy_for_diff ( type(model_state_type), intent(inout)  current_state, integer, intent(in)  halo_depth, logical, intent(in)  involve_corners, type(field_data_wrapper_type), dimension(:), intent(in), optional  source_data  )

private

Does local data copying for diffusion coefficient variable halo swap.

Parameters

current_state	The current model state_mod
source_data	Optional source data which is written into

Definition at line 526 of file diffusion.F90.

    type(model_state_type), intent(inout) :: current_state
    integer, intent(in) :: halo_depth
    logical, intent(in) :: involve_corners
    type(field_data_wrapper_type), dimension(:), intent(in), optional :: source_data
 
    call perform_local_data_copy_for_field(current_state%diff_coefficient%data, current_state%local_grid, &
         current_state%parallel%my_rank, halo_depth, involve_corners)

Here is the caller graph for this function:

perform_q_diffusion()

subroutine diffusion_mod::perform_q_diffusion ( type(model_state_type), intent(inout), target  current_state, integer, intent(in)  local_y, integer, intent(in)  local_x  )

private

Computes the diffusion source terms for each Q field.

Parameters

current_state	The current model state
local_y	Local Y index
local_x	Local X index

Definition at line 433 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
    integer, intent(in) :: local_y, local_x
 
    integer :: n
 
    do n=1, current_state%number_q_fields
      call general_diffusion(current_state, current_state%zq(n), current_state%sq(n), local_y, local_x, q_diffusion(:,n))
    end do

Here is the call graph for this function:

Here is the caller graph for this function:

perform_th_diffusion()

subroutine diffusion_mod::perform_th_diffusion ( type(model_state_type), intent(inout), target  current_state, integer, intent(in)  local_y, integer, intent(in)  local_x  )

private

Computes the diffusion source terms for the theta field.

Parameters

current_state	The current model state
local_y	Local Y index
local_x	Local X index

Definition at line 448 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
    integer, intent(in) :: local_y, local_x
 
    integer :: k
    real(kind=default_precision) :: adjustment
 
    call general_diffusion(current_state, current_state%zth, current_state%sth, local_y, local_x, th_diffusion)
 
    if (current_state%use_anelastic_equations) then
      ! This code only needs to be executed if anelastic, otherwise THREF is constant and the increment calculated here is zero
      do k=2, current_state%local_grid%size(z_index)
        adjustment=(current_state%global_grid%configuration%vertical%cza(k)*&
             current_state%global_grid%configuration%vertical%dthref(k)*&
             current_state%diff_coefficient%data(k, local_y, local_x) - current_state%global_grid%configuration%vertical%czb(k)*&
             current_state%global_grid%configuration%vertical%dthref(k-1)*&
             current_state%diff_coefficient%data(k-1, local_y, local_x))
        current_state%sth%data(k, local_y, local_x)=current_state%sth%data(k, local_y, local_x)+adjustment
        th_diffusion(k)=th_diffusion(k)+adjustment
      end do
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

save_precomponent_tendencies()

subroutine diffusion_mod::save_precomponent_tendencies ( type(model_state_type), intent(in), target  current_state, integer, intent(in)  cxn, integer, intent(in)  cyn, integer, intent(in)  txn, integer, intent(in)  tyn  )

private

Save the 3d tendencies coming into the component.

Parameters

current_state	Current model state
cxn	The current slice, x, index
cyn	The current column, y, index.
txn	target_x_index
tyn	target_y_index

Definition at line 587 of file diffusion.F90.

    type(model_state_type), target, intent(in) :: current_state
    integer, intent(in) ::  cxn, cyn, txn, tyn
 
    ! Save 3d tendency fields upon request (of 3d or profiles) and availability
    if (l_tend_3d_th) then
      tend_3d_th(:,tyn,txn)=current_state%sth%data(:,cyn,cxn)
    endif
    if (l_tend_3d_qv) then
      tend_3d_qv(:,tyn,txn)=current_state%sq(iqv)%data(:,cyn,cxn)
    endif
    if (l_tend_3d_ql) then
      tend_3d_ql(:,tyn,txn)=current_state%sq(iql)%data(:,cyn,cxn)
    endif
    if (l_tend_3d_qi) then
      tend_3d_qi(:,tyn,txn)=current_state%sq(iqi)%data(:,cyn,cxn)
    endif
    if (l_tend_3d_qr) then
      tend_3d_qr(:,tyn,txn)=current_state%sq(iqr)%data(:,cyn,cxn)
    endif
    if (l_tend_3d_qs) then
      tend_3d_qs(:,tyn,txn)=current_state%sq(iqs)%data(:,cyn,cxn)
    endif
    if (l_tend_3d_qg) then
      tend_3d_qg(:,tyn,txn)=current_state%sq(iqg)%data(:,cyn,cxn)
    endif
    if (l_tend_3d_tabs) then
      tend_3d_tabs(:,tyn,txn)=current_state%sth%data(:,cyn,cxn) * current_state%global_grid%configuration%vertical%rprefrcp(:)
    endif
 

Here is the caller graph for this function:

set_published_field_value()

subroutine diffusion_mod::set_published_field_value ( type(component_field_value_type), intent(inout)  field_value, real(kind=default_precision), dimension(:), optional  real_1d_field, real(kind=default_precision), dimension(:,:), optional  real_2d_field, real(kind=default_precision), dimension(:,:,:), optional  real_3d_field  )

private

Sets the published field value from the temporary diagnostic values held by this component.

Parameters

field_value	Populated with the value of the field
real_1d_field	Optional one dimensional real of values to publish
real_2d_field	Optional two dimensional real of values to publish

Definition at line 691 of file diffusion.F90.

    type(component_field_value_type), intent(inout) :: field_value
    real(kind=default_precision), dimension(:), optional :: real_1d_field
    real(kind=default_precision), dimension(:,:), optional :: real_2d_field
    real(kind=default_precision), dimension(:,:,:), optional :: real_3d_field
 
    if (present(real_1d_field)) then
      allocate(field_value%real_1d_array(size(real_1d_field)), source=real_1d_field)
    else if (present(real_2d_field)) then
      allocate(field_value%real_2d_array(size(real_2d_field, 1), size(real_2d_field, 2)), source=real_2d_field)
    else if (present(real_3d_field)) then
      allocate(field_value%real_3d_array(size(real_3d_field, 1), size(real_3d_field, 2), size(real_3d_field, 3)), &
               source=real_3d_field)
    end if

Here is the caller graph for this function:

timestep_callback()

subroutine diffusion_mod::timestep_callback ( type(model_state_type), intent(inout), target  current_state )

private

At each timestep will compute the diffusion source terms for TH and Q fields per column if these fields are active.

Parameters

current_state

The current model state

Definition at line 370 of file diffusion.F90.

    type(model_state_type), target, intent(inout) :: current_state
 
    integer :: local_y, local_x, target_x_index, target_y_index
 
    local_y=current_state%column_local_y
    local_x=current_state%column_local_x
    target_y_index=local_y-current_state%local_grid%halo_size(y_index)
    target_x_index=local_x-current_state%local_grid%halo_size(x_index)
 
    ! Zero profile tendency totals on first instance in the sum
    if (current_state%first_timestep_column) then
      if (l_tend_pr_tot_th) then
        tend_pr_tot_th(:)=0.0_default_precision
      endif
      if (l_tend_pr_tot_qv) then
        tend_pr_tot_qv(:)=0.0_default_precision
      endif
      if (l_tend_pr_tot_ql) then
        tend_pr_tot_ql(:)=0.0_default_precision
      endif
      if (l_tend_pr_tot_qi) then
        tend_pr_tot_qi(:)=0.0_default_precision
      endif
      if (l_tend_pr_tot_qr) then
        tend_pr_tot_qr(:)=0.0_default_precision
      endif
      if (l_tend_pr_tot_qs) then
        tend_pr_tot_qs(:)=0.0_default_precision
      endif
      if (l_tend_pr_tot_qg) then
        tend_pr_tot_qg(:)=0.0_default_precision
      endif
      if (l_tend_pr_tot_tabs) then
        tend_pr_tot_tabs(:)=0.0_default_precision
      endif
    endif  ! zero totals
 
    if (.not. current_state%use_viscosity_and_diffusion .or. current_state%halo_column) return
 
    if (current_state%diffusion_halo_swap_state%swap_in_progress) then
      ! If there is a diffusion halo swap in progress then complete it
      call complete_nonblocking_halo_swap(current_state, current_state%diffusion_halo_swap_state, &
           perform_local_data_copy_for_diff, copy_halo_buffer_to_diff, copy_halo_buffer_to_diff_corners)
    end if
 
    if (mod(current_state%timestep, diagnostic_generation_frequency) == 0) then
      call save_precomponent_tendencies(current_state, local_x, local_y, target_x_index, target_y_index)
    end if
 
    if (current_state%th%active) call perform_th_diffusion(current_state, local_y, local_x)
    if (current_state%number_q_fields .gt. 0) call perform_q_diffusion(current_state, local_y, local_x)
 
    if (mod(current_state%timestep, diagnostic_generation_frequency) == 0) then
      call compute_component_tendencies(current_state, local_x, local_y, target_x_index, target_y_index)
    end if
 

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

diagnostic_generation_frequency

integer diffusion_mod::diagnostic_generation_frequency

private

Definition at line 40 of file diffusion.F90.

  integer :: diagnostic_generation_frequency

iqg

integer diffusion_mod::iqg =0

private

Definition at line 39 of file diffusion.F90.

iqi

integer diffusion_mod::iqi =0

private

Definition at line 39 of file diffusion.F90.

iql

integer diffusion_mod::iql =0

private

Definition at line 39 of file diffusion.F90.

iqr

integer diffusion_mod::iqr =0

private

Definition at line 39 of file diffusion.F90.

iqs

integer diffusion_mod::iqs =0

private

Definition at line 39 of file diffusion.F90.

iqv

integer diffusion_mod::iqv =0

private

Definition at line 39 of file diffusion.F90.

  integer :: iqv=0, iql=0, iqr=0, iqi=0, iqs=0, iqg=0

l_tend_3d_qg

logical diffusion_mod::l_tend_3d_qg

private

Definition at line 31 of file diffusion.F90.

l_tend_3d_qi

logical diffusion_mod::l_tend_3d_qi

private

Definition at line 31 of file diffusion.F90.

l_tend_3d_ql

logical diffusion_mod::l_tend_3d_ql

private

Definition at line 31 of file diffusion.F90.

l_tend_3d_qr

logical diffusion_mod::l_tend_3d_qr

private

Definition at line 31 of file diffusion.F90.

l_tend_3d_qs

logical diffusion_mod::l_tend_3d_qs

private

Definition at line 31 of file diffusion.F90.

l_tend_3d_qv

logical diffusion_mod::l_tend_3d_qv

private

Definition at line 31 of file diffusion.F90.

l_tend_3d_tabs

logical diffusion_mod::l_tend_3d_tabs

private

Definition at line 31 of file diffusion.F90.

l_tend_3d_th

logical diffusion_mod::l_tend_3d_th

private

Definition at line 31 of file diffusion.F90.

  logical :: l_tend_3d_th,l_tend_3d_qv,l_tend_3d_ql,l_tend_3d_qi,l_tend_3d_qr,l_tend_3d_qs,l_tend_3d_qg,l_tend_3d_tabs

l_tend_pr_tot_qg

logical diffusion_mod::l_tend_pr_tot_qg

private

Definition at line 36 of file diffusion.F90.

l_tend_pr_tot_qi

logical diffusion_mod::l_tend_pr_tot_qi

private

Definition at line 36 of file diffusion.F90.

l_tend_pr_tot_ql

logical diffusion_mod::l_tend_pr_tot_ql

private

Definition at line 36 of file diffusion.F90.

l_tend_pr_tot_qr

logical diffusion_mod::l_tend_pr_tot_qr

private

Definition at line 36 of file diffusion.F90.

l_tend_pr_tot_qs

logical diffusion_mod::l_tend_pr_tot_qs

private

Definition at line 36 of file diffusion.F90.

l_tend_pr_tot_qv

logical diffusion_mod::l_tend_pr_tot_qv

private

Definition at line 36 of file diffusion.F90.

l_tend_pr_tot_tabs

logical diffusion_mod::l_tend_pr_tot_tabs

private

Definition at line 36 of file diffusion.F90.

l_tend_pr_tot_th

logical diffusion_mod::l_tend_pr_tot_th

private

Definition at line 36 of file diffusion.F90.

  logical :: l_tend_pr_tot_th,l_tend_pr_tot_qv,l_tend_pr_tot_ql,l_tend_pr_tot_qi,l_tend_pr_tot_qr,  &
             l_tend_pr_tot_qs,l_tend_pr_tot_qg,l_tend_pr_tot_tabs

q_diffusion

real(kind=default_precision), dimension(:,:), allocatable diffusion_mod::q_diffusion

private

Definition at line 25 of file diffusion.F90.

  real(kind=default_precision), dimension(:,:), allocatable :: q_diffusion

tend_3d_qg

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_qg

private

Definition at line 29 of file diffusion.F90.

tend_3d_qi

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_qi

private

Definition at line 29 of file diffusion.F90.

tend_3d_ql

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_ql

private

Definition at line 29 of file diffusion.F90.

tend_3d_qr

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_qr

private

Definition at line 29 of file diffusion.F90.

tend_3d_qs

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_qs

private

Definition at line 29 of file diffusion.F90.

tend_3d_qv

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_qv

private

Definition at line 29 of file diffusion.F90.

tend_3d_tabs

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_tabs

private

Definition at line 29 of file diffusion.F90.

tend_3d_th

real(kind=default_precision), dimension(:,:,:), allocatable diffusion_mod::tend_3d_th

private

Definition at line 29 of file diffusion.F90.

  real(kind=default_precision), dimension(:,:,:), allocatable ::     &
       tend_3d_th,tend_3d_qv,tend_3d_ql,tend_3d_qi,tend_3d_qr,tend_3d_qs,tend_3d_qg,tend_3d_tabs

tend_pr_tot_qg

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_qg

private

Definition at line 33 of file diffusion.F90.

tend_pr_tot_qi

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_qi

private

Definition at line 33 of file diffusion.F90.

tend_pr_tot_ql

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_ql

private

Definition at line 33 of file diffusion.F90.

tend_pr_tot_qr

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_qr

private

Definition at line 33 of file diffusion.F90.

tend_pr_tot_qs

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_qs

private

Definition at line 33 of file diffusion.F90.

tend_pr_tot_qv

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_qv

private

Definition at line 33 of file diffusion.F90.

tend_pr_tot_tabs

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_tabs

private

Definition at line 33 of file diffusion.F90.

tend_pr_tot_th

real(kind=default_precision), dimension(:), allocatable diffusion_mod::tend_pr_tot_th

private

Definition at line 33 of file diffusion.F90.

  real(kind=default_precision), dimension(:), allocatable ::         &
       tend_pr_tot_th,tend_pr_tot_qv,tend_pr_tot_ql,tend_pr_tot_qi,tend_pr_tot_qr,tend_pr_tot_qs,tend_pr_tot_qg,  &
       tend_pr_tot_tabs

th_diffusion

real(kind=default_precision), dimension(:), allocatable diffusion_mod::th_diffusion

private

Definition at line 24 of file diffusion.F90.

  real(kind=default_precision), dimension(:), allocatable :: th_diffusion