Calculate the solar zenith angle for given times and locations
Usage
zenith(sun, lon, lat)
Arguments
sun
list of solar time and declination computed by solar.
lon
vector of longitudes.
lat
vector latitudes.
Details
zenith uses the solar time and declination calculated by
solar to compute the solar zenith angle for given times and
locations, using the same methods as
www.esrl.noaa.gov/gmd/grad/solcalc/. This function does not
adjust for atmospheric refraction see refracted.
Value
A vector of solar zenith angles (degrees) for the given
locations and times.
See Also
solar
Examples
## Approx location of Sydney Harbour Bridge
lon <- 151.211
lat <- -33.852
## Solar zenith angle for noon on the first of May 2000
## at the Sydney Harbour Bridge
s <- solar(as.POSIXct("2000-05-01 12:00:00","EST"))
zenith(s,lon,lat)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(GeoLight)
Loading required package: maps
# maps v3.1: updated 'world': all lakes moved to separate new #
# 'lakes' database. Type '?world' or 'news(package="maps")'. #
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/GeoLight/zenith.Rd_%03d_medium.png", width=480, height=480)
> ### Name: zenith
> ### Title: Solar Zenith Angle
> ### Aliases: zenith
>
> ### ** Examples
>
> ## Approx location of Sydney Harbour Bridge
> lon <- 151.211
> lat <- -33.852
> ## Solar zenith angle for noon on the first of May 2000
> ## at the Sydney Harbour Bridge
> s <- solar(as.POSIXct("2000-05-01 12:00:00","EST"))
> zenith(s,lon,lat)
[1] 133.9373
>
>
>
>
>
> dev.off()
null device
1
>