scale_atmos_solarins Module Reference

module atmosphere / SOLARINS More...

Functions/Subroutines
subroutine, public	atmos_solarins_setup (basepoint_lon, basepoint_lat, iyear)
	setup solar incidence module More...
subroutine, public	atmos_solarins_orbit (iyear)
	setup solar incidence module More...
subroutine, public	atmos_solarins_ecliptic_longitude (Re_factor, sinDEC, cosDEC, hourangle, now_date, offset_year, lambda_out)
	calc factor of Earths solar insolation More...
subroutine	atmos_solarins_insolation_0d (real_lon, real_lat, now_date, offset_year, solins, cosSZA, Re_factor_out)
	calc factor of Earths solar insolation More...
subroutine	atmos_solarins_insolation_2d (IA, IS, IE, JA, JS, JE, real_lon, real_lat, now_date, offset_year, solins, cosSZA)
	calc factor of Earths solar insolation More...

Variables
real(rp), public	atmos_solarins_constant = 1360.250117_RP
logical, public	atmos_solarins_set_ve = .false.
logical, public	atmos_solarins_set_ideal = .false.
real(rp), public	atmos_solarins_obliquity = 23.44_RP
real(rp), public	atmos_solarins_eccentricity = 0.0167_RP
real(rp), public	atmos_solarins_perihelion_lon = 282.94719_RP
integer, dimension(6), public	atmos_solarins_ve_date
real(rp), public	atmos_solarins_diurnal_sec
real(rp), public	atmos_solarins_annual_sec
logical, public	atmos_solarins_fixedlatlon = .false.
real(rp), public	atmos_solarins_lon
real(rp), public	atmos_solarins_lat
logical, public	atmos_solarins_fixeddate = .false.
integer, dimension(6), public	atmos_solarins_date

Detailed Description

module atmosphere / SOLARINS

Description

calculate solar insolation based on berger(1978), berger et al.(1993) The original algorithm is valid to +/- 1,000,000 years from AD1950

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_SOLARINS

  name	type	default value	comment
  ATMOS_SOLARINS_CONSTANT	real(RP)	1360.250117_RP	Solar constant [W/m2]
  ATMOS_SOLARINS_SET_VE	logical	.false.	Set vernal equinox condition?
  ATMOS_SOLARINS_SET_IDEAL	logical	.false.	Set obliquity and eccentricity?
  ATMOS_SOLARINS_OBLIQUITY	real(RP)	23.44_RP	Obliquity [deg] cf. https://nssdc.gsfc.nasa.gov/planetary/factsheet/earthfact.html
  ATMOS_SOLARINS_ECCENTRICITY	real(RP)	0.0167_RP	Eccentricity
  ATMOS_SOLARINS_FIXEDLATLON	logical	.false.	Latitude/Longitude is fixed?
  ATMOS_SOLARINS_LON	real(RP)		Longitude for radiation [rad]
  ATMOS_SOLARINS_LAT	real(RP)		Latitude for radiation [rad]
  ATMOS_SOLARINS_FIXEDDATE	logical	.false.	Date is fixed?
  ATMOS_SOLARINS_DATE	integer, dimension(6)		Date for radiation [Y,M,D,H,M,S]
  ATMOS_SOLARINS_PERIHELION_LON	real(RP)	282.94719_RP	Longitude of perihelion [deg].
  ATMOS_SOLARINS_VE_DATE	integer, dimension(6)		Date of first vernal equinox
  ATMOS_SOLARINS_ANNUAL_SEC	real(RP)		Seconds of the annual period [sec]
  ATMOS_SOLARINS_DIURNAL_SEC	real(RP)		Seconds of the diurnal period [sec]
  DEBUG	logical	.false.

  
  

History Output

No history output

Function/Subroutine Documentation

atmos_solarins_setup()

subroutine, public scale_atmos_solarins::atmos_solarins_setup	(	real(rp), intent(in)	basepoint_lon,
		real(rp), intent(in)	basepoint_lat,
		integer, intent(in)	iyear
	)

setup solar incidence module

