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

subroutine	initialisation_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	total_points

integer	ncl_col

integer	iqv =0

integer	iql =0

integer	iqr =0

integer	iqi =0

integer	iqs =0

integer	iqg =0

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

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

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

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

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

real(kind=default_precision)	qlcrit

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Function/Subroutine Documentation

◆ field_information_retrieval_callback()

subroutine scalar_diagnostics_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 322 of file scalar_diagnostics.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=2
     field_information%dimension_sizes(1)=current_state%local_grid%size(y_index)
     field_information%dimension_sizes(2)=current_state%local_grid%size(x_index)
     field_information%data_type=component_double_data_type
  
     if (name .eq. "senhf") then
       field_information%enabled=current_state%use_surface_boundary_conditions .and. current_state%th%active
     else if (name .eq. "lathf") then
       field_information%enabled=current_state%use_surface_boundary_conditions .and. &
            current_state%water_vapour_mixing_ratio_index .gt. 0 .and. &
            current_state%number_q_fields .ge. current_state%water_vapour_mixing_ratio_index
    else if (name .eq. "qlmax".or. name .eq. "cltop" .or. name .eq. "clbas") then
       field_information%enabled=.not. current_state%passive_q .and. current_state%liquid_water_mixing_ratio_index .gt. 0 &
            .and. current_state%number_q_fields .ge. current_state%liquid_water_mixing_ratio_index
     else if (name .eq. "vwp" .or. name .eq. "lwp") then
       field_information%enabled=current_state%number_q_fields .gt. 0 .and. current_state%water_vapour_mixing_ratio_index .gt. 0 &
            .and. current_state%number_q_fields .ge. current_state%water_vapour_mixing_ratio_index
    else if (name .eq. "rwp" ) then
       field_information%enabled= current_state%rain_water_mixing_ratio_index .gt. 0
    else if (name .eq. "iwp" .or. name .eq. 'tot_iwp') then
       field_information%enabled= current_state%ice_water_mixing_ratio_index .gt. 0
    else if (name .eq. "swp" ) then
       field_information%enabled= current_state%snow_water_mixing_ratio_index .gt. 0
    else if (name .eq. "gwp" ) then
       field_information%enabled= current_state%graupel_water_mixing_ratio_index .gt. 0
    else
       field_information%enabled=.true.
    end if
  

Here is the caller graph for this function:

◆ field_value_retrieval_callback()

subroutine scalar_diagnostics_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 363 of file scalar_diagnostics.F90.

     type(model_state_type), target, intent(inout) :: current_state
     character(len=*), intent(in) :: name
     type(component_field_value_type), intent(out) :: field_value
     
     integer :: i
  
     if (name .eq. "wmax") then
       allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index)))
        field_value%real_2d_array(:,:)=wmax(:,:)
     else if (name .eq. "wmin") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=wmin(:,:)
     else if (name .eq. 'reske') then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index)))
        field_value%real_2d_array(:,:)=reske(:,:)
     else if (name .eq. "qlmax") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=qlmax(:,:)
     !else if (name .eq. "hqlmax") then
     !   allocate(field_value%real_2d_array(current_state%local_grid%size(Y_INDEX), &
     !       current_state%local_grid%size(X_INDEX))) 
     !   field_value%real_2d_array(:,:)=hqlmax(:,:)
     else if (name .eq. "cltop") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=cltop(:,:)
     else if (name .eq. "clbas") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=clbas(:,:)
     else if (name .eq. "vwp") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=vwp(:,:)
     else if (name .eq. "lwp") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=lwp(:,:)
     else if (name .eq. "rwp") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=rwp(:,:)
     else if (name .eq. "iwp") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=iwp(:,:)
      else if (name .eq. "swp") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=swp(:,:)
     else if (name .eq. "gwp") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=gwp(:,:)
     else if (name .eq. "tot_iwp") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=tot_iwp(:,:)
     else if (name .eq. "senhf") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=senhf(:,:)
     else if (name .eq. "lathf") then
        allocate(field_value%real_2d_array(current_state%local_grid%size(y_index), &
            current_state%local_grid%size(x_index))) 
        field_value%real_2d_array(:,:)=lathf(:,:) 
     end if
     

