1 | c----------------------------------------------------------------------- |
---|
2 | c------------------------------------------------------------------------ |
---|
3 | c INCLUDE 'dimradmars.h' |
---|
4 | |
---|
5 | c Declaration and initialisation or radiative transfer calculations |
---|
6 | c------------------------------------------------------------------------ |
---|
7 | c------------------------------------------------------------------------ |
---|
8 | |
---|
9 | c Splitting of horizontal grid (to reduce program size on workstation) |
---|
10 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
11 | c nflev: number of vertical layer |
---|
12 | c ndlon: nombre de points de la grille horizontale |
---|
13 | c NDLO2 et ndomainsz pour le decoupage de l'appel a la physique |
---|
14 | c ATTENTION: Il faut 1 < ndomainsz =< ngridmx |
---|
15 | |
---|
16 | INTEGER NFLEV,NDLON, ndomainsz, NDLO2 |
---|
17 | |
---|
18 | parameter (ndomainsz=ngridmx) |
---|
19 | c parameter (ndomainsz=(ngridmx-1)/20 + 1) |
---|
20 | c parameter (ndomainsz=(ngridmx-1)/5 + 1) |
---|
21 | |
---|
22 | parameter (NFLEV=nlayermx,NDLON=ndomainsz) ! avec decoupage |
---|
23 | parameter (NDLO2=NDLON) |
---|
24 | |
---|
25 | |
---|
26 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
27 | |
---|
28 | REAL longrefir,longrefvis |
---|
29 | REAL long1vis,long2vis,long3vis, long1ir,long2ir |
---|
30 | REAL long1co2,long2co2 |
---|
31 | REAL sunfr(2) |
---|
32 | integer nir, nuco2 |
---|
33 | INTEGER npademx,nabsmx,nt_pademx, NSUN |
---|
34 | |
---|
35 | |
---|
36 | c Definition of spectral intervals at thermal infrared wavelengths (LW) |
---|
37 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
38 | parameter (nir=4) ! Total number of thermal IR bands |
---|
39 | parameter (nuco2=2) ! number of bands in CO2 bands |
---|
40 | PARAMETER (long1ir=5.E-6 , long2ir=200.E-6) |
---|
41 | PARAMETER (long1co2=1.E+0 / 865.E+2 , long2co2=1.E+0 / 500.E+2) |
---|
42 | |
---|
43 | c Warning : the "nir" thermal IR bands are not ordered by wavelength: |
---|
44 | c iir=1 : central 15um CO2 bands \ |
---|
45 | c iir=2 : CO2 band wings [long1co2-long2co2] MINUS central band |
---|
46 | c iir=3 : 9 um band [long1ir - long1co2] |
---|
47 | c iir=4 : Far IR [long2co2 - long2ir] |
---|
48 | |
---|
49 | c Definition of spectral interval at solar wavelengths (SW) |
---|
50 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
51 | PARAMETER (NSUN=2) ! do not change that ! |
---|
52 | c Boundaries of spectral intervals (m) : |
---|
53 | PARAMETER (long1vis=0.1E-6 , long2vis=0.5E-6 , long3vis=5.E-6) |
---|
54 | c Fraction of solar energy in solar band #1 [long1vis-long2vis] |
---|
55 | DATA sunfr(1) / 0.274490 / |
---|
56 | c Fraction of solar energy in solar band #2 [long2vis-long3vis] |
---|
57 | DATA sunfr(2) / 0.725509 / |
---|
58 | |
---|
59 | c Reference wavelengths used to compute reference dust optical depth (m) |
---|
60 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
61 | PARAMETER (longrefir=9.E-6,longrefvis=0.67E-6) |
---|
62 | |
---|
63 | c Number of kind of tracer radiative properties |
---|
64 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
65 | c (ex: naerkind=1 if you use one dust mode without ice ...) |
---|
66 | c (ex: naerkind=2 if you use one dust mode and active ice ...) |
---|
67 | integer naerkind |
---|
68 | parameter (naerkind=1) |
---|
69 | |
---|
70 | c Various initialisation for LW radiative code |
---|
71 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
72 | c npademx : nombre de coef de pade |
---|
73 | c nabsmx : ? |
---|
74 | c nt_pademx : nombre d'intervalles de temperature pour pade |
---|
75 | |
---|
76 | PARAMETER (npademx=4,nabsmx=2,nt_pademx=19) |
---|