R/solarZenithAngle.R

solarZenithAngle <-
  function(lat,hr,i,tr=0.2618,hr0=12) {
    y <- acos(sin(lat)*sin(solarDecl(i))+cos(lat)*cos(solarDecl(i))*cos(tr*(hr-hr0)))
    y
  }

Try the sirad package in your browser

Any scripts or data that you put into this service are public.

sirad documentation built on May 1, 2019, 7:42 p.m.