Here is the caller graph for this function:

◆ initialisation_callback()

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

private

Definition at line 64 of file scalar_diagnostics.F90.

     type(model_state_type), target, intent(inout) :: current_state
  
     integer :: k
     integer :: y_size_local, x_size_local
     
     ! qlcrit declared in the config file
     qlcrit=options_get_real(current_state%options_database, "qlcrit")
  
     y_size_local = current_state%local_grid%size(y_index)
     x_size_local = current_state%local_grid%size(x_index)
     
     ! allocate scalar diagnostics as 2-D fields (horizontal slices) that 
     ! are published to the ioserver so that the ioserver can do the manipulation
     ! to obtain the scalar field
     allocate(wmax(y_size_local, x_size_local), wmin(y_size_local, x_size_local), &
          reske(y_size_local, x_size_local),                                      &
          senhf(y_size_local, x_size_local),lathf(y_size_local, x_size_local))
     ! allocate the 1d velocity arrays for the kinetic energy calc
     allocate(ww_prime_res(current_state%local_grid%size(z_index)), &
          uu_prime_res(current_state%local_grid%size(z_index)),     &
          vv_prime_res(current_state%local_grid%size(z_index)))
     ww_prime_res(:) = 0.0
     uu_prime_res(:) = 0.0
     vv_prime_res(:) = 0.0
  
     if (.not. current_state%passive_q .and. current_state%number_q_fields .gt. 0) then
        iqv = current_state%water_vapour_mixing_ratio_index
        iql = current_state%liquid_water_mixing_ratio_index
        allocate(vwp(y_size_local, x_size_local), lwp(y_size_local, x_size_local), &
             qlmax(y_size_local, x_size_local), hqlmax(y_size_local, x_size_local), &
             cltop(y_size_local, x_size_local), clbas(y_size_local, x_size_local))
        allocate(cloud_content(current_state%local_grid%size(z_index)))
        ! allocate other hydrometeors. Allocation dependent on index being set in
        ! appropriate microphysics scheme (see casim component from example)
        if (current_state%rain_water_mixing_ratio_index > 0) then
           iqr = current_state%rain_water_mixing_ratio_index
           allocate(rwp(y_size_local, x_size_local))
        endif
        if (current_state%ice_water_mixing_ratio_index > 0) then
           iqi = current_state%ice_water_mixing_ratio_index 
           allocate(iwp(y_size_local, x_size_local))
           allocate(tot_iwp(y_size_local, x_size_local))
        endif
        if (current_state%snow_water_mixing_ratio_index > 0) then
           iqs = current_state%snow_water_mixing_ratio_index
           allocate(swp(y_size_local, x_size_local))
        endif
        if (current_state%graupel_water_mixing_ratio_index > 0) then
           iqg = current_state%graupel_water_mixing_ratio_index
           allocate(gwp(y_size_local, x_size_local))
        endif
     endif
     
     allocate(dz_rhon_fac(current_state%local_grid%size(z_index)))    
     do k=2, current_state%local_grid%size(z_index)
        ! used in the water path calculation
        dz_rhon_fac(k)=current_state%global_grid%configuration%vertical%dz(k)*&
             current_state%global_grid%configuration%vertical%rhon(k)
     end do    
  

Here is the caller graph for this function:

◆ scalar_diagnostics_get_descriptor()

type(component_descriptor_type) function, public scalar_diagnostics_mod::scalar_diagnostics_get_descriptor

Provides the component descriptor for the core to register.

Returns: The descriptor describing this component