Definition at line 570 of file scale_atmos_solarins.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       const_d2r,          &
       const_undef
    use scale_calendar, only: &
       calendar_doi,          &
       calendar_hour,         &
       calendar_min,          &
       calendar_sec,          &
       calendar_daysec2date,  &
       calendar_date2daysec,  &
       calendar_date2char
    use scale_time, only: &
       time_offset_year,  &
       time_nowday,       &
       time_nowsec
    implicit none
    integer,  intent(in) :: iyear ! year at setup
    real(RP), intent(in) :: basepoint_lon
    real(RP), intent(in) :: basepoint_lat
 
    namelist / param_atmos_solarins / &
       atmos_solarins_constant,       &
       atmos_solarins_set_ve,         &
       atmos_solarins_set_ideal,      &
       atmos_solarins_obliquity,      &
       atmos_solarins_eccentricity,   &
       atmos_solarins_fixedlatlon,    &
       atmos_solarins_lon,            &
       atmos_solarins_lat,            &
       atmos_solarins_fixeddate,      &
       atmos_solarins_date,           &
       atmos_solarins_perihelion_lon, &
       atmos_solarins_ve_date,        &
       atmos_solarins_annual_sec,     &
       atmos_solarins_diurnal_sec,    &
       debug
 
    real(RP) :: dyear ! delta t [year]
    integer  :: year  ! used year
    integer  :: d
    integer  :: date(6)
    real(RP) :: Re_factor, sinDEC, cosDEC, hourangle
    real(DP) :: lambda, subsec
    character(len=27) :: datechar
 
    integer  :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_SOLARINS_setup",*) 'Setup'
 
    atmos_solarins_lon         = basepoint_lon
    atmos_solarins_lat         = basepoint_lat
    atmos_solarins_date(:)     = -1
    atmos_solarins_diurnal_sec = const_undef
    atmos_solarins_annual_sec  = const_undef
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_solarins,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_SOLARINS_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_SOLARINS_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_SOLARINS. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_solarins)
 
    if ( atmos_solarins_set_ve ) then
       if ( .NOT. atmos_solarins_set_ideal ) then
          atmos_solarins_set_ideal    = .true.
          atmos_solarins_obliquity    = 0.0_rp
          atmos_solarins_eccentricity = 0.0_rp
          atmos_solarins_perihelion_lon  = 0.0_rp
       endif
       if ( .NOT. atmos_solarins_fixedlatlon ) then
          atmos_solarins_fixedlatlon  = .true.
          atmos_solarins_lon          = 0.0_rp
          atmos_solarins_lat          = 0.0_rp
       endif
       if ( .NOT. atmos_solarins_fixeddate ) then
          atmos_solarins_fixeddate    = .true.
          atmos_solarins_date(1)      = ve_date(1)
          atmos_solarins_date(2)      = ve_date(2)
          atmos_solarins_date(3)      = ve_date(3)
          atmos_solarins_date(4)      = 12
          atmos_solarins_date(5)      = 0
          atmos_solarins_date(6)      = 0
       endif
    endif
 
    arcsec2d = 1.0_rp / (60.0_rp*60.0_rp)
    arcsec2r = arcsec2d * const_d2r
 
    year = iyear
 
    if ( atmos_solarins_fixeddate ) then
       if( atmos_solarins_date(1) >= 0 ) year = atmos_solarins_date(1)
    endif
 
    if ( atmos_solarins_set_ideal ) then
       ve_date = atmos_solarins_ve_date
       calyearday = calendar_doi
       caldaysec  = calendar_hour*calendar_min*calendar_sec
 
       if (atmos_solarins_annual_sec == const_undef ) then
         solyearsec = calyearday*caldaysec
         atmos_solarins_annual_sec = solyearsec
       else
         solyearsec = atmos_solarins_annual_sec
         year_calender_synced = .false.
       endif
 
       if (atmos_solarins_diurnal_sec == const_undef ) then
         soldaysec  = caldaysec
         atmos_solarins_diurnal_sec = soldaysec
       else
         soldaysec  = atmos_solarins_diurnal_sec
         day_calender_synced = .false.
       endif
 
    else
       calyearday = 365.2425_dp
       caldaysec  = calendar_hour*calendar_min*calendar_sec
       solyearsec = calyearday*caldaysec
       soldaysec  = caldaysec
    endif
 
    call atmos_solarins_orbit( year )
 
    dyear = real( year-year_ref, kind=rp )
 
    !----- report data -----
    log_newline
    log_info("ATMOS_SOLARINS_setup",'(1x,A)')       'Solar insolation parameters '
    log_info("ATMOS_SOLARINS_setup",'(1x,A,I7)')    'Reference year                        : ', year_ref
    log_info("ATMOS_SOLARINS_setup",'(1x,A,I7)')    'Current   year                        : ', year
    log_info("ATMOS_SOLARINS_setup",'(1x,A,I7)')    'Difference from ref.                  : ', int(dyear)
    log_newline
    log_info("ATMOS_SOLARINS_setup",'(1x,A,F12.7)') 'Solar constant                 [W/m2] : ', atmos_solarins_constant
    log_info("ATMOS_SOLARINS_setup",'(1x,A,F12.7)') 'Obliquity                       [deg] : ', obliquity / const_d2r
    log_info("ATMOS_SOLARINS_setup",'(1x,A,F12.7)') 'Eccentricity                          : ', e
    log_info("ATMOS_SOLARINS_setup",'(1x,A,F12.7)') 'Longitude of perihelion         [deg] : ', omega     / const_d2r
