util.h

#ifndef _UTIL_H
#define _UTIL_H
#include <math.h>
#include <stdio.h>
 
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif
 
struct row
{
  double lat, lon;
  int year, mon, day, hour;
  double temp, rh, ws, rain;
  /* Either need solar radiation to be included, or calculated */
  double solrad;
  /* derived from lat/lon/date/timezone */
  double sunrise, sunset;
  /* grass curing (%) [0-100]*/
  double percent_cured;
  /* grass fuel load (kg/m^2) */
  double grass_fuel_load;
};
 
struct row_daily
{
  double lat, lon;
  int year, mon, day;
  double temp, rh, wind, rain;
};
 
struct row_minmax
{
  double lat, lon;
  int year, mon, day;
  double temp_min, temp_max, rh_min, rh_max, ws_min, ws_max, rain;
};
int read_row(FILE* inp, struct row* r);
int read_row_daily(FILE* inp, struct row_daily* r);
int read_row_minmax(FILE* inp, struct row_minmax* r);
double findQ(double temp, double rh);
double findrh(double q, double temp);
double sun(double lat, double lon, int mon, int day, int hour, double timezone, double* sunrise, double* sunset);
double sun_julian(double lat, double lon, int jd, int hour, double timezone, double* sunrise, double* sunset);
int julian(int mon, int day);
void check_header(FILE* input, const char* header);
void check_inputs(double temp, double rh, double wind, double rain);
 
double seasonal_curing(int julian_date);
 
/* C90 round() only rounds to int but want digits */
double _round(double x, int digits);
 
/* C90 max() also causing problems */
double _max(double x, double y);
 
/* C90 min() also causing problems */
double _min(double x, double y);
#endif