Definition at line 34 of file scalar_diagnostics.F90.

     scalar_diagnostics_get_descriptor%name="scalar_diagnostics"
     scalar_diagnostics_get_descriptor%version=0.1
  
     scalar_diagnostics_get_descriptor%initialisation=>initialisation_callback
     scalar_diagnostics_get_descriptor%timestep=>timestep_callback
  
     scalar_diagnostics_get_descriptor%field_value_retrieval=>field_value_retrieval_callback
     scalar_diagnostics_get_descriptor%field_information_retrieval=>field_information_retrieval_callback
     allocate(scalar_diagnostics_get_descriptor%published_fields(16))
  
     scalar_diagnostics_get_descriptor%published_fields(1)="vwp"
     scalar_diagnostics_get_descriptor%published_fields(2)="lwp"
     scalar_diagnostics_get_descriptor%published_fields(3)="qlmax"
     scalar_diagnostics_get_descriptor%published_fields(4)="hqlmax"
     scalar_diagnostics_get_descriptor%published_fields(5)="cltop"
     scalar_diagnostics_get_descriptor%published_fields(6)="clbas"
     scalar_diagnostics_get_descriptor%published_fields(7)="wmax"
     scalar_diagnostics_get_descriptor%published_fields(8)="wmin"
     scalar_diagnostics_get_descriptor%published_fields(9)="senhf"
     scalar_diagnostics_get_descriptor%published_fields(10)="lathf"
     scalar_diagnostics_get_descriptor%published_fields(11)="rwp"
     scalar_diagnostics_get_descriptor%published_fields(12)="iwp"
     scalar_diagnostics_get_descriptor%published_fields(13)="swp"
     scalar_diagnostics_get_descriptor%published_fields(14)="gwp"
     scalar_diagnostics_get_descriptor%published_fields(15)="tot_iwp"
     scalar_diagnostics_get_descriptor%published_fields(16)="reske"
     

Here is the call graph for this function:

◆ timestep_callback()

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

private

Include ice if present.

