| 1 | /* titan.h: parameters for gptitan.c */ |
|---|
| 2 | |
|---|
| 3 | #include <stdio.h> |
|---|
| 4 | #include <stdlib.h> |
|---|
| 5 | #include <string.h> |
|---|
| 6 | #include <math.h> |
|---|
| 7 | |
|---|
| 8 | #define R0 (double)(2575.0) /* Titan's radius */ |
|---|
| 9 | #define NLEV (int)(55) /* Nbre de niv verticaux - aussi dans titan_for.h */ |
|---|
| 10 | |
|---|
| 11 | /* DEPEND DE LA VERSION CHIMIE: */ |
|---|
| 12 | #define VERCHIM "chimie_simpnit_051006_bis" |
|---|
| 13 | #define NREAC (int)(377) /* nombre de reactions - aussi dans titan_for.h */ |
|---|
| 14 | #define RDISS (int)(54) /* nombre de photodiss - aussi dans titan_for.h */ |
|---|
| 15 | #define NC (int)(44) /* nb de composes - aussi dans titan_for.h */ |
|---|
| 16 | #define ST (int)(NC) /* nb de composes inverses */ |
|---|
| 17 | #define NHC (int)(32) /* nb hydrocarbons */ |
|---|
| 18 | |
|---|
| 19 | #define THETA (double)(0.501) |
|---|
| 20 | #ifndef M_PI |
|---|
| 21 | #define M_PI (double)(3.14159265358979323846e0) |
|---|
| 22 | #endif |
|---|
| 23 | #define RAD (double)(M_PI / 180.0e0) |
|---|
| 24 | #ifndef max |
|---|
| 25 | #define max(a,b) ((a)>(b)?(a):(b)) |
|---|
| 26 | #define min(a,b) ((a)<=(b)?(a):(b)) |
|---|
| 27 | #endif |
|---|
| 28 | |
|---|
| 29 | /* void gptitan_(const int *, float *, float *, float *, float *, float *, |
|---|
| 30 | char (*)[10], float (*)[NLEV], float *, |
|---|
| 31 | float *, float *, float *, int *, float *, |
|---|
| 32 | float *, float (*)[NLEV][RDISS+1][15], float (*)[NLEV], |
|---|
| 33 | int (*)[5], int *, int *, int (*)[200], int (*)[200][2] ); */ |
|---|
| 34 | void chimie_(char (*)[10], double *, double *, double (*)[NLEV], |
|---|
| 35 | int (*)[5], int *, int *, int (*)[200][2], int (*)[200]); |
|---|
| 36 | void comp_(char (*)[10], double *); |
|---|
| 37 | void disso_(double (*)[NLEV][RDISS+1][15], int *); |
|---|
| 38 | void solve( double **, double *, int, int ); |
|---|
| 39 | float *rm1d( int, int ); |
|---|
| 40 | float **rm2d( int, int, int, int ); |
|---|
| 41 | float ***rm3d( int, int, int, int, int, int ); |
|---|
| 42 | float ****rm4d( int, int, int, int, int, int, int, int ); |
|---|
| 43 | double *dm1d( int, int ); |
|---|
| 44 | double **dm2d( int, int, int, int ); |
|---|
| 45 | double ***dm3d( int, int, int, int, int, int ); |
|---|
| 46 | double ****dm4d( int, int, int, int, int, int, int, int ); |
|---|
| 47 | void frm1d( float *, int ); |
|---|
| 48 | void frm2d( float **, int, int, int ); |
|---|
| 49 | void frm3d( float ***, int, int, int, int, int ); |
|---|
| 50 | void fdm1d( double *, int ); |
|---|
| 51 | void fdm2d( double **, int, int, int ); |
|---|
| 52 | void fdm3d( double ***, int, int, int, int, int ); |
|---|
| 53 | int *im1d( int, int ); |
|---|
| 54 | int **im2d( int, int, int, int ); |
|---|
| 55 | int ***im3d( int, int, int, int, int, int ); |
|---|
