source: trunk/mesoscale/TESTS/newphys_tracers/callphys.def.mars11 @ 76

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

LMD_MM_MARS: remise a plat du traitement des traceurs terminee

[en particulier en vue de l'utilisation nouvelle physique]
--> fonctionnement concluant sur un cas Tharsis hydro 61x61x61
--> ne pas utiliser le cas FASTCASE trop instable avec traceurs
--> reste a tester l'effet d'especes radiativement actives

options dans namelist.input :
mars = 0 ---> pas de traceurs
mars = 1 ---> cycle de l'eau : water vapour + ice
mars = 2 ---> cycle poussieres : un dust bin
mars = 3 ---> cycle poussieres : dustq + dustn [NOUVELLE PHYS seulement]
mars = 11 ---> cycle de l'eau + poussieres [1+3] [NOUVELLE PHYS seulement]

NB: pour les deux derniers, reste un petit travail mineur sur qsurf

(voir dans module_lmd_driver.F)

il faut definir conjointement le callphys.def associe et relancer real.exe

(sinon on transporte juste dynamiquement des 'dummy' traceurs)

il n'est necessaire de recompiler que si le nombre total de traceurs change

Fichiers a mettre a jour si l'on ajoute une option "mars" dans le Registry


Attention suivant les inputs GCM, il faut peut etre egalement changer

  1. readmeteo.F90 dans PREP_MARS et 2. METGRID.TBL dans WPS/metgrid

M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
---> definition du scenario et de l'ordre des traceurs dans SCALAR

M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F
---> definition des interpolations verticales des champs de traceurs venant du GCM

M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F
---> definition du traitement a reserver aux bornes
---> --- dans les 4 scenarios precites, on passe aux bornes les champs du GCM

[y compris QH2O_ICE contrairement a ce qui etait par defaut precedemment]

---> --- les lignes relatives a un flux nul aux bornes sont laissees a un cas hypothetique mars>50

M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
---> definition de l'ordre correct des traceurs pour le passage a la physique
---> recuperation des tendances de la physique pour passage a la dynamique
---> NB: c'est dans ce module que sont presents des STOP si mars = 4-10 ou mars > 11

.... il faut donc modifier si l'on ajoute des options

M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/main/real_em.F
---> definition et calcul des champs a appliquer aux bornes
---> generalise desormais, il n'y a plus qu'a ajouter d'eventuelles nouvelles options 'mars'
---> .... des modifications sont necessaires si on passe plus de 4 traceurs aux bornes

M 75 mesoscale/LMD_MM_MARS/SIMU/runmeso
---> definition du bon nombre de traceurs dans la compilation de la physique puis l'execution

[l'option mars est lue par le script dans namelist.input]

Fichiers tests


A 0 mesoscale/TESTS/newphys_tracers/*
---> pour la nouvelle physique (ici seulement les fichiers def)
---> toutes les options precitees ont ete testes avec succes a l'execution [pas de crash]
---> ... plausibilite physique verifiee rapidement, PAS d'ANALYSE APPROFONDIE pour le moment

A 0 mesoscale/TESTS/LMD_MM_MARS_TESTCASE_water.tar.gz
---> pour l'ancienne physique (introduit precedemment mais n'avait pas ete synchronise)

M 75 mars/libf/phymars/meso_dustlift.F
NB: correction mineure, de facon a recuperer alpha_lift de initracer

File size: 3.8 KB
Line 
1##General options
2##~~~~~~~~~~~~~~~
3#Run with or without tracer transport ?
4tracer=.true.
5
6#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
7diurnal=.true.
8
9#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start"
10season = .true.   
11
12#write some more output on the screen ?
13lwrite = .false.   
14
15#Save statistics in file "stats.nc" ?
16callstats =.false.
17
18#Save EOF profiles in file "profiles" for Climate Database?
19calleofdump = .false.
20
21## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
22## ~~~~~~~~~~~~~
23#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
24#  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
25#  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
26#  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
27#  =26 Mars Year 26 from TES assimilation
28iaervar = 26
29# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
30tauvis=0.2
31# Dust vertical distribution:
32# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
33#  =2: Viking scenario; =3 MGS scenario)
34iddist  = 3
35# Dust top altitude (km). (Matters only if iddist=1)
36topdustref = 55.
37
38## Physical Parameterizations :
39## ~~~~~~~~~~~~~~~~~~~~~~~~~~
40# call radiative transfer ?
41callrad   = .true.
42# call NLTE radiative schemes ?   matters only if callrad=T
43callnlte = .false.
44# call CO2 NIR absorption ?   matters only if callrad=T
45callnirco2 = .true.
46# call turbulent vertical diffusion ?
47calldifv  = .true.
48# call convective adjustment ?
49calladj  = .true.
50# call CO2 condensation ?
51callcond =.true.
52# call thermal conduction in the soil ?
53callsoil = .true.
54# call Lott's gravity wave/subgrid topography scheme ?
55calllott  = .false.
56# Impose polar cap surface albedos as observed by TES?
57TESicealbedo = .true.
58
59## Radiative transfer options :
60## ~~~~~~~~~~~~~~~~~~~~~~~~~~
61# the rad.transfer is computed every "iradia" physical timestep
62iradia    = 6
63# Output of the exchange coefficient mattrix ? for diagnostic only
64callg2d  = .false.
65# Rayleigh scattering : (should be .false. for now)
66rayleigh  = .false.
67
68## Tracer (dust water, ice and/or chemical species) options (used if tracer=T):
69## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
70# DUST: Transported dust ? (if >0, use 'dustbin' dust bins)
71dustbin    = 2
72# DUST: Radiatively active dust ? (matters if dustbin>0)
73active  = .false.
74# DUST: use mass and number mixing ratios to predict dust size ?
75doubleq   = .true.
76# DUST: lifted by GCM surface winds ?
77lifting = .true.
78# DUST: lifted by dust devils ?
79callddevil = .false.
80# DUST: Scavenging by CO2 snowfall ?
81scavenging = .false.
82# DUST/WATERICE: Gravitationnal sedimentation ?
83sedimentation = .true.
84# WATERICE: Radiatively active transported atmospheric water ice ?
85activice   = .false.
86# WATER: Compute water cycle
87water = .true.
88# WATER: current permanent caps at both poles. True IS RECOMMENDED
89#        (with .true., North cap is a source of water and South pole
90#         is a cold trap)
91caps  = .false. #.true.
92# PHOTOCHEMISTRY: include chemical species
93photochem  = .false.
94
95## Thermospheric options (relevant if tracer=T) :
96##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
97# call thermosphere ?
98callthermos = .false.
99# WATER: included without cycle (only if water=.false.)
100thermoswater = .false.
101# call thermal conduction ? (only if callthermos=.true.)
102callconduct = .false.
103# call EUV heating ? (only if callthermos=.true.)
104calleuv=.false.
105# call molecular viscosity ? (only if callthermos=.true.)
106callmolvis = .false.
107# call molecular diffusion ? (only if callthermos=.true.)
108callmoldiff = .false.
109# call thermospheric photochemistry ? (only if callthermos=.true.)
110thermochem = .false.
111# date for solar flux calculation: (1985 < date < 2002)
112## (Solar min=1996.4 ave=1993.4 max=1990.6)
113solarcondate = 1993.4
114       
115
116
Note: See TracBrowser for help on using the repository browser.