!    LOG_INFO("ATMOS_SOLARINS_setup",'(1x,A,F12.7)') 'Longitude at the vernal equinox associated with circle orbit [deg] : ', lambda_m0 / CONST_D2R
    log_newline
    if ( year_calender_synced ) then
       log_info("ATMOS_SOLARINS_setup",'(1x,A      )') 'Solar annual cycle is synced with calendar'
    else
       log_info("ATMOS_SOLARINS_setup",'(1x,A      )') 'Solar annual cycle is not synced with calendar'
    endif
    log_info("ATMOS_SOLARINS_setup",'(1x,A,F14.4)') 'Solar annual cycle              [sec] : ', solyearsec
    log_info("ATMOS_SOLARINS_setup",'(1x,A,F14.4)') 'Solar annual cycle              [day] : ', solyearsec/caldaysec
    if ( day_calender_synced ) then
       log_info("ATMOS_SOLARINS_setup",'(1x,A      )') 'Solar diurnal cycle is synced with calendar'
    else
       log_info("ATMOS_SOLARINS_setup",'(1x,A      )') 'Solar diurnal cycle is not synced with calendar'
    endif
    log_info("ATMOS_SOLARINS_setup",'(1x,A,F14.4)') 'Solar diurnal cycle             [sec] : ', soldaysec
    log_newline
    log_info("ATMOS_SOLARINS_setup",*)              'Latitude/Longitude is fixed?          : ', atmos_solarins_fixedlatlon
    if ( atmos_solarins_fixedlatlon ) then
       log_info("ATMOS_SOLARINS_setup",*)           'Longitude                       [deg] : ', atmos_solarins_lon
       log_info("ATMOS_SOLARINS_setup",*)           'Latitude                        [deg] : ', atmos_solarins_lat
       atmos_solarins_lon = atmos_solarins_lon * const_d2r
       atmos_solarins_lat = atmos_solarins_lat * const_d2r
    endif
 
    log_info("ATMOS_SOLARINS_setup",*)              'Date is fixed?                        : ', atmos_solarins_fixeddate
    if ( atmos_solarins_fixeddate ) then
       log_info("ATMOS_SOLARINS_setup",*)           'Date                                  : ', atmos_solarins_date
    endif
 
    log_newline
    log_info("ATMOS_SOLARINS_setup",*)  'Orbit info for next 1 solar year'
    if (io_l) write(io_fid_log, '(A)') 'DATE                         SOLAR_CONST, Scaled dist from star, Longitude'
    do d = time_nowday, time_nowday + int(solyearsec/caldaysec)
       call calendar_daysec2date(date,subsec,d,time_nowsec,time_offset_year)
       call atmos_solarins_ecliptic_longitude(  &
          re_factor,                            & ![OUT]
          sindec,                               & ![OUT]
          cosdec,                               & ![OUT]
          hourangle,                            & ![OUT]
          date,                                 & ![IN]
          time_offset_year,                     & ![IN]
          lambda_out = lambda                   & ![optional OUT]
       )
       call calendar_date2char(datechar, date, subsec)
       if (lambda < 0.0_rp) then
          lambda = lambda/const_d2r + 360.0_rp
       else
          lambda = lambda/const_d2r
       endif
       if (io_l) then
          write(io_fid_log,'(2A,F11.4,F23.6,F11.4)') datechar, "  ", re_factor*atmos_solarins_constant, 1.0_dp/sqrt(re_factor), lambda
       endif
    enddo
 
    return

References atmos_solarins_annual_sec, atmos_solarins_constant, atmos_solarins_date, atmos_solarins_diurnal_sec, atmos_solarins_eccentricity, atmos_solarins_ecliptic_longitude(), atmos_solarins_fixeddate, atmos_solarins_fixedlatlon, atmos_solarins_lat, atmos_solarins_lon, atmos_solarins_obliquity, atmos_solarins_orbit(), atmos_solarins_perihelion_lon, atmos_solarins_set_ideal, atmos_solarins_set_ve, atmos_solarins_ve_date, scale_calendar::calendar_date2char(), scale_calendar::calendar_date2daysec(), scale_calendar::calendar_daysec2date(), scale_calendar::calendar_doi, scale_calendar::calendar_hour, scale_calendar::calendar_min, scale_calendar::calendar_sec, scale_const::const_d2r, scale_const::const_undef, scale_io::io_fid_conf, scale_io::io_fid_log, scale_io::io_l, scale_prc::prc_abort(), scale_precision::rp, scale_time::time_nowday, scale_time::time_nowsec, and scale_time::time_offset_year.

