| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149 |
- #include <stdio.h>
- #define FILL 0x1
- #define DIR 0x2
- #define BELEV 0x4
- #define TOPIDX 0x8
- #define IDXSTATS 0x10
- #define OUTPUT 0x20
- #define BUFSIZE 1024
- #define ZERO 0.0000001
- #define TOLERANCE 0.00001
- #define MAXITER 20
- #define NTERMS 10
- /* file_io.c */
- void get_line(FILE * fp, char *buffer);
- void read_input(void);
- void write_output(void);
- /* topmodel.c */
- void create_topidxstats(char *topidx, int ntopidxclasses, char *outtopidxstats);
- double calculate_lambda(void);
- void initialize(void);
- void calculate_flows(void);
- double calculate_efficiency(void);
- void calculate_others(void);
- void run_topmodel(void);
- /* infiltration.c */
- double calculate_infiltration(int timestep, double R);
- /* Topographic index statistics file */
- struct topidxstats
- {
- /* misc.ntopidxclasses */
- double *atb;
- double *Aatb_r;
- };
- /* Parameters file */
- struct params
- {
- char *name;
- double A;
- double qs0;
- double lnTe;
- double m;
- double Sr0;
- double Srmax;
- double td;
- double vch;
- double vr;
- int infex;
- double K0;
- double psi;
- double dtheta;
- int nch;
- /* params.nch's */
- double *d;
- double *Ad_r;
- };
- /* Input file */
- struct input
- {
- int ntimesteps;
- double dt;
- /* input.ntimestep's */
- double *R;
- double *Ep;
- };
- /* File names */
- struct file
- {
- char *params;
- char *topidxstats;
- char *input;
- char *output;
- char *obsflow;
- };
- /* Miscellaneous TOPMODEL variables */
- struct misc
- {
- /* Number of non-null cells */
- int ncells;
- /* Number of topographic index classes */
- int ntopidxclasses;
- /* Model efficiency */
- double Em;
- int ndelays;
- int nreaches;
- double lnTe;
- double vch;
- double vr;
- double lambda;
- double qss;
- double qs0;
- double Qobs_peak;
- double Qt_peak;
- double Qobs_mean;
- double Qt_mean;
- int tobs_peak;
- int tt_peak;
- /* params.nch's */
- double *tch;
- /* misc.nreach's */
- double *Ad;
- /* input.ntimestep's */
- double *Qobs;
- double *Qt;
- double *qs; /* spatially constant? */
- double *S_mean;
- double *f;
- double *fex;
- /* input.ntimestep * (misc.ntopidxclasses + 1)'s */
- double **qt;
- double **qo;
- double **qv;
- /* input.ntimestep * misc.ntopidxclassess */
- double **Srz;
- double **Suz;
- double **S;
- double **Ea;
- double **ex;
- /* Miscellaneous variables */
- int timestep;
- int topidxclass;
- };
- #ifdef _MAIN_C_
- #define GLOBAL
- #else
- #define GLOBAL extern
- #endif
- GLOBAL struct params params;
- GLOBAL struct topidxstats topidxstats;
- GLOBAL struct input input;
- GLOBAL struct file file;
- GLOBAL struct misc misc;
- /* Miscellaneous variables */
- GLOBAL char buf[BUFSIZE];
|
