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