CalcSZA.c

Go to the documentation of this file.

#include <math.h>
#include <string.h>
#include <stdio.h>

#include "solpos.h"     /* <-- There is the 'include' I was talking about */
/*============================================================================
*
*    NAME:  soltest.c
*
*    PURPOSE:  Exercises the solar position algorithms in 'solpos.c'.
*
*        S_solpos
*            INPUTS:     year, daynum, hour, minute, second, latitude,
*                        longitude, timezone
*
*            OPTIONAL:   press   DEFAULT 1013.0 (standard pressure)
*                        temp    DEFAULT   10.0 (standard temperature)
*                        tilt    DEFAULT    0.0 (horizontal panel)
*                        aspect  DEFAULT  180.0 (South-facing panel)
*                        month   (if the S_DOY function is turned off)
*                        day     ( "             "             "     )
*
*            OUTPUTS:    amass, ampress, azim, cosinc, coszen, day, daynum,
*                        elevref, etr, etrn, etrtilt, month, prime,
*                        sbcf, sretr, ssetr, unprime, zenref
*
*       S_init        (optional initialization for all input parameters in
*                      the posdata struct)
*           INPUTS:     struct posdata*
*           OUTPUTS:    struct posdata*
*
*                     (Note: initializes the required S_solpos INPUTS above
*                      to out-of-bounds conditions, forcing the user to
*                      supply the parameters; initializes the OPTIONAL
*                      S_solpos inputs above to nominal values.)
*
*      S_decode       (optional utility for decoding the S_solpos return code)
*           INPUTS:     long int S_solpos return value, struct posdata*
*           OUTPUTS:    text to stderr
*
*
*        All variables are defined as members of the struct posdata
*        in 'solpos.h'.
*
*    Usage:
*         In calling program, along with other 'includes', insert:
*
*              #include "solpos.h"
*
*    Martin Rymes
*    National Renewable Energy Laboratory
*    25 March 1998
*
*    28 March 2001 REVISION:  SMW changed benchmark numbers to reflect the
*                             February 2001 changes to solpos00.c
*
*----------------------------------------------------------------------------*/


int SZACalc(double d_lat, double d_long, double d_tzone, int d_year, int d_daynum,
            int d_hour, int d_minute, int d_second, float d_temp, float d_press,
            float d_tilt, float d_aspect)
{
       //Coordinate TNB
 //   S_longitude = 164.123;  /* Note that latitude and longitude are  */
 //   S_latitude  =  -74.695;  /*   in DECIMAL DEGREES, not Deg/Min/Sec */
 //   S_timezone  =  11.0;   /* Eastern time zone, even though longitude would
    //              suggest Central.  We use what they use.
    //              DO NOT ADJUST FOR DAYLIGHT SAVINGS TIME. */
    //Coordinate StaraZagora BG
    //  S_longitude = 25.633;
    //  S_latitude  = 42.416;
    //  S_timezone  = 2;

    //Coordinate Bologna
    //  S_longitude = 11.333;
    //  S_latitude  = 44.483;
    //  S_timezone  = 1;

    //Coordinate GreenWitch
    //  S_longitude = 0;
    //  S_latitude  = 51.29;
    //  S_timezone  = 0;

    //Coordinate 0,90,0
    //  S_longitude = 0;
    //  S_latitude  = 90;
    //  S_timezone  = 0;

    //Coordinate Mt Cimone
 //     S_longitude = 10.42;
 //     S_latitude  = 44.11;
 //     S_timezone  = 1;

      S_longitude = d_long;
      S_latitude  = d_lat;
      S_timezone  = d_tzone;
      S_year      = d_year;
      S_daynum    = d_daynum;

    S_date ( );

    S_hour = d_hour;
    S_minute = d_minute;
    S_second = d_second;

    S_temp      =   d_temp;
    S_press     = d_press;

    S_tilt      = S_latitude;  /* Tilted at latitude */
    S_aspect    = d_aspect;       /* 135 deg. = SE */

    S_solpos ( );

    return 0;
}