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

Last change on this file since 54 was 54, checked in by aslmd, 15 years ago

LMD_LES_MARS: Tests with new physics in unified frame
M 53 mesoscale/LMD_LES_MARS/modif_mars/module_model_constants.F
M 53 mesoscale/LMD_LES_MARS/modif_mars/module_first_rk_step_part1.F
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
Modifications made to cope with new soil scheme

M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/meso_physiq.F
Commented aeronomars in old physics so that the model can compile with lots of grid points


LMD_MM_MARS: Implement tracers with the new physics. Test with JBM runs with 5 tracers

including radiatively active dust and water ice

TODO: propagate reading "traceur.def" and few initialisations from

testphys1d.F into module_lmd_driver.F

M 53 mars/libf/phymars/dimradmars.h
M 53 mars/libf/phymars/callradite.F
Change it back to both dust and ice radiatively active

M 53 mesoscale/LMD_MM_MARS/makemeso
Allow to "fresh start" with the inclusion of "clean -a" in the script

M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/main/real_em.F
M 53 mesoscale/LMD_MM_MARS/SIMU/runmeso
Added the option mars==11 with 5 tracers: CO2 H2Ovap H2Oice DUST DUSTN
-- Note that the order matters for H2Ovap and H2Oice

LMD_MM_MARS: Create a folder to run LMD GCM with exact same physics as mesoscale

A 0 mesoscale/LMDZ.MARS/libf/phymars/physiq.F
All files in mesoscale/LMDZ.MARS/libf are links except this one
--> Because specific WRITEDIAGFI commands are set to output what is needed PREP_MARS
--> Something must be done so that qCO2, qdust, qdustN can be propagated too !!!

A 0 mesoscale/LMDZ.MARS/libf_gcm
A 0 mesoscale/LMDZ.MARS/in_lmdz_mars_newphys
Two important links that must not be broken

A 0 mesoscale/LMDZ.MARS/myGCM
This folder is for runs. Only links here -- except temp files which are not synchronized.

M 53 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/compile
A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/run.def
A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/traceur.def
M 53 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/callphys.def
A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/callphys.def.csttau
Those files are needed to run the GCM in order to prepare initial & boundary

conditions for the mesoscale... but basically it is the exact same GCM

--> compile is a convenient script

A 0 mesoscale/LMDZ.MARS/makegcm
A 0 mesoscale/LMDZ.MARS/makegcm_g95
A 0 mesoscale/LMDZ.MARS/create_make_gcm
This is supposed to be merged one day with files in trunk/mars/

A 0 mesoscale/LMDZ.MARS/myGCM/DEFS_JB
Necessary files to restart a GCM run corresponding to new water cycle in JB's thesis

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.