Referenced by mod_atmos_driver::atmos_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_solarins_orbit()

subroutine, public scale_atmos_solarins::atmos_solarins_orbit

(

integer, intent(in)

iyear

)

setup solar incidence module

Definition at line 777 of file scale_atmos_solarins.F90.

    use scale_const, only: &
       pi  => const_pi, &
       d2r => const_d2r,&
       const_undef
    implicit none
 
    integer, intent(in) :: iyear ! year at setup
 
    real(RP) :: dyear       ! delta t [year]
    real(RP) :: EsinOMG     ! e * sin(w) : ecliptic parameter
    real(RP) :: EcosOMG     ! e * cos(w) : ecliptic parameter
    real(RP) :: Perih_pi    ! pi         : longitude of fixed perihelion [rad]
    real(RP) :: Perih_psi   ! psi        : general precession [degree]
    real(RP) :: Perih_omega ! omega_bar  : longitude of perihelion [degree]
 
    real(RP) :: temp, beta, EcosOMG_mod
    integer  :: i
    !---------------------------------------------------------------------------
 
    if ( .not. atmos_solarins_set_ideal ) then ! real Earth orbit
       ! time from reference year(1950.0 AD)
       dyear = real( iyear-year_ref, kind=rp )
 
       ! obliquity
       temp = 0.0_rp
       do i = 1, nobliq
          temp = temp + obliq_amp(i)*arcsec2d * cos( obliq_rate(i)*arcsec2r*dyear + obliq_phase(i)*d2r )
       enddo
       obliquity = ( obliquity_ref + temp ) * d2r
 
       ! eccentricity
       esinomg = 0.0_rp
       ecosomg = 0.0_rp
       do i = 1, neclip
          esinomg = esinomg + eclip_amp(i) * sin( eclip_rate(i)*arcsec2r*dyear + eclip_phase(i)*d2r )
          ecosomg = ecosomg + eclip_amp(i) * cos( eclip_rate(i)*arcsec2r*dyear + eclip_phase(i)*d2r )
       enddo
       e = sqrt( esinomg*esinomg + ecosomg*ecosomg )
 
       ! longitude of fixed perihelion
       if ( ecosomg == 0.0_rp ) then
          ecosomg_mod  = 1.e-30_rp
       else
          ecosomg_mod  = ecosomg
       endif
       perih_pi = atan(esinomg/ecosomg_mod) + pi * ( 0.5_rp - sign(0.5_rp, ecosomg) )
 
      ! if( abs(EcosOMG) < 1.E-8_RP ) then
      !    if    ( EsinOMG == 0.0_RP ) then
      !       Perih_pi = 0.0_RP
      !    elseif( EsinOMG <  0.0_RP ) then
      !       Perih_pi = 1.5_RP * PI
      !    elseif( EsinOMG >  0.0_RP ) then
      !       Perih_pi = 0.5_RP * PI
      !    endif
 
      ! elseif( EcosOMG < 0.0_RP ) then
 
      !    Perih_pi = atan(EsinOMG/EcosOMG) + PI
 
      ! elseif( EcosOMG > 0.0_RP ) then
 
      !    if ( EsinOMG <  0.0_RP ) then
      !       Perih_pi = atan(EsinOMG/EcosOMG) + 2.0_RP * PI
      !    else
      !       Perih_pi = atan(EsinOMG/EcosOMG)
      !    endif
 
      ! endif
        perih_pi = perih_pi / d2r ! [rad]->[degree]
 
       ! general precession
       temp = 0.0_rp
       do i = 1, nprece
          temp = temp + prece_amp(i)*arcsec2d * sin( prece_rate(i)*arcsec2r*dyear + prece_phase(i)*d2r )
       enddo
       perih_psi = psi_bar * arcsec2d * dyear + zeta + temp
 
       ! longitude of perihelion
       perih_omega = perih_pi + perih_psi
 
       do i = 1, 1000
          if    ( perih_omega + 180.0_rp <    0.0_rp ) then
             perih_omega = perih_omega + 360.0_rp
          elseif( perih_omega + 180.0_rp >= 360.0_rp ) then
             perih_omega = perih_omega - 360.0_rp
          else
             exit
          endif
       enddo
 
       ! longitude of perihelion (see berger et al.(1993))
       omega = ( perih_omega + 180.0_rp ) * d2r
 
       atmos_solarins_eccentricity    = e
       atmos_solarins_obliquity       = obliquity / d2r
       atmos_solarins_perihelion_lon  = omega / d2r
 
    else ! ideal orbit
 
      obliquity = atmos_solarins_obliquity * d2r
      e         = atmos_solarins_eccentricity
      omega     = atmos_solarins_perihelion_lon * d2r
 
    endif
 
    beta = sqrt( 1.0_rp - e*e )
 
    ! longitude at the vernal equinox associated with circle orbit
    lambda_m0 = 2.0_rp * ( ( 1.0_rp/2.0_rp*e + 1.0_rp/8.0_rp*e*e*e ) * ( 1.0_rp        + beta ) * sin(       omega) &
                         - (                   1.0_rp/4.0_rp*e*e   ) * ( 1.0_rp/2.0_rp + beta ) * sin(2.0_rp*omega) &
                         + (                   1.0_rp/8.0_rp*e*e*e ) * ( 1.0_rp/3.0_rp + beta ) * sin(3.0_rp*omega) )
 
    return

References atmos_solarins_eccentricity, atmos_solarins_obliquity, atmos_solarins_perihelion_lon, atmos_solarins_set_ideal, scale_const::const_d2r, scale_const::const_pi, scale_const::const_undef, and scale_precision::rp.

Referenced by atmos_solarins_setup().

Here is the caller graph for this function:

atmos_solarins_ecliptic_longitude()

subroutine, public scale_atmos_solarins::atmos_solarins_ecliptic_longitude	(	real(rp), intent(out)	Re_factor,
		real(rp), intent(out)	sinDEC,
		real(rp), intent(out)	cosDEC,
		real(rp), intent(out)	hourangle,
		integer, dimension(6), intent(in)	now_date,
		integer, intent(in)	offset_year,
		real(dp), intent(out), optional	lambda_out
	)

calc factor of Earths solar insolation

Definition at line 904 of file scale_atmos_solarins.F90.

    use scale_const, only: &
       pi => const_pi
    use scale_calendar, only: &
       calendar_ymd2absday,    &
       calendar_hms2abssec,    &
       i_year, i_month, i_day, &
       i_hour, i_min, i_sec
    implicit none
 
    real(RP), intent(out) :: Re_factor   ! factor of the distance of Earth from the sun (1/rho2)
    real(RP), intent(out) :: sinDEC      ! sin/cos(solar declination)
    real(RP), intent(out) :: cosDEC      ! sin/cos(solar declination)
    real(RP), intent(out) :: hourangle   ! hour angle: relative longitude of subsolar point (noon=0deg)
    integer,  intent(in)  :: now_date(6) ! date(yyyy,mm,dd,hh,mm,ss)
    integer,  intent(in)  :: offset_year ! year offset
    real(DP), intent(out), optional :: lambda_out   ! actual longitude from vernal equinox (for output)
 
 
    integer  :: date(6)
    integer  :: oyear
    integer  :: absday, absday_ve
    real(DP) :: abssec, abssec_ve
    real(DP) :: diffday, diffyr, diffsec
 
    real(DP) :: lambda_m ! mean   longitude from vernal equinox
    real(DP) :: lambda   ! actual longitude from vernal equinox
    real(DP) :: nu
 
    !---------------------------------------------------------------------------
 
    date(:) = now_date(:)
    oyear   = offset_year
 
    if ( atmos_solarins_fixeddate ) then
       if( atmos_solarins_date(1) >= 0 ) oyear   = 0
       if( atmos_solarins_date(1) >= 0 ) date(1) = atmos_solarins_date(1)
       if( atmos_solarins_date(2) >= 1 ) date(2) = atmos_solarins_date(2)
       if( atmos_solarins_date(3) >= 1 ) date(3) = atmos_solarins_date(3)
       if( atmos_solarins_date(4) >= 0 ) date(4) = atmos_solarins_date(4)
       if( atmos_solarins_date(5) >= 0 ) date(5) = atmos_solarins_date(5)
       if( atmos_solarins_date(6) >= 0 ) date(6) = atmos_solarins_date(6)
    endif
 
    call calendar_ymd2absday( absday,           & ! [OUT]
                              date(i_year),     & ! [IN]
                              date(i_month),    & ! [IN]
                              date(i_day),      & ! [IN]
                              oyear             ) ! [IN]
 
    call calendar_hms2abssec( abssec,           & ! [OUT]
                              date(i_hour),     & ! [IN]
                              date(i_min),      & ! [IN]
                              date(i_sec),      & ! [IN]
                              0.0_dp            ) ! [IN]
 
    call calendar_ymd2absday( absday_ve,        & ! [OUT]
                              ve_date(i_year),  & ! [IN]
                              ve_date(i_month), & ! [IN]
                              ve_date(i_day),   & ! [IN]
                              oyear             ) ! [IN]
 
    call calendar_hms2abssec( abssec_ve,        & ! [OUT]
                              ve_date(i_hour),  & ! [IN]
                              ve_date(i_min),   & ! [IN]
                              ve_date(i_sec),   & ! [IN]
                              0.0_dp            ) ! [IN]
 
    if ( year_calender_synced ) then
       diffday = real(absday-absday_ve,kind=dp) + (abssec-abssec_ve) / caldaysec
       diffyr  = diffday / calyearday - real( nint( diffday / calyearday ),kind=dp)
    else
       diffsec = real(absday-absday_ve,kind=dp)*caldaysec + (abssec-abssec_ve)
       diffyr  = diffsec / solyearsec - real( nint( diffsec / solyearsec ),kind=dp)
    endif
 
    lambda_m = lambda_m0 + 2.0_rp * pi * diffyr
 
    nu = lambda_m - omega
 
    ! 1 / (rho*rho)
    re_factor = 1.0_dp                          &
              + 2.0_dp*e     * cos(        nu ) &
              + 0.5_dp*e*e   * cos( 2.0_dp*nu ) &
              + 2.5_dp*e*e                      &
              + 4.0_dp*e*e*e * cos( 3.0_dp*nu )
 
    ! actual longitude from vernal equinox
    lambda = lambda_m &
           + ( 2.0_dp*e -  1.0_dp/ 4.0_dp*e*e*e ) * sin(        nu ) &
           + (             5.0_dp/ 4.0_dp*e*e   ) * sin( 2.0_dp*nu ) &
           + (            13.0_dp/12.0_dp*e*e*e ) * sin( 3.0_dp*nu )
 
    ! solar declination
    sindec = sin(lambda) * sin(obliquity)
    cosdec = sqrt( 1.0_rp - sindec*sindec )
 
    ! hour angle
    if ( day_calender_synced ) then
       hourangle = pi + 2.0_rp * pi * abssec / caldaysec
    else
       hourangle = pi + 2.0_rp * pi                                                    &
                      * ( (absday*caldaysec + abssec) / atmos_solarins_diurnal_sec     &
                        - real( nint((absday*caldaysec + abssec)/ soldaysec ),kind=dp) )
    endif
 
    if ( debug ) then
       log_info("ATMOS_SOLARINS_ecliptic_longitude",*) lambda_m, nu, lambda, re_factor
    endif
 
    if ( present(lambda_out) ) lambda_out = lambda
 
    return

