[3919] | 1 | # paramÚtres isotopiques principaux |
---|
| 2 | |
---|
| 3 | # lambda: supersaturation: valeurs typiques: 0.002. Range possible: 0 (pas de sursat) Ã 0.004 (forte sursat) |
---|
| 4 | lambda=0.004 |
---|
| 5 | |
---|
| 6 | # 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) |
---|
| 7 | thumxt1=0.9 |
---|
| 8 | |
---|
| 9 | # nombre d'iterations pour calcul fractionnement lors revap precip: |
---|
| 10 | #default: 20. |
---|
| 11 | # pour O17: 40. |
---|
| 12 | ntot= 40 |
---|
| 13 | |
---|
| 14 | # 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. |
---|
| 15 | #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) |
---|
| 16 | h_land_ice=20. |
---|
| 17 | |
---|
| 18 | # 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. |
---|
| 19 | P_veg=1.0 |
---|
| 20 | |
---|
| 21 | # 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 |
---|
| 22 | bidouille_anti_divergence=y |
---|
| 23 | |
---|
| 24 | # 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... |
---|
| 25 | essai_convergence=n |
---|
| 26 | |
---|
| 27 | # initialisation_iso: la façon dont les isos sont initialisés. |
---|
| 28 | # - Lors du lancement de gcm: |
---|
| 29 | # Pour la physique: |
---|
| 30 | # 0: lecture d un fichier, souvent de start (phyiso_etat0_fichier.F) |
---|
| 31 | # >0: initialisation des isotopes académique, donnée dans phyiso_etat0_dur.F |
---|
| 32 | # Pour la dynamique: |
---|
| 33 | # 0: lecture d un fichier (dyniso_etat0_fichier.F) |
---|
| 34 | # >0: initialisation académique: |
---|
| 35 | # 1: isos à 0 sauf H2016=eau normale |
---|
| 36 | # 2: distill de Rayleigh |
---|
| 37 | # 3: tout est uniforme au SMOW |
---|
| 38 | # En général, on met initialisation_iso=0 pour utiliser les fichiers starts |
---|
| 39 | # - Lors du lancement de create_etat0_limit: |
---|
| 40 | # Pour la physique: initialisation en dure, toujours |
---|
| 41 | # Pour la dynamique: initialisation académique toujours: |
---|
| 42 | # 1: isos à 0 sauf H2016=eau normale |
---|
| 43 | # 2 ou 0: distill de Rayleigh |
---|
| 44 | # 3: tout est uniforme au SMOW |
---|
| 45 | # En général, on laisse à initialisation_iso=0 pour initialiser avec distill Rayleigh |
---|
| 46 | initialisation_iso=2 |
---|
| 47 | |
---|
| 48 | # modif_sst: 0 par defaut, 1 si on veut modifier la sst de maniÚre uniforme |
---|
| 49 | # 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) |
---|
| 50 | modif_sst=0 |
---|
| 51 | |
---|
| 52 | # deltaTtest: sensibilité à une variation uniforme de la SST. Utilse seulement si modif_sst=1 |
---|
| 53 | deltaTtest=0.0 |
---|
| 54 | |
---|
| 55 | deltaTtestpoles=-15.0 |
---|
| 56 | sstlatcrit=30.0 |
---|
| 57 | |
---|
| 58 | # 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. |
---|
| 59 | deltaO18_oce=0.0 |
---|
| 60 | |
---|
| 61 | # 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) |
---|
| 62 | deltaP_BL=10.0 |
---|
| 63 | |
---|
| 64 | # 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é. |
---|
| 65 | ruissellement_pluie=1 |
---|
| 66 | |
---|
| 67 | # alphak_stewart: si 1: alphak=(D/Diso)^nsol. si 0: alphak=1/(1-kcin(vsurf)) |
---|
| 68 | alphak_stewart=1 |
---|
| 69 | |
---|
| 70 | # 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. |
---|
| 71 | tdifexp_sol=0.67 |
---|
| 72 | |
---|
| 73 | # 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). |
---|
| 74 | # mettre 1 quand on guide, 0 sinon. |
---|
| 75 | calendrier_guide=1 |
---|
| 76 | |
---|
| 77 | # 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 |
---|
| 78 | cste_surf_cond=0 |
---|
| 79 | |
---|
| 80 | # lev_histiso |
---|
| 81 | niveaux de sorties pour les isotopes et aussi sans les isotopes. |
---|
| 82 | #Pas encore utilisé |
---|
| 83 | lev_histiso=1 |
---|
| 84 | |
---|
| 85 | # 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 |
---|
| 86 | # par défaut: evap_cont_cste=0 |
---|
| 87 | # si evap_cont_cste=1, préciser deltaO18 et dexcess de l''évap |
---|
| 88 | evap_cont_cste=0 |
---|
| 89 | deltaO18_evap_cont=-5.0 |
---|
| 90 | d_evap_cont=14.0 |
---|
| 91 | |
---|
| 92 | lim_tag20=35.0 |
---|
| 93 | |
---|
| 94 | # quand adv 14, saturation limiter tel que on advec min (q,A_satlim*qs) |
---|
| 95 | # par defaut, A_satlim=1: c'est le code original de Francis Codron |
---|
| 96 | # ce param n'a d'effet que quand adv 14. |
---|
| 97 | # si ok_restrict_A_satlim=1, alors on n'active adv 14 que dans les extra-tropiques |
---|
| 98 | A_satlim=0.8 |
---|
| 99 | ok_restrict_A_satlim=0 |
---|
| 100 | |
---|
| 101 | # 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 |
---|
| 102 | # ceci n'a un impact que si adv 14 dans traceurs.def |
---|
| 103 | # par défaut, samepoles=0: le schéma 14 std |
---|
| 104 | samepoles=1 |
---|
| 105 | |
---|
| 106 | #on multiplie par 3 ratqs aux extraropiques si modif_ratqs=1. |
---|
| 107 | modif_ratqs=0 |
---|
| 108 | |
---|
| 109 | # deltaD aberrant |
---|
| 110 | deltalim=2000.0 |
---|
| 111 | deltaDmin=-990.0 |
---|
| 112 | cond_temp_env = y |
---|
| 113 | ok_iso_verif=n |
---|