Definition at line 127 of file scalar_diagnostics.F90.

     type(model_state_type), target, intent(inout) :: current_state
  
     integer :: k, i
     integer :: current_y_index, current_x_index, target_x_index, target_y_index
  
     current_y_index=current_state%column_local_y
     current_x_index=current_state%column_local_x
     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 (current_state%first_timestep_column) then
        ! maximum vertical velocity for each column
        wmax(:,:)=0.0
        ! minimum vertical velocity for each column
        wmin(:,:)=0.0
        ! resolved ke
        reske(:,:) = 0.0
        
        if (.not. current_state%passive_q .and. current_state%number_q_fields .gt. 0) then
           ! maximum liquid water content in a column
           qlmax(:,:)=0.0
           ! the height of the maximum liquid water content in a column
           hqlmax(:,:)=0.0
           ! cloud top height where liqud water content is greater than qlcrit
           cltop(:,:)=0.0
           ! minimum cloud base where liquid water content is greater than qlcrit
           clbas(:,:)=0.0
           ! water vapour path for each column
           vwp(:,:)=0.0
           ! liquid water path for each column
           lwp(:,:)=0.0
           ! rain water path for each column
           if (current_state%rain_water_mixing_ratio_index > 0) rwp(:,:)=0.0
           ! ice water path for each column
           if (current_state%ice_water_mixing_ratio_index > 0) then
              iwp(:,:)=0.0
              ! total ice water path (iwp + swp + gwp) for each column
              tot_iwp(:,:)=0.0
           endif
           ! snow water path for each column
           if (current_state%snow_water_mixing_ratio_index > 0) swp(:,:)=0.0
           ! graupel water path for each column
           if (current_state%graupel_water_mixing_ratio_index > 0) gwp(:,:)=0.0
        endif
      
        ! surface sensible heat flux
        senhf(:,:)=0.0
        ! surface latent heat flux
        lathf(:,:)=0.0
     end if
     
     if (.not. current_state%halo_column) then
        ! maximum and minimum vertical velocity in each column
        wmax(target_y_index, target_x_index)=maxval(current_state%w%data(:, current_y_index, current_x_index))
        wmin(target_y_index, target_x_index)=minval(current_state%w%data(:, current_y_index, current_x_index))
                
        ! work out the column resolved ww, uu, vv
        ww_prime_res(:) = &
             (current_state%w%data(:,current_state%column_local_y,current_state%column_local_x)**2.)
        if (allocated(current_state%global_grid%configuration%vertical%olubar)) then
           uu_prime_res(:) = &
                ((current_state%u%data(:,current_state%column_local_y,current_state%column_local_x) &
                - (current_state%global_grid%configuration%vertical%olubar(:) - current_state%ugal))**2.)
        else
           uu_prime_res(:) = 0.0
        endif
        if (allocated(current_state%global_grid%configuration%vertical%olvbar)) then
           vv_prime_res(:) = &
                ((current_state%v%data(:,current_state%column_local_y,current_state%column_local_x) &
                - (current_state%global_grid%configuration%vertical%olvbar(:) - current_state%vgal))**2.)
        else
           vv_prime_res(:) = 0.0
        endif
        ! use column resolved ww, uu, vv to derive total resolved KE for column
        do k=2, current_state%local_grid%size(z_index)-1
           reske(target_y_index, target_x_index) = reske(target_y_index, target_x_index) + &
                (uu_prime_res(k) + uu_prime_res(k+1)                                       &
                + vv_prime_res(k) + vv_prime_res(k+1)                                      &
                + 2_default_precision * ww_prime_res(k))*0.25_default_precision*           &
                current_state%global_grid%configuration%vertical%dzn(k+1)*                 &
                current_state%global_grid%configuration%vertical%rho(k)
        enddo
        ! divide reske by altitude to make it column mean (as in the LEM)
        reske(target_y_index, target_x_index) = reske(target_y_index, target_x_index)/ &
             current_state%global_grid%configuration%vertical%z(current_state%local_grid%size(z_index))
        ! subke is derive in the subgrid_profile_diagnostics component
        
        if (.not. current_state%passive_q .and. current_state%number_q_fields .gt. 0) then
           !
           ! calculate the lwc maximum and height of lwc max for each column
           !
           if (current_state%liquid_water_mixing_ratio_index .gt. 0 .and. &
                current_state%number_q_fields .ge. current_state%liquid_water_mixing_ratio_index) then
              qlmax(target_y_index, target_x_index) = &
                   maxval(current_state%q(current_state%liquid_water_mixing_ratio_index)%data &
                   (:,current_y_index, current_x_index))
              !hqlmax(current_y_index, current_x_index) = &
              !     current_state%global_grid%configuration%vertical% &
              !     zn(maxloc(current_state%q(current_state%liquid_water_mixing_ratio_index)%data &
              !     (:,current_y_index, current_x_index)))
              
  
              ! calculate the cloud top maximum and minimum for each column 
              !
              cloud_content(:) = current_state%q(current_state%liquid_water_mixing_ratio_index)%data( &
                                                                      :,current_y_index, current_x_index)
              if (current_state%ice_water_mixing_ratio_index .gt. 0)  &
                  cloud_content(:) = cloud_content(:) + current_state%q(current_state%ice_water_mixing_ratio_index)%data( &
                                                                                       :,current_y_index, current_x_index)
              do k = 2, current_state%local_grid%size(z_index)
                if (cloud_content(k) .gt. qlcrit) then
                  cltop(target_y_index, target_x_index) = &
                        current_state%global_grid%configuration%vertical%zn(k)
                endif
                if (cloud_content(current_state%local_grid%size(z_index)+1-k) .gt. qlcrit) then
                   clbas(target_y_index, target_x_index) = &
                        current_state%global_grid%configuration%vertical%zn(current_state%local_grid%size(z_index)+1-k)
                end if
              enddo ! k loop over height
           endif
           !
           ! calculate the vapour and liquid water path
           !
           if (current_state%water_vapour_mixing_ratio_index .gt. 0 .and. &
                current_state%number_q_fields .ge. current_state%water_vapour_mixing_ratio_index) then
              do k = 2, current_state%local_grid%size(z_index)
                 vwp(target_y_index, target_x_index)=vwp(target_y_index, target_x_index) &
                      +dz_rhon_fac(k)*current_state%q(current_state%water_vapour_mixing_ratio_index)%data(k, &
                      current_y_index, current_x_index)
                 lwp(target_y_index, target_x_index)=lwp(target_y_index, target_x_index) &
                      +dz_rhon_fac(k)*current_state%q(current_state%liquid_water_mixing_ratio_index)%data(k, &
                      current_y_index, current_x_index)
              enddo
              if (current_state%rain_water_mixing_ratio_index > 0) then
                 do k = 2, current_state%local_grid%size(z_index)
                    rwp(target_y_index, target_x_index)=rwp(target_y_index, target_x_index) &
                         +dz_rhon_fac(k)*current_state%q(iqr)%data(k,current_y_index, current_x_index)
                 enddo
              endif
              if (current_state%ice_water_mixing_ratio_index > 0) then
                 do k = 2, current_state%local_grid%size(z_index)
                    iwp(target_y_index, target_x_index)=iwp(target_y_index, target_x_index) &
                         +dz_rhon_fac(k)*current_state%q(iqi)%data(k,current_y_index, current_x_index)
                 enddo
                 tot_iwp(target_y_index, target_x_index)=tot_iwp(target_y_index, target_x_index)+ &
                      iwp(target_y_index, target_x_index)
              endif
              if (current_state%snow_water_mixing_ratio_index > 0) then
                 do k = 2, current_state%local_grid%size(z_index)
                    swp(target_y_index, target_x_index)=swp(target_y_index, target_x_index) &
                         +dz_rhon_fac(k)*current_state%q(iqs)%data(k,current_y_index, current_x_index)
                 enddo
                 tot_iwp(target_y_index, target_x_index)=tot_iwp(target_y_index, target_x_index)+ &
                      swp(target_y_index, target_x_index)
              endif
              if (current_state%graupel_water_mixing_ratio_index > 0) then
                 do k = 2, current_state%local_grid%size(z_index)
                    gwp(target_y_index, target_x_index)=gwp(target_y_index, target_x_index) &
                         +dz_rhon_fac(k)*current_state%q(iqg)%data(k,current_y_index, current_x_index)
                 enddo
                 tot_iwp(target_y_index, target_x_index)=tot_iwp(target_y_index, target_x_index)+ &
                      gwp(target_y_index, target_x_index)
              endif
           end if
        end if
  
        ! surface flux diagnostics
        if (current_state%use_surface_boundary_conditions) then
           if (current_state%water_vapour_mixing_ratio_index .gt. 0 .and. &
                current_state%number_q_fields .ge. current_state%water_vapour_mixing_ratio_index) then 
              lathf(target_y_index, target_x_index)= (current_state%diff_coefficient%data(1, current_y_index, current_x_index) *   &                   
                   current_state%global_grid%configuration%vertical%rdzn(2) *  & 
                   (current_state%q(current_state%water_vapour_mixing_ratio_index)%data(1,current_y_index,current_x_index) - & 
                   current_state%q(current_state%water_vapour_mixing_ratio_index)%data(2,current_y_index,current_x_index))) &
                   * rlvap * current_state%global_grid%configuration%vertical%rhon(1)
           endif
  
       if (current_state%th%active) then      
               senhf(target_y_index, target_x_index)=(current_state%diff_coefficient%data(1, current_y_index, current_x_index)  &            
                    * current_state%global_grid%configuration%vertical%rdzn(2)  &            
                    * (current_state%th%data(1, current_y_index, current_x_index) &            
                    - current_state%th%data(2, current_y_index, current_x_index) &            
                    - current_state%global_grid%configuration%vertical%dthref(1))) &
                    * current_state%global_grid%configuration%vertical%rhon(1)*cp
           endif
        endif      
     end if

