Implimentation of CASIM microphysics. More...

Functions/Subroutines
type(component_descriptor_type) function, public	casim_get_descriptor ()
	Provides the descriptor back to the caller and is used in component registration. More...

subroutine	initialisation_callback (current_state)
	The initialisation callback sets up the microphysics. More...

subroutine	timestep_callback (current_state)
	Called for each column per timestep this will calculate the microphysical tendencies. More...

subroutine	read_configuration (current_state)

subroutine	field_information_retrieval_callback (current_state, name, field_information)

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
real(wp), dimension(:,:,:), allocatable	theta

real(wp), dimension(:,:,:), allocatable	pressure

real(wp), dimension(:,:,:), allocatable	z_half

real(wp), dimension(:,:,:), allocatable	z_centre

real(wp), dimension(:,:,:), allocatable	dz

real(wp), dimension(:,:,:), allocatable	qv

real(wp), dimension(:,:,:), allocatable	qc

real(wp), dimension(:,:,:), allocatable	nc

real(wp), dimension(:,:,:), allocatable	qr

real(wp), dimension(:,:,:), allocatable	nr

real(wp), dimension(:,:,:), allocatable	m3r

real(wp), dimension(:,:,:), allocatable	rho

real(wp), dimension(:,:,:), allocatable	exner

real(wp), dimension(:,:,:), allocatable	w

real(wp), dimension(:,:,:), allocatable	tke

real(wp), dimension(:,:,:), allocatable	qi

real(wp), dimension(:,:,:), allocatable	ni

real(wp), dimension(:,:,:), allocatable	qs

real(wp), dimension(:,:,:), allocatable	ns

real(wp), dimension(:,:,:), allocatable	m3s

real(wp), dimension(:,:,:), allocatable	qg

real(wp), dimension(:,:,:), allocatable	ng

real(wp), dimension(:,:,:), allocatable	m3g

real(wp), dimension(:,:,:), allocatable	accumsolmass

real(wp), dimension(:,:,:), allocatable	accumsolnumber

real(wp), dimension(:,:,:), allocatable	activesolliquid

real(wp), dimension(:,:,:), allocatable	aitkensolmass

real(wp), dimension(:,:,:), allocatable	aitkensolnumber

real(wp), dimension(:,:,:), allocatable	coarsesolmass

real(wp), dimension(:,:,:), allocatable	coarsesolnumber

real(wp), dimension(:,:,:), allocatable	activesolrain

real(wp), dimension(:,:,:), allocatable	coarsedustmass

real(wp), dimension(:,:,:), allocatable	coarsedustnumber

real(wp), dimension(:,:,:), allocatable	activeinsolice

real(wp), dimension(:,:,:), allocatable	activesolice

real(wp), dimension(:,:,:), allocatable	activeinsolliquid

real(wp), dimension(:,:,:), allocatable	accuminsolmass

real(wp), dimension(:,:,:), allocatable	accuminsolnumber

real(wp), dimension(:,:,:), allocatable	activesolnumber

real(wp), dimension(:,:,:), allocatable	activeinsolnumber

real(wp), dimension(:,:,:), allocatable	aitkensolbk

real(wp), dimension(:,:,:), allocatable	accumsolbk

real(wp), dimension(:,:,:), allocatable	coarsesolbk

real(wp), dimension(:,:,:), allocatable	dqv

real(wp), dimension(:,:,:), allocatable	dth

real(wp), dimension(:,:,:), allocatable	dqc

real(wp), dimension(:,:,:), allocatable	dnc

real(wp), dimension(:,:,:), allocatable	dqr

real(wp), dimension(:,:,:), allocatable	dnr

real(wp), dimension(:,:,:), allocatable	dm3r

real(wp), dimension(:,:,:), allocatable	dqi

real(wp), dimension(:,:,:), allocatable	dni

real(wp), dimension(:,:,:), allocatable	dqs

real(wp), dimension(:,:,:), allocatable	dns

real(wp), dimension(:,:,:), allocatable	dm3s

real(wp), dimension(:,:,:), allocatable	dqg

real(wp), dimension(:,:,:), allocatable	dng

real(wp), dimension(:,:,:), allocatable	dm3g

real(wp), dimension(:,:,:), allocatable	daccumsolmass

real(wp), dimension(:,:,:), allocatable	daccumsolnumber

real(wp), dimension(:,:,:), allocatable	dactivesolliquid

real(wp), dimension(:,:,:), allocatable	daitkensolmass

real(wp), dimension(:,:,:), allocatable	daitkensolnumber

real(wp), dimension(:,:,:), allocatable	dcoarsesolmass

real(wp), dimension(:,:,:), allocatable	dcoarsesolnumber

real(wp), dimension(:,:,:), allocatable	dactivesolrain

real(wp), dimension(:,:,:), allocatable	dcoarsedustmass

real(wp), dimension(:,:,:), allocatable	dcoarsedustnumber

real(wp), dimension(:,:,:), allocatable	dactiveinsolice

real(wp), dimension(:,:,:), allocatable	dactivesolice

real(wp), dimension(:,:,:), allocatable	dactiveinsolliquid

real(wp), dimension(:,:,:), allocatable	daccuminsolmass

real(wp), dimension(:,:,:), allocatable	daccuminsolnumber

real(wp), dimension(:,:,:), allocatable	dactivesolnumber

real(wp), dimension(:,:,:), allocatable	dactiveinsolnumber

real(wp), dimension(:,:), allocatable	surface_precip

integer	ils

integer	ile

integer	jls

integer	jle

integer	kls

integer	kle

integer	its

integer	ite

integer	jts

integer	jte

integer	kts

integer	kte

integer	iqv =0

integer	iql =0

integer	iqr =0

integer	iqi =0

integer	iqs =0

integer	iqg =0

integer	inl =0

integer	inr =0

integer	ini =0

integer	ins =0

integer	ing =0

integer	i3mr =0

integer	i3ms =0

integer	i3mg =0

integer	i_accumsolmass =0

integer	i_accumsolnumber =0

integer	i_activesolliquid =0

integer	i_aitkensolmass =0

integer	i_aitkensolnumber =0

integer	i_coarsesolmass =0

integer	i_coarsesolnumber =0

integer	i_activesolrain =0

integer	i_coarsedustmass =0

integer	i_coarsedustnumber =0

integer	i_activeinsolice =0

integer	i_activesolice =0

integer	i_activeinsolliquid =0

integer	i_accuminsolmass =0

integer	i_accuminsolnumber =0

integer	i_activesolnumber =0

integer	i_activeinsolnumber =0

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Detailed Description

Implimentation of CASIM microphysics.

Dummy stub when not compiling with CASIM microphysics.

Function/Subroutine Documentation

◆ casim_get_descriptor()

type(component_descriptor_type) function, public casim_mod::casim_get_descriptor

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

Returns: The termination check component descriptor

Definition at line 178 of file casim.F90.

     casim_get_descriptor%name="casim"
     casim_get_descriptor%version=0.1
     casim_get_descriptor%initialisation=>initialisation_callback
     casim_get_descriptor%timestep=>timestep_callback
  
     ! Set up fields to be published for diagnostics. These are made available to
     ! IO server in the field_value_retrieval_callback at the end of this module
     !
     casim_get_descriptor%field_value_retrieval=>field_value_retrieval_callback
     casim_get_descriptor%field_information_retrieval=>field_information_retrieval_callback
     
     allocate(casim_get_descriptor%published_fields(27))
     
     casim_get_descriptor%published_fields(1)="surface_precip"    
     casim_get_descriptor%published_fields(2)="homogeneous_freezing_rate"
     casim_get_descriptor%published_fields(3)="ice_nucleations_rate"
     casim_get_descriptor%published_fields(4)="ice_deposition_rate"
     casim_get_descriptor%published_fields(5)="snow_deposition_rate"
     casim_get_descriptor%published_fields(6)="ice_acc_cloud_rate"
     casim_get_descriptor%published_fields(7)="snow_acc_cloud_rate"
     casim_get_descriptor%published_fields(8)="snow_acc_rain_rate"
     casim_get_descriptor%published_fields(9)="ice_sublime_rate"
     casim_get_descriptor%published_fields(10)="snow_sublime_rate"
     casim_get_descriptor%published_fields(11)="ice_melt_rate"
     casim_get_descriptor%published_fields(12)="snow_melt_rate"
     casim_get_descriptor%published_fields(13)="snow_autoconversion_rate"
     casim_get_descriptor%published_fields(14)="snow_acc_ice_rate"
     casim_get_descriptor%published_fields(15)="rain_autoconversion_rate"
     casim_get_descriptor%published_fields(16)="rain_acc_cloud_rate"
     casim_get_descriptor%published_fields(17)="rain_evap_rate"
     casim_get_descriptor%published_fields(18)="graup_acc_cloud_rate"
     casim_get_descriptor%published_fields(19)="graup_acc_snow_rate"
     casim_get_descriptor%published_fields(20)="graup_melt_rate"
     casim_get_descriptor%published_fields(21)="graup_sublime_rate"
     casim_get_descriptor%published_fields(22)="ice_sed_rate"
     casim_get_descriptor%published_fields(23)="snow_sed_rate"
     casim_get_descriptor%published_fields(24)="rain_sed_rate"
     casim_get_descriptor%published_fields(25)="graup_sed_rate"
     casim_get_descriptor%published_fields(26)="cloud_sed_rate"
     casim_get_descriptor%published_fields(27)="condensation_rate"
     

