set_bound.F90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	set_bound (control, isolir, atm, dimen, spectrum, bound, socrates_derived_fields)

Function/Subroutine Documentation

set_bound()

subroutine set_bound	(	type(strctrl), intent(in)	control,
		integer, intent(in)	isolir,
		type(stratm), intent(in)	atm,
		type(strdim), intent(in)	dimen,
		type (strspecdata), intent(in)	spectrum,
		type(strbound), intent(out)	bound,
		type (str_socrates_derived_fields), intent(in)	socrates_derived_fields
	)

Definition at line 10 of file set_bound.F90.

 
USE rad_pcf
USE def_control, ONLY: strctrl
USE def_atm,     ONLY: stratm 
USE def_spectrum, ONLY: strspecdata
USE def_dimen,    ONLY: strdim
USE def_bound,    ONLY: strbound, allocate_bound
use def_socrates_derived_fields, only: str_socrates_derived_fields
 
IMPLICIT NONE
 
! Control options:
TYPE(StrCtrl),      INTENT(IN)    :: control
 
!Spectral region
INTEGER,            INTENT(IN)    :: isolir
 
! Atmospheric properties:
TYPE(StrAtm),       INTENT(IN) :: atm
 
! Dimensions:
TYPE(StrDim),       INTENT(IN)  :: dimen
 
! Spectral data:
TYPE (StrSpecData), INTENT(IN)  :: spectrum
 
! Boundary conditions:
TYPE(StrBound),     INTENT(OUT) :: bound
 
! solar variables from MONC solar_position_angle_mod and timestep_callback
type (str_socrates_derived_fields), intent(in) :: socrates_derived_fields
 
call allocate_bound(bound, dimen, spectrum)
 
!     Set the surface basis functions for a Lambertian surface.
bound%n_brdf_basis_fnc = 1
!     By setting F_{1,0,0,0} equal to 4 we can set rho_alb equal to
!     the diffuse albedo
bound%f_brdf(1,0,0,0)  = 4.0 
if ( isolir == ip_solar ) then
!     Need to set this initially using input from config
  bound%rho_alb(1:atm%n_profile, ip_surf_alb_diff, 1:spectrum%basic%n_band) = &
       socrates_derived_fields%albedoin1
  bound%rho_alb(1:atm%n_profile, ip_surf_alb_dir,  1:spectrum%basic%n_band) = & 
       socrates_derived_fields%albedoin1
else if ( isolir == ip_infra_red ) then
  bound%rho_alb(1:atm%n_profile, ip_surf_alb_diff, 1:spectrum%basic%n_band) = &
       0.0
  bound%rho_alb(1:atm%n_profile, ip_surf_alb_dir,  1:spectrum%basic%n_band) = & 
       0.0
end if
!     Holds the secant of the zenith angle for the fluxes
bound%zen_0(1:atm%n_profile) = socrates_derived_fields%sec_out
! sol_const = default_solar_const * scs (calced in timestep_callback)
! fractional_lit = lit in UM (is 0 or 1)
bound%solar_irrad(:) = socrates_derived_fields%sol_const * socrates_derived_fields%fraction_lit
!     pass from the MONC surface temperature
bound%t_ground(1:atm%n_profile) = socrates_derived_fields%srf_temperature ! Needs sorting
 

Here is the caller graph for this function: