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	atmos_solarins_ecliptic_longitude (Re_factor, sinDEC, cosDEC, hourangle, now_date, offset_year)
	calc factor of Earths solar insolation More...
subroutine	atmos_solarins_insolation_0d (real_lon, real_lat, now_date, offset_year, solins, cosSZA)
	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 = 0.0_RP
real(rp), public	atmos_solarins_eccentricity = 0.0_RP
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)	0.0_RP	Obliquity [deg]
  ATMOS_SOLARINS_ECCENTRICITY	real(RP)	0.0_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]
  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 556 of file scale_atmos_solarins.F90.

References atmos_solarins_constant, atmos_solarins_date, atmos_solarins_eccentricity, atmos_solarins_fixeddate, atmos_solarins_fixedlatlon, atmos_solarins_lat, atmos_solarins_lon, atmos_solarins_obliquity, atmos_solarins_orbit(), atmos_solarins_set_ideal, atmos_solarins_set_ve, scale_const::const_d2r, scale_io::io_fid_conf, and scale_prc::prc_abort().

Referenced by mod_atmos_driver::atmos_driver_setup().

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       const_d2r
    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,         &
       debug

    real(RP) :: dyear ! delta t [year]
    integer  :: year  ! used year

    integer  :: ierr
    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_SOLARINS_setup",*) 'Setup'

    atmos_solarins_lon     = basepoint_lon
    atmos_solarins_lat     = basepoint_lat
    atmos_solarins_date(:) = -1

    !--- 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
       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

    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 [deg] : ', lambda_m0 / const_d2r
    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

    return

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 672 of file scale_atmos_solarins.F90.

References atmos_solarins_eccentricity, atmos_solarins_obliquity, atmos_solarins_set_ideal, scale_const::const_d2r, and scale_const::const_pi.

Referenced by atmos_solarins_setup().

    use scale_const, only: &
       pi  => const_pi, &
       d2r => const_d2r
    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
    !---------------------------------------------------------------------------

    ! 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 )

    if ( atmos_solarins_set_ideal ) then
       obliquity = atmos_solarins_obliquity * d2r
       e         = atmos_solarins_eccentricity
    endif

    ! 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 perigee (see berger et al.(1993))
    omega = ( perih_omega + 180.0_rp ) * d2r

    beta = sqrt( 1.0_rp - e*e )

    ! longitude at the vernal equinox
    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

Here is the caller graph for this function:

atmos_solarins_ecliptic_longitude()

subroutine 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
	)

calc factor of Earths solar insolation

Definition at line 788 of file scale_atmos_solarins.F90.

References atmos_solarins_date, atmos_solarins_fixeddate, scale_calendar::calendar_getdayofyear(), scale_calendar::calendar_hms2abssec(), scale_calendar::calendar_ymd2absday(), scale_const::const_pi, 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(), and atmos_solarins_insolation_2d().

    use scale_const, only: &
       pi => const_pi
    use scale_calendar, only: &
       calendar_getdayofyear,  &
       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
    integer,  intent(in)  :: now_date(6) ! date(yyyy,mm,dd,hh,mm,ss)
    integer,  intent(in)  :: offset_year ! year offset

    integer  :: date(6)
    integer  :: oyear
    integer  :: absday, absday_ve
    real(DP) :: dayofyear, abssec

    real(RP) :: lambda_m ! mean   longitude from vernal equinox
    real(RP) :: lambda   ! actual longitude from vernal equinox
    real(RP) :: 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_getdayofyear( dayofyear, date(i_year) )

    call calendar_ymd2absday( absday,           & ! [OUT]
                              date(i_year),     & ! [IN]
                              date(i_month),    & ! [IN]
                              date(i_day),      & ! [IN]
                              oyear             ) ! [IN]

    call calendar_ymd2absday( absday_ve,        & ! [OUT]
                              date(i_year),  & ! [IN]
                              ve_date(i_month), & ! [IN]
                              ve_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]

    lambda_m = lambda_m0 + 2.0_rp * pi * real(absday-absday_ve,kind=RP) / dayofyear

    nu = lambda_m - omega

    ! 1 / (rho*rho)
    re_factor = 1.0_rp                          &
              + 2.0_rp*e     * cos(        nu ) &
              + 0.5_rp*e*e   * cos( 2.0_rp*nu ) &
              + 2.5_rp*e*e                      &
              + 4.0_rp*e*e*e * cos( 3.0_rp*nu )

    ! actual longitude from vernal equinox
    lambda = lambda_m &
           + ( 2.0_rp*e -  1.0_rp/ 4.0_rp*e*e*e ) * sin(        nu ) &
           + (             5.0_rp/ 4.0_rp*e*e   ) * sin( 2.0_rp*nu ) &
           + (            13.0_rp/12.0_rp*e*e*e ) * sin( 3.0_rp*nu )

    ! solar declination
    sindec = sin(lambda) * sin(obliquity)
    cosdec = sqrt( 1.0_rp - sindec*sindec )

    ! hour angle
    hourangle = 2.0_rp * pi * abssec / (24.0_rp*60.0_rp*60.0_rp)

    if ( debug ) then
       log_info("ATMOS_SOLARINS_ecliptic_longitude",*) lambda_m, nu, lambda
    endif

    return

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
	)

