simplecloud_mod Module Reference

A very simple saturation adjustment scheme without any microphysics. More...

Functions/Subroutines
type(component_descriptor_type) function, public	simplecloud_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 moisture fields. More...
subroutine	timestep_callback (current_state)
	Called for each column per timestep this will apply a forcing term to the aerosol fields. More...

Variables
integer	iqv
integer	iql
integer	k_cloudmax
real(kind=default_precision)	max_height_cloud

Detailed Description

A very simple saturation adjustment scheme without any microphysics.

Function/Subroutine Documentation

initialisation_callback()

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

private

The initialisation callback sets up the moisture fields.

Parameters

current_state

The current model state

Definition at line 41 of file simplecloud.F90.

    type(model_state_type), target, intent(inout) :: current_state

    integer :: k ! look counter

    if (is_component_enabled(current_state%options_database, "casim")) then
      call log_master_log(log_error, "Casim and Simplecloud are enabled, this does not work yet. Please disable one")
    end if 

    iqv=get_q_index(standard_q_names%VAPOUR, 'simplecloud')
    iql=get_q_index(standard_q_names%CLOUD_LIQUID_MASS, 'simplecloud')

    ! set buoyancy coefficient (value for vapour should be set
    ! elsewhere for a moist model
    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
    current_state%cq(iql) = -1.0 

    max_height_cloud=options_get_real(current_state%options_database, "max_height_cloud")
    do k=2, current_state%local_grid%size(z_index)-1
      if (current_state%global_grid%configuration%vertical%zn(k) > max_height_cloud) exit
    end do
    k_cloudmax=k-1

Here is the caller graph for this function:

simplecloud_get_descriptor()

type(component_descriptor_type) function, public simplecloud_mod::simplecloud_get_descriptor

(

)

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

Returns

The termination check component descriptor

Definition at line 32 of file simplecloud.F90.

    simplecloud_get_descriptor%name="simplecloud"
    simplecloud_get_descriptor%version=0.1
    simplecloud_get_descriptor%initialisation=>initialisation_callback
    simplecloud_get_descriptor%timestep=>timestep_callback

Here is the call graph for this function:

timestep_callback()

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

private

Called for each column per timestep this will apply a forcing term to the aerosol fields.

Parameters

current_state

The current model state

Definition at line 72 of file simplecloud.F90.

    type(model_state_type), target, intent(inout) :: current_state

    real(DEFAULT_PRECISION) :: tdegk   ! Temperature in Kelvin
    real(DEFAULT_PRECISION) :: pmb     ! Pressure in mb
    real(DEFAULT_PRECISION) :: exner   ! Exner pressure
    real(DEFAULT_PRECISION) :: one_over_exner ! Reciprocal of Exner pressure
    real(DEFAULT_PRECISION) :: qv,qc   ! Shorthand for vapour and cloud mass mixing ratio
    real(DEFAULT_PRECISION) :: qs      ! Saturation mixing ratio
    real(DEFAULT_PRECISION) :: dqsdt   ! Rate of change of qs with temperature
    real(DEFAULT_PRECISION) :: qsatfac ! Multiplicative factor
    real(DEFAULT_PRECISION) :: dmass   ! Mass transfer mixing ratio
    
    integer :: k          ! Loop counter
    integer :: icol, jcol ! Shorthand column indices

    real(DEFAULT_PRECISION) :: dtm  ! Local timestep variable

    if (current_state%halo_column) return

    dtm = current_state%dtm*2.0
    if (current_state%field_stepping == forward_stepping) dtm=current_state%dtm! Should this be revised to scalar_stepping

    icol=current_state%column_local_x
    jcol=current_state%column_local_y

    do k=2,k_cloudmax

      exner = current_state%global_grid%configuration%vertical%rprefrcp(k)
      one_over_exner = current_state%global_grid%configuration%vertical%prefrcp(k)
      pmb   = (current_state%global_grid%configuration%vertical%prefn(k)/100.) 

      if (current_state%field_stepping == forward_stepping) then  ! Should this be revised to scalar_stepping
        qv    = current_state%q(iqv)%data(k, jcol, icol) + current_state%sq(iqv)%data(k, jcol, icol)*dtm
        qc    = current_state%q(iql)%data(k, jcol, icol) + current_state%sq(iql)%data(k, jcol, icol)*dtm
        tdegk = (current_state%th%data(k, jcol, icol) + current_state%sth%data(k, jcol, icol)*dtm &
           + current_state%global_grid%configuration%vertical%thref(k))*exner
      else
        qv    = current_state%zq(iqv)%data(k, jcol, icol) + current_state%sq(iqv)%data(k, jcol, icol)*dtm
        qc    = current_state%zq(iql)%data(k, jcol, icol) + current_state%sq(iql)%data(k, jcol, icol)*dtm
        tdegk = (current_state%zth%data(k, jcol, icol) + current_state%sth%data(k, jcol, icol)*dtm &
           + current_state%global_grid%configuration%vertical%thref(k))*exner
      end if
      ! Calculate the cloud/vapour increments

      qs = qsaturation(tdegk, pmb)

      if (qv > qs .or. qc >0.0)then
        dqsdt = dqwsatdt(qs, tdegk)
        
        qsatfac = 1.0/(1.0 + rlvap_over_cp*dqsdt)
      
        dmass = max(-qc,(qv-qs)*qsatfac)/dtm
      
        current_state%sq(iqv)%data(k, jcol, icol) = current_state%sq(iqv)%data(k, jcol, icol) - dmass
        current_state%sq(iql)%data(k, jcol, icol) = current_state%sq(iql)%data(k, jcol, icol) + dmass

        current_state%sth%data(k, jcol, icol) = current_state%sth%data(k, jcol, icol) &
           + rlvap_over_cp*dmass*one_over_exner

      end if

    end do

    ! If there's any cloud above then evaporate it
    do k=k_cloudmax+1, current_state%local_grid%size(z_index)
      if (current_state%scalar_stepping == forward_stepping) then
        qv    = current_state%q(iqv)%data(k, jcol, icol) + current_state%sq(iqv)%data(k, jcol, icol)*dtm
        qc    = current_state%q(iql)%data(k, jcol, icol) + current_state%sq(iql)%data(k, jcol, icol)*dtm
        tdegk = (current_state%th%data(k, jcol, icol) + current_state%sth%data(k, jcol, icol)*dtm &
           + current_state%global_grid%configuration%vertical%thref(k))*exner
      else
        qv    = current_state%zq(iqv)%data(k, jcol, icol) + current_state%sq(iqv)%data(k, jcol, icol)*dtm
        qc    = current_state%zq(iql)%data(k, jcol, icol) + current_state%sq(iql)%data(k, jcol, icol)*dtm
        tdegk = (current_state%zth%data(k, jcol, icol) + current_state%sth%data(k, jcol, icol)*dtm &
           + current_state%global_grid%configuration%vertical%thref(k))*exner
      end if
      if (qc >0.0)then
        dmass = -qc/dtm
      
        current_state%sq(iqv)%data(k, jcol, icol) = current_state%sq(iqv)%data(k, jcol, icol) - dmass
        current_state%sq(iql)%data(k, jcol, icol) = current_state%sq(iql)%data(k, jcol, icol) + dmass

        current_state%sth%data(k, jcol, icol) = current_state%sth%data(k, jcol, icol) &
           + rlvap_over_cp*dmass*one_over_exner

      end if
    end do

Here is the caller graph for this function:

Variable Documentation

iql

integer simplecloud_mod::iql

private

Definition at line 21 of file simplecloud.F90.

iqv

integer simplecloud_mod::iqv

private

Definition at line 21 of file simplecloud.F90.

    integer :: iqv, iql

k_cloudmax

integer simplecloud_mod::k_cloudmax

private

Definition at line 23 of file simplecloud.F90.

    integer :: k_cloudmax ! max k index for height

max_height_cloud

real(kind=default_precision) simplecloud_mod::max_height_cloud

private

Definition at line 24 of file simplecloud.F90.

    real(kind=DEFAULT_PRECISION) :: max_height_cloud