Here is the caller graph for this function:

Variable Documentation

◆ clbas

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::clbas

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ cloud_content

real(kind=default_precision), dimension(:), allocatable scalar_diagnostics_mod::cloud_content

private

Definition at line 21 of file scalar_diagnostics.F90.

◆ cltop

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::cltop

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ dz_rhon_fac

real(kind=default_precision), dimension(:), allocatable scalar_diagnostics_mod::dz_rhon_fac

private

Definition at line 20 of file scalar_diagnostics.F90.

20 real(kind=default_precision), dimension(:), allocatable :: dz_rhon_fac

◆ gwp

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::gwp

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ hqlmax

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::hqlmax

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ iqg

integer scalar_diagnostics_mod::iqg =0

private

Definition at line 18 of file scalar_diagnostics.F90.

◆ iqi

integer scalar_diagnostics_mod::iqi =0

private

Definition at line 18 of file scalar_diagnostics.F90.

◆ iql

integer scalar_diagnostics_mod::iql =0

private

Definition at line 18 of file scalar_diagnostics.F90.

◆ iqr

integer scalar_diagnostics_mod::iqr =0

private

Definition at line 18 of file scalar_diagnostics.F90.

◆ iqs

integer scalar_diagnostics_mod::iqs =0

private

Definition at line 18 of file scalar_diagnostics.F90.

