# paramètres isotopiques principaux

# lambda: supersaturation: valeurs typiques: 0.002. Range possible: 0 (pas de sursat) à 0.004 (forte sursat)
lambda=0.004

# 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)
thumxt1=0.9

# nombre d'iterations pour calcul fractionnement lors revap precip: 
#default: 20.
# pour O17: 40.
ntot=  40

# 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.
#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)
h_land_ice=20.

# 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.
P_veg=1.0

# 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
bidouille_anti_divergence=y

# 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...
essai_convergence=n

# initialisation_iso: la façon dont les isos sont initialisés.
#   - Lors du lancement de gcm:
#	Pour la physique:
#	0: lecture d un fichier, souvent de start (phyiso_etat0_fichier.F)
#	>0: initialisation des isotopes académique, donnée dans phyiso_etat0_dur.F
#	Pour la dynamique:
#	0: lecture d un fichier (dyniso_etat0_fichier.F)
#	>0: initialisation académique:
#		1: isos à 0 sauf H2016=eau normale
#		2: distill de Rayleigh
#		3: tout est uniforme au SMOW
#	En général, on met initialisation_iso=0 pour utiliser les fichiers starts
#   - Lors du lancement de create_etat0_limit:
#	Pour la physique: initialisation en dure, toujours
#	Pour la dynamique: initialisation académique toujours:
#		1: isos à 0 sauf H2016=eau normale
#		2 ou 0: distill de Rayleigh
#		3: tout est uniforme au SMOW
#	En général, on laisse à initialisation_iso=0 pour initialiser avec distill Rayleigh
initialisation_iso=2

# modif_sst: 0 par defaut, 1 si on veut modifier la sst de manière uniforme
# 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)
modif_sst=0

# deltaTtest: sensibilité à une variation uniforme de la SST. Utilse seulement si modif_sst=1
deltaTtest=0.0

deltaTtestpoles=-15.0
sstlatcrit=30.0

# 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.
deltaO18_oce=0.0

# 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)
deltaP_BL=10.0 

# 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é.
ruissellement_pluie=1

# alphak_stewart: si 1: alphak=(D/Diso)^nsol. si 0: alphak=1/(1-kcin(vsurf))
alphak_stewart=1

# 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.
tdifexp_sol=0.67

# 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).
# mettre 1 quand on guide, 0 sinon.
calendrier_guide=1

# 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=2
cste_surf_cond=0

# lev_histiso
niveaux de sorties pour les isotopes et aussi sans les isotopes.
#Pas encore utilisé
lev_histiso=1

# pour que les continents évaporent une compositon constante, pour voir si les cycles de dexce        ss en Amazonie sont liés à la compo de l''évap sur continent
# par défaut: evap_cont_cste=0
# si evap_cont_cste=1, préciser deltaO18 et dexcess de l''évap
evap_cont_cste=0
deltaO18_evap_cont=-5.0
d_evap_cont=14.0

lim_tag20=35.0

# quand adv 14, saturation limiter tel que on advec min (q,A_satlim*qs)
# par defaut, A_satlim=1: c'est le code original de Francis Codron
# ce param n'a d'effet que quand adv 14.
# si ok_restrict_A_satlim=1, alors on n'active adv 14 que dans les extra-tropiques
A_satlim=0.8
ok_restrict_A_satlim=0

# quand adv 14, le traitement des poles est légèrement différent -> pour isoler la différence entre adv 10 et adv 14, on prend le même traitement des poles que pour adv 10
# ceci n'a un impact que si adv 14 dans traceurs.def
# par défaut, samepoles=0: le schéma 14 std
samepoles=1

#on multiplie par 3 ratqs aux extraropiques si modif_ratqs=1.
modif_ratqs=0

# deltaD aberrant
deltalim=2000.0
deltaDmin=-990.0
cond_temp_env = y
ok_iso_verif=n