---

WindAzimuth.c

Go to the documentation of this file.

#include "WindAzimuth.h"

int GetWindAzimuthFIXEDFromProps(ChHashTable * proptbl, int month, int day, int hour,
                                                unsigned int row, unsigned int col, double * waz)   {
    /* static variables used to store state across function calls */
    static int smonth               = 0;
    static int sday                 = 0;
    static int shour                = 0;
    static double swaz              = 0.0;
    static DblTwoDArray * swaz_tbl  = NULL;
    /* stack variables */
    KeyVal * entry                  = NULL;             /* key/val instances from table */
    FILE * fstream                  = NULL;             /* file stream */
    int rec;

    if ( (smonth != month) || (sday != day) || (shour != hour) )    {
        /* new wind azimuth needed */
        if (swaz_tbl == NULL )  {
            *waz = swaz;        
            /* table of fixed values only needs to be initialized once */
            if ( ChHashTableRetrieve(proptbl, GetFireProp(PROP_WAZFFILE), (void *)&entry) ) {
                ERR_ERROR("Unable to retrieve WIND_AZIMUTH_FIXED_FILE property. \n", ERR_EINVAL);
                }
            /* open fixed azimuth file */
            if ( (fstream = fopen((char *) entry->val, "r")) == NULL )  {
                ERR_ERROR("Unable to open file containing fixed Wind Azimuths. \n", ERR_EIOFAIL);
                }
            /* create the table of values */
            if ( (swaz_tbl = GetDblTwoDArrayTableFStreamIO(fstream, 
                    WIND_AZIMUTH_SEP_CHARS, WIND_AZIMUTH_WAZ_COMMENT_CHAR)) == NULL )   {
                fclose(fstream);
                ERR_ERROR("Unable to initialize fixed Wind Azimuth table. \n", ERR_EFAILED);
                }
            fclose(fstream);
            }
        /* find record in table */
        row = DBLTWODARRAY_SIZE_ROW(swaz_tbl);
        for(rec = 0; rec < row; rec++)  {
            /* find matching month and day */
            if ( DBLTWODARRAY_GET_DATA(swaz_tbl, rec, WIND_AZIMUTH_MO_WAZ_TBL_INDEX) == month 
                 && DBLTWODARRAY_GET_DATA(swaz_tbl, rec, WIND_AZIMUTH_DY_WAZ_TBL_INDEX) == day )    {
                /* retrieve waz on that month and day */
                col = WIND_AZIMUTH_HR_TO_WAZ_TBL_INDEX(hour);
                swaz = DBLTWODARRAY_GET_DATA(swaz_tbl, rec, col);
                break;
                }
            }
        /* set {month, day, hour} for future calls */
        smonth = month;
        sday = day;
        shour = hour;
        }

    *waz = swaz;

    return ERR_SUCCESS; 
    }

int GetWindAzimuthRANDUFromProps(ChHashTable * proptbl, int month, int day, int hour,
                                                unsigned int row, unsigned int col, double * waz)   {
    /* static variables used to store state across function calls */
    static int smonth               = 0;
    static int sday                 = 0;
    static int shour                = 0;
    static double swaz              = 0.0;
    
    /* WA: stop warnings from unused vars */
    row = col = 0;
    if ( proptbl == NULL )  {
        ERR_ERROR("Simulation properties table not initialized. \n", ERR_EINVAL);
        }

    if ( (smonth != month) || (sday != day) || (shour != hour) )    {   
        /* new wind azimuth needed */
        swaz = randu(WIND_AZIMUTH_RANDU_MIN_AZ, WIND_AZIMUTH_RANDU_MAX_AZ);
        /* set {month, day, hour} for future calls */
        smonth = month;
        sday = day;
        shour = hour;       
        }

    *waz = swaz;

    return ERR_SUCCESS;
    }

int GetWindAzimuthRANDHFromProps(ChHashTable * proptbl, int month, int day, int hour,
                                                unsigned int row, unsigned int col, double * waz)   {
    /* static variables used to store state across function calls */
    static int smonth               = 0;
    static int sday                 = 0;
    static int shour                = 0;
    static double swaz              = 0.0;
    static DblTwoDArray * swaz_tbl  = NULL;
    /* stack variables */
    KeyVal * entry                  = NULL;             /* key/val instances from table */
    FILE * fstream                  = NULL;             /* file stream */

    if ( (smonth != month) || (sday != day) || (shour != hour) )    {
        /* new wind azimuth needed */
        if (swaz_tbl == NULL )  {
            *waz = swaz;        
            /* table of fixed values only needs to be initialized once */
            if ( ChHashTableRetrieve(proptbl, GetFireProp(PROP_WAZHFILE), (void *)&entry) ) {
                ERR_ERROR("Unable to retrieve WIND_AZIMUTH_HISTORICAL_FILE property. \n", ERR_EINVAL);
                }
            /* open fixed azimuth file */
            if ( (fstream = fopen((char *) entry->val, "r")) == NULL )  {
                ERR_ERROR("Unable to open file containing historical Wind Azimuths. \n", ERR_EIOFAIL);
                }
            /* create the table of values */
            if ( (swaz_tbl = GetDblTwoDArrayTableFStreamIO(fstream, 
                    WIND_AZIMUTH_SEP_CHARS, WIND_AZIMUTH_WAZ_COMMENT_CHAR)) == NULL )   {
                fclose(fstream);
                ERR_ERROR("Unable to initialize historical Wind Azimuth table. \n", ERR_EFAILED);
                }
            fclose(fstream);                    
            }
        /* find wind azimuth from random record in table that is not NO DATA */
        do  {
            /* retrieve a record at random from table within range 0 to num_recs */
            /* convert u.r.n. range 0-max to range a-b = a + (((b-a)*rand)/randmax) */
            /*row = 0 + (((DBLTWODARRAY_SIZE_ROW(swaz_tbl) - 0) * randi(0)) / LONG_MAX);*/
            /* appears to be flawed bc always returns 0 */
            row = randi(0) % DBLTWODARRAY_SIZE_ROW(swaz_tbl);
            col = WIND_AZIMUTH_HR_TO_WAZ_TBL_INDEX(hour);
            /* retrieve waz from random recno using current hour */
            swaz = DBLTWODARRAY_GET_DATA(swaz_tbl, row, col);       
            } while( swaz == WIND_AZIMUTH_WAZ_NO_DATA_VALUE );
        /* set {month, day, hour} for future calls */
        smonth = month;
        sday = day;
        shour = hour;
        }

    *waz = swaz;

    return ERR_SUCCESS;         
    }

int GetWindAzimuthSPATIALFromProps(ChHashTable * proptbl, int month, int day, int hour, 
                                                unsigned int row, unsigned int col, double * waz)   {
    /* WA: stop warnings from unused vars */
    month = day = hour = row = col = *waz =0;
    if ( proptbl == NULL )  {
        ERR_ERROR("Simulation properties table not initialized. \n", ERR_EINVAL);
        }
    
    ERR_ERROR("SPATIAL option for WIND_AZIMUTH_TYPE not yet implemented. \n", ERR_EUNIMPL);
    }
 
/* end of WindAzimuth.c */

---