| 1 | ## $Id$ |
|---|
| 2 | # |
|---|
| 3 | ### type_ocean = force / slab /couple |
|---|
| 4 | type_ocean=couple |
|---|
| 5 | ### version_ocean = nemo / opa8 |
|---|
| 6 | version_ocean=nemo |
|---|
| 7 | ### VEGET= y si ORCHIDEE, =n si bucket |
|---|
| 8 | VEGET=y |
|---|
| 9 | # |
|---|
| 10 | # PARAMETRES ANCIENNEMENT DANS gcm.def |
|---|
| 11 | ## Cycle diurne ou non |
|---|
| 12 | cycle_diurne=y |
|---|
| 13 | ## Soil Model ou non |
|---|
| 14 | soil_model=y |
|---|
| 15 | ## Choix ou non de New oliq |
|---|
| 16 | new_oliq=y |
|---|
| 17 | ## Activation ou non de la parametrisation de Hines pour la strato |
|---|
| 18 | ok_hines=n |
|---|
| 19 | ## Orodr ou non pour l orographie |
|---|
| 20 | ok_orodr=y |
|---|
| 21 | ## Orolf ou non pour l orographie |
|---|
| 22 | ok_orolf=y |
|---|
| 23 | ## Si = .T. , lecture du fichier limit avec la bonne annee |
|---|
| 24 | ok_limitvrai=n |
|---|
| 25 | ## Nombre d'appels des routines de rayonnements ( par jour) |
|---|
| 26 | nbapp_rad=12 |
|---|
| 27 | ## Flag pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect |
|---|
| 28 | iflag_con=30 |
|---|
| 29 | ## Facteur multiplication des precip convectives dans KE |
|---|
| 30 | cvl_corr=1.0 |
|---|
| 31 | ## Facteur additif pour l'albedo |
|---|
| 32 | pmagic=0.008 |
|---|
| 33 | # |
|---|
| 34 | # |
|---|
| 35 | # |
|---|
| 36 | # Parametres fichiers de sortie |
|---|
| 37 | # |
|---|
| 38 | ### type_run = type run par rapport aux fichiers et variables de sortie |
|---|
| 39 | # - type_run = CLIM/ENSP (=1) |
|---|
| 40 | # - type_run = AMIP/CFMI (=2) |
|---|
| 41 | type_run=CLIM |
|---|
| 42 | ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc |
|---|
| 43 | OK_journe=n |
|---|
| 44 | ### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc |
|---|
| 45 | OK_mensuel=y |
|---|
| 46 | ### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la physique) |
|---|
| 47 | OK_instan=n |
|---|
| 48 | ### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc |
|---|
| 49 | ok_hf=n |
|---|
| 50 | # |
|---|
| 51 | # Parametres niveau de sorties differents fichiers |
|---|
| 52 | # |
|---|
| 53 | ### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" |
|---|
| 54 | # - lev_histhf=0 => pas de sorties histhf.nc |
|---|
| 55 | # - lev_histhf=2 => defaut |
|---|
| 56 | # - lev_histhf=3 => variables sur niveaux standards |
|---|
| 57 | # - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc |
|---|
| 58 | lev_histhf=2 |
|---|
| 59 | ### lev_histday=0-5, niveau de sortie fichier "histday.nc" |
|---|
| 60 | # - lev_histday=0 => pas de sorties lev_histday.nc |
|---|
| 61 | # - lev_histday=2 => defaut |
|---|
| 62 | # - lev_histday=3 => + champs 3D => F. Lott |
|---|
| 63 | # - lev_histday=4 => + champs sous-surfaces |
|---|
| 64 | # - lev_histday=5 => + champs F. Aires |
|---|
| 65 | lev_histday=2 |
|---|
| 66 | ### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" |
|---|
| 67 | # - lev_histmth=0 => pas de sorties lev_histmth.nc |
|---|
| 68 | # - lev_histmth=2 => defaut |
|---|
| 69 | # - lev_histmth=3 => albedo, rugosite sous-surfaces |
|---|
| 70 | # - lev_histmth=4 => champs tendances 3d |
|---|
| 71 | lev_histmth=2 |
|---|
| 72 | ### ecrit_hf = frequence ecriture fichier histhf.nc en jours |
|---|
| 73 | ecrit_hf=0.250 |
|---|
| 74 | ### ecrit_day = frequence ecriture fichier histday.nc en jours |
|---|
| 75 | ecrit_day=1. |
|---|
| 76 | ### ecrit_mth = frequence ecriture fichier histmth.nc en jours |
|---|
| 77 | ecrit_mth=_ecrit-mth_ |
|---|
| 78 | ### freqin_isccp = frequence input en secondes du simulateur ISCCP |
|---|
| 79 | freq_ISCCP=10800. |
|---|
| 80 | ### freqout_isccp = frequence output en jours du simulateur ISCCP |
|---|
| 81 | ecrit_ISCCP=_ecrit-ISCCP_ |
|---|
| 82 | ### niveau du diagnostique de conservation d energie |
|---|
| 83 | if_ebil=0 |
|---|
| 84 | # |
|---|
| 85 | # parametres KE |
|---|
| 86 | # |
|---|
| 87 | ### epmax = Efficacite precipitation maximale |
|---|
| 88 | epmax = .99 |
|---|
| 89 | ### ok_adj_ema = ?? pas utilise |
|---|
| 90 | ok_adj_ema = n |
|---|
| 91 | ### iflag_clw Flag calcul eau liquide |
|---|
| 92 | # - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l) |
|---|
| 93 | # - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l) |
|---|
| 94 | # - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip |
|---|
| 95 | iflag_clw = 1 |
|---|
| 96 | # |
|---|
| 97 | # parametres nuages |
|---|
| 98 | # |
|---|
| 99 | ### cld_lc_lsc contenu en eau liquide des nuages large-scale (fisrtilp) |
|---|
| 100 | cld_lc_lsc = 2.6e-4 |
|---|
| 101 | ### cld_lc_con contenu en eau liquide des nuages convectifs (fisrtilp) |
|---|
| 102 | cld_lc_con = 2.6e-4 |
|---|
| 103 | ### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) |
|---|
| 104 | cld_tau_lsc = 3600. |
|---|
| 105 | ### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) |
|---|
| 106 | cld_tau_con = 3600. |
|---|
| 107 | ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) |
|---|
| 108 | ffallv_lsc = 1. |
|---|
| 109 | ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) |
|---|
| 110 | ffallv_con = 1. |
|---|
| 111 | ### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?) |
|---|
| 112 | coef_eva = 2.e-5 |
|---|
| 113 | ### reevap_ice reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) |
|---|
| 114 | reevap_ice = y |
|---|
| 115 | ### iflag_cldcon flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) |
|---|
| 116 | # - iflag_cldcon<=-1 diag. rain_Tiedtke |
|---|
| 117 | # - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k) |
|---|
| 118 | # - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc) |
|---|
| 119 | # - iflag_cldcon=3, ratqs=ratqss |
|---|
| 120 | iflag_cldcon = 3 |
|---|
| 121 | ### iflag_pdf : flag calcul distribution sous-maille de l eau et des nuages |
|---|
| 122 | # - iflag_pdf=0, version ratqs, |
|---|
| 123 | # - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir |
|---|
| 124 | # - des PDFS de Sandrine Bony |
|---|
| 125 | iflag_pdf = 1 |
|---|
| 126 | ### fact_cldcon constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) |
|---|
| 127 | fact_cldcon = 1. |
|---|
| 128 | ### facttemps= facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3) |
|---|
| 129 | facttemps = 1.e-4 |
|---|
| 130 | ## ok_newmicro =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) |
|---|
| 131 | ok_newmicro = y |
|---|
| 132 | ### iflag_ratqs=0 correspond a la version IPCC AR4 |
|---|
| 133 | iflag_ratqs=0 |
|---|
| 134 | ### ratqsbas ratqs en bas si iflag_cldcon=1 |
|---|
| 135 | ratqsbas = 0.005 |
|---|
| 136 | ### ratqshaut ratqs en haut pour ratqss "stables" |
|---|
| 137 | ratqshaut = 0.33 |
|---|
| 138 | ### rad_froid rayon cristaux des nuages de glace "froids" |
|---|
| 139 | rad_froid = 35 |
|---|
| 140 | ### rad_chau1 rayon goutelettes d eau chauds", en haut: k=4-klev |
|---|
| 141 | rad_chau1=12 |
|---|
| 142 | ### rad_chau2 rayon goutelettes d eau chauds", en bas: k=1-3 |
|---|
| 143 | rad_chau2=11 |
|---|
| 144 | # |
|---|
| 145 | # Coefficient et parametres sur les drags |
|---|
| 146 | # |
|---|
| 147 | f_cdrag_ter=1. |
|---|
| 148 | f_cdrag_stable=1. |
|---|
| 149 | f_cdrag_oce=0.8 |
|---|
| 150 | f_rugoro=0. |
|---|
| 151 | ### ksta_ter coef.diffusion minimale sur terre/sic/lic |
|---|
| 152 | ksta_ter=1.e-7 |
|---|
| 153 | ### cdmmax = cdrag maximum pour le moment |
|---|
| 154 | cdmmax = 2.5E-3 |
|---|
| 155 | ### cdhmax = cdrag maximum pour l energie |
|---|
| 156 | cdhmax = 2.0E-3 |
|---|
| 157 | # |
|---|
| 158 | # Parametres "orbitaux/ ere geologique" |
|---|
| 159 | # |
|---|
| 160 | ### R_ecc = Excentricite |
|---|
| 161 | R_ecc = 0.016715 |
|---|
| 162 | ### R_peri = Equinoxe |
|---|
| 163 | R_peri = 102.7 |
|---|
| 164 | ### R_incl = Inclinaison |
|---|
| 165 | R_incl = 23.441 |
|---|
| 166 | ### solaire = Constante solaire |
|---|
| 167 | solaire = 1361. |
|---|
| 168 | # |
|---|
| 169 | # Taux gaz a effet de serre |
|---|
| 170 | # |
|---|
| 171 | ### co2_ppm = taux CO2 en ppm |
|---|
| 172 | co2_ppm = 367. |
|---|
| 173 | ### CH4_ppb = taux CH4 en ppb |
|---|
| 174 | CH4_ppb = 1760. |
|---|
| 175 | ### N2O_ppb = taux N2O en ppb |
|---|
| 176 | N2O_ppb = 316. |
|---|
| 177 | ### CFC11_ppt = taux CFC11 en ppt |
|---|
| 178 | CFC11_ppt = 741.2 |
|---|
| 179 | ### CFC12_ppt = taux CFC12 en ppt |
|---|
| 180 | CFC12_ppt = 535. |
|---|
| 181 | # |
|---|
| 182 | # Parametres effets directs/indirects des "aerosols" |
|---|
| 183 | # |
|---|
| 184 | ### ok_ade=y/n flag Aerosol direct effect |
|---|
| 185 | ok_ade=y |
|---|
| 186 | ### ok_aie=y/n flag Aerosol indirect effect |
|---|
| 187 | ok_aie=y |
|---|
| 188 | ### aer_type = Aerosol variation type : actuel / preind / scenario |
|---|
| 189 | aer_type=actuel |
|---|
| 190 | ### type of coupled aerosol =1 (default) =2 => bc only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol |
|---|
| 191 | flag_aerosol=1 |
|---|
| 192 | ### bl95_b0 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
|---|
| 193 | bl95_b0=1.7 |
|---|
| 194 | ### bl95_b1 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
|---|
| 195 | bl95_b1=0.2 |
|---|
| 196 | # |
|---|
| 197 | # Parametres simulateur ISCCP |
|---|
| 198 | # |
|---|
| 199 | ### ok_isccp=y/n flag simulateur ISCCP |
|---|
| 200 | ok_isccp=n |
|---|
| 201 | ### top_height = flag choix calcul nuages par le simulateur en utilisant |
|---|
| 202 | # - les donnees IR et/ou VIS et l algorithme ISCCP-D1 |
|---|
| 203 | # - top_height = 1 -> algo IR-VIS |
|---|
| 204 | # - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" |
|---|
| 205 | # - top_height = 3 -> algo IR |
|---|
| 206 | top_height = 1 |
|---|
| 207 | ### overlap = Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP |
|---|
| 208 | # - overlap=1 Max overlap |
|---|
| 209 | # - overlap=2 Random overlap |
|---|
| 210 | # - overlap=3 Max/Random overlap |
|---|
| 211 | overlap = 3 |
|---|