calc factor of Earths solar insolation

Definition at line 885 of file scale_atmos_solarins.F90.

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

    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) :: 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

    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) )

    return

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 932 of file scale_atmos_solarins.F90.

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

    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

    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]

    if ( atmos_solarins_fixedlatlon ) then
       !$omp parallel do OMP_SCHEDULE_ collapse(2)
       do j = js, je
       do i = is, ie
          lon(i,j) = atmos_solarins_lon
          lat(i,j) = atmos_solarins_lat
       enddo
       enddo
    else
       !$omp parallel do OMP_SCHEDULE_ collapse(2)
       do j = js, je
       do i = is, ie
          lon(i,j) = real_lon(i,j)
          lat(i,j) = real_lat(i,j)
       enddo
       enddo
    endif

    !$omp parallel do OMP_SCHEDULE_ collapse(2)
    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

    return

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 42 of file scale_atmos_solarins.F90.

Referenced by atmos_solarins_insolation_0d(), atmos_solarins_insolation_2d(), atmos_solarins_setup(), and scale_atmos_phy_rd_mm5sw::swrad().

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

atmos_solarins_set_ve

logical, public scale_atmos_solarins::atmos_solarins_set_ve = .false.

Definition at line 44 of file scale_atmos_solarins.F90.

Referenced by atmos_solarins_setup().

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

atmos_solarins_set_ideal

logical, public scale_atmos_solarins::atmos_solarins_set_ideal = .false.

Definition at line 46 of file scale_atmos_solarins.F90.

Referenced by atmos_solarins_orbit(), and atmos_solarins_setup().

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

atmos_solarins_obliquity

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

Definition at line 47 of file scale_atmos_solarins.F90.

Referenced by atmos_solarins_orbit(), and atmos_solarins_setup().

  real(RP), public :: atmos_solarins_obliquity    = 0.0_rp         ! Obliquity [deg]

atmos_solarins_eccentricity

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

Definition at line 48 of file scale_atmos_solarins.F90.

Referenced by atmos_solarins_orbit(), and atmos_solarins_setup().

  real(RP), public :: atmos_solarins_eccentricity = 0.0_rp         ! Eccentricity

atmos_solarins_fixedlatlon

logical, public scale_atmos_solarins::atmos_solarins_fixedlatlon = .false.

Definition at line 50 of file scale_atmos_solarins.F90.

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

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

atmos_solarins_lon

real(rp), public scale_atmos_solarins::atmos_solarins_lon

Definition at line 51 of file scale_atmos_solarins.F90.

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

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

atmos_solarins_lat

real(rp), public scale_atmos_solarins::atmos_solarins_lat

Definition at line 52 of file scale_atmos_solarins.F90.

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

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

atmos_solarins_fixeddate

logical, public scale_atmos_solarins::atmos_solarins_fixeddate = .false.

Definition at line 54 of file scale_atmos_solarins.F90.

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

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

atmos_solarins_date

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

Definition at line 55 of file scale_atmos_solarins.F90.

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

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