References atmos_solarins_date, atmos_solarins_diurnal_sec, atmos_solarins_fixeddate, scale_calendar::calendar_hms2abssec(), scale_calendar::calendar_ymd2absday(), scale_const::const_pi, scale_precision::dp, scale_calendar::i_day, scale_calendar::i_hour, scale_calendar::i_min, scale_calendar::i_month, scale_calendar::i_sec, and scale_calendar::i_year.

Referenced by atmos_solarins_insolation_0d(), atmos_solarins_insolation_2d(), and atmos_solarins_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_solarins_insolation_0d()

subroutine scale_atmos_solarins::atmos_solarins_insolation_0d	(	real(rp), intent(in)	real_lon,
		real(rp), intent(in)	real_lat,
		integer, dimension(6), intent(in)	now_date,
		integer, intent(in)	offset_year,
		real(rp), intent(out)	solins,
		real(rp), intent(out)	cosSZA,
		real(rp), intent(out), optional	Re_factor_out
	)

calc factor of Earths solar insolation

Definition at line 1025 of file scale_atmos_solarins.F90.

    use scale_const, only: &
       eps => const_eps
    implicit none
    real(RP), intent(in)  :: real_lon    ! longitude
    real(RP), intent(in)  :: real_lat    ! latitude
    integer,  intent(in)  :: now_date(6) ! date(yyyy,mm,dd,hh,mm,ss)
    integer,  intent(in)  :: offset_year ! year offset
 
    real(RP), intent(out) :: solins      ! solar insolation
    real(RP), intent(out) :: cosSZA      ! cos(Solar Zenith Angle)
    real(RP), intent(out), optional :: Re_factor_out ! factor of the distance of Earth from the sun (1/rho2) [for output]
 
    real(RP) :: Re_factor      ! factor of the distance of Earth from the sun (1/rho2)
    real(RP) :: sinDEC, cosDEC ! sin/cos(solar declination)
    real(RP) :: hourangle      ! hour angle: relative longitude of subsolar point (noon=0deg)
 
    real(RP) :: lon
    real(RP) :: lat
    !---------------------------------------------------------------------------
 
    call atmos_solarins_ecliptic_longitude( re_factor,   & ! [OUT]
                                            sindec,      & ! [OUT]
                                            cosdec,      & ! [OUT]
                                            hourangle,   & ! [OUT]
                                            now_date(:), & ! [IN]
                                            offset_year  ) ! [IN]
 
    if ( atmos_solarins_fixedlatlon ) then
       lon = atmos_solarins_lon
       lat = atmos_solarins_lat
    else
       lon = real_lon
       lat = real_lat
    endif
 
    cossza = sin(lat)*sindec + cos(lat)*cosdec*cos(lon+hourangle)
    solins = atmos_solarins_constant * re_factor * ( 0.5_rp + sign(0.5_rp,cossza-eps) )
 
    if ( present( re_factor_out ) ) then
       re_factor_out = re_factor
    endif
 
    return

