lw_control_default_mod.F90

Go to the documentation of this file.

! *****************************COPYRIGHT*******************************
! (C) Crown copyright Met Office. All rights reserved.
! For further details please refer to the file COPYRIGHT.txt
! which you should have received as part of this distribution.
! *****************************COPYRIGHT*******************************
 
! Subroutine to set the default values of the control structure.
 
! Description:
!   Performs suitable default initializations of the controlling 
!   structure for LW radiation.
 
! Code Owner: Please refer to the UM file CodeOwners.txt
! This file belongs in section: Radiation Control
 
MODULE lw_control_default_mod
 
CONTAINS
 
SUBROUTINE lw_control_default(lw_control)
 
USE def_control, ONLY: strctrl
USE yomhook,     ONLY: lhook, dr_hook
USE parkind1,    ONLY: jprb, jpim
USE rad_pcf,     ONLY: ip_infra_red, ip_two_stream, ip_elsasser, &
                       ip_solver_homogen_direct, ip_trunc_azim_sym, &
                       ip_sph_direct, ip_sph_mode_rad,ip_max_rand, &
                       ip_homogeneous, ip_overlap_k_eqv_scl, &
                       ip_cloud_ice_water, ip_scatter_full
 
IMPLICIT NONE
 
TYPE (StrCtrl), INTENT(INOUT) :: lw_control
!   The block of controlling options for the code
 
INTEGER(KIND=jpim), PARAMETER :: zhook_in  = 0
INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1
REAL(KIND=jprb)               :: zhook_handle
 
 
IF (lhook) CALL dr_hook('LW_CONTROL_DEFAULT',zhook_in,zhook_handle)
 
! Spectral region and bands
lw_control%isolir     = ip_infra_red
lw_control%first_band = 1
 
! Physical processes
lw_control%l_gas          = .true.
lw_control%l_continuum    = .true.
lw_control%l_cloud        = .true.
lw_control%l_drop         = .true.
lw_control%l_ice          = .true.
lw_control%l_aerosol      = .false.
lw_control%l_aerosol_ccn  = .false.
 
! Properties of clouds
lw_control%l_local_cnv_partition = .false.
lw_control%l_global_cloud_top    = .true.
 
! Angular integration (including algorithmic options)
lw_control%n_channel              = 1
lw_control%i_angular_integration  = ip_two_stream
lw_control%i_2stream              = ip_elsasser
lw_control%i_solver_clear         = ip_solver_homogen_direct
lw_control%n_order_gauss          = 0
lw_control%i_truncation           = ip_trunc_azim_sym
lw_control%i_sph_algorithm        = ip_sph_direct
lw_control%n_order_phase_solar    = 1
lw_control%ls_global_trunc        = 9
lw_control%ms_min                 = 0
lw_control%ms_max                 = 0
lw_control%ls_brdf_trunc          = 0
lw_control%accuracy_adaptive      = 1.0e-04
lw_control%l_rescale              = .true.
lw_control%l_henyey_greenstein_pf = .true.
lw_control%i_sph_mode             = ip_sph_mode_rad
lw_control%i_solar_src            = 3
lw_control%l_ir_source_quad       = .false.
 
! Satellite data
lw_control%sat_hgt      = 0.0
lw_control%sat_lon      = 0.0
lw_control%sat_lat      = 0.0
lw_control%max_view_lon = 0.0
lw_control%min_view_lon = 0.0
lw_control%max_view_lat = 0.0
lw_control%min_view_lat = 0.0
 
! Originally in lw_rad_iput, defaults do not need to change for MONC
lw_control%l_solar_tail_flux      = .false.
lw_control%i_gas_overlap          = 2 ! set to 2 (random overlap, slow) 1 = water vapour only 
!                                     ! 5 = fast gas overlap
lw_control%i_gas                  = 2 ! forces gas to be vapour only when i_gas_overlap = 1 
lw_control%i_cloud_representation = ip_cloud_ice_water
lw_control%i_overlap              = ip_max_rand
lw_control%i_inhom                = ip_homogeneous
lw_control%l_n2o                  = .true.
lw_control%l_ch4                  = .true.
lw_control%l_o3                   = .true.
lw_control%l_h2o                  = .true.
lw_control%l_co2                  = .true.
lw_control%l_cfc11                = .true.
lw_control%l_cfc12                = .true.
lw_control%l_cfc113               = .true.
lw_control%l_cfc114               = .true.
lw_control%l_hcfc22               = .true.
lw_control%l_hfc125               = .true.
lw_control%l_hfc134a              = .true.
lw_control%l_co                   = .false.
lw_control%l_nh3                  = .false.
lw_control%l_tio                  = .false.
lw_control%l_vo                   = .false.
lw_control%l_h2                   = .false.
lw_control%l_he                   = .false.
lw_control%l_na                   = .false.
lw_control%l_k                    = .false.
lw_control%i_st_water             = 5
lw_control%i_cnv_water            = 5
lw_control%i_st_ice               = 8
lw_control%i_cnv_ice              = 8
lw_control%i_scatter_method       = ip_scatter_full
lw_control%l_microphysics         = .true.
 
IF (lhook) CALL dr_hook('LW_CONTROL_DEFAULT',zhook_out,zhook_handle)
 
END SUBROUTINE lw_control_default
 
END MODULE lw_control_default_mod