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

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

LMD_MM_MARS avec nouvelle physique entierement interfacee avec LMDZ.MARS


M mars/libf/phymars/dimradmars.h
M mars/libf/phymars/callradite.F

simples changements pour faire tourner un cas test

  • dust conrath
  • un seul traceur

un petit domaine splitting de plus

A mesoscale/TESTS/LMD_MM_MARS_FASTCASEnewphys.tar.gz

cas test du mesochelle avec la nouvelle physique
essai concluant !

M mesoscale/NOTES.txt
M mesoscale/LMD_MM_MARS/prepare_ini
M chantiers/commit_importants.log

MAJ notes

M mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F

correctif pour le cas ou on tourne avec d'anciens diagfi avec 10 niveaux
dans le sol seulement. si tsoil < 20 K (disons) dans un niveau donne, alors
tsoil est regle a la temperature du niveau superieur.

D mesoscale/COMMON_GCM/gr_fi_dyn.F
D mesoscale/COMMON_GCM/diff_oldgcm_oldmeso_checked
D mesoscale/COMMON_GCM/meso_dustlift.txt
D mesoscale/COMMON_GCM/dimradmars.txt
D mesoscale/COMMON_GCM/nocompile.txt
D mesoscale/COMMON_GCM/meso_callkeys.txt
D mesoscale/COMMON_GCM/dimradmars.h

menage

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.