source: trunk/mars/libf/phymars/dimradmars.h @ 119

Last change on this file since 119 was 118, checked in by aslmd, 14 years ago

LMD_MM_MARS: ajout mars_lmd_new_storm avec routines modifiees par J. Faure (stage M2). modification makemeso et runmeso avec ajout d une option -s qui permet de preciser un scenario base sur une physique legerement differente comme par exemple storm. ceci cree un nouveau dossier de compilation dont le nom commence par le nom du scenario.

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
32! Reference wavelengths used to compute reference optical depth (m)
33! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
34
35      REAL longrefir(naerkind),longrefvis(naerkind)
36
37      REAL long1vis,long2vis,long3vis, long1ir,long2ir
38      REAL long1co2,long2co2
39      REAL sunfr(2)
40      integer nir, nuco2
41      INTEGER npademx,nabsmx,nt_pademx, NSUN
42
43
44! Definition of spectral intervals at thermal infrared wavelengths (LW)
45! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
46      parameter (nir=4) ! Total number of thermal IR bands
47      parameter (nuco2=2) ! number of bands in CO2 bands
48      PARAMETER (long1ir=5.E-6 , long2ir=200.E-6)
49      PARAMETER (long1co2=1.E+0 / 865.E+2 , long2co2=1.E+0 / 500.E+2)
50
51!  Warning : the "nir" thermal IR bands are not ordered by wavelength:
52!      iir=1 : central 15um CO2 bands \   
53!      iir=2 : CO2 band wings    [long1co2-long2co2] MINUS central band
54!      iir=3 : 9 um band [long1ir - long1co2]
55!      iir=4 : Far IR    [long2co2 - long2ir]
56   
57!  Definition of spectral interval at solar wavelengths (SW)
58!  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
59      PARAMETER (NSUN=2)   ! do not change that !
60!  Boundaries of spectral intervals (m) : 
61      PARAMETER (long1vis=0.1E-6 , long2vis=0.5E-6 , long3vis=5.E-6)
62!  Fraction of solar energy in solar band #1 [long1vis-long2vis]
63      DATA sunfr(1) / 0.274490 / 
64!  Fraction of solar energy in solar band #2 [long2vis-long3vis]
65      DATA sunfr(2) / 0.725509 /
66
67! Maximum number of grain size classes
68! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
69! This parameter has to be set to the maximum number of particle
70!   sizes contained in the optical parameter database; For example,
71!   if only one grain size is used to describe dust, and 30 are used
72!   to describe water-ice crystals in the visible and 15 in the IR,
73!   nsizemax has to be set to 30.
74! If only one grain size is considered for all the aerosols, set
75!   this parameter to 1 and convolution will be turned off during
76!   the radiative calculations.
77
78      INTEGER, PARAMETER :: nsizemax = 60
79!      INTEGER, PARAMETER :: nsizemax = 1
80
81! Various initialisation for LW radiative code
82! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
83! npademx : nombre de coef de pade
84! nabsmx : ?
85! nt_pademx : nombre d'intervalles de temperature pour pade
86
87      PARAMETER (npademx=4,nabsmx=2,nt_pademx=19)
Note: See TracBrowser for help on using the repository browser.