root/trunk/EBaporate/ebsurf.h

/***************************************
  $Header$

  This header file contains definitions, structures and function prototypes for
  ebsurf.c and shorts.c.
***************************************/


#ifndef EBSURF_H
#define EBSURF_H    /*+ To stop multiple inclusions. +*/

#include <stdio.h>
#include <stdlib.h> /* For getenv() */
#include <math.h>
#define MALLOC(a,b) (a *)malloc((b)*sizeof(a))
#define MIN(x,y) (((x)>(y))?(y):(x))
#ifndef M_PI
#define M_PI 3.141592653589793
#endif
#ifndef _ABS
#define _ABS(x) (((x)<0)?-(x):(x))
#endif
#define GR 8.314e7               /* Gas constant in cgs units    */
#define TORR2CGSP (1013000./760.)/* Torr to dyn/cm^2 conversion  */
#define SIGMA 5.67e-5            /* Radiation constant in cgs    */
#define MAXCYC 30                /* Max number of cycles in auto */

/* Messages passed to output() */
#define GENZ 0 /* Graph beam heat gen vs. depth z */
#define GENT 5 /* Graph beam heat gen vs. time t */
#define TMPS 1 /* Graph temperatures vs. z and t (3-D) */
#define SURF 2 /* Graph surface temperatures vs. t */
#define ONEC 3 /* Graph one cycle's temperatures vs. z and t */
#define FLUX 4 /* Graph evaporative flux from surface vs. t */
#define PATT 6 /* Graph x-y pattern */
#define PSDIR     "/home/hazelsct/cvs/EBaporate/" /* Dir with PS templates */
#define OUTPUTDIR "/home/hazelsct/cvs/EBaporate/" /* Dir to put PS output in */
#define OUTPUTURLDIR OUTPUTDIR                    /* Dir for graph links */

/* Flags to pass to cycles() */
#define GREET 1
#define OUTPUT 2
#define TSUR 4   /* T0 = surface temperature, not layer bottom temperature */

struct param {
  unsigned long zs,ts,cyc,ptyp,xw,yw;
  double (*gtime)(); /* gen(time) function pointer */
  double (*xpat)(), (*ypat)(); /* x(t) and y(t) function pointers */
  double al,k,   /* Material props: alpha(=k/rho cp), therm conductivity */
    rho,mu,dsig, /* Dens, visc and Marangoni dsigma/dt to test regime III/IV */
    Ma,eps,Tenv, /* Molar mass, radiative emissivity, env temp for radiation */
    Q0,S,        /* Heat of vaporiz, modeled as Q0 exp (ST) */
    A,B,C,D,E;   /* Clausius-Clapeyron evap coefficients, Langmuir coeff */
  double pow,volt,dp,   /* Run parameters: power, voltage, penetration depth */
    tfin,xfrq,yfrq,plen,pwid, /* period, x and y freqs, pattern length * wid */
    vb,spx,spy,               /* beam veloc, x and y spot size */
    x,y,tres; /* x and y point where our model runs, beam dwell time */
  double zr,tr,b,ap,bp,w,th,it,ct,Tin,Tsurf; /* Sim params */
  char *gout,*gtou,*tout,*sout,*lout,*fout,*pout; }; /* Output file names */

/* Function prototypes from ebsurf.c */
int cycles (double *X, double *t, double *T, double T0, double *big,
            struct param *p, int m, double *ou, FILE *fule, int pt);
double timestep (struct param p, double *X, double *t, double *T, double *big,
                 double Gto, int s, int m);
char *paraminp (struct param *pars, char *qstring);
void initcond (double *X, double *t, double *T, double T0, double TS,
               double *big, struct param *p);
void output (struct param *p, double *X, double *t, double *T, int msg);
void tridag (double *a, double *b, double *c, double *r, double *u, int n);
void errer (char *strong);

/* Function prototypes from shorts.c */
double radloss (double T, struct param *p);
double loss (double T, struct param *p);
double vloss (double T, double E, double A, double C, double D);
double vloss2 (double T, double E, double A, double C, double D);
double gtnorm (double t, struct param *p);
double gtliss (double t, struct param *p);
double gtdisk (double t, struct param *p);
double sino(double t, double freq);
double triang(double t, double freq);
double sawt(double t, double freq);

#endif /* EBSURF_H */