Here is the call graph for this function:

◆ field_information_retrieval_callback()

subroutine casim_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
	)

Definition at line 953 of file casim.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%data_type=component_double_data_type
     if (name .eq. "surface_precip") then
        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)
     !else if (name .eq. "pcond_total" .or. name .eq. "psedl_total") then
     !   field_information%number_dimensions=1
     !   field_information%dimension_sizes(1)=current_state%local_grid%size(Z_INDEX)
     else
        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)
     endif
        
     field_information%enabled=.true.    
  

Here is the caller graph for this function:

◆ field_value_retrieval_callback()

subroutine casim_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 982 of file casim.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. "surface_precip") 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(:,:)= surface_precip(:,:)
     else if (name .eq. "condensation_rate") then
        allocate(field_value%real_3d_array(current_state%local_grid%size(z_index),  &
             current_state%local_grid%size(y_index),                                &
             current_state%local_grid%size(x_index)))
        field_value%real_3d_array(:,:,:) = casim_monc_dgs % pcond(:,:,:)
 !!$    else if (name .eq. "pcond_total") then
 !!$       allocate(field_value%real_1d_array(current_state%local_grid%size(Z_INDEX)))
 !!$       field_value%real_1d_array(:)=pcond_tot(:)
     end if
     

Here is the caller graph for this function:

◆ initialisation_callback()

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

private

The initialisation callback sets up the microphysics.

Parameters

current_state The current model state

Set up and allocate the local arrays

Definition at line 224 of file casim.F90.

     type(model_state_type), target, intent(inout) :: current_state
  
     integer :: y_size_local, x_size_local
  
     if (is_component_enabled(current_state%options_database, "simplecloud")) then
       call log_master_log(log_error, "Casim and Simplecloud are enabled, this does not work yet. Please disable one")
    end if
  
    !allocate(psedl_tot(current_state%local_grid%size(Z_INDEX)), &
    !      pcond_tot(current_state%local_grid%size(Z_INDEX)))
     
     y_size_local = current_state%local_grid%size(y_index)
     x_size_local = current_state%local_grid%size(x_index)
     
     call read_configuration(current_state)   
  
     ils=1
     ile=1
     jls=1
     jle=1
     kls=2
     kle=current_state%local_grid%size(z_index)
     its=1
     ite=1
     jts=1
     jte=1
     kts=1
     kte=current_state%local_grid%size(z_index)
  
  
     allocate(pressure(kte,1,1))
     allocate(z_half(0:kte,1,1))
     allocate(z_centre(kte,1,1))
     allocate(dz(kte,1,1))
     allocate(rho(kte,1,1))
     allocate(exner(kte,1,1))
     allocate(w(kte,1,1))
     allocate(tke(kte,1,1))
  
     allocate(theta(kte,1,1))
     allocate(qv(kte,1,1))
     allocate(qc(kte,1,1))
     allocate(nc(kte,1,1))
     allocate(qr(kte,1,1))
     allocate(nr(kte,1,1))
     allocate(m3r(kte,1,1))
     allocate(qi(kte,1,1))
     allocate(ni(kte,1,1))
     allocate(qs(kte,1,1))
     allocate(ns(kte,1,1))
     allocate(m3s(kte,1,1))
     allocate(qg(kte,1,1))
     allocate(ng(kte,1,1))
     allocate(m3g(kte,1,1))
  
     allocate(accumsolmass(kte,1,1))
     allocate(accumsolnumber(kte,1,1))
     allocate(activesolliquid(kte,1,1))
     allocate(aitkensolmass(kte,1,1))
     allocate(aitkensolnumber(kte,1,1))
     allocate(coarsesolmass(kte,1,1))
     allocate(coarsesolnumber(kte,1,1))
     allocate(activesolrain(kte,1,1))
     allocate(coarsedustmass(kte,1,1))
     allocate(coarsedustnumber(kte,1,1))
     allocate(activeinsolice(kte,1,1))
     allocate(activesolice(kte,1,1))
     allocate(activeinsolliquid(kte,1,1))
     allocate(accuminsolmass(kte,1,1))
     allocate(accuminsolnumber(kte,1,1))
     allocate(activesolnumber(kte,1,1))
     allocate(activeinsolnumber(kte,1,1))
     ! allocate the hygoscopicity arrays
     allocate(aitkensolbk(kte,1,1))
     allocate(accumsolbk(kte,1,1))
     allocate(coarsesolbk(kte,1,1))
  
     allocate(dth(kte,1,1))
     allocate(dqv(kte,1,1))
     allocate(dqc(kte,1,1))
     allocate(dnc(kte,1,1))
     allocate(dqr(kte,1,1))
     allocate(dnr(kte,1,1))
     allocate(dm3r(kte,1,1))
     allocate(dqi(kte,1,1))
     allocate(dni(kte,1,1))
     allocate(dqs(kte,1,1))
     allocate(dns(kte,1,1))
     allocate(dm3s(kte,1,1))
     allocate(dqg(kte,1,1))
     allocate(dng(kte,1,1))
     allocate(dm3g(kte,1,1))
  
     allocate(daccumsolmass(kte,1,1))
     allocate(daccumsolnumber(kte,1,1))
     allocate(dactivesolliquid(kte,1,1))
     allocate(daitkensolmass(kte,1,1))
     allocate(daitkensolnumber(kte,1,1))
     allocate(dcoarsesolmass(kte,1,1))
     allocate(dcoarsesolnumber(kte,1,1))
     allocate(dactivesolrain(kte,1,1))
     allocate(dcoarsedustmass(kte,1,1))
     allocate(dcoarsedustnumber(kte,1,1))
     allocate(dactiveinsolice(kte,1,1))
     allocate(dactivesolice(kte,1,1))
     allocate(dactiveinsolliquid(kte,1,1))
     allocate(daccuminsolmass(kte,1,1))
     allocate(daccuminsolnumber(kte,1,1))
     allocate(dactivesolnumber(kte,1,1))
     allocate(dactiveinsolnumber(kte,1,1))
  
     call set_mphys_switches(option,aerosol_option)
     call mphys_init(its, ite, jts, jte, kts, kte, ils, ile, jls, jle, kls, kle, l_tendency=.true.)
  
     ! Need to allocate the appropriate indices, e.g. iqv, iql...
     ! This needs to be compatible with the rest of the model
     ! This essentially reproduces the switching in the main microphysics
     ! code already done above (set_mphys_switches), so could be combined 
     ! once the MONC method has been finalized.
     ! Note the numbers assigned here may be different from those assigned 
     ! in the microphysics since we share the q array with other components. 
  
     if (.not. allocated(current_state%cq))then
        allocate(current_state%cq(current_state%number_q_fields))
        current_state%cq=0.0_default_precision
     end if
  
     ! Mass
     iqv = get_q_index(standard_q_names%VAPOUR, 'casim')
     if (nq_l>0)then
        iql = get_q_index(standard_q_names%CLOUD_LIQUID_MASS, 'casim')
        current_state%cq(iql) = -1.0
     end if
     if (nq_r>0)then
        iqr = get_q_index(standard_q_names%RAIN_MASS, 'casim')
        current_state%rain_water_mixing_ratio_index=iqr
       current_state%cq(iqr) = -1.0
     end if
     if (.not. l_warm)then
       if (nq_i>0)then
          iqi = get_q_index(standard_q_names%ICE_MASS, 'casim')
          current_state%ice_water_mixing_ratio_index=iqi
         current_state%cq(iqi) = -1.0
       end if
       if (nq_s>0)then
          iqs = get_q_index(standard_q_names%SNOW_MASS, 'casim')
          current_state%snow_water_mixing_ratio_index=iqs
         current_state%cq(iqs) = -1.0
       end if
       if (nq_g>0)then
          iqg = get_q_index(standard_q_names%GRAUPEL_MASS, 'casim')
          current_state%graupel_water_mixing_ratio_index=iqg
         current_state%cq(iqg) = -1.0
       end if
     end if
  
     ! Number
     if (l_2mc)inl = get_q_index(standard_q_names%CLOUD_LIQUID_NUMBER, 'casim')
     if (l_2mr)inr = get_q_index(standard_q_names%RAIN_NUMBER, 'casim')
     if (.not. l_warm)then
       if (l_2mi)ini = get_q_index(standard_q_names%ICE_NUMBER, 'casim')
       if (l_2ms)ins = get_q_index(standard_q_names%SNOW_NUMBER, 'casim')
       if (l_2mg)ing = get_q_index(standard_q_names%GRAUPEL_NUMBER, 'casim')
     end if
  
     ! Third moments
     if (l_3mr)i3mr = get_q_index(standard_q_names%RAIN_THIRD_MOMENT, 'casim')
     if (.not. l_warm)then
       if (l_3ms)i3ms = get_q_index(standard_q_names%SNOW_THIRD_MOMENT, 'casim')
       if (l_3mg)i3mg = get_q_index(standard_q_names%GRAUPEL_THIRD_MOMENT, 'casim')
     end if
  
     ! Aerosol
     if (soluble_modes(1) > 1)   i_aitkensolmass     = &
        get_q_index(standard_q_names%AITKEN_SOL_MASS,     'casim') 
     if (soluble_modes(1) > 0)   i_aitkensolnumber   = &
        get_q_index(standard_q_names%AITKEN_SOL_NUMBER,   'casim') 
     if (soluble_modes(2) > 1)   i_accumsolmass      = &
        get_q_index(standard_q_names%ACCUM_SOL_MASS,      'casim') 
     if (soluble_modes(2) > 0)   i_accumsolnumber    = &
        get_q_index(standard_q_names%ACCUM_SOL_NUMBER,    'casim') 
     if (soluble_modes(3) > 1)   i_coarsesolmass     = &
        get_q_index(standard_q_names%COARSE_SOL_MASS,     'casim') 
     if (soluble_modes(3) > 0)   i_coarsesolnumber   = &
        get_q_index(standard_q_names%COARSE_SOL_NUMBER,   'casim') 
     if (active_cloud(isol))     i_activesolliquid   = &
        get_q_index(standard_q_names%ACTIVE_SOL_LIQUID,   'casim') 
     if (active_rain(isol))      i_activesolrain     = &
        get_q_index(standard_q_names%ACTIVE_SOL_RAIN,     'casim') 
     if (insoluble_modes(2) > 1) i_coarsedustmass    = &
        get_q_index(standard_q_names%COARSE_DUST_MASS,    'casim') 
     if (insoluble_modes(2) > 0) i_coarsedustnumber  = &
        get_q_index(standard_q_names%COARSE_DUST_NUMBER,  'casim') 
     if (active_ice(iinsol))     i_activeinsolice    = &
        get_q_index(standard_q_names%ACTIVE_INSOL_ICE,    'casim') 
     if (active_ice(isol))       i_activesolice      = &
        get_q_index(standard_q_names%ACTIVE_SOL_ICE,      'casim') 
     if (active_cloud(iinsol))   i_activeinsolliquid = &
        get_q_index(standard_q_names%ACTIVE_INSOL_LIQUID, 'casim') 
     if (insoluble_modes(1) > 1) i_accuminsolmass    = &
        get_q_index(standard_q_names%ACCUM_INSOL_MASS,    'casim') 
     if (insoluble_modes(1) > 0) i_accuminsolnumber  = &
        get_q_index(standard_q_names%ACCUM_INSOL_NUMBER,  'casim') 
     if (active_number(isol))    i_activesolnumber   = &
        get_q_index(standard_q_names%ACTIVE_SOL_NUMBER,   'casim') 
     if (active_number(iinsol))  i_activeinsolnumber = &
          get_q_index(standard_q_names%ACTIVE_INSOL_NUMBER, 'casim')
  
     ! set logicals for the microphysics diagnostics: process rates
     casdiags % l_dth = .true.
     casdiags % l_dqv = .true.
     if ( nq_l>0 ) casdiags % l_dqc = .true.
     if ( nq_r>0 ) casdiags % l_dqr = .true.
     if ( l_pcond ) casdiags % l_pcond = .true.
     if ( l_psedl ) then
        casdiags % l_psedl = .true.
        casdiags % l_surface_rain = .true.
        casdiags % l_precip = .true.
     endif
     if ( l_praut ) casdiags % l_praut = .true.
     if ( l_pracw ) casdiags % l_pracw = .true.
     if ( l_prevp ) casdiags % l_prevp = .true.
     if ( l_psedr ) then
        casdiags % l_psedr = .true.
        casdiags % l_surface_rain = .true.
        casdiags % l_precip = .true.
     endif
     if (.not. l_warm) then
        if ( nq_i>0 ) casdiags % l_dqi = .true.
        if ( nq_s>0 ) casdiags % l_dqs = .true.
        if ( nq_g>0 ) casdiags % l_dqg = .true.
        if ( l_phomc ) casdiags % l_phomc = .true.
        if ( l_pinuc ) casdiags % l_pinuc = .true.
        if ( l_pidep ) casdiags % l_pidep = .true.
        if ( l_piacw ) casdiags % l_piacw = .true.
        if ( l_pisub ) casdiags % l_pisub = .true.
        if ( l_pimlt ) casdiags % l_pimlt = .true.
        if ( l_psedi ) then
           casdiags % l_psedi = .true.
           casdiags % l_surface_snow = .true.
           casdiags % l_precip = .true.
        endif
        if ( l_psmlt ) casdiags % l_psmlt = .true.
        if ( l_psaut ) casdiags % l_psaut = .true.
        if ( l_psaci ) casdiags % l_psaci = .true.
        if ( l_psacw ) casdiags % l_psacw = .true.
        if ( l_psacr ) casdiags % l_psacr = .true.
        if ( l_pssub ) casdiags % l_pssub = .true.
        if ( l_psdep ) casdiags % l_psdep = .true.
        if ( l_pseds ) then
           casdiags % l_pseds = .true.
           casdiags % l_surface_snow = .true.
           casdiags % l_precip = .true.
        endif
        if ( l_pgacw ) casdiags % l_pgacw = .true.
        if ( l_pgacs ) casdiags % l_pgacs = .true.
        if ( l_pgmlt ) casdiags % l_pgmlt = .true.
        if ( l_pgsub ) casdiags % l_pgsub = .true.
        if ( l_psedg ) then
           casdiags % l_psedg = .true.
           casdiags % l_surface_graup = .true.
           casdiags % l_precip = .true.
        endif
     endif
  
     ! allocate diagnostic space in casdiags depending on the logicals defined above
     CALL allocate_diagnostic_space(its, ite, jts, jte, kts, kte)
     ! this is no longer needed since can use cas_monc_dgs structure but keep for now
     allocate(surface_precip(y_size_local, x_size_local))
     ! allocate diagnostic space for MONC fields to export to IO server
     call allocate_casim_monc_dgs_space(current_state, casdiags)
  

