stepfields_mod Module Reference

Does the field stepping Stepping is called at the end of processing a column and steps the x-2 column. More...

Functions/Subroutines
type(component_descriptor_type) function, public	stepfields_get_descriptor ()
	Provides the descriptor back to the caller and is used in component registration. More...
subroutine	initialisation_callback (current_state)
	Initialisation callback. More...
subroutine	finalisation_callback (current_state)
	Finalisation callback. More...
subroutine	timestep_callback (current_state)
	Called at each timestep and will perform swapping and smoothing as required. More...
subroutine	step_all_fields (current_state)
	Steps all fields. More...
subroutine	determine_local_flow_minmax (current_state, local_y, local_x)
	Determines the minimum and maximum values for the local flow field. These are before the stepping, and are all reduced later on in the cfl test. More...
subroutine	reset_local_minmax_values (current_state)
	Resets the local min and max values for the flow fields. More...
subroutine	remove_negative_rounding_errors_for_single_field (x_local_index, y_local_index, x_prev, y_prev, field, local_grid)
	Removes the negative rounding errors from a specific single field. This works two columns behind and then catches up on the last column. More...
subroutine	remove_negative_rounding_errors_in_slice (y_local_index, x_prev, y_prev, field, local_grid)
	Removes the negative rounding errors from a slice of a single field. This works two columns behind and then catches up on the last column. More...
subroutine	step_single_field (x_local_index, y_local_index, x_prev, y_prev, field, zfield, sfield, local_grid, flow_field, direction, dtm, gal, c1, c2, do_timesmoothing, sav)
	Steps a single specific field. This will step on the yth column of the x-2 slice and x-1 and x if this is the last slice. More...
subroutine	perform_timesmooth_for_field (field, zfield, local_grid, x_index, y_index, c1, c2)
	Performs initial timesmoothing for a theta or Q field using Robert filter. This is finished off in swapsmooth. More...
subroutine	step_column_in_slice (y_local_index, x_prev, y_prev, field, zfield, sfield, local_grid, flow_field, direction, dtm, gal, c1, c2, do_timesmoothing, sav)
	Will step a column in a specific slice. If y_prev is large enough then will step the y-1 column and if this is the last column of the slice then will also step the current column. More...
subroutine	step_field (x_local_index, y_local_index, field, zfield, sfield, local_grid, flow_field, direction, dtm, gal, sav)
	Will do the actual field stepping. More...
subroutine	compute_component_tendencies (current_state, cxn, cyn, txn, tyn)
	Computation of component tendencies. More...
subroutine	field_information_retrieval_callback (current_state, name, field_information)
	Field information retrieval callback, this returns information for a specific component's 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	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
logical	determine_flow_minmax =.false.
logical	cfl_is_enabled
real(kind=default_precision), dimension(:), allocatable	resetq_min
logical	l_nonconservative_positive_q =.true.
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

Does the field stepping Stepping is called at the end of processing a column and steps the x-2 column.

Function/Subroutine Documentation

compute_component_tendencies()

subroutine stepfields_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 628 of file stepfields.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)
    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
 
   ! 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:

determine_local_flow_minmax()

subroutine stepfields_mod::determine_local_flow_minmax ( type(model_state_type), intent(inout), target  current_state, integer, intent(in)  local_y, integer, intent(in)  local_x  )

private

Determines the minimum and maximum values for the local flow field. These are before the stepping, and are all reduced later on in the cfl test.

Parameters

current_state	The current model state
local_y	The local y index
local_x	The local x index

Definition at line 365 of file stepfields.F90.

    type(model_state_type), target, intent(inout) :: current_state
    integer, intent(in) :: local_x, local_y
    
    integer :: k
    
    do k=2, current_state%local_grid%local_domain_end_index(z_index)
#ifdef U_ACTIVE
      current_state%local_zumax = max(current_state%local_zumax, current_state%zu%data(k,local_y,local_x))
      current_state%local_zumin = min(current_state%local_zumin, current_state%zu%data(k,local_y,local_x))
#endif
#ifdef V_ACTIVE
      current_state%local_zvmax = max(current_state%local_zvmax, current_state%zv%data(k,local_y,local_x))
      current_state%local_zvmin = min(current_state%local_zvmin, current_state%zv%data(k,local_y,local_x))
