[57] | 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) |
---|