Here is the call graph for this function:

Here is the caller graph for this function:

◆ read_configuration()

subroutine casim_mod::read_configuration ( type(model_state_type), intent(inout), target current_state )

private

Definition at line 856 of file casim.F90.

  
  
     Use mphys_parameters, only: p1, p2, p3, sp1, sp2, sp3   
  
     type(model_state_type), target, intent(inout) :: current_state
  
     integer :: ierr
  
     option          = options_get_integer(current_state%options_database, 'option')
     diag_mu_option  = options_get_integer(current_state%options_database, 'diag_mu_option')
     iopt_act        = options_get_integer(current_state%options_database, 'iopt_act')
     iopt_inuc       = options_get_integer(current_state%options_database, 'iopt_inuc')
     process_level   = options_get_integer(current_state%options_database, 'process_level')
     aerosol_option  = options_get_integer(current_state%options_database, 'aerosol_option')
     max_step_length = options_get_real(current_state%options_database, 'max_step_length')
     max_sed_length  = options_get_real(current_state%options_database, 'max_sed_length')
     p1              = options_get_real(current_state%options_database, 'p1')
     p2              = options_get_real(current_state%options_database, 'p2')
     p3              = options_get_real(current_state%options_database, 'p3')
     sp1             = options_get_real(current_state%options_database, 'sp1')
     sp2             = options_get_real(current_state%options_database, 'sp2')
     sp3             = options_get_real(current_state%options_database, 'sp3')
     max_mu          = options_get_real(current_state%options_database, 'max_mu')
     fix_mu          = options_get_real(current_state%options_database, 'fix_mu')
  
     l_aaut          = options_get_logical(current_state%options_database, 'l_aaut')
     l_aacc          = options_get_logical(current_state%options_database, 'l_aacc')
     l_aevp          = options_get_logical(current_state%options_database, 'l_aevp')
     l_ased          = options_get_logical(current_state%options_database, 'l_ased')
     l_warm          = options_get_logical(current_state%options_database, 'l_warm')
     l_inuc          = options_get_logical(current_state%options_database, 'l_inuc')
     l_iaut          = options_get_logical(current_state%options_database, 'l_iaut')
     l_idep          = options_get_logical(current_state%options_database, 'l_idep')
     l_iacw          = options_get_logical(current_state%options_database, 'l_iacw')
     l_active_inarg2000 = options_get_logical(current_state%options_database, 'l_active_inarg2000')
     l_separate_rain = options_get_logical(current_state%options_database, 'l_separate_rain')
     l_sg            = options_get_logical(current_state%options_database, 'l_sg')
     l_g             = options_get_logical(current_state%options_database, 'l_g')
     l_passive       = options_get_logical(current_state%options_database, 'l_passive')
     l_passive3m     = options_get_logical(current_state%options_database, 'l_passive3m')
     l_limit_psd     = options_get_logical(current_state%options_database, 'l_limit_psd')
     l_override_checks = options_get_logical(current_state%options_database, 'l_override_checks')
     l_raci_g        = options_get_logical(current_state%options_database, 'l_raci_g')
     l_onlycollect   = options_get_logical(current_state%options_database, 'l_onlycollect')
     l_abelshipway   = options_get_logical(current_state%options_database, 'l_abelshipway')
     l_cons          = options_get_logical(current_state%options_database, 'l_cons')
     l_rain          = options_get_logical(current_state%options_database, 'l_rain')
     l_sed_3mdiff    = options_get_logical(current_state%options_database, 'l_sed_3mdiff')
     l_sed_icecloud_as_1m = options_get_logical(current_state%options_database, 'l_sed_icecloud_as_1m')
     l_tidy_conserve_e = options_get_logical(current_state%options_database, 'l_tidy_conserve_E')
     l_tidy_conserve_q = options_get_logical(current_state%options_database, 'l_tidy_conserve_q')
  
     l_inhom_revp    = options_get_logical(current_state%options_database, 'l_inhom_revp')
     l_pcond         = options_get_logical(current_state%options_database, 'l_pcond')
     l_praut         = options_get_logical(current_state%options_database, 'l_praut')
     l_pracw         = options_get_logical(current_state%options_database, 'l_pracw')
     l_pracr         = options_get_logical(current_state%options_database, 'l_pracr')
     l_prevp         = options_get_logical(current_state%options_database, 'l_prevp')
     l_psedl         = options_get_logical(current_state%options_database, 'l_psedl')
     l_psedr         = options_get_logical(current_state%options_database, 'l_psedr')
     l_ptidy         = options_get_logical(current_state%options_database, 'l_ptidy')
     l_ptidy2        = options_get_logical(current_state%options_database, 'l_ptidy2')
     l_pinuc         = options_get_logical(current_state%options_database, 'l_pinuc')
     l_pidep         = options_get_logical(current_state%options_database, 'l_pidep')
     l_piacw         = options_get_logical(current_state%options_database, 'l_piacw')
     l_psaut         = options_get_logical(current_state%options_database, 'l_psaut')
     l_psdep         = options_get_logical(current_state%options_database, 'l_psdep')
     l_psacw         = options_get_logical(current_state%options_database, 'l_psacw')
     l_pgdep         = options_get_logical(current_state%options_database, 'l_pgdep')
     l_pseds         = options_get_logical(current_state%options_database, 'l_pseds')
     l_psedi         = options_get_logical(current_state%options_database, 'l_psedi')
     l_psedg         = options_get_logical(current_state%options_database, 'l_psedg')
     l_psaci         = options_get_logical(current_state%options_database, 'l_psaci')
     l_praci         = options_get_logical(current_state%options_database, 'l_praci')
     l_psacr         = options_get_logical(current_state%options_database, 'l_psacr')
     l_pgacr         = options_get_logical(current_state%options_database, 'l_pgacr')
     l_pgacw         = options_get_logical(current_state%options_database, 'l_pgacw')
     l_pgaci         = options_get_logical(current_state%options_database, 'l_pgaci')
     l_pgacs         = options_get_logical(current_state%options_database, 'l_pgacs')
     l_piagg         = options_get_logical(current_state%options_database, 'l_piagg')
     l_psagg         = options_get_logical(current_state%options_database, 'l_psagg')
     l_pgagg         = options_get_logical(current_state%options_database, 'l_pgagg')
     l_psbrk         = options_get_logical(current_state%options_database, 'l_psbrk')
     l_pgshd         = options_get_logical(current_state%options_database, 'l_pgshd')
     l_pihal         = options_get_logical(current_state%options_database, 'l_pihal')
     l_psmlt         = options_get_logical(current_state%options_database, 'l_psmlt')
     l_pgmlt         = options_get_logical(current_state%options_database, 'l_pgmlt')
     l_phomr         = options_get_logical(current_state%options_database, 'l_phomr')
     l_phomc         = options_get_logical(current_state%options_database, 'l_phomc')
     l_pssub         = options_get_logical(current_state%options_database, 'l_pssub')
     l_pgsub         = options_get_logical(current_state%options_database, 'l_pgsub')
     l_pisub         = options_get_logical(current_state%options_database, 'l_pisub')
     l_pimlt         = options_get_logical(current_state%options_database, 'l_pimlt')
  