#endif
#ifdef W_ACTIVE
      if (k .lt. current_state%local_grid%local_domain_end_index(z_index)) then
        current_state%abswmax(k) = max(current_state%abswmax(k), abs(current_state%zw%data(k,local_y,local_x)))
      end if
#endif
    end do    

Here is the caller graph for this function:

field_information_retrieval_callback()

subroutine stepfields_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 component's 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
strcomp	Starting index within 1st argument string that matches substring (2nd argument); 0 if not a match.

Definition at line 693 of file stepfields.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 information for 3d
    strcomp=index(name, "_total_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_total_3d_local") then
        field_information%enabled=l_tend_3d_th
      else if (name .eq. "tend_qv_total_3d_local") then
        field_information%enabled=l_tend_3d_qv
      else if (name .eq. "tend_ql_total_3d_local") then
        field_information%enabled=l_tend_3d_ql
      else if (name .eq. "tend_qi_total_3d_local") then
        field_information%enabled=l_tend_3d_qi
      else if (name .eq. "tend_qr_total_3d_local") then
        field_information%enabled=l_tend_3d_qr
      else if (name .eq. "tend_qs_total_3d_local") then
        field_information%enabled=l_tend_3d_qs
      else if (name .eq. "tend_qg_total_3d_local") then
        field_information%enabled=l_tend_3d_qg
      else if (name .eq. "tend_tabs_total_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, "_total_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_total_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_th
      else if (name .eq. "tend_qv_total_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qv
      else if (name .eq. "tend_ql_total_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_ql
      else if (name .eq. "tend_qi_total_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qi
      else if (name .eq. "tend_qr_total_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qr
      else if (name .eq. "tend_qs_total_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qs
      else if (name .eq. "tend_qg_total_profile_total_local") then
        field_information%enabled=l_tend_pr_tot_qg
      else if (name .eq. "tend_tabs_total_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 stepfields_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 768 of file stepfields.F90.

    type(model_state_type), target, intent(inout) :: current_state
    character(len=*), intent(in) :: name
    type(component_field_value_type), intent(out) :: field_value
 
    ! 3d Tendency Fields
    if      (name .eq. "tend_th_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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_total_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 stepfields_mod::finalisation_callback ( type(model_state_type), intent(inout), target  current_state )

private

Finalisation callback.

Parameters

current_state

The current model state

Definition at line 209 of file stepfields.F90.

    type(model_state_type), target, intent(inout) :: current_state
 
    deallocate(resetq_min)
 
    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:

initialisation_callback()

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

private

Initialisation callback.

Parameters

current_state

The current model state

Definition at line 89 of file stepfields.F90.

    type(model_state_type), target, intent(inout) :: current_state
    logical :: l_qdiag
 
#ifdef U_ACTIVE
    allocate(current_state%global_grid%configuration%vertical%savolubar(current_state%local_grid%size(z_index)))
#endif
#ifdef V_ACTIVE
    allocate(current_state%global_grid%configuration%vertical%savolvbar(current_state%local_grid%size(z_index)))
#endif
    
    allocate(resetq_min(current_state%number_q_fields))
    cfl_is_enabled=is_component_enabled(current_state%options_database, "cfltest") 
    if (cfl_is_enabled) call reset_local_minmax_values(current_state)
 
    ! 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, 'stepfields')
      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, 'stepfields')
      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, 'stepfields')
      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, 'stepfields')
      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, 'stepfields')
      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, 'stepfields')
      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 call graph for this function:

Here is the caller graph for this function:

perform_timesmooth_for_field()

subroutine stepfields_mod::perform_timesmooth_for_field ( type(prognostic_field_type), intent(inout)  field, type(prognostic_field_type), intent(inout)  zfield, type(local_grid_type), intent(inout)  local_grid, integer, intent(in)  x_index, integer, intent(in)  y_index, real(kind=default_precision), intent(in)  c1, real(kind=default_precision), intent(in)  c2  )

private

Performs initial timesmoothing for a theta or Q field using Robert filter. This is finished off in swapsmooth.

Parameters

field	The field to smooth
zfield	The zfield to use in smoothing
local_grid	Description of the local grid
x_index	The X index to work on
y_index	The Y index to work on
c1	Constant to use in smoothing
c2	Constant to use in smoothing

Definition at line 515 of file stepfields.F90.

    type(prognostic_field_type), intent(inout) :: field, zfield
    type(local_grid_type), intent(inout) :: local_grid
    integer, intent(in) :: x_index, y_index
    real(kind=default_precision), intent(in) :: c1, c2
 
    integer :: k
 
    do k=1,local_grid%size(z_index)
      field%data(k, y_index, x_index)=c1*field%data(k, y_index, x_index)+c2*zfield%data(k, y_index, x_index)
    end do    

Here is the caller graph for this function:

remove_negative_rounding_errors_for_single_field()

subroutine stepfields_mod::remove_negative_rounding_errors_for_single_field ( integer, intent(in)  x_local_index, integer, intent(in)  y_local_index, integer, intent(in)  x_prev, integer, intent(in)  y_prev, type(prognostic_field_type), intent(inout)  field, type(local_grid_type), intent(inout)  local_grid  )

private

Removes the negative rounding errors from a specific single field. This works two columns behind and then catches up on the last column.

Parameters

x_local_index	The current local x index
y_local_index	The current local y index
x_prev	The previous x index to step
y_prev	The previous y index to step
field	The prognostic field
local_grid	Description of the local grid

Definition at line 409 of file stepfields.F90.

    integer, intent(in) :: x_local_index, y_local_index, x_prev, y_prev
    type(local_grid_type), intent(inout) :: local_grid
    type(prognostic_field_type), intent(inout) :: field
 
    if (x_prev .ge. local_grid%local_domain_start_index(x_index)) then
      call remove_negative_rounding_errors_in_slice(y_local_index, x_prev, y_prev, field, local_grid)
    end if
    if (x_local_index == local_grid%local_domain_end_index(x_index)) then
      if (x_local_index .gt. 1) then
        call remove_negative_rounding_errors_in_slice(y_local_index, x_local_index-1, y_prev, field, local_grid)
      end if
      call remove_negative_rounding_errors_in_slice(y_local_index, x_local_index, y_prev, field, local_grid)
    end if    

Here is the call graph for this function:

Here is the caller graph for this function:

remove_negative_rounding_errors_in_slice()

subroutine stepfields_mod::remove_negative_rounding_errors_in_slice ( integer, intent(in)  y_local_index, integer, intent(in)  x_prev, integer, intent(in)  y_prev, type(prognostic_field_type), intent(inout)  field, type(local_grid_type), intent(inout)  local_grid  )

private

Removes the negative rounding errors from a slice of a single field. This works two columns behind and then catches up on the last column.

Parameters

y_local_index	The current local y index
x_prev	The previous x index to step
y_prev	The previous y index to step
field	The prognostic field
local_grid	Description of the local grid

Definition at line 432 of file stepfields.F90.

    integer, intent(in) :: y_local_index, x_prev, y_prev
    type(local_grid_type), intent(inout) :: local_grid
    type(prognostic_field_type), intent(inout) :: field
 
    if (y_prev .ge. local_grid%local_domain_start_index(y_index)) then
      where (field%data(:, y_prev,  x_prev) < 0.0_default_precision)
        field%data(:, y_prev,  x_prev)=0.0_default_precision
      end where
    end if
    if (y_local_index == local_grid%local_domain_end_index(y_index)) then
      where (field%data(:, y_local_index,  x_prev) < 0.0_default_precision)
        field%data(:, y_local_index,  x_prev)=0.0_default_precision
      end where
    end if    

Here is the caller graph for this function:

reset_local_minmax_values()

subroutine stepfields_mod::reset_local_minmax_values ( type(model_state_type), intent(inout), target  current_state )

private

Resets the local min and max values for the flow fields.

Parameters

current_state

The current model state

Definition at line 390 of file stepfields.F90.

    type(model_state_type), intent(inout), target :: current_state
 
    ! Reset the local values for the next timestep
    current_state%local_zumin=rlargep
    current_state%local_zumax=-rlargep
    current_state%local_zvmin=rlargep
    current_state%local_zvmax=-rlargep
    current_state%abswmax=-rlargep

Here is the caller graph for this function:

set_published_field_value()

subroutine stepfields_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 817 of file stepfields.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:

step_all_fields()

subroutine stepfields_mod::step_all_fields ( type(model_state_type), intent(inout), target  current_state )

private

Steps all fields.

Parameters

current_state

The current model state_mod

Definition at line 315 of file stepfields.F90.

    type(model_state_type), target, intent(inout) :: current_state
 
    integer :: x_prev, y_prev, i, k
    real(kind=default_precision) :: c1, c2
 
    x_prev = current_state%column_local_x-2
    y_prev = current_state%column_local_y-1
 
    c1 = 1.0_default_precision - 2.0_default_precision*current_state%tsmth                                               
    c2 = current_state%tsmth
 
#ifdef U_ACTIVE
    current_state%global_grid%configuration%vertical%savolubar=current_state%global_grid%configuration%vertical%olzubar
    call step_single_field(current_state%column_local_x,  current_state%column_local_y, &
         x_prev, y_prev, current_state%u, current_state%zu, current_state%su, current_state%local_grid, .true., &
         current_state%field_stepping, current_state%dtm, current_state%ugal, c1, c2, .false., current_state%savu)
#endif
#ifdef V_ACTIVE
    current_state%global_grid%configuration%vertical%savolvbar=current_state%global_grid%configuration%vertical%olzvbar
    call step_single_field(current_state%column_local_x,  current_state%column_local_y, &
         x_prev, y_prev, current_state%v, current_state%zv, current_state%sv, current_state%local_grid, .true., &
         current_state%field_stepping, current_state%dtm, current_state%vgal, c1, c2, .false., current_state%savv)
#endif
#ifdef W_ACTIVE
    call step_single_field(current_state%column_local_x,  current_state%column_local_y, &
         x_prev, y_prev, current_state%w, current_state%zw, current_state%sw, current_state%local_grid, .false., &
         current_state%field_stepping, current_state%dtm, real(0., kind=default_precision), c1, c2, .false., current_state%savw)
#endif
    if (current_state%th%active) then       
       call step_single_field(current_state%column_local_x,  current_state%column_local_y, &
         x_prev, y_prev, current_state%th, current_state%zth, current_state%sth, current_state%local_grid, .false., &
         current_state%field_stepping, current_state%dtm, real(0., kind=default_precision), c1, c2, &
         current_state%field_stepping == centred_stepping)
    endif
    do i=1,current_state%number_q_fields
      if (current_state%q(i)%active) then         
        call step_single_field(current_state%column_local_x,  current_state%column_local_y, x_prev, y_prev, &
             current_state%q(i), current_state%zq(i), current_state%sq(i), current_state%local_grid, .false., &
             current_state%field_stepping, current_state%dtm, real(0., kind=default_precision), c1, c2, &
             current_state%field_stepping == centred_stepping)
      end if
    end do

Here is the call graph for this function:

Here is the caller graph for this function:

step_column_in_slice()

subroutine stepfields_mod::step_column_in_slice ( integer, intent(in)  y_local_index, integer, intent(in)  x_prev, integer, intent(in)  y_prev, type(prognostic_field_type), intent(inout)  field, type(prognostic_field_type), intent(inout)  zfield, type(prognostic_field_type), intent(inout)  sfield, type(local_grid_type), intent(inout)  local_grid, logical, intent(in)  flow_field, integer, intent(in)  direction, real(kind=default_precision), intent(in)  dtm, real(kind=default_precision), intent(in)  gal, real(kind=default_precision), intent(in)  c1, real(kind=default_precision), intent(in)  c2, logical, intent(in)  do_timesmoothing, type(prognostic_field_type), intent(inout), optional  sav  )

private

Will step a column in a specific slice. If y_prev is large enough then will step the y-1 column and if this is the last column of the slice then will also step the current column.

Parameters

x_local_index	The current local x index
y_local_index	The current local y index
x_prev	The previous x index to step
y_prev	The previous y index to step
field	The prognostic field
zfield	Z prognostic field
sfield	Source terms for the prognostic field
local_grid	Description of the local grid
flow_field	Whether or not this is a flow field
direction	The stepping direction (centred or forward)
dtm	The delta time per timestep
gal	Galilean transformation
c1	Constant to use in smoothing
c2	Constant to use in smoothing
do_timesmoothing	Whether timesmoothing using Robert filter should be done on the field
sav	Optional sav field

Definition at line 546 of file stepfields.F90.

    integer, intent(in) :: y_local_index, x_prev, y_prev, direction
    real(kind=default_precision), intent(in) :: dtm, gal
    logical, intent(in) :: flow_field, do_timesmoothing
    type(local_grid_type), intent(inout) :: local_grid
    type(prognostic_field_type), intent(inout) :: field, zfield, sfield
    real(kind=default_precision), intent(in) :: c1, c2
    type(prognostic_field_type), optional, intent(inout) :: sav
 
    if (y_prev .ge. local_grid%local_domain_start_index(y_index)) then
      if (do_timesmoothing) then
        call perform_timesmooth_for_field(field, zfield, local_grid, x_prev, y_prev, c1, c2)
      end if      
      if (present(sav)) then
        call step_field(x_prev, y_prev, field, zfield, sfield, local_grid, flow_field, direction, dtm, gal, sav)
      else
        call step_field(x_prev, y_prev, field, zfield, sfield, local_grid, flow_field, direction, dtm, gal)
      end if
    end if
 
    if (y_local_index == local_grid%local_domain_end_index(y_index)) then
      if (do_timesmoothing) then
        call perform_timesmooth_for_field(field, zfield, local_grid, x_prev, y_local_index, c1, c2)
      end if 
      if (present(sav)) then
        call step_field(x_prev, y_local_index, field, zfield, sfield, local_grid, flow_field, direction, dtm, gal, sav)
      else
        call step_field(x_prev, y_local_index, field, zfield, sfield, local_grid, flow_field, direction, dtm, gal)
      end if
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

step_field()

subroutine stepfields_mod::step_field ( integer, intent(in)  x_local_index, integer, intent(in)  y_local_index, type(prognostic_field_type), intent(inout)  field, type(prognostic_field_type), intent(inout)  zfield, type(prognostic_field_type), intent(inout)  sfield, type(local_grid_type), intent(inout)  local_grid, logical, intent(in)  flow_field, integer, intent(in)  direction, real(kind=default_precision), intent(in)  dtm, real(kind=default_precision), intent(in)  gal, type(prognostic_field_type), intent(inout), optional  sav  )

private

Will do the actual field stepping.

Parameters

flow_field	Whether or not we are stepping a flow field
direction	1=forward, 0=centred
x_index	The local X slice index
y_index	The local Y column index
kkp	Points in the vertical column
dtm	The model timestep
field	The prognostic field
zfield	The prognostic z field (which goes to timestep t+1)
xfield	The tendency of the field
gal	The galilean transformation

Definition at line 590 of file stepfields.F90.

    integer, intent(in) :: x_local_index, y_local_index, direction
    real(kind=default_precision), intent(in) :: dtm, gal
    logical, intent(in) :: flow_field
    type(local_grid_type), intent(inout) :: local_grid
    type(prognostic_field_type), intent(inout) :: field, zfield, sfield
    type(prognostic_field_type), optional, intent(inout) :: sav
 
    integer :: k
    real(kind=default_precision) :: actual_gal, dtm_x2
 
    dtm_x2 = 2.0_default_precision * dtm
 
    actual_gal = merge(gal, real(0.0_default_precision, kind=default_precision), flow_field)
 
    sfield%data(1,y_local_index, x_local_index)=0.0_default_precision
 
    do k=1,local_grid%size(z_index)
      ! Save the Z field which is used in the Robert filter
      if (present(sav) .and. direction .eq. centred_stepping) &
           sav%data(k,y_local_index, x_local_index) = zfield%data(k, y_local_index, x_local_index) + actual_gal
      if (flow_field) field%data(k, y_local_index, x_local_index) = actual_gal + field%data(k, y_local_index, x_local_index)
      if (direction == forward_stepping) then
        zfield%data(k, y_local_index, x_local_index) = field%data(k, y_local_index, x_local_index) + dtm * &
             sfield%data(k, y_local_index, x_local_index)
      else
        zfield%data(k, y_local_index, x_local_index) = actual_gal+zfield%data(k, y_local_index, x_local_index)+dtm_x2*&
             sfield%data(k, y_local_index, x_local_index)
      end if
    end do

Here is the caller graph for this function:

step_single_field()

subroutine stepfields_mod::step_single_field ( integer, intent(in)  x_local_index, integer, intent(in)  y_local_index, integer, intent(in)  x_prev, integer, intent(in)  y_prev, type(prognostic_field_type), intent(inout)  field, type(prognostic_field_type), intent(inout)  zfield, type(prognostic_field_type), intent(inout)  sfield, type(local_grid_type), intent(inout)  local_grid, logical, intent(in)  flow_field, integer, intent(in)  direction, real(kind=default_precision), intent(in)  dtm, real(kind=default_precision), intent(in)  gal, real(kind=default_precision), intent(in)  c1, real(kind=default_precision), intent(in)  c2, logical, intent(in)  do_timesmoothing, type(prognostic_field_type), intent(inout), optional  sav  )

private

Steps a single specific field. This will step on the yth column of the x-2 slice and x-1 and x if this is the last slice.

Parameters

x_local_index	The current local x index
y_local_index	The current local y index
x_prev	The previous x index to step
y_prev	The previous y index to step
field	The prognostic field
zfield	Z prognostic field
sfield	Source terms for the prognostic field
local_grid	Description of the local grid
flow_field	Whether or not this is a flow field
direction	The stepping direction (centred or forward)
dtm	The delta time per timestep
gal	Galilean transformation
c1	Constant to use in smoothing
c2	Constant to use in smoothing
do_timesmoothing	Whether timesmoothing using Robert filter should be done on the field
sav	Optional sav field

Definition at line 466 of file stepfields.F90.

    integer, intent(in) :: x_local_index, y_local_index, x_prev, y_prev, direction
    real(kind=default_precision), intent(in) :: dtm, gal
    logical, intent(in) :: flow_field, do_timesmoothing
    type(local_grid_type), intent(inout) :: local_grid
    type(prognostic_field_type), intent(inout) :: field, zfield, sfield
    real(kind=default_precision), intent(in) :: c1, c2
    type(prognostic_field_type), optional, intent(inout) :: sav
 
    if (x_prev .ge. local_grid%local_domain_start_index(x_index)) then
      if (present(sav)) then
        call step_column_in_slice(y_local_index, x_prev, y_prev, field, zfield, sfield, local_grid, &
             flow_field, direction, dtm, gal, c1, c2, do_timesmoothing, sav)
      else               
        call step_column_in_slice(y_local_index, x_prev, y_prev, field, zfield, sfield, local_grid, &
             flow_field, direction, dtm, gal, c1, c2, do_timesmoothing)
      end if
    end if
 
    if (x_local_index == local_grid%local_domain_end_index(x_index)) then
      ! If this is the last slice then process x-1 (if applicable) and x too
      if (x_local_index .gt. 1) then
        if (present(sav)) then
          call step_column_in_slice(y_local_index, x_local_index-1, y_prev, field, zfield, sfield, local_grid, &
               flow_field, direction, dtm, gal, c1, c2, do_timesmoothing, sav)
        else
          call step_column_in_slice(y_local_index, x_local_index-1, y_prev, field, zfield, sfield, local_grid, &
               flow_field, direction, dtm, gal, c1, c2, do_timesmoothing)
        end if
      end if
      if (present(sav)) then
        call step_column_in_slice(y_local_index, x_local_index, y_prev, field, zfield, sfield, local_grid, &
             flow_field, direction, dtm, gal, c1, c2, do_timesmoothing, sav)     
      else
        call step_column_in_slice(y_local_index, x_local_index, y_prev, field, zfield, sfield, local_grid, &
             flow_field, direction, dtm, gal, c1, c2, do_timesmoothing)     
      end if
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

stepfields_get_descriptor()

type(component_descriptor_type) function, public stepfields_mod::stepfields_get_descriptor

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

Returns

The KidReader component descriptor

Definition at line 56 of file stepfields.F90.

    stepfields_get_descriptor%name="stepfields"
    stepfields_get_descriptor%version=0.1
    stepfields_get_descriptor%initialisation=>initialisation_callback
    stepfields_get_descriptor%timestep=>timestep_callback
    stepfields_get_descriptor%finalisation=>finalisation_callback
 
    stepfields_get_descriptor%field_value_retrieval=>field_value_retrieval_callback
    stepfields_get_descriptor%field_information_retrieval=>field_information_retrieval_callback
    allocate(stepfields_get_descriptor%published_fields(8+8))
 
    stepfields_get_descriptor%published_fields(1)="tend_th_total_3d_local"
    stepfields_get_descriptor%published_fields(2)="tend_qv_total_3d_local"
    stepfields_get_descriptor%published_fields(3)="tend_ql_total_3d_local"
    stepfields_get_descriptor%published_fields(4)="tend_qi_total_3d_local"
    stepfields_get_descriptor%published_fields(5)="tend_qr_total_3d_local"
    stepfields_get_descriptor%published_fields(6)="tend_qs_total_3d_local"
    stepfields_get_descriptor%published_fields(7)="tend_qg_total_3d_local"
    stepfields_get_descriptor%published_fields(8)="tend_tabs_total_3d_local"
 
    stepfields_get_descriptor%published_fields(8+1)="tend_th_total_profile_total_local"
    stepfields_get_descriptor%published_fields(8+2)="tend_qv_total_profile_total_local"
    stepfields_get_descriptor%published_fields(8+3)="tend_ql_total_profile_total_local"
    stepfields_get_descriptor%published_fields(8+4)="tend_qi_total_profile_total_local"
    stepfields_get_descriptor%published_fields(8+5)="tend_qr_total_profile_total_local"
    stepfields_get_descriptor%published_fields(8+6)="tend_qs_total_profile_total_local"
    stepfields_get_descriptor%published_fields(8+7)="tend_qg_total_profile_total_local"
    stepfields_get_descriptor%published_fields(8+8)="tend_tabs_total_profile_total_local"
 

Here is the call graph for this function:

timestep_callback()

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

private

Called at each timestep and will perform swapping and smoothing as required.

Parameters

current_state

The current model state_mod

Definition at line 236 of file stepfields.F90.

    type(model_state_type), target, intent(inout) :: current_state
 
    integer :: iq
    integer :: current_x_index, current_y_index, target_x_index, target_y_index
 
    current_x_index=current_state%column_local_x
    current_y_index=current_state%column_local_y
    target_y_index=current_y_index-current_state%local_grid%halo_size(y_index)
    target_x_index=current_x_index-current_state%local_grid%halo_size(x_index)
 
 
    if (cfl_is_enabled .and. current_state%first_timestep_column) then
      if (mod(current_state%timestep, current_state%cfl_frequency) == 1 .or. &
         current_state%timestep-current_state%start_timestep .le. current_state%cfl_frequency) then
        determine_flow_minmax=.true.
        call reset_local_minmax_values(current_state)
      else
        determine_flow_minmax=.false.
      end if
    end if
 
    if (.not. current_state%halo_column) then
      if (determine_flow_minmax .and. cfl_is_enabled) &
         call determine_local_flow_minmax(current_state, current_state%column_local_y,  current_state%column_local_x)
      call step_all_fields(current_state)
    end if
 
    ! Remove negative rounding errors
    if (l_nonconservative_positive_q)then
      do iq=1,current_state%number_q_fields
        if (current_state%first_timestep_column)then
          resetq_min(iq)=minval(current_state%zq(iq)%data(:,current_state%column_local_y,  current_state%column_local_x))
        else
          resetq_min(iq)=min(resetq_min(iq),&
             minval(current_state%zq(iq)%data(:,current_state%column_local_y,  current_state%column_local_x)))
        end if
        call remove_negative_rounding_errors_for_single_field(current_state%column_local_x, current_state%column_local_y, &
             current_state%column_local_x-2, current_state%column_local_y-1, current_state%zq(iq), current_state%local_grid)
      end do
    end if
 
 
    ! 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 (mod(current_state%timestep, diagnostic_generation_frequency) == 0 .and. .not. current_state%halo_column) then
      call compute_component_tendencies(current_state, current_x_index, current_y_index, 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

cfl_is_enabled

logical stepfields_mod::cfl_is_enabled

private

Definition at line 22 of file stepfields.F90.

determine_flow_minmax

logical stepfields_mod::determine_flow_minmax =.false.

private

Definition at line 22 of file stepfields.F90.

  logical :: determine_flow_minmax=.false., cfl_is_enabled

diagnostic_generation_frequency

integer stepfields_mod::diagnostic_generation_frequency

private

Definition at line 47 of file stepfields.F90.

  integer :: diagnostic_generation_frequency

iqg

integer stepfields_mod::iqg =0

private

Definition at line 46 of file stepfields.F90.

iqi

integer stepfields_mod::iqi =0

private

Definition at line 46 of file stepfields.F90.

iql

integer stepfields_mod::iql =0

private

Definition at line 46 of file stepfields.F90.

iqr

integer stepfields_mod::iqr =0

private

Definition at line 46 of file stepfields.F90.

iqs

integer stepfields_mod::iqs =0

private

Definition at line 46 of file stepfields.F90.

iqv

integer stepfields_mod::iqv =0

private

Definition at line 46 of file stepfields.F90.

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

l_nonconservative_positive_q

logical stepfields_mod::l_nonconservative_positive_q =.true.

private

Definition at line 26 of file stepfields.F90.

  logical :: l_nonconservative_positive_q=.true.

l_tend_3d_qg

logical stepfields_mod::l_tend_3d_qg

private

Definition at line 34 of file stepfields.F90.

l_tend_3d_qi

logical stepfields_mod::l_tend_3d_qi

private

Definition at line 34 of file stepfields.F90.

l_tend_3d_ql

logical stepfields_mod::l_tend_3d_ql

private

Definition at line 34 of file stepfields.F90.

l_tend_3d_qr

logical stepfields_mod::l_tend_3d_qr

private

Definition at line 34 of file stepfields.F90.

l_tend_3d_qs

logical stepfields_mod::l_tend_3d_qs

private

Definition at line 34 of file stepfields.F90.

l_tend_3d_qv

logical stepfields_mod::l_tend_3d_qv

private

Definition at line 34 of file stepfields.F90.

l_tend_3d_tabs

logical stepfields_mod::l_tend_3d_tabs

private

Definition at line 34 of file stepfields.F90.

l_tend_3d_th

logical stepfields_mod::l_tend_3d_th

private

Definition at line 34 of file stepfields.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 stepfields_mod::l_tend_pr_tot_qg

private

Definition at line 42 of file stepfields.F90.

l_tend_pr_tot_qi

logical stepfields_mod::l_tend_pr_tot_qi

private

Definition at line 42 of file stepfields.F90.

l_tend_pr_tot_ql

logical stepfields_mod::l_tend_pr_tot_ql

private

Definition at line 42 of file stepfields.F90.

l_tend_pr_tot_qr

logical stepfields_mod::l_tend_pr_tot_qr

private

Definition at line 42 of file stepfields.F90.

l_tend_pr_tot_qs

logical stepfields_mod::l_tend_pr_tot_qs

private

Definition at line 42 of file stepfields.F90.

l_tend_pr_tot_qv

logical stepfields_mod::l_tend_pr_tot_qv

private

Definition at line 42 of file stepfields.F90.

l_tend_pr_tot_tabs

logical stepfields_mod::l_tend_pr_tot_tabs

private

Definition at line 42 of file stepfields.F90.

l_tend_pr_tot_th

logical stepfields_mod::l_tend_pr_tot_th

private

Definition at line 42 of file stepfields.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

resetq_min

real(kind=default_precision), dimension(:), allocatable stepfields_mod::resetq_min

private

Definition at line 25 of file stepfields.F90.

  real(kind=default_precision), allocatable :: resetq_min(:)

tend_3d_qg

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

private

Definition at line 30 of file stepfields.F90.

tend_3d_qi

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

private

Definition at line 30 of file stepfields.F90.

tend_3d_ql

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

private

Definition at line 30 of file stepfields.F90.

tend_3d_qr

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

private

Definition at line 30 of file stepfields.F90.

tend_3d_qs

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

private

Definition at line 30 of file stepfields.F90.

tend_3d_qv

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

private

Definition at line 30 of file stepfields.F90.

tend_3d_tabs

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

private

Definition at line 30 of file stepfields.F90.

tend_3d_th

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

private

Definition at line 30 of file stepfields.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 stepfields_mod::tend_pr_tot_qg

private

Definition at line 38 of file stepfields.F90.

tend_pr_tot_qi

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

private

Definition at line 38 of file stepfields.F90.

tend_pr_tot_ql

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

private

Definition at line 38 of file stepfields.F90.

tend_pr_tot_qr

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

private

Definition at line 38 of file stepfields.F90.

tend_pr_tot_qs

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

private

Definition at line 38 of file stepfields.F90.

tend_pr_tot_qv

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

private

Definition at line 38 of file stepfields.F90.

tend_pr_tot_tabs

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

private

Definition at line 38 of file stepfields.F90.

tend_pr_tot_th

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

private

Definition at line 38 of file stepfields.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