Conditionally averaged diagnostics, Part 1 of 2. More...

Functions/Subroutines
type(component_descriptor_type) function, public	conditional_diagnostics_column_get_descriptor ()
	Provides the component descriptor for the core to register. More...

subroutine	initialisation_callback (current_state)

subroutine	finalisation_callback (current_state)

subroutine	timestep_callback (current_state)

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...

Variables
integer	iqv

integer	iql

integer	iqr

integer	iqi

integer	iqs

integer	iqg

real(kind=default_precision), dimension(:,:,:), allocatable, public	conddiags_tot

real(kind=default_precision)	qlcrit

real(kind=default_precision)	qicrit

real(kind=default_precision)	qpptcrit

real(kind=default_precision)	vpptcrit

real(kind=default_precision)	thvprcrit

real(kind=default_precision)	wsdwncrit

real(kind=default_precision)	wsupcrit

real(kind=default_precision)	wupcrit

real(kind=default_precision)	wdwncrit

integer, public	ncond

integer, public	ndiag

integer	x_size

integer	y_size

real(kind=default_precision), public	gpts_total

logical	thv_from_th_with_liqice

logical	l_qi_qr_qs_qg

logical	l_do_near

character(len=string_length), dimension(23)	master_conditions_list

character(len=string_length), dimension(31)	master_diagnostics_list

character(len=string_length), dimension(:), allocatable, public	cond_request

character(len=string_length), dimension(:), allocatable, public	cond_long

character(len=string_length), dimension(:), allocatable, public	diag_request

character(len=string_length), dimension(:), allocatable, public	diag_long

integer, dimension(:), allocatable	diag_locations

integer, dimension(:), allocatable	cond_locations

integer, public	requested_area

integer	diagnostic_generation_frequency

Detailed Description

Conditionally averaged diagnostics, Part 1 of 2.

Function/Subroutine Documentation

◆ conditional_diagnostics_column_get_descriptor()

type(component_descriptor_type) function, public conditional_diagnostics_column_mod::conditional_diagnostics_column_get_descriptor

Provides the component descriptor for the core to register.

Returns: The descriptor describing this component

Definition at line 68 of file conditional_diagnostics_column.F90.

     conditional_diagnostics_column_get_descriptor%name="conditional_diagnostics_column"
     conditional_diagnostics_column_get_descriptor%version=0.1
     conditional_diagnostics_column_get_descriptor%initialisation=>initialisation_callback
     conditional_diagnostics_column_get_descriptor%timestep=>timestep_callback
     conditional_diagnostics_column_get_descriptor%finalisation=>finalisation_callback
  
     conditional_diagnostics_column_get_descriptor%field_value_retrieval=>field_value_retrieval_callback
     conditional_diagnostics_column_get_descriptor%field_information_retrieval=>field_information_retrieval_callback
     allocate(conditional_diagnostics_column_get_descriptor%published_fields(1))  
  
     conditional_diagnostics_column_get_descriptor%published_fields(1)="CondDiags_tot"

Here is the call graph for this function:

◆ field_information_retrieval_callback()

subroutine conditional_diagnostics_column_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 768 of file conditional_diagnostics_column.F90.

     type(model_state_type), target, intent(inout) :: current_state
     character(len=*), intent(in) :: name
     type(component_field_information_type), intent(out) :: field_information
  
  
     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) = ncond * 2 
     field_information%dimension_sizes(3) = ndiag
  
     field_information%data_type = component_double_data_type
  
     if (name .eq. "CondDiags_tot") then 
       field_information%enabled = allocated(conddiags_tot) 
     end if

Here is the caller graph for this function:

◆ field_value_retrieval_callback()

subroutine conditional_diagnostics_column_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 793 of file conditional_diagnostics_column.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. "CondDiags_tot") then  
        allocate(field_value%real_3d_array(current_state%local_grid%size(z_index),   &
                                           ncond * 2,   &
                                           ndiag)) 
        field_value%real_3d_array(:,:,:) = conddiags_tot(:,:,:) 
     end if

Here is the caller graph for this function:

◆ finalisation_callback()

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

private

Definition at line 360 of file conditional_diagnostics_column.F90.

     type(model_state_type), target, intent(inout) :: current_state
  
     if (allocated(conddiags_tot)) deallocate(conddiags_tot)
     if (allocated(cond_request)) deallocate(cond_request)
     if (allocated(cond_long)) deallocate(cond_long)
     if (allocated(cond_locations)) deallocate(cond_locations)
     if (allocated(diag_request)) deallocate(diag_request)
     if (allocated(diag_long)) deallocate(diag_long)
     if (allocated(diag_locations)) deallocate(diag_locations)
  

Here is the caller graph for this function:

◆ initialisation_callback()

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

private

This is section contains the full list of available conditions and

Check for requested diagnostics.

Record which conditions and diagnostics have been requested

