[3] | 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 */ |
---|
[1965] | 9 | #define NLEV (int)(133) /* Nbre de niv verticaux -> Need to be coherent with the vertical grid used !! */ |
---|
[1908] | 10 | #define NLD (int)(40) /* Nbre de niv verticaux faits sans diff -> Need to be coherent with the vertical grid used !! */ |
---|
[2326] | 11 | //#define NLD (int)(0) /* -> This is for 1D, diffusion on all column */ |
---|
[1126] | 12 | #define NLRT (int)(650) /* Nbre de niv verticaux dans table fmoy - aussi dans common_mod */ |
---|
[3] | 13 | |
---|
[1126] | 14 | /* fluxes at 1300 km : upward is +, downward is - */ |
---|
| 15 | #define top_H (double)(+1.1e4) |
---|
| 16 | #define top_H2 (double)(+3.7e3) |
---|
| 17 | #define top_N4S (double)(-1.1e8) /* = -2.5e8/2.27 ... */ |
---|
| 18 | |
---|
[3] | 19 | /* DEPEND DE LA VERSION CHIMIE: */ |
---|
[2326] | 20 | //#define VERCHIM "chimie_simpnit_051006_bis" -> If you want to use this one make sure you use an old version of disso.c |
---|
| 21 | //#define NREAC (int)(377) /* nombre de reactions - aussi dans common_mod */ |
---|
| 22 | //#define RDISS (int)(54) /* nombre de photodiss - aussi dans common_mod */ |
---|
| 23 | //#define NC (int)(44) /* nb de composes - aussi dans common_mod */ |
---|
| 24 | //#define ST (int)(NC) /* nb de composes inverses */ |
---|
| 25 | //#define NHC (int)(32) /* nb hydrocarbons */ |
---|
| 26 | #define VERCHIM "chimie_2019_43spc_vuitton" |
---|
| 27 | #define NREAC (int)(330) /* nombre de reactions - aussi dans common_mod */ |
---|
| 28 | #define RDISS (int)(55) /* nombre de photodiss - aussi dans common_mod */ |
---|
| 29 | #define NC (int)(43) /* nb de composes - aussi dans common_mod */ |
---|
[3] | 30 | #define ST (int)(NC) /* nb de composes inverses */ |
---|
[2326] | 31 | #define NHC (int)(31) /* nb hydrocarbons */ |
---|
[3] | 32 | |
---|
| 33 | #define THETA (double)(0.501) |
---|
| 34 | #ifndef M_PI |
---|
| 35 | #define M_PI (double)(3.14159265358979323846e0) |
---|
| 36 | #endif |
---|
| 37 | #define RAD (double)(M_PI / 180.0e0) |
---|
| 38 | #ifndef max |
---|
| 39 | #define max(a,b) ((a)>(b)?(a):(b)) |
---|
| 40 | #define min(a,b) ((a)<=(b)?(a):(b)) |
---|
| 41 | #endif |
---|
| 42 | |
---|
| 43 | void chimie_(char (*)[10], double *, double *, double (*)[NLEV], |
---|
| 44 | int (*)[5], int *, int *, int (*)[200][2], int (*)[200]); |
---|
[1126] | 45 | void comp_(char (*)[10], double *, double *, double *, double (*)[NLEV]); |
---|
[1950] | 46 | void disso_(double (*)[15][NLRT][RDISS+1], int *); |
---|
[1126] | 47 | double omega( double, double, double ); |
---|
| 48 | void solve( double ***, int, int, int ); |
---|
[1965] | 49 | void solve_lapack( double ***, int, int, int ); |
---|
[1126] | 50 | void solve_b( double ***, double **, int, int, int ); |
---|
[3] | 51 | float *rm1d( int, int ); |
---|
| 52 | float **rm2d( int, int, int, int ); |
---|
| 53 | float ***rm3d( int, int, int, int, int, int ); |
---|
| 54 | float ****rm4d( int, int, int, int, int, int, int, int ); |
---|
| 55 | double *dm1d( int, int ); |
---|
| 56 | double **dm2d( int, int, int, int ); |
---|
| 57 | double ***dm3d( int, int, int, int, int, int ); |
---|
| 58 | double ****dm4d( int, int, int, int, int, int, int, int ); |
---|
| 59 | void frm1d( float *, int ); |
---|
| 60 | void frm2d( float **, int, int, int ); |
---|
| 61 | void frm3d( float ***, int, int, int, int, int ); |
---|
| 62 | void fdm1d( double *, int ); |
---|
| 63 | void fdm2d( double **, int, int, int ); |
---|
| 64 | void fdm3d( double ***, int, int, int, int, int ); |
---|
| 65 | int *im1d( int, int ); |
---|
| 66 | int **im2d( int, int, int, int ); |
---|
| 67 | int ***im3d( int, int, int, int, int, int ); |
---|