◆ iqv

integer scalar_diagnostics_mod::iqv =0

private

Definition at line 18 of file scalar_diagnostics.F90.

◆ iwp

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::iwp

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ lathf

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::lathf

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ lwp

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::lwp

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ ncl_col

integer scalar_diagnostics_mod::ncl_col

private

Definition at line 18 of file scalar_diagnostics.F90.

◆ qlcrit

real(kind=default_precision) scalar_diagnostics_mod::qlcrit

private

Definition at line 23 of file scalar_diagnostics.F90.

23 real(kind=default_precision) :: qlcrit

◆ qlmax

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::qlmax

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ reske

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::reske

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ rwp

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::rwp

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ senhf

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::senhf

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ swp

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::swp

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ tot_iwp

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::tot_iwp

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ total_points

integer scalar_diagnostics_mod::total_points

private

Definition at line 18 of file scalar_diagnostics.F90.

18 integer :: total_points, ncl_col, iqv=0, iql=0, iqr=0, iqi=0, iqs=0, &

19 iqg=0

◆ uu_prime_res

real(kind=default_precision), dimension(:), allocatable scalar_diagnostics_mod::uu_prime_res

private

Definition at line 21 of file scalar_diagnostics.F90.

◆ vv_prime_res

real(kind=default_precision), dimension(:), allocatable scalar_diagnostics_mod::vv_prime_res

private

Definition at line 21 of file scalar_diagnostics.F90.

◆ vwp

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::vwp

private

Definition at line 24 of file scalar_diagnostics.F90.

   real(kind=default_precision), dimension(:,:), allocatable :: vwp, lwp, wmax, wmin, &
        qlmax, hqlmax, cltop, clbas,  senhf, lathf, rwp, iwp, swp, gwp, tot_iwp,      &
        reske

◆ wmax

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::wmax

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ wmin

real(kind=default_precision), dimension(:,:), allocatable scalar_diagnostics_mod::wmin

private

Definition at line 24 of file scalar_diagnostics.F90.

◆ ww_prime_res

real(kind=default_precision), dimension(:), allocatable scalar_diagnostics_mod::ww_prime_res

private

Definition at line 21 of file scalar_diagnostics.F90.

21 real(kind=default_precision), dimension(:), allocatable :: ww_prime_res, uu_prime_res, &

22 vv_prime_res, cloud_content

Functions/Subroutines

Variables

Function/Subroutine Documentation

◆ field_information_retrieval_callback()

◆ field_value_retrieval_callback()

◆ initialisation_callback()

◆ scalar_diagnostics_get_descriptor()

◆ timestep_callback()

Variable Documentation

◆ clbas

◆ cloud_content

◆ cltop

◆ dz_rhon_fac

◆ gwp

◆ hqlmax

◆ iqg

◆ iqi

◆ iql

◆ iqr

◆ iqs

◆ iqv

◆ iwp

◆ lathf

◆ lwp

◆ ncl_col

◆ qlcrit

◆ qlmax

◆ reske

◆ rwp

◆ senhf

◆ swp

◆ tot_iwp

◆ total_points

◆ uu_prime_res

◆ vv_prime_res

◆ vwp

◆ wmax

◆ wmin

◆ ww_prime_res