References atmos_solarins_constant, atmos_solarins_ecliptic_longitude(), atmos_solarins_fixedlatlon, atmos_solarins_lat, atmos_solarins_lon, and scale_const::const_eps.

Here is the call graph for this function:

atmos_solarins_insolation_2d()

subroutine scale_atmos_solarins::atmos_solarins_insolation_2d	(	integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension(ia,ja), intent(in)	real_lon,
		real(rp), dimension(ia,ja), intent(in)	real_lat,
		integer, dimension(6), intent(in)	now_date,
		integer, intent(in)	offset_year,
		real(rp), dimension (ia,ja), intent(out)	solins,
		real(rp), dimension (ia,ja), intent(out)	cosSZA
	)

calc factor of Earths solar insolation

Definition at line 1077 of file scale_atmos_solarins.F90.

    use scale_const, only: &
       eps => const_eps
    implicit none
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(in)  :: real_lon(IA,JA) ! longitude [rad]
    real(RP), intent(in)  :: real_lat(IA,JA) ! latitude  [rad]
    integer,  intent(in)  :: now_date(6)     ! date(yyyy,mm,dd,hh,mm,ss)
    integer,  intent(in)  :: offset_year     ! year offset
 
    real(RP), intent(out) :: solins  (IA,JA) ! solar insolation
    real(RP), intent(out) :: cosSZA  (IA,JA) ! cos(Solar Zenith Angle)
 
    real(RP) :: Re_factor      ! factor of the distance of Earth from the sun (1/rho2)
    real(RP) :: sinDEC, cosDEC ! sin/cos(solar declination)
    real(RP) :: hourangle      ! hour angle: relative longitude of subsolar point (noon=0deg)
 
    real(RP) :: lon(IA,JA)
    real(RP) :: lat(IA,JA)
 
    integer  :: i, j
    !---------------------------------------------------------------------------
 
    call atmos_solarins_ecliptic_longitude( re_factor,   & ! [OUT]
                                            sindec,      & ! [OUT]
                                            cosdec,      & ! [OUT]
                                            hourangle,   & ! [OUT]
                                            now_date(:), & ! [IN]
                                            offset_year  ) ! [IN]
 
    !$acc data create(lon, lat)
 
    if ( atmos_solarins_fixedlatlon ) then
       !$omp parallel do OMP_SCHEDULE_ collapse(2)
       !$acc kernels
       do j = js, je
       do i = is, ie
          lon(i,j) = atmos_solarins_lon
          lat(i,j) = atmos_solarins_lat
       enddo
       enddo
       !$acc end kernels
    else
       !$omp parallel do OMP_SCHEDULE_ collapse(2)
       !$acc kernels
       do j = js, je
       do i = is, ie
          lon(i,j) = real_lon(i,j)
          lat(i,j) = real_lat(i,j)
       enddo
       enddo
       !$acc end kernels
    endif
 
    !$omp parallel do OMP_SCHEDULE_ collapse(2)
    !$acc kernels
    do j = js, je
    do i = is, ie
       cossza(i,j) = sin(lat(i,j))*sindec + cos(lat(i,j))*cosdec*cos(lon(i,j)+hourangle)
       solins(i,j) = atmos_solarins_constant * re_factor * ( 0.5_rp + sign(0.5_rp,cossza(i,j)-eps) )
    enddo
    enddo
    !$acc end kernels
 
    !$acc end data
 
    return

References atmos_solarins_constant, atmos_solarins_ecliptic_longitude(), atmos_solarins_fixedlatlon, atmos_solarins_lat, atmos_solarins_lon, and scale_const::const_eps.

Here is the call graph for this function:

Variable Documentation

atmos_solarins_constant

real(rp), public scale_atmos_solarins::atmos_solarins_constant = 1360.250117_RP

Definition at line 43 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_constant     = 1360.250117_rp ! Solar constant [W/m2]

Referenced by atmos_solarins_insolation_0d(), atmos_solarins_insolation_2d(), atmos_solarins_setup(), mod_rm_driver::rm_driver(), mod_rm_prep::rm_prep(), and scale_atmos_phy_rd_mm5sw::swrad().

atmos_solarins_set_ve

logical, public scale_atmos_solarins::atmos_solarins_set_ve = .false.

