source: BOL/LMDZ_Setup_amaury/DEF/iso.def_iso @ 5182

Last change on this file since 5182 was 4615, checked in by asima, 18 months ago

Initial import of LMDZ_Setup environment
(former TEST_PROD which was distributed as tutorial_prod.tar).

See also related doc LMDZ_Setup_HowTo (former tutorial_prod_FAQ) :
https://docs.google.com/document/d/1OLZG6e-86NiXuv5-aALxKIh-QPkp4BdCwWtiBFot-6c

This initial import corresponds to TEST_PROD from
https://lmdz.lmd.jussieu.fr/pub/Training/Archive/tutorial_prod.tar_v20230707.

Previous versions of tutorial_prod.tar are available at
https://lmdz.lmd.jussieu.fr/pub/Training/Archive

File size: 7.3 KB
Line 
1# paramètres isotopiques principaux
2
3ok_iso_verif=n
4
5# lambda: supersaturation: valeurs typiques: 0.002. Range possible: 0 (pas de sursat) à 0.004 (forte sursat)
6lambda=0.002
7
8# thumxt1: gère le fractionnement cinétique lors de la réévaporation des gouttes. Valeur typique: 0.9. Range possible: 0 (rh au contact des gouttes = rh autour des gouttes, fort fractionnement cinétique) à 1 (rh=1 au contact des gouttes, pas de fractionnement cinétique)
9thumxt1=0.9
10
11# h_land_ice: hauteur caractéristique de la glace de terre: la glace de terre à t+1 est égale à h_land_ice de glace de terre à t + la precip neigeuse en t.
12#Valeur typique: 20m. Range possible: proche de 0 (-> la glace de terre prend la composition de la dernière précip neigeuse) à très grand (-> la compo de la glace de terre est constante)
13h_land_ice=20.
14
15# P_veg: fraction de l evaporation continentale provenant de la transpiration sans fractionnement de la végétation, au lieu de l évaporation avec fractionnement de sol nu ou d eau libre. Valeur standard: 1 (comme autres GCMs). Range possible: 0 à 1.
16P_veg=1.0
17
18# bidouille_anti_divergence: si T: on rappel régulièrement l isotope H2O16 vers l eau normale. Si F: pas de rappel. En général, on met T pour être sur, mais on a verifie que ça a peut d importance
19bidouille_anti_divergence=y
20
21# essai_convergence: si T, on fait tout exactement comme pour l eau normale, pour que les 2 versions convergent. Si F, on modifie un peu: ex: initialisations supplémentaires de variables, on zappe l homogénéisation sous le nuage dans cv3_yield...
22essai_convergence=n
23
24# initialisation_iso: la façon dont les isos sont initialisés.
25#   - Lors du lancement de gcm:
26#       Pour la physique:
27#       0: lecture d un fichier, souvent de start (phyiso_etat0_fichier.F)
28#       >0: initialisation des isotopes académique, donnée dans phyiso_etat0_dur.F
29#       Pour la dynamique:
30#       0: lecture d un fichier (dyniso_etat0_fichier.F)
31#       >0: initialisation académique:
32#               1: isos à 0 sauf H2016=eau normale
33#               2: distill de Rayleigh
34#               3: tout est uniforme au SMOW
35#       En général, on met initialisation_iso=0
36#   - Lors du lancement de create_etat0_limit:
37#       Pour la physique: initialisation en dure, toujours
38#       Pour la dynamique: initialisation académique toujours:
39#               1: isos à 0 sauf H2016=eau normale
40#               2 ou 0: distill de Rayleigh
41#               3: tout est uniforme au SMOW
42#       En général, on laisse à initialisation_iso=0
43initialisation_iso=2
44
45# initialisation_isotrac: comment initialiser les traceurs d'eau?
46# 1: idéalisé: on met toute l'eau dans le tag izone_init, 0 ailleurs
47# 0: en lisant dans un fichier. Possible seulement si initialisation_iso=0
48initialisation_isotrac=1
49
50# modif_sst: 0 par defaut, 1 si on veut modifier la sst de manière uniforme
51# si 2, on impose un profil de changement de SST, avec 0 entre sstlatcrit et - sstlatcrit en latitude. Ailleurs: +deltaTtestpoles (rlat(ii)**2-sstlatcrit**2) /(90.0**2-sstlatcrit**2)
52# si 2, on impose changement de SST linéaire par morceau: 0 entre sstlatcrit et - sstlatcrit, deltaTtestpoles poleward of sstlatcrit+dsstlatcrit, et linéaire entre les 2.
53modif_sst=0
54
55# deltaTtest: sensibilité à une variation uniforme de la SST. Utilise seulement si modif_sst=1
56deltaTtest=0.0
57
58deltaTtestpoles=-15.0
59sstlatcrit=30.0
60dsstlatcrit=20.0
61
62# deltaO18_oce: deltaO18 de l'océan: ex: +1.1 au LGM, 0 actuellement. On suppose que le dexcess et O17excess océan sont nuls.
63deltaO18_oce=0.0
64
65# deltaP_BL:  epaisseur de la couche limite avec laquelle l'eau du sol tend à se rééquilibrer. Par défaut: 10mb (optim dans ORCHIDEE_LMDZ)
66deltaP_BL=10.0
67
68# ruissellement_pluie: qu'est-ce qui ruisselle? si 1: c'est la pluie qui ruisselle. elle ne s'infiltre. elle ne s'infiltre donc jamais dans un sol saturé. si 0: c'est le sol qui ruisselle. La pluie s'inglitre donc dans le sol saturé.
69ruissellement_pluie=1
70
71# alphak_stewart: si 1: alphak=(D/Diso)^nsol. si 0: alphak=1/(1-kcin(vsurf))
72alphak_stewart=1
73
74# tdifexp_sol: tdifexp_sol est l'exposant de D/Diso. Il paramétrise la turbulence. D'abitude, il est de 0.58. Mais d'après Mathieu et Bariac, il est entre 0.67 et 1: 0.67 pour les sols secs et 1 pour les sols saturés.
75tdifexp_sol=0.67
76
77# calendrier_guide: si 1, on fait en sorte que l'année puisse avoir 365 ou 366 jours. Si 0, l'année à 360 jours (c'est LMDZ par défaut).
78# mettre 1 quand on guide, 0 sinon.
79calendrier_guide=1
80
81# cste_surf_cond: si 2, la témprature et la rh de surface sont constantes: rhs=0.7 et Ts=15°=288K. si 1: la rh de surface est conste égale à 0.7. Par défaut, cste_surf_cond=0: tout est explicit.
82cste_surf_cond=0
83
84# lev_histiso
85# niveaux de sorties pour les isotopes et aussi sans les isotopes.
86# Pas encore utilisé
87lev_histiso=1
88
89# lev_histiso_diag
90# ne s applique que si DIAGISO. 1 par defaut. si 0, seul xtvp est outputé
91lev_histiso_diag=1
92
93# deltaD aberrant: si on le dépasse, on plante
94deltalim=1000000.0
95
96# deltaD aberrant pour les traceurs
97deltalimtrac=1000000.0
98dexcess_max=2000
99dexcess_min=-1000
100
101# est-ce qu'on verifie l'O17?
102O17_verif=y
103o17excess_bas=-1000.0
104o17excess_haut=900.0
105nlevmaxO17=18
106
107# pour que les continents évaporent une compositon constante, pour voir si les cycles de dexcess en Amazonie sont liés à la compo de l''évap sur continent
108# par défaut: evap_cont_cste=0
109# si evap_cont_cste=1, préciser deltaO18 et dexcess de l''évap
110evap_cont_cste=0
111deltaO18_evap_cont=-5.0
112d_evap_cont=14.0
113
114# quand on utilise tag 17, on taggue la vapeur résiduelle quand la fraction condensée est supérieure à seuil_tag_tmin, c'est à dire quand l'évènement de condensation affecte suffisemment la compo de la vapeur résiduelle. Si seuil_tag_tmin=0, le moindre évènement de condensation même très faibel agit sur les tags, et les proportions de sfc, condensat et rev sont presques nulles.
115# par defaut: seuil_tag_tmin=0.01: 1% de condensation <-> 0.1 permil d'effet sur la vapeur résiduelle.
116seuil_tag_tmin=0.01
117
118# idem pour condensation dans LS: par defaut égal à seuil_tag_tmin
119seuil_tag_tmin_ls=0.01
120
121# option_seuil_tag_tmin: si 1, on recolorise la vapeur résiduelle dès que cond/qt>seuil_tag_tmin. Ca taggue bien la vap résiduelle mais le problème est que plus on humidifie par la microphysique, plus le seuil est facilement attent -> plus on retaggue -> au final, on n'a pas plus de tag mcrophysique alors que c'est la sourec d'humdification.
122# dans l'option 2: on taggue si (cond-qmicro)/((qt-qmicro)>seuil
123option_seuil_tag_tmin=2
124
125# desactivier post-condensational exchanges comme dans Field, Jones et Brown à réviewer
126no_pce=0
127
128# latitude de taggage des extra-tropiques quand option de taggage n° 20
129# 35.0 par défaut
130lim_tag20=35.0
131
132# quand tag9, recolorise t'on que le condensat convectif (option 2) ou aussi le condensat starti (otion 1)
133option_cond=2
134
135# quand tag 22, on recolorise la vapeur dans les zones ou precip>lim_precip_tag22
136# 10.0 par défaut
137lim_precip_tag22=10.0
138
139# quand adv 14, saturation limiter tel que on advec min (q,A_satlim*qs)
140# par defaut, A_satlim=1: c'est le code original de Francis Codron
141# ce param n'a d'effet que quand adv 14.
142A_satlim=0.8
143
144# si ok_restrict_A_satlim=1, alors on n'active adv 14 que dans les extra-tropiques
145# si ok_restrict_A_satlim=2, alors on n'active adv 14 que dans les tropiques
146# si ok_restrict_A_satlim=3, alors on n'active adv 14 nulle part
147ok_restrict_A_satlim=0
148
149slope_limiter=2.0
150# defaut: 2; pour adv 10: 0
151
152#on multiplie par 3 ratqs aux extraropiques si modif_ratqs=1.
153modif_ratqs=0
154
Note: See TracBrowser for help on using the repository browser.