Definition at line 83 of file conditional_diagnostics_column.F90.

     type(model_state_type), target, intent(inout) :: current_state
  
     integer :: total_points, inc
  
     !----------------------------------------------------------------------
     !   diagnostics.  
     !
     !   If you add more, you will need to update some of these:
     !         1. master_conditions_list (and its declared size)
     !         2. master_diagnostics_list (and its declared size)
     !         3. cond_locations CASE checks
     !         4. diag_locations CASE checks
     !         5. local_diag value
     !         6. l_cond logical test
     !         7. add short name to global_config default list
     !
     !   TAKE CARE TO ENSURE that indices correspond directly between:
     !         1. master_conditions_list / cond_locations
     !         2. master_diagnostics_list / diag_locations / local_diag
     !
     !----------------------------------------------------------------------
  
     master_conditions_list(1) = "all points"
     master_conditions_list(2) = "buoyant updraft"
     master_conditions_list(3) = "buoyant cloudy updraft"
     master_conditions_list(4) = "near buoyant cloudy updraft"
     master_conditions_list(5) = "cloudy"
     master_conditions_list(6) = "cloudy updraft"
     master_conditions_list(7) = "cloudy downdraft"
     master_conditions_list(8) = "strong updraft"
     master_conditions_list(9) = "strong downdraft"
     master_conditions_list(10) = "updraft"
     master_conditions_list(11) = "downdraft"
     master_conditions_list(12) = "liquid cloudy updraft"
     master_conditions_list(13) = "liquid cloudy downdraft"
     master_conditions_list(14) = "hydrometeors"
     master_conditions_list(15) = "cloud liquid water"
     master_conditions_list(16) = "cloud ice"
     master_conditions_list(17) = "precipitating downdraft"
     master_conditions_list(18) = "large precipitating downdraft"
     master_conditions_list(19) = "precipitating strong downdraft"
     master_conditions_list(20) = "saturated updraft"
     master_conditions_list(21) = "buoyant saturated updraft"
     master_conditions_list(22) = "cloud-free precipitating downdraft"
     master_conditions_list(23) = "cloud-free large precipitating downdraft"
  
     master_diagnostics_list(1) = "Fractional area"
     master_diagnostics_list(2) = "Vertical velocity"
     master_diagnostics_list(3) = "Vertical velocity variance"
     master_diagnostics_list(4) = "Potential temperature"
     master_diagnostics_list(5) = "w * th"
     master_diagnostics_list(6) = "Potential temperature anomalies"
     master_diagnostics_list(7) = "Flux of potential temperature"
     master_diagnostics_list(8) = "Virtual potential temperature anomalies"
     master_diagnostics_list(9) = "Flux of virtual potential temperature"
     master_diagnostics_list(10) = "Variance of potential temperature"
     master_diagnostics_list(11) = "Diff_coeff * dth/dz"
     master_diagnostics_list(12) = "w**3"
     master_diagnostics_list(13) = "Relative humidity"
     master_diagnostics_list(14) = "Zonal velocity"
     master_diagnostics_list(15) = "Meridional velovity"
     master_diagnostics_list(16) = "w * u"
     master_diagnostics_list(17) = "w * v"
     master_diagnostics_list(18) = "Vis_coeff * du/dz"
     master_diagnostics_list(19) = "Vis_coeff * dv/dz"
     master_diagnostics_list(20) = "Temperature"
     master_diagnostics_list(21) = "th+lvap*qv/cp"
     master_diagnostics_list(22) = "Anomalies of THL"
     master_diagnostics_list(23) = "Variance of THL"
     master_diagnostics_list(24) = "qv+ql+qi"
     master_diagnostics_list(25) = "Anomalies of QVLI"
     master_diagnostics_list(26) = "Variance of QVLI"
     master_diagnostics_list(27) = "qr+qs+qg"
     master_diagnostics_list(28) = "Anomalies of QRSG"
     master_diagnostics_list(29) = "Variance of QRSG"
     master_diagnostics_list(30) = "Flux of QVLI"
     master_diagnostics_list(31) = "Flux of QRSG"
  
  
     total_points=(current_state%local_grid%size(y_index) * current_state%local_grid%size(x_index))
  
     if (.not. current_state%passive_q .and. current_state%number_q_fields .gt. 0) then                        
       iqv       =get_q_index(standard_q_names%VAPOUR, 'profile_diags')                                                 
       iql       =get_q_index(standard_q_names%CLOUD_LIQUID_MASS, 'profile_diags')
       iqi       =get_q_index(standard_q_names%ICE_MASS, 'profile_diags')
       iqr       =get_q_index(standard_q_names%RAIN_MASS, 'profile_diags')
       iqs       =get_q_index(standard_q_names%SNOW_MASS, 'profile_diags')
       iqg       =get_q_index(standard_q_names%GRAUPEL_MASS, 'profile_diags')
                                                                                                                      
       qlcrit    =options_get_real(current_state%options_database, "qlcrit")  
       qicrit    =options_get_real(current_state%options_database, "qicrit")  
       qpptcrit  =options_get_real(current_state%options_database, "qpptcrit") 
       vpptcrit  =options_get_real(current_state%options_database, "vpptcrit")                                        
     endif
  
     !calculate qi, qr, qs and qg if and only if number_q_fields=11
     l_qi_qr_qs_qg =  (.not. current_state%passive_q .and. current_state%number_q_fields .ge. 11)
  
     x_size       =options_get_integer(current_state%options_database, "x_size")
     y_size       =options_get_integer(current_state%options_database, "y_size")
     thvprcrit    =options_get_real(current_state%options_database, "thvprcrit")  
     wsdwncrit    =options_get_real(current_state%options_database, "wSdwncrit")  
     wsupcrit     =options_get_real(current_state%options_database, "wSupcrit")  
     wupcrit      =options_get_real(current_state%options_database, "wupcrit")  
     wdwncrit     =options_get_real(current_state%options_database, "wdwncrit")  
     ncond        =options_get_array_size(current_state%options_database, "cond_request")
     ndiag        =options_get_array_size(current_state%options_database, "diag_request")
     thv_from_th_with_liqice =options_get_logical(current_state%options_database, "thv_from_th_with_liqice") 
  
     gpts_total   = 1.0_default_precision * x_size * y_size
  
     if (ncond .lt. 1 .or. ndiag .lt. 1) call log_master_log(log_error, &
          "When conditional_diagnostics_column are enabled, condition and diagnostic request lists must be provided.")
  
     if (current_state%th%active) then
       allocate(conddiags_tot(current_state%local_grid%size(z_index),  &
                              ncond * 2 ,  &
                              ndiag))
     endif
  
  
     l_do_near = .false.
     requested_area = 0
     allocate(cond_request(ncond))
     allocate(cond_long(ncond))
     allocate(cond_locations(ncond))
     allocate(diag_request(ndiag))
     allocate(diag_long(ndiag))
     allocate(diag_locations(ndiag))
  
     do inc = 1, ncond
       cond_request(inc) = trim(options_get_string(current_state%options_database, "cond_request", inc))
  
       SELECT CASE ( trim(cond_request(inc)) )
         CASE DEFAULT
           call log_master_log(log_error, &
               "Condition '"//trim(cond_request(inc))//"' has not been set up.")
         CASE ("ALL")
           cond_locations(inc) = 1
         CASE ("BYu")
           cond_locations(inc) = 2
         CASE ("BCu")
           cond_locations(inc) = 3
         CASE ("NrBCu")
           cond_locations(inc) = 4
           l_do_near = .true.
         CASE ("AC")
           cond_locations(inc) = 5
         CASE ("ACu")
           cond_locations(inc) = 6
         CASE ("ACd")
           cond_locations(inc) = 7
         CASE ("WG1")
           cond_locations(inc) = 8
         CASE ("WL1")
           cond_locations(inc) = 9
         CASE ("ALu")
           cond_locations(inc) = 10
         CASE ("ALd")
           cond_locations(inc) = 11
         CASE ("CLu")
           cond_locations(inc) = 12
         CASE ("CLd")
           cond_locations(inc) = 13
         CASE ("AH")
           cond_locations(inc) = 14
         CASE ("AL")
           cond_locations(inc) = 15
         CASE ("AI")
           cond_locations(inc) = 16
         CASE ("PPd")
           cond_locations(inc) = 17
         CASE ("VPd")
           cond_locations(inc) = 18
         CASE ("PVd")
           cond_locations(inc) = 19
         CASE ("MO")
           cond_locations(inc) = 20
         CASE ("BM")
           cond_locations(inc) = 21
         CASE ("AA")
           cond_locations(inc) = 22
         CASE ("AV")
           cond_locations(inc) = 23
       END SELECT
     end do ! loop over requested conditions
  
  
     do inc = 1, ndiag
       diag_request(inc) = trim(options_get_string(current_state%options_database, "diag_request", inc))
       if ( trim(diag_request(inc)) .eq. "area" ) requested_area = inc
  
       SELECT CASE ( trim(diag_request(inc)) )
         CASE DEFAULT
           call log_master_log(log_error, &
                      "Diagnostic '"//trim(diag_request(inc))//"' has not been set up.")
         CASE ("area")
           diag_locations(inc) = 1
         CASE ("W")
           diag_locations(inc) = 2
         CASE ("W2")
           diag_locations(inc) = 3
         CASE ("TH")
           diag_locations(inc) = 4
         CASE ("WTH")
           diag_locations(inc) = 5
         CASE ("THP")
           diag_locations(inc) = 6
         CASE ("WTHP")
           diag_locations(inc) = 7
         CASE ("THVP")
           diag_locations(inc) = 8
         CASE ("WTHVP")
           diag_locations(inc) = 9
         CASE ("THP2")
           diag_locations(inc) = 10
         CASE ("WTHSG")
           diag_locations(inc) = 11
         CASE ("W3")
           diag_locations(inc) = 12
         CASE ("RH")
           diag_locations(inc) = 13
         CASE ("U")
           diag_locations(inc) = 14
         CASE ("V")
           diag_locations(inc) = 15
         CASE ("WU")
           diag_locations(inc) = 16
         CASE ("WV")
           diag_locations(inc) = 17
         CASE ("WUSG")
           diag_locations(inc) = 18
         CASE ("WVSG")
           diag_locations(inc) = 19
         CASE ("TEMP")
           diag_locations(inc) = 20
         CASE ("THL")
           diag_locations(inc) = 21
         CASE ("THLP")
           diag_locations(inc) = 22
         CASE ("THLP2")
           diag_locations(inc) = 23
         CASE ("QVLI")
           diag_locations(inc) = 24
         CASE ("QVLIP")
           diag_locations(inc) = 25
         CASE ("QVLIP2")
           diag_locations(inc) = 26
         CASE ("QRSG")
           diag_locations(inc) = 27
         CASE ("QRSGP")
           diag_locations(inc) = 28
         CASE ("QRSGP2")
           diag_locations(inc) = 29
         CASE ("WQVLIP")
           diag_locations(inc) = 30
         CASE ("WQRSGP")
           diag_locations(inc) = 31
       END SELECT
     end do ! loop over requested diagnostics
  
     if (requested_area == 0) call log_master_log(log_error, &
           "The diagnostic 'area' must be provided to complete the conditional diagnostics process.")
  
     cond_long=master_conditions_list(cond_locations)
     diag_long=master_diagnostics_list(diag_locations)
  
     ! 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:

◆ timestep_callback()

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

private

Reset sums on this process

Decide if conditions are appropriate to proceed with calculations

Begin the calculations Loop over levels

work out all potential terms at this location

end of potential term calculations

Store the diagnostics

Determine status of requested conditions

Near it, but not it.

Record the value for this diagnostic meeting this condition and its complement by adding to cumulative sum

Definition at line 374 of file conditional_diagnostics_column.F90.

     type(model_state_type), target, intent(inout) :: current_state
  
     logical ::  l_cond                                     ! condition at local position (k,j,i)
     real(kind=default_precision), dimension(size(master_diagnostics_list)) ::  local_diag ! diagnostics at local position (k,j,i)
  
     real(kind=default_precision) ::  qi, ql,qv,qli,qsat,qppt,qvl
     real(kind=default_precision) ::  qi_pr, ql_pr,qv_pr,qli_pr,qvli, qvli_pr, qppt_pr
     real(kind=default_precision) ::  thv_pr,thv,exner,tdegk,th_pr,th_pr2
     real(kind=default_precision) ::  tmp_th !=thref+th
     real(kind=default_precision) ::  pottemp !=thref+olzthbar
     real(kind=default_precision) ::  relhum,pmb,wth,wthsg, wthv_pr, wth_pr
     real(kind=default_precision) ::  w_zn ! vertical velocity at theta levels
     real(kind=default_precision) ::  tmp_u,tmp_v ! horizontal veleocities
     real(kind=default_precision) ::  w_zn2, w_zn3, wu, wv, wusg, wvsg, wqvli_pr,wqppt_pr
     real(kind=default_precision) ::  th_h,th_h_pr1,th_h_pr2 
     real(kind=default_precision) ::  qvli_pr2, qppt_pr2
  
     real(kind=default_precision) ::  qv_jip1,qv_jim1,qv_jp1i,qv_jm1i,qv_pr_jip1,qv_pr_jim1,qv_pr_jp1i,qv_pr_jm1i
     real(kind=default_precision) ::  qi_jip1,qi_jim1,qi_jp1i,qi_jm1i
     real(kind=default_precision) ::  ql_jip1,ql_jim1,ql_jp1i,ql_jm1i
     real(kind=default_precision) ::  qli_jip1,qli_jim1,qli_jp1i,qli_jm1i,qli_pr_jip1,qli_pr_jim1,qli_pr_jp1i,qli_pr_jm1i
     real(kind=default_precision) ::  w_zn_jip1,w_zn_jim1,w_zn_jp1i,w_zn_jm1i
     real(kind=default_precision) ::  thv_pr_jip1,thv_pr_jim1,thv_pr_jp1i,thv_pr_jm1i
     real(kind=default_precision) ::  tmp_th_jip1,tmp_th_jim1,tmp_th_jp1i,tmp_th_jm1i
     real(kind=default_precision) ::  th_pr_jip1,th_pr_jim1,th_pr_jp1i,th_pr_jm1i
  
     integer :: k, j, i
     integer :: inc     ! loop increment variable
     integer :: local_y, locaL_x, target_x_index, target_y_index
  
     j=current_state%column_local_y
     i=current_state%column_local_x
  
     if (current_state%first_timestep_column) then
       conddiags_tot(:,:,:)    = 0.0_default_precision  
     end if
  
     if (current_state%halo_column) return
     if ( .not. (current_state%th%active .and.       &
                 .not. current_state%passive_q .and. &
                 current_state%number_q_fields .gt. 0) ) return
     if (.not. mod(current_state%timestep, diagnostic_generation_frequency) == 0) return
  
     !  Due to above/below level solution dependency, no data will exist at first/last levels,
     !  but k-dimension size kept equal to size(Z_INDEX) for convenience
     do k = 2, current_state%local_grid%size(z_index)-1
  
       qv      = current_state%zq(iqv)%data(k,j,i)
       ql      = current_state%zq(iql)%data(k,j,i)
       qvl     = qv + ql
       qv_pr   = qv - current_state%global_grid%configuration%vertical%olzqbar(k,iqv)
       ql_pr   = ql - current_state%global_grid%configuration%vertical%olzqbar(k,iql)
  
       if (.not. l_qi_qr_qs_qg) then
         qi      = 0.0
         qli     = ql
         qvli    = qvl
         qppt    = 0.0
         qi_pr   = 0.0
         qli_pr  = ql_pr
         qvli_pr = qv_pr + ql_pr
         qppt_pr = 0.0
       end if
  
       if (l_qi_qr_qs_qg) then
         qi      = current_state%zq(iqi)%data(k,j,i)  
         qli     = ql + qi
         qvli    = qv + qli
         qppt    = current_state%zq(iqr)%data(k,j,i) + current_state%zq(iqs)%data(k,j,i)+current_state%zq(iqg)%data(k,j,i)
         qi_pr   = qli - current_state%global_grid%configuration%vertical%olzqbar(k,iqi)
         qli_pr  = ql_pr + qi_pr
         qvli_pr = qv_pr + qli_pr
         qppt_pr = current_state%zq(iqr)%data(k,j,i) - current_state%global_grid%configuration%vertical%olzqbar(k,iqr)+   &
                   current_state%zq(iqs)%data(k,j,i) - current_state%global_grid%configuration%vertical%olzqbar(k,iqs)+   &
                   current_state%zq(iqg)%data(k,j,i) - current_state%global_grid%configuration%vertical%olzqbar(k,iqg)
       end if
  
       w_zn     = 0.5 * (current_state%zw%data(k,j,i) + current_state%zw%data(k-1,j,i)) !w at theta levels
       w_zn2    = w_zn * w_zn
       w_zn3    = w_zn * w_zn * w_zn
  
       qvli_pr2 = qvli_pr * qvli_pr
       qppt_pr2 = qppt_pr * qppt_pr
       wqvli_pr = w_zn * qvli_pr
       wqppt_pr = w_zn * qppt_pr
  
       tmp_u    = 0.5 * (current_state%zu%data(k,j,i-1) + current_state%zu%data(k,j,i)) + current_state%ugal !u at theta points
       tmp_v    = 0.5 * (current_state%zv%data(k,j-1,i) + current_state%zv%data(k,j,i)) + current_state%vgal !v at theta points
       wu       = w_zn * tmp_u
       wv       = w_zn * tmp_v
       tmp_th   = current_state%global_grid%configuration%vertical%thref(k) + current_state%zth%data(k,j,i)
       wth      = w_zn * tmp_th
       wthsg    = -0.5 * current_state%diff_coefficient%data(k,j,i) * (current_state%zth%data(k+1,j,i) +          &
                  current_state%global_grid%configuration%vertical%thref(k+1) - current_state%zth%data(k,j,i) -   &
                  current_state%global_grid%configuration%vertical%thref(k)) *                                    &
                  current_state%global_grid%configuration%vertical%rdzn(k+1)                                      &
                  -0.5 * current_state%diff_coefficient%data(k-1,j,i) * (current_state%zth%data(k,j,i) +          &
                  current_state%global_grid%configuration%vertical%thref(k) - current_state%zth%data(k-1,j,i) -   &
                  current_state%global_grid%configuration%vertical%thref(k-1)) *                                  &
                  current_state%global_grid%configuration%vertical%rdzn(k)
  
       wusg     = -0.25 * current_state%vis_coefficient%data(k  ,j,i-1) *                      &
                  (current_state%zu%data(k+1,j,i-1) - current_state%zu%data(k  ,j,i-1)) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k+1)                   &
                  -0.25 * current_state%vis_coefficient%data(k  ,j,i  ) *                      &
                  (current_state%zu%data(k+1,j,i  ) - current_state%zu%data(k  ,j,i  )) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k+1)                   &
                  -0.25 * current_state%vis_coefficient%data(k-1,j,i-1) *                      &
                  (current_state%zu%data(k  ,j,i-1) - current_state%zu%data(k-1,j,i-1)) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k)                     &
                  -0.25 * current_state%vis_coefficient%data(k-1,j,i  ) *                      &
                  (current_state%zu%data(k  ,j,i  ) - current_state%zu%data(k-1,j,i  )) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k)  
  
       wvsg     = -0.25 * current_state%vis_coefficient%data(k  ,j-1,i) *                      &
                  (current_state%zv%data(k+1,j-1,i) - current_state%zv%data(k  ,j-1,i)) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k+1)                   &
                  -0.25 * current_state%vis_coefficient%data(k  ,j  ,i) *                      &
                  (current_state%zv%data(k+1,j  ,i) - current_state%zv%data(k  ,j  ,i)) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k+1)                   &
                  -0.25 * current_state%vis_coefficient%data(k-1,j-1,i) *                      &
                  (current_state%zv%data(k  ,j-1,i) - current_state%zv%data(k-1,j-1,i)) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k)                     &
                  -0.25 * current_state%vis_coefficient%data(k-1,j  ,i) *                      &
                  (current_state%zv%data(k  ,j  ,i) - current_state%zv%data(k-1,j  ,i)) *      &
                  current_state%global_grid%configuration%vertical%rdzn(k)
  
       pottemp  = current_state%global_grid%configuration%vertical%thref(k) +    &
                  current_state%global_grid%configuration%vertical%olzthbar(k)
  
       th_pr    = current_state%zth%data(k,j,i) - current_state%global_grid%configuration%vertical%olzthbar(k)
       wth_pr   = w_zn * th_pr
       th_pr2   = th_pr * th_pr
       exner    = current_state%global_grid%configuration%vertical%rprefrcp(k)
       tdegk    = (current_state%global_grid%configuration%vertical%thref(k) + current_state%zth%data(k,j,i)) * exner
       pmb      = current_state%global_grid%configuration%vertical%prefn(k) / 100.
       qsat     = qsaturation(tdegk,pmb)
       relhum   = qv / qsat
       th_h     = tmp_th + rlvap * qv / cp
       th_h_pr1 = th_pr + rlvap * qv_pr / cp
       th_h_pr2 = th_h_pr1 * th_h_pr1
  
       !calculate values at nearby points
       if (l_do_near) then
         qv_jip1      = current_state%zq(iqv)%data(k,j,i+1) !qv at the position j and i+1
         qv_jim1      = current_state%zq(iqv)%data(k,j,i-1) !qv at the position j and i-1
         qv_jp1i      = current_state%zq(iqv)%data(k,j+1,i) !qv at the position j+1 and i
         qv_jm1i      = current_state%zq(iqv)%data(k,j-1,i) !qv at the position j-1 and i
         ql_jip1      = current_state%zq(iql)%data(k,j,i+1)
         ql_jim1      = current_state%zq(iql)%data(k,j,i-1)
         ql_jp1i      = current_state%zq(iql)%data(k,j+1,i)
         ql_jm1i      = current_state%zq(iql)%data(k,j-1,i)
         qv_pr_jip1   = current_state%zq(iqv)%data(k,j,i+1)-current_state%global_grid%configuration%vertical%olzqbar(k,iqv)
         qv_pr_jim1   = current_state%zq(iqv)%data(k,j,i-1)-current_state%global_grid%configuration%vertical%olzqbar(k,iqv)
         qv_pr_jp1i   = current_state%zq(iqv)%data(k,j+1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iqv)
         qv_pr_jm1i   = current_state%zq(iqv)%data(k,j-1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iqv)
  
         if (.not. l_qi_qr_qs_qg) then
           qi_jip1      = 0.0
           qi_jim1      = 0.0
           qi_jp1i      = 0.0
           qi_jm1i      = 0.0
           qli_jip1     = current_state%zq(iql)%data(k,j,i+1)+0.0
           qli_jim1     = current_state%zq(iql)%data(k,j,i-1)+0.0
           qli_jp1i     = current_state%zq(iql)%data(k,j+1,i)+0.0
           qli_jm1i     = current_state%zq(iql)%data(k,j-1,i)+0.0
           qli_pr_jip1  = current_state%zq(iql)%data(k,j,i+1)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+0.0
           qli_pr_jim1  = current_state%zq(iql)%data(k,j,i-1)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+0.0
           qli_pr_jp1i  = current_state%zq(iql)%data(k,j+1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+0.0
           qli_pr_jm1i  = current_state%zq(iql)%data(k,j-1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+0.0
         end if
  
         if (l_qi_qr_qs_qg) then
           qi_jip1      = current_state%zq(iqi)%data(k,j,i+1)
           qi_jim1      = current_state%zq(iqi)%data(k,j,i-1)
           qi_jp1i      = current_state%zq(iqi)%data(k,j+1,i)
           qi_jm1i      = current_state%zq(iqi)%data(k,j-1,i)
           qli_jip1     = current_state%zq(iql)%data(k,j,i+1)+current_state%zq(iqi)%data(k,j,i+1)
           qli_jim1     = current_state%zq(iql)%data(k,j,i-1)+current_state%zq(iqi)%data(k,j,i-1)
           qli_jp1i     = current_state%zq(iql)%data(k,j+1,i)+current_state%zq(iqi)%data(k,j+1,i)
           qli_jm1i     = current_state%zq(iql)%data(k,j-1,i)+current_state%zq(iqi)%data(k,j-1,i)
           qli_pr_jip1  = current_state%zq(iql)%data(k,j,i+1)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+ &
                          current_state%zq(iqi)%data(k,j,i+1)-current_state%global_grid%configuration%vertical%olzqbar(k,iqi)
           qli_pr_jim1  = current_state%zq(iql)%data(k,j,i-1)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+ &
                          current_state%zq(iqi)%data(k,j,i-1)-current_state%global_grid%configuration%vertical%olzqbar(k,iqi)
           qli_pr_jp1i  = current_state%zq(iql)%data(k,j+1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+ &
                          current_state%zq(iqi)%data(k,j+1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iqi)
           qli_pr_jm1i  = current_state%zq(iql)%data(k,j-1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iql)+ &
                          current_state%zq(iqi)%data(k,j-1,i)-current_state%global_grid%configuration%vertical%olzqbar(k,iqi)
         end if
  
         tmp_th_jip1  = current_state%global_grid%configuration%vertical%thref(k) + current_state%zth%data(k,j,i+1) 
         tmp_th_jim1  = current_state%global_grid%configuration%vertical%thref(k) + current_state%zth%data(k,j,i-1) 
         tmp_th_jp1i  = current_state%global_grid%configuration%vertical%thref(k) + current_state%zth%data(k,j+1,i) 
         tmp_th_jm1i  = current_state%global_grid%configuration%vertical%thref(k) + current_state%zth%data(k,j-1,i) 
         th_pr_jip1   = current_state%zth%data(k,j,i+1)-current_state%global_grid%configuration%vertical%olzthbar(k)
         th_pr_jim1   = current_state%zth%data(k,j,i-1)-current_state%global_grid%configuration%vertical%olzthbar(k)
         th_pr_jp1i   = current_state%zth%data(k,j-1,i)-current_state%global_grid%configuration%vertical%olzthbar(k)
         th_pr_jm1i   = current_state%zth%data(k,j+1,i)-current_state%global_grid%configuration%vertical%olzthbar(k)
         w_zn_jip1    = 0.5*(current_state%zw%data(k,j,i+1)+current_state%zw%data(k-1,j,i+1)) 
         w_zn_jim1    = 0.5*(current_state%zw%data(k,j,i-1)+current_state%zw%data(k-1,j,i-1)) 
         w_zn_jp1i    = 0.5*(current_state%zw%data(k,j+1,i)+current_state%zw%data(k-1,j+1,i)) 
         w_zn_jm1i    = 0.5*(current_state%zw%data(k,j-1,i)+current_state%zw%data(k-1,j-1,i)) 
  
         if (thv_from_th_with_liqice) then  
           thv       = tmp_th*(1.0+0.61*qv-qli)
           !the term th_pr*(0.61*qv-qli) is very small and can be neglected
           thv_pr    = th_pr+tmp_th*(0.61*qv_pr-qli_pr)+th_pr*(0.61*qv-qli)
           wthv_pr   = w_zn*thv_pr
  
           thv_pr_jip1   =  th_pr_jip1+tmp_th_jip1*(0.61*qv_pr_jip1-qli_pr_jip1)+th_pr_jip1*(0.61*qv_jip1-qli_jip1)
           thv_pr_jim1   =  th_pr_jim1+tmp_th_jim1*(0.61*qv_pr_jim1-qli_pr_jim1)+th_pr_jim1*(0.61*qv_jim1-qli_jim1)
           thv_pr_jp1i   =  th_pr_jp1i+tmp_th_jp1i*(0.61*qv_pr_jp1i-qli_pr_jp1i)+th_pr_jp1i*(0.61*qv_jp1i-qli_jp1i)
           thv_pr_jm1i   =  th_pr_jm1i+tmp_th_jm1i*(0.61*qv_pr_jm1i-qli_pr_jm1i)+th_pr_jm1i*(0.61*qv_jm1i-qli_jm1i)
         end if
  
         if (.not. thv_from_th_with_liqice) then  
           thv      = tmp_th*(1.0+0.61*qv)
           ! the 0.61*qv*th_pr is very small and can be neglected  
           thv_pr   = th_pr + 0.61 * tmp_th * qv_pr + 0.61 * qv * th_pr 
           wthv_pr  = w_zn * thv_pr
           thv_pr_jip1   = th_pr_jip1+0.61*tmp_th_jip1*qv_pr_jip1+0.61*qv_jip1*th_pr_jip1 
           thv_pr_jim1   = th_pr_jim1+0.61*tmp_th_jim1*qv_pr_jim1+0.61*qv_jim1*th_pr_jim1 
           thv_pr_jp1i   = th_pr_jp1i+0.61*tmp_th_jp1i*qv_pr_jp1i+0.61*qv_jp1i*th_pr_jp1i 
           thv_pr_jm1i   = th_pr_jm1i+0.61*tmp_th_jm1i*qv_pr_jm1i+0.61*qv_jm1i*th_pr_jm1i 
         end if 
  
       end if ! l_do_near
  
  
       local_diag(1) = 1.0_default_precision
       local_diag(2) = w_zn
       local_diag(3) = w_zn2
       local_diag(4) = tmp_th
       local_diag(5) = wth
       local_diag(6) = th_pr
       local_diag(7) = wth_pr
       local_diag(8) = thv_pr
       local_diag(9) = wthv_pr
       local_diag(10) = th_pr2
       local_diag(11) = wthsg
       local_diag(12) = w_zn3
       local_diag(13) = relhum
       local_diag(14) = tmp_u
       local_diag(15) = tmp_v
       local_diag(16) = wu
       local_diag(17) = wv
       local_diag(18) = wusg
       local_diag(19) = wvsg
       local_diag(20) = tdegk
       local_diag(21) = th_h
       local_diag(22) = th_h_pr1
       local_diag(23) = th_h_pr2
       local_diag(24) = qvli
       local_diag(25) = qvli_pr
       local_diag(26) = qvli_pr2
       local_diag(27) = qppt
       local_diag(28) = qppt_pr
       local_diag(29) = qppt_pr2
       local_diag(30) = wqvli_pr
       local_diag(31) = wqppt_pr
  
       do inc = 1, ncond
  
         l_cond = .false.  ! initialization
  
         !Condition processing on requested condition
         SELECT CASE ( trim(cond_request(inc)) )
  
           CASE DEFAULT
             call log_master_log(log_error, &
                   "Condition '"//trim(cond_request(inc))//"' has not been set up.")
  
           CASE ("ALL")
             l_cond = .true.
  
           CASE ("BYu")
             l_cond = ( (w_zn .gt. wupcrit) .and. (thv_pr .gt. thvprcrit) )
  
           CASE ("BCu")
             l_cond = ( (w_zn .gt. wupcrit .and. thv_pr .gt. thvprcrit ) .and.      &
                        ( ql .gt. qlcrit .or. qi .gt. qicrit ) )
  
           CASE ("NrBCu")
             if ( .not. ( (w_zn .gt. wupcrit .and. thv_pr .gt. thvprcrit ) .and.      &
                        ( ql .gt. qlcrit .or. qi .gt. qicrit ) )    ) then
  
               l_cond = (  ( (w_zn_jip1 .gt. wupcrit .and. thv_pr_jip1 .gt. thvprcrit) .and.   &
                             (ql_jip1 .gt. qlcrit .or. qi_jip1 .gt. qicrit) )                  &
                           .or.                                                                &
                           ( (w_zn_jim1 .gt. wupcrit .and. thv_pr_jim1 .gt. thvprcrit) .and.   &
                             (ql_jim1 .gt. qlcrit .or. qi_jim1 .gt. qicrit) )                  &
                           .or.                                                                &
                           ( (w_zn_jp1i .gt. wupcrit .and. thv_pr_jp1i .gt. thvprcrit) .and.   &
                             (ql_jp1i .gt. qlcrit .or. qi_jp1i .gt. qicrit) )                  &
                           .or.                                                                &
                           ( (w_zn_jm1i .gt. wupcrit .and. thv_pr_jm1i .gt. thvprcrit) .and.   &
                             (ql_jm1i .gt. qlcrit .or. qi_jm1i .gt. qicrit) )  )
             
             end if
  
           CASE ("AC")
             l_cond = ( (ql .gt. qlcrit) .or. (qi .gt. qicrit) )
  
           CASE ("ACu")
             l_cond = ( (w_zn .gt. wupcrit) .and.          &
                        ( (ql .gt. qlcrit) .or. (qi .gt. qicrit) ) )
  
           CASE ("ACd")
             l_cond = ( (w_zn .lt. wdwncrit) .and.         &
                        ( (ql .gt. qlcrit) .or. (qi .gt. qicrit) ) )
  
           CASE ("WG1")
             l_cond = ( w_zn .gt. wsupcrit )
  
           CASE ("WL1")
             l_cond = ( w_zn .lt. wsdwncrit )
  
           CASE ("ALu")
             l_cond = ( w_zn .gt. wupcrit )
  
           CASE ("ALd")
             l_cond = ( w_zn .lt. wdwncrit )
  
           CASE ("CLu")
             l_cond = ( (w_zn .gt. wupcrit) .and. (ql .gt. qlcrit) )
  
           CASE ("CLd")
             l_cond = ( (w_zn .lt. wdwncrit) .and. (ql .gt. qlcrit) )
  
           CASE ("AH")
             l_cond = ( (ql .gt. qlcrit) .or. (qi .gt. qicrit) .or. (qppt .gt. qpptcrit) )
  
           CASE ("AL")
             l_cond = ( ql .gt. qlcrit )
  
           CASE ("AI")
             l_cond = ( qi .gt. qicrit )
  
           CASE ("PPd")
             l_cond = ( (w_zn .lt. wdwncrit) .and. (qppt .gt. qpptcrit) )
  
           CASE ("VPd")
             l_cond = ( (w_zn .lt. wdwncrit) .and. (qppt .gt. vpptcrit) )
  
           CASE ("PVd")
             l_cond = ( (w_zn .lt. wsdwncrit) .and. (qppt .gt. qpptcrit) )
  
           CASE ("MO")
             l_cond = ( (w_zn .gt. wupcrit) .and. (qvl .gt. qsat) )
  
           CASE ("BM")
             l_cond = ( (w_zn .gt. wupcrit) .and. (thv_pr .gt. thvprcrit) .and. (qvl .gt. qsat) )
  
           CASE ("AA")
             l_cond = ( (ql .lt. qlcrit) .and. (qi .lt. qicrit) .and.   &
                        (w_zn .lt. wdwncrit) .and. (qppt .gt. qpptcrit)  )
  
           CASE ("AV")
             l_cond = ( (ql .lt. qlcrit) .and. (qi .lt. qicrit) .and.   &
                        (w_zn .lt. wdwncrit) .and. (qppt .gt. vpptcrit)  )
  
         END SELECT
  
             
         if (l_cond) then
           conddiags_tot(k,inc,:) = conddiags_tot(k,inc,:) + local_diag(diag_locations)
         else
           conddiags_tot(k,ncond+inc,:) = conddiags_tot(k,ncond+inc,:) + local_diag(diag_locations)
         end if
  
       end do ! inc loop over requested conditions
  
     enddo ! k loop over vertical levels
  

Here is the caller graph for this function:

Variable Documentation

◆ cond_locations

integer, dimension(:), allocatable conditional_diagnostics_column_mod::cond_locations

private

Definition at line 55 of file conditional_diagnostics_column.F90.

◆ cond_long

character(len=string_length), dimension(:), allocatable, public conditional_diagnostics_column_mod::cond_long

Definition at line 53 of file conditional_diagnostics_column.F90.

◆ cond_request

character(len=string_length), dimension(:), allocatable, public conditional_diagnostics_column_mod::cond_request

Definition at line 53 of file conditional_diagnostics_column.F90.

53 character(len=STRING_LENGTH), dimension(:), allocatable :: cond_request, cond_long

◆ conddiags_tot

real(kind=default_precision), dimension(:,:,:), allocatable, public conditional_diagnostics_column_mod::conddiags_tot

Definition at line 27 of file conditional_diagnostics_column.F90.

27 real(kind=default_precision), dimension(:,:,:), allocatable :: conddiags_tot

◆ diag_locations

integer, dimension(:), allocatable conditional_diagnostics_column_mod::diag_locations

private

Definition at line 55 of file conditional_diagnostics_column.F90.

55 integer, dimension(:), allocatable :: diag_locations, cond_locations

◆ diag_long

character(len=string_length), dimension(:), allocatable, public conditional_diagnostics_column_mod::diag_long

Definition at line 54 of file conditional_diagnostics_column.F90.

◆ diag_request

character(len=string_length), dimension(:), allocatable, public conditional_diagnostics_column_mod::diag_request

Definition at line 54 of file conditional_diagnostics_column.F90.

54 character(len=STRING_LENGTH), dimension(:), allocatable :: diag_request, diag_long

◆ diagnostic_generation_frequency

integer conditional_diagnostics_column_mod::diagnostic_generation_frequency

private

Definition at line 58 of file conditional_diagnostics_column.F90.

58 integer :: diagnostic_generation_frequency

◆ gpts_total

real(kind=default_precision), public conditional_diagnostics_column_mod::gpts_total

Definition at line 44 of file conditional_diagnostics_column.F90.

44 real(kind=default_precision) :: gpts_total

◆ iqg

integer conditional_diagnostics_column_mod::iqg

private

Definition at line 25 of file conditional_diagnostics_column.F90.

◆ iqi

integer conditional_diagnostics_column_mod::iqi

private

Definition at line 25 of file conditional_diagnostics_column.F90.

◆ iql

integer conditional_diagnostics_column_mod::iql

private

Definition at line 25 of file conditional_diagnostics_column.F90.

◆ iqr

integer conditional_diagnostics_column_mod::iqr

private

Definition at line 25 of file conditional_diagnostics_column.F90.

◆ iqs

integer conditional_diagnostics_column_mod::iqs

private

Definition at line 25 of file conditional_diagnostics_column.F90.

◆ iqv

integer conditional_diagnostics_column_mod::iqv

private

Definition at line 25 of file conditional_diagnostics_column.F90.

25 integer :: iqv, iql, iqr, iqi, iqs, iqg

◆ l_do_near

logical conditional_diagnostics_column_mod::l_do_near

private

Definition at line 48 of file conditional_diagnostics_column.F90.

48 logical :: l_do_near

◆ l_qi_qr_qs_qg

logical conditional_diagnostics_column_mod::l_qi_qr_qs_qg

private

Definition at line 47 of file conditional_diagnostics_column.F90.

47 logical :: l_qi_qr_qs_qg

◆ master_conditions_list

character(len=string_length), dimension(23) conditional_diagnostics_column_mod::master_conditions_list

private

Definition at line 50 of file conditional_diagnostics_column.F90.

50 character(len=STRING_LENGTH), dimension(23) :: master_conditions_list

◆ master_diagnostics_list

character(len=string_length), dimension(31) conditional_diagnostics_column_mod::master_diagnostics_list

private

Definition at line 51 of file conditional_diagnostics_column.F90.

51 character(len=STRING_LENGTH), dimension(31) :: master_diagnostics_list

◆ ncond

integer, public conditional_diagnostics_column_mod::ncond

Definition at line 39 of file conditional_diagnostics_column.F90.

39 integer :: ncond

◆ ndiag

integer, public conditional_diagnostics_column_mod::ndiag

Definition at line 40 of file conditional_diagnostics_column.F90.

40 integer :: ndiag

◆ qicrit

real(kind=default_precision) conditional_diagnostics_column_mod::qicrit

private

Definition at line 30 of file conditional_diagnostics_column.F90.

30 real(kind=default_precision) :: qicrit

◆ qlcrit

real(kind=default_precision) conditional_diagnostics_column_mod::qlcrit

private

Definition at line 29 of file conditional_diagnostics_column.F90.

29 real(kind=default_precision) :: qlcrit

◆ qpptcrit

real(kind=default_precision) conditional_diagnostics_column_mod::qpptcrit

private

Definition at line 31 of file conditional_diagnostics_column.F90.

31 real(kind=default_precision) :: qpptcrit

◆ requested_area

integer, public conditional_diagnostics_column_mod::requested_area

Definition at line 56 of file conditional_diagnostics_column.F90.

56 integer :: requested_area

◆ thv_from_th_with_liqice

logical conditional_diagnostics_column_mod::thv_from_th_with_liqice

private

Definition at line 46 of file conditional_diagnostics_column.F90.

46 logical :: thv_from_th_with_liqice

◆ thvprcrit

real(kind=default_precision) conditional_diagnostics_column_mod::thvprcrit

private

Definition at line 33 of file conditional_diagnostics_column.F90.

33 real(kind=default_precision) :: thvprcrit

◆ vpptcrit

real(kind=default_precision) conditional_diagnostics_column_mod::vpptcrit

private

Definition at line 32 of file conditional_diagnostics_column.F90.

32 real(kind=default_precision) :: vpptcrit

◆ wdwncrit

real(kind=default_precision) conditional_diagnostics_column_mod::wdwncrit

private

Definition at line 37 of file conditional_diagnostics_column.F90.

37 real(kind=default_precision) :: wdwncrit

◆ wsdwncrit

real(kind=default_precision) conditional_diagnostics_column_mod::wsdwncrit

private

Definition at line 34 of file conditional_diagnostics_column.F90.

34 real(kind=default_precision) :: wsdwncrit

◆ wsupcrit

real(kind=default_precision) conditional_diagnostics_column_mod::wsupcrit

private

Definition at line 35 of file conditional_diagnostics_column.F90.

35 real(kind=default_precision) :: wsupcrit

◆ wupcrit

real(kind=default_precision) conditional_diagnostics_column_mod::wupcrit

private

Definition at line 36 of file conditional_diagnostics_column.F90.

36 real(kind=default_precision) :: wupcrit

◆ x_size

integer conditional_diagnostics_column_mod::x_size

private

Definition at line 41 of file conditional_diagnostics_column.F90.

41 integer :: x_size

◆ y_size

integer conditional_diagnostics_column_mod::y_size

private

Definition at line 42 of file conditional_diagnostics_column.F90.

42 integer :: y_size

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ conditional_diagnostics_column_get_descriptor()

◆ field_information_retrieval_callback()

◆ field_value_retrieval_callback()

◆ finalisation_callback()

◆ initialisation_callback()

◆ timestep_callback()

Variable Documentation

◆ cond_locations

◆ cond_long

◆ cond_request

◆ conddiags_tot

◆ diag_locations

◆ diag_long

◆ diag_request

◆ diagnostic_generation_frequency

◆ gpts_total

◆ iqg

◆ iqi

◆ iql

◆ iqr

◆ iqs

◆ iqv

◆ l_do_near

◆ l_qi_qr_qs_qg

◆ master_conditions_list

◆ master_diagnostics_list

◆ ncond

◆ ndiag

◆ qicrit

◆ qlcrit

◆ qpptcrit

◆ requested_area

◆ thv_from_th_with_liqice

◆ thvprcrit

◆ vpptcrit

◆ wdwncrit

◆ wsdwncrit

◆ wsupcrit

◆ wupcrit

◆ x_size

◆ y_size