Definition at line 45 of file scale_atmos_solarins.F90.

  logical,  public :: ATMOS_SOLARINS_set_ve       = .false.        ! Set vernal equinox condition?

Referenced by atmos_solarins_setup().

atmos_solarins_set_ideal

logical, public scale_atmos_solarins::atmos_solarins_set_ideal = .false.

Definition at line 47 of file scale_atmos_solarins.F90.

  logical,  public :: ATMOS_SOLARINS_set_ideal    = .false.        ! Set obliquity and eccentricity?

Referenced by atmos_solarins_orbit(), and atmos_solarins_setup().

atmos_solarins_obliquity

real(rp), public scale_atmos_solarins::atmos_solarins_obliquity = 23.44_RP

Definition at line 48 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_obliquity    = 23.44_rp       ! Obliquity [deg] cf. https://nssdc.gsfc.nasa.gov/planetary/factsheet/earthfact.html

Referenced by atmos_solarins_orbit(), and atmos_solarins_setup().

atmos_solarins_eccentricity

real(rp), public scale_atmos_solarins::atmos_solarins_eccentricity = 0.0167_RP

Definition at line 49 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_eccentricity = 0.0167_rp      ! Eccentricity

Referenced by atmos_solarins_orbit(), and atmos_solarins_setup().

atmos_solarins_perihelion_lon

real(rp), public scale_atmos_solarins::atmos_solarins_perihelion_lon = 282.94719_RP

Definition at line 50 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_perihelion_lon = 282.94719_rp ! Longitude of perihelion [deg].

Referenced by atmos_solarins_orbit(), and atmos_solarins_setup().

atmos_solarins_ve_date

integer, dimension(6), public scale_atmos_solarins::atmos_solarins_ve_date

Definition at line 51 of file scale_atmos_solarins.F90.

  integer,  public :: ATMOS_SOLARINS_ve_date(6)                    ! Date of first vernal equinox

Referenced by atmos_solarins_setup().

atmos_solarins_diurnal_sec

real(rp), public scale_atmos_solarins::atmos_solarins_diurnal_sec

Definition at line 53 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_diurnal_sec                   ! Seconds of the diurnal period [sec]

Referenced by atmos_solarins_ecliptic_longitude(), and atmos_solarins_setup().

atmos_solarins_annual_sec

real(rp), public scale_atmos_solarins::atmos_solarins_annual_sec

Definition at line 54 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_annual_sec                    ! Seconds of the annual period [sec]

Referenced by atmos_solarins_setup().

atmos_solarins_fixedlatlon

logical, public scale_atmos_solarins::atmos_solarins_fixedlatlon = .false.

Definition at line 56 of file scale_atmos_solarins.F90.

  logical,  public :: ATMOS_SOLARINS_fixedlatlon  = .false.        ! Latitude/Longitude is fixed?

Referenced by scale_atmos_phy_rd_profile::atmos_phy_rd_profile_read(), atmos_solarins_insolation_0d(), atmos_solarins_insolation_2d(), and atmos_solarins_setup().

atmos_solarins_lon

real(rp), public scale_atmos_solarins::atmos_solarins_lon

Definition at line 57 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_lon                           ! Longitude for radiation [rad]

Referenced by atmos_solarins_insolation_0d(), atmos_solarins_insolation_2d(), and atmos_solarins_setup().

atmos_solarins_lat

real(rp), public scale_atmos_solarins::atmos_solarins_lat

Definition at line 58 of file scale_atmos_solarins.F90.

  real(RP), public :: ATMOS_SOLARINS_lat                           ! Latitude  for radiation [rad]

Referenced by scale_atmos_phy_rd_profile::atmos_phy_rd_profile_read(), atmos_solarins_insolation_0d(), atmos_solarins_insolation_2d(), and atmos_solarins_setup().

atmos_solarins_fixeddate

logical, public scale_atmos_solarins::atmos_solarins_fixeddate = .false.

Definition at line 60 of file scale_atmos_solarins.F90.

  logical,  public :: ATMOS_SOLARINS_fixeddate    = .false.        ! Date is fixed?

Referenced by scale_atmos_phy_rd_profile::atmos_phy_rd_profile_read(), atmos_solarins_ecliptic_longitude(), and atmos_solarins_setup().

atmos_solarins_date

integer, dimension(6), public scale_atmos_solarins::atmos_solarins_date

Definition at line 61 of file scale_atmos_solarins.F90.

  integer,  public :: ATMOS_SOLARINS_date(6)                       ! Date      for radiation [Y,M,D,H,M,S]

Referenced by scale_atmos_phy_rd_profile::atmos_phy_rd_profile_read(), atmos_solarins_ecliptic_longitude(), and atmos_solarins_setup().