Here is the caller graph for this function:

◆ timestep_callback()

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

private

Called for each column per timestep this will calculate the microphysical tendencies.

Parameters

current_state The current model state

Definition at line 502 of file casim.F90.

     type(model_state_type), target, intent(inout) :: current_state
  
     REAL(wp) :: dtwp
     INTEGER :: icol, jcol, iqx, target_x_index, target_y_index
  
     icol=current_state%column_local_x
     jcol=current_state%column_local_y
     target_y_index=jcol-current_state%local_grid%halo_size(y_index)
     target_x_index=icol-current_state%local_grid%halo_size(x_index)
  
     !if (current_state%first_timestep_column) then
     !   psedl_tot(:)= 0.0_DEFAULT_PRECISION
     !   pcond_tot(:)= 0.0_DEFAULT_PRECISION
     !endif
     
     if (current_state%halo_column .or. current_state%timestep < 2) return
  
     if (current_state%field_stepping == forward_stepping)then
       call log_master_log(log_error, 'Currently, CASIM assumes CENTERED_STEPPING')
       dtwp = current_state%dtm
     else
       dtwp = 2.0*current_state%dtm
     end if
  
     ! Initialize aerosol fields to zero...
     aitkensolmass = 0.0
     daitkensolmass = 0.0
     aitkensolnumber = 0.0
     daitkensolnumber = 0.0
     accumsolmass = 0.0
     daccumsolmass = 0.0
     accumsolnumber = 0.0
     daccumsolnumber = 0.0
     coarsesolmass = 0.0
     dcoarsesolmass = 0.0
     coarsesolnumber = 0.0
     dcoarsesolnumber = 0.0
     activesolliquid = 0.0
     dactivesolliquid = 0.0
     coarsedustmass = 0.0
     dcoarsedustmass = 0.0
     coarsedustnumber = 0.0
     dcoarsedustnumber = 0.0
     activeinsolice = 0.0
     dactiveinsolice = 0.0
     activesolice = 0.0
     dactivesolice = 0.0
     activeinsolliquid = 0.0
     dactiveinsolliquid = 0.0
     accuminsolmass = 0.0
     daccuminsolmass = 0.0
     accuminsolnumber = 0.0
     daccuminsolnumber = 0.0
     activesolnumber = 0.0
     dactivesolnumber = 0.0
     activeinsolnumber = 0.0
     dactiveinsolnumber = 0.0
     aitkensolbk = 0.0
     accumsolbk = 0.0 
     coarsesolbk = 0.0
  
     theta(:,1,1) = current_state%zth%data(:, jcol, icol) + current_state%global_grid%configuration%vertical%thref(:)
     dth(:,1,1) = current_state%sth%data(:, jcol, icol)
     exner(:,1,1) = current_state%global_grid%configuration%vertical%rprefrcp(:)
     pressure(:,1,1) = current_state%global_grid%configuration%vertical%prefn(:)
     z_centre(:,1,1) = current_state%global_grid%configuration%vertical%zn(:)
     dz(:,1,1) = current_state%global_grid%configuration%vertical%dz(:)
     z_half(:kte-1,1,1) = current_state%global_grid%configuration%vertical%z(:)
     rho(:,1,1) = current_state%global_grid%configuration%vertical%rhon(:)
     w(:,1,1) = current_state%zw%data(:, jcol, icol)
     tke(:,1,1) = 0.1 ! Test value
  
     iqx = iqv
     qv(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
     dqv(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
  
     ! Warm microphysical fields
     IF (nq_l > 0)then
       iqx = iql
       qc(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dqc(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_r > 0)then
       iqx = iqr
       qr(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dqr(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_l > 1)then
       iqx = inl
       nc(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dnc(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_r > 1)then
       iqx = inr
       nr(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dnr(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_r > 2)then
       iqx = i3mr
       m3r(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dm3r(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
  
     ! Ice microphysical fields
     IF (nq_i > 0)then
       iqx = iqi
       qi(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dqi(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_s > 0)then
       iqx = iqs
       qs(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dqs(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_g > 0)then
       iqx = iqg
       qg(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dqg(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_i > 1)then
       iqx = ini
       ni(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dni(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_s > 1)then
       iqx = ins
       ns(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dns(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_g > 1)then
       iqx = ing
       ng(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dng(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_s > 2)then
       iqx = i3ms
       m3s(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dm3s(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
     IF (nq_g > 2)then
       iqx = i3mg
       m3g(:,1,1) = current_state%zq(iqx)%data(:,jcol,icol)
       dm3g(:,1,1) = current_state%sq(iqx)%data(:,jcol,icol)
     end IF
  
     ! Aerosol fields
  
     if (i_aitkensolmass>0)     aitkensolmass(:,1,1)      = current_state%zq(i_aitkensolmass)%data(:,jcol,icol)
     if (i_aitkensolmass>0)    daitkensolmass(:,1,1)     = current_state%sq(i_aitkensolmass)%data(:,jcol,icol)
     if (i_aitkensolnumber>0)   aitkensolnumber(:,1,1)    = current_state%zq(i_aitkensolnumber)%data(:,jcol,icol)
     if (i_aitkensolnumber>0)  daitkensolnumber(:,1,1)   = current_state%sq(i_aitkensolnumber)%data(:,jcol,icol)
     if (i_accumsolmass>0)      accumsolmass(:,1,1)       = current_state%zq(i_accumsolmass)%data(:,jcol,icol)
     if (i_accumsolmass>0)     daccumsolmass(:,1,1)      = current_state%sq(i_accumsolmass)%data(:,jcol,icol)
     if (i_accumsolnumber>0)    accumsolnumber(:,1,1)     = current_state%zq(i_accumsolnumber)%data(:,jcol,icol)
     if (i_accumsolnumber>0)   daccumsolnumber(:,1,1)    = current_state%sq(i_accumsolnumber)%data(:,jcol,icol)
     if (i_coarsesolmass>0)     coarsesolmass(:,1,1)      = current_state%zq(i_coarsesolmass)%data(:,jcol,icol)
     if (i_coarsesolmass>0)    dcoarsesolmass(:,1,1)     = current_state%sq(i_coarsesolmass)%data(:,jcol,icol)
     if (i_coarsesolnumber>0)   coarsesolnumber(:,1,1)    = current_state%zq(i_coarsesolnumber)%data(:,jcol,icol)
     if (i_coarsesolnumber>0)  dcoarsesolnumber(:,1,1)   = current_state%sq(i_coarsesolnumber)%data(:,jcol,icol)
     if (i_activesolliquid>0)   activesolliquid(:,1,1)    = current_state%zq(i_activesolliquid)%data(:,jcol,icol)
     if (i_activesolliquid>0)  dactivesolliquid(:,1,1)   = current_state%sq(i_activesolliquid)%data(:,jcol,icol)
     if (i_coarsedustmass>0)    coarsedustmass(:,1,1)     = current_state%zq(i_coarsedustmass)%data(:,jcol,icol)
     if (i_coarsedustmass>0)   dcoarsedustmass(:,1,1)    = current_state%sq(i_coarsedustmass)%data(:,jcol,icol)
     if (i_coarsedustnumber>0)  coarsedustnumber(:,1,1)   = current_state%zq(i_coarsedustnumber)%data(:,jcol,icol)
     if (i_coarsedustnumber>0) dcoarsedustnumber(:,1,1)  = current_state%sq(i_coarsedustnumber)%data(:,jcol,icol)
     if (i_activeinsolice>0)    activeinsolice(:,1,1)     = current_state%zq(i_activeinsolice)%data(:,jcol,icol)
     if (i_activeinsolice>0)   dactiveinsolice(:,1,1)    = current_state%sq(i_activeinsolice)%data(:,jcol,icol)
     if (i_activesolice>0)      activesolice(:,1,1)       = current_state%zq(i_activesolice)%data(:,jcol,icol)
     if (i_activesolice>0)     dactivesolice(:,1,1)      = current_state%sq(i_activesolice)%data(:,jcol,icol)
     if (i_activeinsolliquid>0) activeinsolliquid(:,1,1)  = current_state%zq(i_activeinsolliquid)%data(:,jcol,icol)
     if (i_activeinsolliquid>0)dactiveinsolliquid(:,1,1) = current_state%sq(i_activeinsolliquid)%data(:,jcol,icol)
     if (i_accuminsolmass>0)    accuminsolmass(:,1,1)     = current_state%zq(i_accuminsolmass)%data(:,jcol,icol)
     if (i_accuminsolmass>0)   daccuminsolmass(:,1,1)    = current_state%sq(i_accuminsolmass)%data(:,jcol,icol)
     if (i_accuminsolnumber>0)  accuminsolnumber(:,1,1)   = current_state%zq(i_accuminsolnumber)%data(:,jcol,icol)
     if (i_accuminsolnumber>0) daccuminsolnumber(:,1,1)  = current_state%sq(i_accuminsolnumber)%data(:,jcol,icol)
     if (i_activesolnumber>0)   activesolnumber(:,1,1)    = current_state%zq(i_activesolnumber)%data(:,jcol,icol)
     if (i_activesolnumber>0)  dactivesolnumber(:,1,1)   = current_state%sq(i_activesolnumber)%data(:,jcol,icol)
     if (i_activeinsolnumber>0) activeinsolnumber(:,1,1)  = current_state%zq(i_activeinsolnumber)%data(:,jcol,icol)
     if (i_activeinsolnumber>0)dactiveinsolnumber(:,1,1) = current_state%sq(i_activeinsolnumber)%data(:,jcol,icol)
  
     CALL shipway_microphysics(                     &
                                 ! in
        its, ite,                                   &
        jts, jte,                                   &
        kts, kte,                                   &
        dtwp,                                       &
        qv, qc, qr,                                 &
        nc, nr, m3r,                                &
        qi, qs, qg,                                 &
        ni, ns, ng,                                 &
        m3s, m3g,                                   &
        theta,                                      &
        aitkensolmass, aitkensolnumber,             &
        accumsolmass, accumsolnumber,               &
        coarsesolmass, coarsesolnumber,             &
        activesolliquid,                            &
        activesolrain,                              &
        coarsedustmass, coarsedustnumber,           &
        activeinsolice,                             &
        activesolice,                               &
        activeinsolliquid,                          &
        accuminsolmass,                             &
        accuminsolnumber,                           &
        activesolnumber,                            &
        activeinsolnumber,                          &
        aitkensolbk,                                &
        accumsolbk,                                 & 
        coarsesolbk,                                &
        exner,                                      &
        pressure, rho,                              &
        w, tke,                                     &
        z_half, z_centre,                           &
        dz,                                         &
                                 ! in/out
        dqv, dqc, dqr, dnc, dnr, dm3r,              &
        dqi, dqs, dqg, dni, dns, dng, dm3s, dm3g,   &
        dth,                                        &
        daitkensolmass, daitkensolnumber,           &
        daccumsolmass, daccumsolnumber,             &
        dcoarsesolmass, dcoarsesolnumber,           &
        dactivesolliquid,                           &
        dactivesolrain,                             &
        dcoarsedustmass, dcoarsedustnumber,         &
        dactiveinsolice,                            &
        dactivesolice,                              &
        dactiveinsolliquid,                         &
        daccuminsolmass,                            &
        daccuminsolnumber,                          &
        dactivesolnumber,                           &
        dactiveinsolnumber,                         &
        ils, ile,                                   &
        jls, jle,                                   &
        kls, kle,                                   &
        l_tendency=.true.                           &
        )
  
     ! write back the tendencies
     current_state%sth%data(:,jcol,icol) = current_state%sth%data(:,jcol,icol) + dth(:,1,1)
  
     iqx = iqv
     current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dqv(:,1,1)
  
     ! Warm microphysical fields
     IF (nq_l > 0)then
       iqx = iql
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dqc(:,1,1)
     end IF
  
     IF (nq_r > 0)then
       iqx = iqr
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dqr(:,1,1)
     end IF
     IF (nq_l > 1)then
       iqx = inl
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dnc(:,1,1)
     end IF
     IF (nq_r > 1)then
       iqx = inr
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dnr(:,1,1)
     end IF
     IF (nq_r > 2)then
       iqx = i3mr
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dm3r(:,1,1)
     end IF
  
     ! Ice microphysical fields
     IF (nq_i > 0)then
       iqx = iqi
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dqi(:,1,1)
     end IF
     IF (nq_s > 0)then
       iqx = iqs
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dqs(:,1,1)
     end IF
     IF (nq_g > 0)then
       iqx = iqg
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dqg(:,1,1)
     end IF
     IF (nq_i > 1)then
       iqx = ini
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dni(:,1,1)
     end IF
     IF (nq_s > 1)then
       iqx = ins
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dns(:,1,1)
     end IF
     IF (nq_g > 1)then
       iqx = ing
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dng(:,1,1)
     end IF
     IF (nq_s > 2)then
       iqx = i3ms
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dm3s(:,1,1)
     end IF
     IF (nq_g > 2)then
       iqx = i3mg
       current_state%sq(iqx)%data(:,jcol,icol) = current_state%sq(iqx)%data(:,jcol,icol) + dm3g(:,1,1)
     end IF
  
     ! Aerosol fields
  
     if (i_aitkensolmass>0)     current_state%sq(i_aitkensolmass)%data(:,jcol,icol) &
        = current_state%sq(i_aitkensolmass)%data(:,jcol,icol) + daitkensolmass(:,1,1)
     if (i_aitkensolnumber>0)   current_state%sq(i_aitkensolnumber)%data(:,jcol,icol) &
        =  current_state%sq(i_aitkensolnumber)%data(:,jcol,icol) + daitkensolnumber(:,1,1)
     if (i_accumsolmass>0)      current_state%sq(i_accumsolmass)%data(:,jcol,icol) &
        =  current_state%sq(i_accumsolmass)%data(:,jcol,icol) + daccumsolmass(:,1,1)
     if (i_accumsolnumber>0)    current_state%sq(i_accumsolnumber)%data(:,jcol,icol) &
        =  current_state%sq(i_accumsolnumber)%data(:,jcol,icol) + daccumsolnumber(:,1,1)
     if (i_coarsesolmass>0)     current_state%sq(i_coarsesolmass)%data(:,jcol,icol) &
        =  current_state%sq(i_coarsesolmass)%data(:,jcol,icol) + dcoarsesolmass(:,1,1)
     if (i_coarsesolnumber>0)   current_state%sq(i_coarsesolnumber)%data(:,jcol,icol) &
        =  current_state%sq(i_coarsesolnumber)%data(:,jcol,icol) + dcoarsesolnumber(:,1,1)
     if (i_activesolliquid>0)   current_state%sq(i_activesolliquid)%data(:,jcol,icol) &
        =  current_state%sq(i_activesolliquid)%data(:,jcol,icol) + dactivesolliquid(:,1,1)
     if (i_coarsedustmass>0)    current_state%sq(i_coarsedustmass)%data(:,jcol,icol) &
        =  current_state%sq(i_coarsedustmass)%data(:,jcol,icol) + dcoarsedustmass(:,1,1)
     if (i_coarsedustnumber>0)  current_state%sq(i_coarsedustnumber)%data(:,jcol,icol) &
        =  current_state%sq(i_coarsedustnumber)%data(:,jcol,icol) + dcoarsedustnumber(:,1,1)
     if (i_activeinsolice>0)    current_state%sq(i_activeinsolice)%data(:,jcol,icol) &
        =  current_state%sq(i_activeinsolice)%data(:,jcol,icol) + dactiveinsolice(:,1,1)
     if (i_activesolice>0)      current_state%sq(i_activesolice)%data(:,jcol,icol) &
        =  current_state%sq(i_activesolice)%data(:,jcol,icol) + dactivesolice(:,1,1)
     if (i_activeinsolliquid>0) current_state%sq(i_activeinsolliquid)%data(:,jcol,icol) &
        =  current_state%sq(i_activeinsolliquid)%data(:,jcol,icol) + dactiveinsolliquid(:,1,1)
     if (i_accuminsolmass>0)    current_state%sq(i_accuminsolmass)%data(:,jcol,icol) &
        =  current_state%sq(i_accuminsolmass)%data(:,jcol,icol) + daccuminsolmass(:,1,1)
     if (i_accuminsolnumber>0)  current_state%sq(i_accuminsolnumber)%data(:,jcol,icol) &
        =  current_state%sq(i_accuminsolnumber)%data(:,jcol,icol) + daccuminsolnumber(:,1,1)
     if (i_activesolnumber>0)   current_state%sq(i_activesolnumber)%data(:,jcol,icol) &
        =  current_state%sq(i_activesolnumber)%data(:,jcol,icol) + dactivesolnumber(:,1,1)
     if (i_activeinsolnumber>0) current_state%sq(i_activeinsolnumber)%data(:,jcol,icol) &
        =  current_state%sq(i_activeinsolnumber)%data(:,jcol,icol) + dactiveinsolnumber(:,1,1)
  
     ! for total surface precipitation, sum the surface rain rate (cloud + rain which is precip_r)
     ! and surface
     ! snow rate (precip_s), which is the sum of ice, snow and graupel (See micromain.F90 in casim for
     ! calculation).
     if (l_warm) then
        surface_precip(target_y_index,target_x_index) = &
             casdiags % SurfaceRainR(1,1)
     else
        surface_precip(target_y_index,target_x_index) = &
             casdiags % SurfaceRainR(1,1) + casdiags % SurfaceSnowR(1,1)
     endif
     call populate_casim_monc_dg(current_state, casdiags)
  
     

Here is the caller graph for this function:

Variable Documentation

◆ accuminsolmass

real(wp), dimension(:,:,:), allocatable casim_mod::accuminsolmass

private

Definition at line 107 of file casim.F90.

107 REAL(wp), allocatable :: AccumInsolMass(:,:,:) ! Accum mode dust mass

◆ accuminsolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::accuminsolnumber

private

Definition at line 108 of file casim.F90.

108 REAL(wp), allocatable :: AccumInsolNumber(:,:,:) ! Accum mode dust number

◆ accumsolbk

real(wp), dimension(:,:,:), allocatable casim_mod::accumsolbk

private

Definition at line 116 of file casim.F90.

116 REAL(wp), allocatable :: AccumSolBk(:,:,:)

◆ accumsolmass

real(wp), dimension(:,:,:), allocatable casim_mod::accumsolmass

private

Definition at line 98 of file casim.F90.

98 REAL(wp), allocatable :: AccumSolMass(:,:,:), AccumSolNumber(:,:,:) ! Accumulation mode aerosol

◆ accumsolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::accumsolnumber

private

Definition at line 98 of file casim.F90.

◆ activeinsolice

real(wp), dimension(:,:,:), allocatable casim_mod::activeinsolice

private

Definition at line 104 of file casim.F90.

104 REAL(wp), allocatable :: ActiveInsolIce(:,:,:) ! Activeated dust

◆ activeinsolliquid

real(wp), dimension(:,:,:), allocatable casim_mod::activeinsolliquid

private

Definition at line 106 of file casim.F90.

106 REAL(wp), allocatable :: ActiveInsolLiquid(:,:,:) ! Activeated dust in cloud

◆ activeinsolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::activeinsolnumber

private

Definition at line 110 of file casim.F90.

110 REAL(wp), allocatable :: ActiveInsolNumber(:,:,:) ! Activated insoluble number (if we need a tracer)

◆ activesolice

real(wp), dimension(:,:,:), allocatable casim_mod::activesolice

private

Definition at line 105 of file casim.F90.

105 REAL(wp), allocatable :: ActiveSolIce(:,:,:) ! Activeated aerosol in ice

◆ activesolliquid

real(wp), dimension(:,:,:), allocatable casim_mod::activesolliquid

private

Definition at line 99 of file casim.F90.

99 REAL(wp), allocatable :: ActiveSolLiquid(:,:,:) ! Activated aerosol

◆ activesolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::activesolnumber

private

Definition at line 109 of file casim.F90.

109 REAL(wp), allocatable :: ActiveSolNumber(:,:,:) ! Activated soluble number (if we need a tracer)

◆ activesolrain

real(wp), dimension(:,:,:), allocatable casim_mod::activesolrain

private

Definition at line 102 of file casim.F90.

102 REAL(wp), allocatable :: ActiveSolRain(:,:,:) ! Activeated aerosol in rain

◆ aitkensolbk

real(wp), dimension(:,:,:), allocatable casim_mod::aitkensolbk

private

Definition at line 115 of file casim.F90.

115 REAL(wp), allocatable :: AitkenSolBk(:,:,:)

◆ aitkensolmass

real(wp), dimension(:,:,:), allocatable casim_mod::aitkensolmass

private

Definition at line 100 of file casim.F90.

100 REAL(wp), allocatable :: AitkenSolMass(:,:,:), AitkenSolNumber(:,:,:) ! Aitken mode aerosol

◆ aitkensolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::aitkensolnumber

private

Definition at line 100 of file casim.F90.

◆ coarsedustmass

real(wp), dimension(:,:,:), allocatable casim_mod::coarsedustmass

private

Definition at line 103 of file casim.F90.

103 REAL(wp), allocatable :: CoarseDustMass(:,:,:), CoarseDustNumber(:,:,:) ! Coarse Dust

◆ coarsedustnumber

real(wp), dimension(:,:,:), allocatable casim_mod::coarsedustnumber

private

Definition at line 103 of file casim.F90.

◆ coarsesolbk

real(wp), dimension(:,:,:), allocatable casim_mod::coarsesolbk

private

Definition at line 117 of file casim.F90.

117 REAL(wp), allocatable :: CoarseSolBk(:,:,:)

◆ coarsesolmass

real(wp), dimension(:,:,:), allocatable casim_mod::coarsesolmass

private

Definition at line 101 of file casim.F90.

101 REAL(wp), allocatable :: CoarseSolMass(:,:,:), CoarseSolNumber(:,:,:) ! Course mode aerosol

◆ coarsesolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::coarsesolnumber

private

Definition at line 101 of file casim.F90.

◆ daccuminsolmass

real(wp), dimension(:,:,:), allocatable casim_mod::daccuminsolmass

private

Definition at line 135 of file casim.F90.

135 REAL(wp), allocatable :: dAccumInsolMass(:,:,:) ! Accum mode dust mass

◆ daccuminsolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::daccuminsolnumber

private

Definition at line 136 of file casim.F90.

136 REAL(wp), allocatable :: dAccumInsolNumber(:,:,:) ! Accum mode dust number

◆ daccumsolmass

real(wp), dimension(:,:,:), allocatable casim_mod::daccumsolmass

private

Definition at line 126 of file casim.F90.

126 REAL(wp), allocatable :: dAccumSolMass(:,:,:), dAccumSolNumber(:,:,:) ! Accumulation mode aerosol

◆ daccumsolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::daccumsolnumber

private

Definition at line 126 of file casim.F90.

◆ dactiveinsolice

real(wp), dimension(:,:,:), allocatable casim_mod::dactiveinsolice

private

Definition at line 132 of file casim.F90.

132 REAL(wp), allocatable :: dActiveInsolIce(:,:,:) ! Activeated dust

◆ dactiveinsolliquid

real(wp), dimension(:,:,:), allocatable casim_mod::dactiveinsolliquid

private

Definition at line 134 of file casim.F90.

134 REAL(wp), allocatable :: dActiveInsolLiquid(:,:,:) ! Activeated dust in cloud

◆ dactiveinsolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::dactiveinsolnumber

private

Definition at line 138 of file casim.F90.

138 REAL(wp), allocatable :: dActiveInsolNumber(:,:,:) ! Activated insoluble number (if we need a tracer)

◆ dactivesolice

real(wp), dimension(:,:,:), allocatable casim_mod::dactivesolice

private

Definition at line 133 of file casim.F90.

133 REAL(wp), allocatable :: dActiveSolIce(:,:,:) ! Activeated aerosol in ice

◆ dactivesolliquid

real(wp), dimension(:,:,:), allocatable casim_mod::dactivesolliquid

private

Definition at line 127 of file casim.F90.

127 REAL(wp), allocatable :: dActiveSolLiquid(:,:,:) ! Activated aerosol

◆ dactivesolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::dactivesolnumber

private

Definition at line 137 of file casim.F90.

137 REAL(wp), allocatable :: dActiveSolNumber(:,:,:) ! Activated soluble number (if we need a tracer)

◆ dactivesolrain

real(wp), dimension(:,:,:), allocatable casim_mod::dactivesolrain

private

Definition at line 130 of file casim.F90.

130 REAL(wp), allocatable :: dActiveSolRain(:,:,:) ! Activeated aerosol in rain

◆ daitkensolmass

real(wp), dimension(:,:,:), allocatable casim_mod::daitkensolmass

private

Definition at line 128 of file casim.F90.

128 REAL(wp), allocatable :: dAitkenSolMass(:,:,:), dAitkenSolNumber(:,:,:) ! Aitken mode aerosol

◆ daitkensolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::daitkensolnumber

private

Definition at line 128 of file casim.F90.

◆ dcoarsedustmass

real(wp), dimension(:,:,:), allocatable casim_mod::dcoarsedustmass

private

Definition at line 131 of file casim.F90.

131 REAL(wp), allocatable :: dCoarseDustMass(:,:,:), dCoarseDustNumber(:,:,:) ! Dust

◆ dcoarsedustnumber

real(wp), dimension(:,:,:), allocatable casim_mod::dcoarsedustnumber

private

Definition at line 131 of file casim.F90.

◆ dcoarsesolmass

real(wp), dimension(:,:,:), allocatable casim_mod::dcoarsesolmass

private

Definition at line 129 of file casim.F90.

129 REAL(wp), allocatable :: dCoarseSolMass(:,:,:), dCoarseSolNumber(:,:,:) ! Course mode aerosol

◆ dcoarsesolnumber

real(wp), dimension(:,:,:), allocatable casim_mod::dcoarsesolnumber

private

Definition at line 129 of file casim.F90.

◆ dm3g

real(wp), dimension(:,:,:), allocatable casim_mod::dm3g

private

Definition at line 121 of file casim.F90.

◆ dm3r

real(wp), dimension(:,:,:), allocatable casim_mod::dm3r

private

Definition at line 121 of file casim.F90.

◆ dm3s

real(wp), dimension(:,:,:), allocatable casim_mod::dm3s

private

Definition at line 121 of file casim.F90.

◆ dnc

real(wp), dimension(:,:,:), allocatable casim_mod::dnc

private

Definition at line 121 of file casim.F90.

◆ dng

real(wp), dimension(:,:,:), allocatable casim_mod::dng

private

Definition at line 121 of file casim.F90.

◆ dni

real(wp), dimension(:,:,:), allocatable casim_mod::dni

private

Definition at line 121 of file casim.F90.

◆ dnr

real(wp), dimension(:,:,:), allocatable casim_mod::dnr

private

Definition at line 121 of file casim.F90.

◆ dns

real(wp), dimension(:,:,:), allocatable casim_mod::dns

private

Definition at line 121 of file casim.F90.

◆ dqc

real(wp), dimension(:,:,:), allocatable casim_mod::dqc

private

Definition at line 121 of file casim.F90.

◆ dqg

real(wp), dimension(:,:,:), allocatable casim_mod::dqg

private

Definition at line 121 of file casim.F90.

◆ dqi

real(wp), dimension(:,:,:), allocatable casim_mod::dqi

private

Definition at line 121 of file casim.F90.

◆ dqr

real(wp), dimension(:,:,:), allocatable casim_mod::dqr

private

Definition at line 121 of file casim.F90.

◆ dqs

real(wp), dimension(:,:,:), allocatable casim_mod::dqs

private

Definition at line 121 of file casim.F90.

◆ dqv

real(wp), dimension(:,:,:), allocatable casim_mod::dqv

private

Definition at line 121 of file casim.F90.

   REAL(wp), allocatable :: dqv(:,:,:), dth(:,:,:), dqc(:,:,:), dnc(:,:,:)              &
      , dqr(:,:,:), dnr(:,:,:), dm3r(:,:,:)                                  &
      , dqi(:,:,:), dni(:,:,:), dqs(:,:,:), dns(:,:,:), dm3s(:,:,:)      & 
      , dqg(:,:,:), dng(:,:,:), dm3g(:,:,:) 

◆ dth

real(wp), dimension(:,:,:), allocatable casim_mod::dth

private

Definition at line 121 of file casim.F90.

◆ dz

real(wp), dimension(:,:,:), allocatable casim_mod::dz

private

Definition at line 91 of file casim.F90.

◆ exner

real(wp), dimension(:,:,:), allocatable casim_mod::exner

private

Definition at line 91 of file casim.F90.

◆ i3mg

integer casim_mod::i3mg =0

private

Definition at line 147 of file casim.F90.

◆ i3mr

integer casim_mod::i3mr =0

private

Definition at line 147 of file casim.F90.

147 INTEGER :: i3mr=0, i3ms=0, i3mg=0

◆ i3ms

integer casim_mod::i3ms =0

private

Definition at line 147 of file casim.F90.

◆ i_accuminsolmass

integer casim_mod::i_accuminsolmass =0

private

Definition at line 148 of file casim.F90.

◆ i_accuminsolnumber

integer casim_mod::i_accuminsolnumber =0

private

Definition at line 148 of file casim.F90.

◆ i_accumsolmass

integer casim_mod::i_accumsolmass =0

private

Definition at line 148 of file casim.F90.

   INTEGER ::                &
      i_AccumSolMass=0,      &
      i_accumsolnumber=0,    &
      i_activesolliquid=0,   &
      i_aitkensolmass=0,     &
      i_aitkensolnumber=0,   &
      i_coarsesolmass=0,     &
      i_coarsesolnumber=0,   &
      i_activesolrain=0,     &
      i_coarsedustmass=0,    &
      i_coarsedustnumber=0,  &
      i_activeinsolice=0,    &
      i_activesolice=0,      &
      i_activeinsolliquid=0, &
      i_accuminsolmass=0,    &
      i_accuminsolnumber=0,  &
      i_activesolnumber=0,   &
      i_activeinsolnumber=0

◆ i_accumsolnumber

integer casim_mod::i_accumsolnumber =0

private

Definition at line 148 of file casim.F90.

◆ i_activeinsolice

integer casim_mod::i_activeinsolice =0

private

Definition at line 148 of file casim.F90.

◆ i_activeinsolliquid

integer casim_mod::i_activeinsolliquid =0

private

Definition at line 148 of file casim.F90.

◆ i_activeinsolnumber

integer casim_mod::i_activeinsolnumber =0

private

Definition at line 148 of file casim.F90.

◆ i_activesolice

integer casim_mod::i_activesolice =0

private

Definition at line 148 of file casim.F90.

◆ i_activesolliquid

integer casim_mod::i_activesolliquid =0

private

Definition at line 148 of file casim.F90.

◆ i_activesolnumber

integer casim_mod::i_activesolnumber =0

private

Definition at line 148 of file casim.F90.

◆ i_activesolrain

integer casim_mod::i_activesolrain =0

private

Definition at line 148 of file casim.F90.

◆ i_aitkensolmass

integer casim_mod::i_aitkensolmass =0

private

Definition at line 148 of file casim.F90.

◆ i_aitkensolnumber

integer casim_mod::i_aitkensolnumber =0

private

Definition at line 148 of file casim.F90.

◆ i_coarsedustmass

integer casim_mod::i_coarsedustmass =0

private

Definition at line 148 of file casim.F90.

◆ i_coarsedustnumber

integer casim_mod::i_coarsedustnumber =0

private

Definition at line 148 of file casim.F90.

◆ i_coarsesolmass

integer casim_mod::i_coarsesolmass =0

private

Definition at line 148 of file casim.F90.

◆ i_coarsesolnumber

integer casim_mod::i_coarsesolnumber =0

private

Definition at line 148 of file casim.F90.

◆ ile

integer casim_mod::ile

private

Definition at line 142 of file casim.F90.

◆ ils

integer casim_mod::ils

private

Definition at line 142 of file casim.F90.

142 INTEGER :: ils,ile, jls,jle, kls,kle, &

143 its,ite, jts,jte, kts,kte

◆ ing

integer casim_mod::ing =0

private

Definition at line 146 of file casim.F90.

◆ ini

integer casim_mod::ini =0

private

Definition at line 146 of file casim.F90.

◆ inl

integer casim_mod::inl =0

private

Definition at line 146 of file casim.F90.

146 INTEGER :: inl=0, inr=0, ini=0, ins=0, ing=0

◆ inr

integer casim_mod::inr =0

private

Definition at line 146 of file casim.F90.

◆ ins

integer casim_mod::ins =0

private

Definition at line 146 of file casim.F90.

◆ iqg

integer casim_mod::iqg =0

private

Definition at line 145 of file casim.F90.

◆ iqi

integer casim_mod::iqi =0

private

Definition at line 145 of file casim.F90.

◆ iql

integer casim_mod::iql =0

private

Definition at line 145 of file casim.F90.

◆ iqr

integer casim_mod::iqr =0

private

Definition at line 145 of file casim.F90.

◆ iqs

integer casim_mod::iqs =0

private

Definition at line 145 of file casim.F90.

◆ iqv

integer casim_mod::iqv =0

private

Definition at line 145 of file casim.F90.

145 INTEGER :: iqv=0, iql=0, iqr=0, iqi=0, iqs=0, iqg=0

◆ ite

integer casim_mod::ite

private

Definition at line 142 of file casim.F90.

◆ its

integer casim_mod::its

private

Definition at line 142 of file casim.F90.

◆ jle

integer casim_mod::jle

private

Definition at line 142 of file casim.F90.

◆ jls

integer casim_mod::jls

private

Definition at line 142 of file casim.F90.

◆ jte

integer casim_mod::jte

private

Definition at line 142 of file casim.F90.

◆ jts

integer casim_mod::jts

private

Definition at line 142 of file casim.F90.

◆ kle

integer casim_mod::kle

private

Definition at line 142 of file casim.F90.

◆ kls

integer casim_mod::kls

private

Definition at line 142 of file casim.F90.

◆ kte

integer casim_mod::kte

private

Definition at line 142 of file casim.F90.

◆ kts

integer casim_mod::kts

private

Definition at line 142 of file casim.F90.

◆ m3g

real(wp), dimension(:,:,:), allocatable casim_mod::m3g

private

Definition at line 91 of file casim.F90.

◆ m3r

real(wp), dimension(:,:,:), allocatable casim_mod::m3r

private

Definition at line 91 of file casim.F90.

◆ m3s

real(wp), dimension(:,:,:), allocatable casim_mod::m3s

private

Definition at line 91 of file casim.F90.

◆ nc

real(wp), dimension(:,:,:), allocatable casim_mod::nc

private

Definition at line 91 of file casim.F90.

◆ ng

real(wp), dimension(:,:,:), allocatable casim_mod::ng

private

Definition at line 91 of file casim.F90.

◆ ni

real(wp), dimension(:,:,:), allocatable casim_mod::ni

private

Definition at line 91 of file casim.F90.

◆ nr

real(wp), dimension(:,:,:), allocatable casim_mod::nr

private

Definition at line 91 of file casim.F90.

◆ ns

real(wp), dimension(:,:,:), allocatable casim_mod::ns

private

Definition at line 91 of file casim.F90.

◆ pcond_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pcond_tot

private

Definition at line 167 of file casim.F90.

◆ pgacs_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pgacs_tot

private

Definition at line 167 of file casim.F90.

◆ pgacw_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pgacw_tot

private

Definition at line 167 of file casim.F90.

◆ pgmlt_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pgmlt_tot

private

Definition at line 167 of file casim.F90.

◆ pgsub_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pgsub_tot

private

Definition at line 167 of file casim.F90.

◆ phomc_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::phomc_tot

private

Definition at line 167 of file casim.F90.

   real(kind=default_precision), dimension(:), allocatable ::     &
        phomc_tot, pinuc_tot, pidep_tot, psdep_tot, piacw_tot, psacw_tot, psacr_tot, pisub_tot,   &
        pssub_tot, pimlt_tot, psmlt_tot, psaut_tot, psaci_tot, praut_tot, pracw_tot, prevp_tot,   &
        pgacw_tot, pgacs_tot, pgmlt_tot, pgsub_tot, psedi_tot, pseds_tot, psedr_tot, psedg_tot,   &
        psedl_tot, pcond_tot

◆ piacw_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::piacw_tot

private

Definition at line 167 of file casim.F90.

◆ pidep_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pidep_tot

private

Definition at line 167 of file casim.F90.

◆ pimlt_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pimlt_tot

private

Definition at line 167 of file casim.F90.

◆ pinuc_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pinuc_tot

private

Definition at line 167 of file casim.F90.

◆ pisub_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pisub_tot

private

Definition at line 167 of file casim.F90.

◆ pracw_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pracw_tot

private

Definition at line 167 of file casim.F90.

◆ praut_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::praut_tot

private

Definition at line 167 of file casim.F90.

◆ pressure

real(wp), dimension(:,:,:), allocatable casim_mod::pressure

private

Definition at line 91 of file casim.F90.

◆ prevp_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::prevp_tot

private

Definition at line 167 of file casim.F90.

◆ psaci_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psaci_tot

private

Definition at line 167 of file casim.F90.

◆ psacr_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psacr_tot

private

Definition at line 167 of file casim.F90.

◆ psacw_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psacw_tot

private

Definition at line 167 of file casim.F90.

◆ psaut_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psaut_tot

private

Definition at line 167 of file casim.F90.

◆ psdep_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psdep_tot

private

Definition at line 167 of file casim.F90.

◆ psedg_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psedg_tot

private

Definition at line 167 of file casim.F90.

◆ psedi_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psedi_tot

private

Definition at line 167 of file casim.F90.

◆ psedl_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psedl_tot

private

Definition at line 167 of file casim.F90.

◆ psedr_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psedr_tot

private

Definition at line 167 of file casim.F90.

◆ pseds_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pseds_tot

private

Definition at line 167 of file casim.F90.

◆ psmlt_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::psmlt_tot

private

Definition at line 167 of file casim.F90.

◆ pssub_tot

real(kind=default_precision), dimension(:), allocatable casim_mod::pssub_tot

private

Definition at line 167 of file casim.F90.

◆ qc

real(wp), dimension(:,:,:), allocatable casim_mod::qc

private

Definition at line 91 of file casim.F90.

◆ qg

real(wp), dimension(:,:,:), allocatable casim_mod::qg

private

Definition at line 91 of file casim.F90.

◆ qi

real(wp), dimension(:,:,:), allocatable casim_mod::qi

private

Definition at line 91 of file casim.F90.

◆ qr

real(wp), dimension(:,:,:), allocatable casim_mod::qr

private

Definition at line 91 of file casim.F90.

◆ qs

real(wp), dimension(:,:,:), allocatable casim_mod::qs

private

Definition at line 91 of file casim.F90.

◆ qv

real(wp), dimension(:,:,:), allocatable casim_mod::qv

private

Definition at line 91 of file casim.F90.

◆ rho

real(wp), dimension(:,:,:), allocatable casim_mod::rho

private

Definition at line 91 of file casim.F90.

◆ surface_precip

real(wp), dimension(:,:), allocatable casim_mod::surface_precip

private

Definition at line 140 of file casim.F90.

140 REAL(wp), allocatable :: surface_precip(:,:)

◆ theta

real(wp), dimension(:,:,:), allocatable casim_mod::theta

private

Definition at line 91 of file casim.F90.

   REAL(wp), allocatable :: theta(:,:,:), pressure(:,:,:),  &
      z_half(:,:,:), z_centre(:,:,:), dz(:,:,:), qv(:,:,:),qc(:,:,:) &
      , nc(:,:,:), qr(:,:,:), nr(:,:,:), m3r(:,:,:),rho(:,:,:) &
      , exner(:,:,:), w(:,:,:), tke(:,:,:)                               &
      , qi(:,:,:), ni(:,:,:), qs(:,:,:), ns(:,:,:), m3s(:,:,:) &
      , qg(:,:,:), ng(:,:,:), m3g(:,:,:) 

◆ tke

real(wp), dimension(:,:,:), allocatable casim_mod::tke

private

Definition at line 91 of file casim.F90.

◆ w

real(wp), dimension(:,:,:), allocatable casim_mod::w

private

Definition at line 91 of file casim.F90.

◆ z_centre

real(wp), dimension(:,:,:), allocatable casim_mod::z_centre

private

Definition at line 91 of file casim.F90.

◆ z_half

real(wp), dimension(:,:,:), allocatable casim_mod::z_half

private

Definition at line 91 of file casim.F90.

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ casim_get_descriptor()

◆ field_information_retrieval_callback()

◆ field_value_retrieval_callback()

◆ initialisation_callback()

◆ read_configuration()

◆ timestep_callback()

Variable Documentation

◆ accuminsolmass

◆ accuminsolnumber

◆ accumsolbk

◆ accumsolmass

◆ accumsolnumber

◆ activeinsolice

◆ activeinsolliquid

◆ activeinsolnumber

◆ activesolice

◆ activesolliquid

◆ activesolnumber

◆ activesolrain

◆ aitkensolbk

◆ aitkensolmass

◆ aitkensolnumber

◆ coarsedustmass

◆ coarsedustnumber

◆ coarsesolbk

◆ coarsesolmass

◆ coarsesolnumber

◆ daccuminsolmass

◆ daccuminsolnumber

◆ daccumsolmass

◆ daccumsolnumber

◆ dactiveinsolice

◆ dactiveinsolliquid

◆ dactiveinsolnumber

◆ dactivesolice

◆ dactivesolliquid

◆ dactivesolnumber

◆ dactivesolrain

◆ daitkensolmass

◆ daitkensolnumber

◆ dcoarsedustmass

◆ dcoarsedustnumber

◆ dcoarsesolmass

◆ dcoarsesolnumber

◆ dm3g

◆ dm3r

◆ dm3s

◆ dnc

◆ dng

◆ dni

◆ dnr

◆ dns

◆ dqc

◆ dqg

◆ dqi

◆ dqr

◆ dqs

◆ dqv

◆ dth

◆ dz

◆ exner

◆ i3mg

◆ i3mr

◆ i3ms

◆ i_accuminsolmass

◆ i_accuminsolnumber

◆ i_accumsolmass

◆ i_accumsolnumber

◆ i_activeinsolice

◆ i_activeinsolliquid

◆ i_activeinsolnumber

◆ i_activesolice

◆ i_activesolliquid

◆ i_activesolnumber

◆ i_activesolrain

◆ i_aitkensolmass