scale_atm_phy_rd_solarins_simple Module Reference

module FElib / Atmosphere / Physics radiation / Solar insolation / Simple gray-radiation scheme More...

Functions/Subroutines
subroutine, public	atm_phy_rd_solarins_simple_setup ()
	Setup the simplified solar insolation module.
subroutine, public	atm_phy_rd_solarins_simple_update ()
subroutine, public	atm_phy_rd_solarins_simple_get (solins, cossza, lat, np)
	Get solar insolation and cosine of the solar zenith angle for the given latitude array.

Detailed Description

module FElib / Atmosphere / Physics radiation / Solar insolation / Simple gray-radiation scheme

Description

A module for idealized solar insolation scheme

Author

Yuta Kawai, Team SCALE

Reference

• Frierson, D.M.W., Held, I.M., and Zurita-Gotor, P. 2006: A gray-radiation aquaplanet moist GCM. Part I: Static stability and eddy scale. J. Atmos. Sci., 63, 2548–2566.

NAMELIST

• PARAM_ATMOS_SOLARINS_SIMPLE

  name	type	default value	comment
  SOLARINS_TYPE	character(len=H_SHORT)	'ANNUAL_MEAN'
  CONST_FLUX	real(RP)	340.0_RP	Constant solar flux value [W/m2]
  ORBIT_REFERENCE_YEAR	integer	2000	Reference year for orbital calculations
  ANNUAL_MEAN_YEAR	integer	2000	Year for computing annual mean insolation
  ANNUAL_SAMPLES_PER_DAY	integer	1	Number of samples per day for annual mean calculation
  CACHE_ANNUAL_MEAN	logical	.true.	Flag to cache annual mean insolation
  FRIERSON2006_DELTA_S	real(RP)	1.4E0_RP	Parameter for latitudinal variation (Delta_s) in Frierson et al. (2006) insolation scheme

  
  

History Output

No history output

Function/Subroutine Documentation

atm_phy_rd_solarins_simple_setup()

subroutine, public scale_atm_phy_rd_solarins_simple::atm_phy_rd_solarins_simple_setup

Setup the simplified solar insolation module.

Definition at line 84 of file scale_atm_phy_rd_solarins_simple.F90.

    use scale_atmos_solarins, only: &
        atmos_solarins_setup    
    implicit none
    character(len=H_SHORT) :: SOLARINS_TYPE = 'ANNUAL_MEAN'
    namelist / param_atmos_solarins_simple / &
        solarins_type,                    &
        const_flux,                       &
        orbit_reference_year,             &
        annual_mean_year,                 &
        annual_samples_per_day,           &
        cache_annual_mean,                &
        frierson2006_delta_s
 
    integer :: ierr
    !------------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_RD_SOLARINS_SIMPLE_setup",*) 'Setup'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_solarins_simple,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_RD_SOLARINS_SIMPLE_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_RD_SOLARINS_SIMPLE_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_SOLARINS_SIMPLE. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_solarins_simple)
 
    select case (solarins_type)
    case ('CONST')
       solarins_type_id = solarins_simple_type_id_const
    case ('ANNUAL_MEAN')
       solarins_type_id = solarins_simple_type_id_annual_mean
    case ('FRIERSON2006')
       solarins_type_id = solarins_simple_type_id_frierson2006
    case default
       log_error("ATMOS_PHY_RD_SOLARINS_SIMPLE_setup",*) 'Not appropriate SOLARINS_TYPE. Check!'
       call prc_abort
    end select
 
    !- Initializes orbital parameters, the vernal-equinox reference date, and related internal state.
    !  Here, longitude and latitude are irrelevant because this module uses only the orbital state returned by ecliptic_longitude.
    call atmos_solarins_setup( 0.0_rp, 0.0_rp, orbit_reference_year )
 
    return

atm_phy_rd_solarins_simple_update()

subroutine, public scale_atm_phy_rd_solarins_simple::atm_phy_rd_solarins_simple_update

Definition at line 135 of file scale_atm_phy_rd_solarins_simple.F90.

atm_phy_rd_solarins_simple_get()

subroutine, public scale_atm_phy_rd_solarins_simple::atm_phy_rd_solarins_simple_get	(	real(rp), dimension(np), intent(out)	solins,
		real(rp), dimension(np), intent(out)	cossza,
		real(rp), dimension(np), intent(in)	lat,
		integer, intent(in)	np )

Get solar insolation and cosine of the solar zenith angle for the given latitude array.

Parameters

[in]	np	Number of points
[out]	solins	Solar insolation [W/m2] for each latitude point
[out]	cossza	Cosine of the solar zenith angle for each latitude point
[in]	lat	Array storing latitude values [radians] for each point

Definition at line 140 of file scale_atm_phy_rd_solarins_simple.F90.

    use scale_atmos_solarins, only: &
        atmos_solarins_constant
    implicit none
    integer, intent(in) :: Np           !< Number of points
    real(RP), intent(out) :: solins(Np) !< Solar insolation [W/m2] for each latitude point
    real(RP), intent(out) :: cosSZA(Np) !< Cosine of the solar zenith angle for each latitude point
    real(RP), intent(in) :: lat(Np)     !< Array storing latitude values [radians] for each point
    !------------------------------------------------------------------------------
 
    select case (solarins_type_id)
    case (solarins_simple_type_id_const)
      solins(:) = const_flux
      cossza(:) = undef
    case (solarins_simple_type_id_annual_mean)
      call annual_mean_insol( solins, & ! (out)
        lat, np )                       ! (in)
      cossza(:) = undef
    case (solarins_simple_type_id_frierson2006)
      solins(:) = 0.25_rp * atmos_solarins_constant &
                * ( 1.0_rp + 0.25_rp * frierson2006_delta_s * ( 1.0_rp - 3.0_rp * sin(lat(:))**2 ) )
      cossza(:) = undef
    end select
    return