/////////////////////////////////////////////////////////////////////////////////////////// // Calculate all SPA parameters and put into structure // Note: All inputs values (listed in header file) must already be in structure /////////////////////////////////////////////////////////////////////////////////////////// void spa_calculate(spa_data *spa) { spa->jd = julian_day (spa->year, spa->month, spa->day, spa->hour, spa->minute, spa->second, spa->timezone); calculate_geocentric_sun_right_ascension_and_declination(spa); spa->h = observer_hour_angle(spa->nu, spa->longitude, spa->alpha); spa->xi = sun_equatorial_horizontal_parallax(spa->r); sun_right_ascension_parallax_and_topocentric_dec(spa->latitude, spa->elevation, spa->xi, spa->h, spa->delta, &(spa->del_alpha), &(spa->delta_prime)); spa->alpha_prime = topocentric_sun_right_ascension(spa->alpha, spa->del_alpha); spa->h_prime = topocentric_local_hour_angle(spa->h, spa->del_alpha); spa->e0 = topocentric_elevation_angle(spa->latitude, spa->delta_prime, spa->h_prime); spa->del_e = atmospheric_refraction_correction(spa->pressure, spa->temperature, spa->e0); spa->e = topocentric_elevation_angle_corrected(spa->e0, spa->del_e); spa->zenith = topocentric_zenith_angle(spa->e); spa->azimuth180 = topocentric_azimuth_angle_neg180_180(spa->h_prime, spa->latitude, spa->delta_prime); spa->azimuth = topocentric_azimuth_angle_zero_360(spa->azimuth180); spa->incidence = surface_incidence_angle(spa->zenith, spa->azimuth180, spa->azm_rotation, spa->slope); spa->altitude = 90.0 - spa->incidence; //**KMB calculate_eot_and_sun_rise_transit_set(spa); }
/////////////////////////////////////////////////////////////////////////////////////////// // Calculate all SPA parameters and put into structure // Note: All inputs values (listed in header file) must already be in structure /////////////////////////////////////////////////////////////////////////////////////////// int spa_calculate(spa_data *spa) { int result; result = validate_inputs(spa); if (result == 0) { spa->jd = julian_day (spa->year, spa->month, spa->day, spa->hour, spa->minute, spa->second, spa->timezone); calculate_geocentric_sun_right_ascension_and_declination(spa); spa->h = observer_hour_angle(spa->nu, spa->longitude, spa->alpha); spa->xi = sun_equatorial_horizontal_parallax(spa->r); sun_right_ascension_parallax_and_topocentric_dec(spa->latitude, spa->elevation, spa->xi, spa->h, spa->delta, &(spa->del_alpha), &(spa->delta_prime)); spa->alpha_prime = topocentric_sun_right_ascension(spa->alpha, spa->del_alpha); spa->h_prime = topocentric_local_hour_angle(spa->h, spa->del_alpha); spa->e0 = topocentric_elevation_angle(spa->latitude, spa->delta_prime, spa->h_prime); spa->del_e = atmospheric_refraction_correction(spa->pressure, spa->temperature, spa->atmos_refract, spa->e0); spa->e = topocentric_elevation_angle_corrected(spa->e0, spa->del_e); spa->zenith = topocentric_zenith_angle(spa->e); spa->azimuth180 = topocentric_azimuth_angle_neg180_180(spa->h_prime, spa->latitude, spa->delta_prime); spa->azimuth = topocentric_azimuth_angle_zero_360(spa->azimuth180); if ((spa->function == SPA_ZA_INC) || (spa->function == SPA_ALL)) spa->incidence = surface_incidence_angle(spa->zenith, spa->azimuth180, spa->azm_rotation, spa->slope); if ((spa->function == SPA_ZA_RTS) || (spa->function == SPA_ALL)) calculate_eot_and_sun_rise_transit_set(spa); } return result; }
void calculate_eot_and_sun_rise_transit_set(spa_data *spa) { spa_data sun_rts; double nu, m, h0, n; double alpha[JD_COUNT], delta[JD_COUNT]; double m_rts[SUN_COUNT], nu_rts[SUN_COUNT], h_rts[SUN_COUNT]; double alpha_prime[SUN_COUNT], delta_prime[SUN_COUNT], h_prime[SUN_COUNT]; double h0_prime = -1*(SUN_RADIUS + spa->atmos_refract); int i; sun_rts = *spa; m = sun_mean_longitude(spa->jme); spa->eot = eot(m, spa->alpha, spa->del_psi, spa->epsilon); sun_rts.hour = sun_rts.minute = sun_rts.second = 0; sun_rts.timezone = 0.0; sun_rts.jd = julian_day (sun_rts.year, sun_rts.month, sun_rts.day, sun_rts.hour, sun_rts.minute, sun_rts.second, sun_rts.timezone); calculate_geocentric_sun_right_ascension_and_declination(&sun_rts); nu = sun_rts.nu; sun_rts.delta_t = 0; sun_rts.jd--; for (i = 0; i < JD_COUNT; i++) { calculate_geocentric_sun_right_ascension_and_declination(&sun_rts); alpha[i] = sun_rts.alpha; delta[i] = sun_rts.delta; sun_rts.jd++; } m_rts[SUN_TRANSIT] = approx_sun_transit_time(alpha[JD_ZERO], spa->longitude, nu); h0 = sun_hour_angle_at_rise_set(spa->latitude, delta[JD_ZERO], h0_prime); if (h0 >= 0) { approx_sun_rise_and_set(m_rts, h0); for (i = 0; i < SUN_COUNT; i++) { nu_rts[i] = nu + 360.985647*m_rts[i]; n = m_rts[i] + spa->delta_t/86400.0; alpha_prime[i] = rts_alpha_delta_prime(alpha, n); delta_prime[i] = rts_alpha_delta_prime(delta, n); h_prime[i] = limit_degrees180pm(nu_rts[i] + spa->longitude - alpha_prime[i]); h_rts[i] = rts_sun_altitude(spa->latitude, delta_prime[i], h_prime[i]); } spa->srha = h_prime[SUN_RISE]; spa->ssha = h_prime[SUN_SET]; spa->sta = h_rts[SUN_TRANSIT]; spa->suntransit = dayfrac_to_local_hr(m_rts[SUN_TRANSIT] - h_prime[SUN_TRANSIT] / 360.0, spa->timezone); spa->sunrise = dayfrac_to_local_hr(sun_rise_and_set(m_rts, h_rts, delta_prime, spa->latitude, h_prime, h0_prime, SUN_RISE), spa->timezone); spa->sunset = dayfrac_to_local_hr(sun_rise_and_set(m_rts, h_rts, delta_prime, spa->latitude, h_prime, h0_prime, SUN_SET), spa->timezone); } else spa->srha= spa->ssha= spa->sta= spa->suntransit= spa->sunrise= spa->sunset= -99999; }