source: trunk/LMDZ.MARS/libf/phymars/dimradmars.h @ 146

Last change on this file since 146 was 132, checked in by aslmd, 14 years ago

mars: default naerkind=1 nsizemax=1 in dimradmars.h [so that GCM+meso+LES works out-of-the-box. this would change in the future when RAC physics would become the reference version]

File size: 3.5 KB
Line 
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)
Note: See TracBrowser for help on using the repository browser.