[650] | 1 | ## $Id$ |
---|
| 2 | # |
---|
[686] | 3 | # PARAMETRES ANCIENNEMENT DANS gcm.def |
---|
| 4 | ## Cycle diurne ou non |
---|
| 5 | cycle_diurne=y |
---|
| 6 | ## Soil Model ou non |
---|
| 7 | soil_model=y |
---|
| 8 | ## Choix ou non de New oliq |
---|
| 9 | new_oliq=y |
---|
| 10 | ## Activation ou non de la parametrisation de Hines pour la strato |
---|
| 11 | ok_hines=y |
---|
| 12 | ## Orodr ou non pour l orographie |
---|
| 13 | ok_orodr=y |
---|
| 14 | ## Orolf ou non pour l orographie |
---|
| 15 | ok_orolf=y |
---|
| 16 | ## Si = .T. , lecture du fichier limit avec la bonne annee |
---|
| 17 | ok_limitvrai=n |
---|
| 18 | ## Nombre d'appels des routines de rayonnements ( par jour) |
---|
| 19 | nbapp_rad=24 |
---|
| 20 | ## Flag pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect |
---|
| 21 | iflag_con=30 |
---|
[762] | 22 | ## Facteur multiplication des precip convectives dans KE |
---|
| 23 | cvl_corr=1.0 |
---|
| 24 | ## Facteur additif pour l'albedo |
---|
[984] | 25 | pmagic=0.0 |
---|
[686] | 26 | # |
---|
[762] | 27 | # |
---|
| 28 | # |
---|
[650] | 29 | # Parametres fichiers de sortie |
---|
| 30 | # |
---|
| 31 | ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc |
---|
| 32 | OK_journe=n |
---|
| 33 | ### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc |
---|
| 34 | OK_mensuel=y |
---|
| 35 | ### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la physique) |
---|
| 36 | OK_instan=n |
---|
| 37 | ### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc |
---|
| 38 | ok_hf=n |
---|
[988] | 39 | # Noms des fichiers |
---|
[1507] | 40 | phys_out_filenames= histmth histday histhf histhf3h histhf3hm histstn |
---|
[988] | 41 | # Sortir ou non les fichiers |
---|
[1507] | 42 | phys_out_filekeys= y y y y y n |
---|
[988] | 43 | # Niveaux de sorties |
---|
[1507] | 44 | phys_out_filelevels= 5 5 5 5 5 5 |
---|
| 45 | ### Type de fichier : global (n) ou stations (y) |
---|
| 46 | phys_out_filestations = n n n n n y |
---|
| 47 | ### Frequences des sorties |
---|
| 48 | phys_out_filetimesteps = 1.mth, 1.day, 0.25day, 0.125day, 0.125day, 1800.s |
---|
| 49 | ### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott |
---|
| 50 | flag_t100 = 1, 7, 7, 5, 10, 10 |
---|
| 51 | flag_u100 = 1, 7, 7, 5, 10, 10 |
---|
| 52 | flag_v100 = 1, 7, 7, 5, 10, 10 |
---|
| 53 | flag_w100 = 1, 7, 7, 5, 10, 10 |
---|
| 54 | flag_t50 = 1, 7, 7, 5, 10, 10 |
---|
| 55 | flag_u50 = 1, 7, 7, 5, 10, 10 |
---|
| 56 | flag_v50 = 1, 7, 7, 5, 10, 10 |
---|
| 57 | flag_w50 = 1, 7, 7, 5, 10, 10 |
---|
| 58 | flag_t10 = 1, 7, 7, 5, 10, 10 |
---|
| 59 | flag_u10 = 1, 7, 7, 5, 10, 10 |
---|
| 60 | flag_v10 = 1, 7, 7, 5, 10, 10 |
---|
| 61 | flag_w10 = 1, 7, 7, 5, 10, 10 |
---|
| 62 | # more variables for level 5 histday/month: vitw and pres |
---|
[1065] | 63 | # Vertical wind |
---|
[1507] | 64 | flag_vitw = 2, 3, 7, 6, 6, 6 |
---|
[1065] | 65 | # Air pressure |
---|
[1507] | 66 | flag_pres = 2, 3, 7, 6, 6, 6 |
---|
[1283] | 67 | # w500 for histday |
---|
[1507] | 68 | flag_w500 = 1, 5, 3, 10, 10, 10 |
---|
[1283] | 69 | # t700 for histday |
---|
[1507] | 70 | flag_t700 = 1, 5, 3, 10, 10, 10 |
---|
[1382] | 71 | # Cloud fraction |
---|
[1507] | 72 | flag_rneb = 2, 5, 10, 5, 10, 10 |
---|
| 73 | # Aga sortie avec niveau 2 |
---|
| 74 | flag_Aga = 2, 10, 10, 10, 10, 10 |
---|
| 75 | # niveaux maximum 5 mensuel pour variables CFMIP |
---|
| 76 | flag_rlu = 5, 10, 10, 10, 10, 10 |
---|
| 77 | flag_rsu = 5, 10, 10, 10, 10, 10 |
---|
| 78 | flag_rld = 5, 10, 10, 10, 10, 10 |
---|
| 79 | flag_rsd = 5, 10, 10, 10, 10, 10 |
---|
| 80 | flag_rlucs = 5, 10, 10, 10, 10, 10 |
---|
| 81 | flag_rsucs = 5, 10, 10, 10, 10, 10 |
---|
| 82 | flag_rldcs = 5, 10, 10, 10, 10, 10 |
---|
| 83 | flag_rsdcs = 5, 10, 10, 10, 10, 10 |
---|
| 84 | flag_tnt = 5, 10, 10, 10, 10, 10 |
---|
| 85 | flag_dtdyn = 4, 10, 10, 10, 10, 10 |
---|
| 86 | flag_dtphy = 4, 10, 10, 10, 10, 10 |
---|
| 87 | flag_dtlschr = 5, 10, 10, 10, 10, 10 |
---|
| 88 | flag_tntr = 5, 10, 10, 10, 10, 10 |
---|
| 89 | flag_tntc = 5, 10, 10, 10, 10, 10 |
---|
| 90 | flag_tntscpbl = 5, 10, 10, 10, 10, 10 |
---|
| 91 | flag_tnhus = 5, 10, 10, 10, 10, 10 |
---|
| 92 | flag_dqdyn = 4, 10, 10, 10, 10, 10 |
---|
| 93 | flag_tnhusc = 5, 10, 10, 10, 10, 10 |
---|
| 94 | flag_dqvdf = 5, 10, 10, 10, 10, 10 |
---|
| 95 | flag_tnhusscpbl = 5, 10, 10, 10, 10, 10 |
---|
| 96 | flag_dqphy = 2, 10, 10, 10, 10, 10 |
---|
| 97 | flag_evu = 5, 10, 10, 10, 10, 10 |
---|
| 98 | flag_kz = 4, 10, 10, 10, 10, 10 |
---|
| 99 | ### |
---|
| 100 | # more variables for t tendency monthly |
---|
| 101 | ### |
---|
| 102 | flag_dtvdf = 5, 10, 10, 10, 10, 10 |
---|
| 103 | flag_dteva = 5, 10, 10, 10, 10, 10 |
---|
| 104 | flag_dtcon = 4, 10, 10, 10, 10, 10 |
---|
| 105 | flag_dtajs = 4, 10, 10, 10, 10, 10 |
---|
| 106 | flag_dtlsc = 4, 10, 10, 10, 10, 10 |
---|
| 107 | flag_dtoro = 4, 10, 10, 10, 10, 10 |
---|
| 108 | flag_dtlif = 4, 10, 10, 10, 10, 10 |
---|
| 109 | flag_dthin = 4, 10, 10, 10, 10, 10 |
---|
| 110 | flag_dtec = 4, 10, 10, 10, 10, 10 |
---|
| 111 | ### more variables for q tendency monthly |
---|
| 112 | flag_dqeva = 5, 10, 10, 10, 10, 10 |
---|
| 113 | flag_dqcon = 4, 10, 10, 10, 10, 10 |
---|
| 114 | flag_dqajs = 4, 10, 10, 10, 10, 10 |
---|
| 115 | flag_dqlsc = 5, 10, 10, 10, 10, 10 |
---|
| 116 | ### nouvelle physique monthly |
---|
| 117 | flag_dtthe = 4, 10, 10, 10, 10, 10 |
---|
| 118 | flag_dtwak = 4, 5, 10, 10, 10, 10 |
---|
| 119 | flag_dqthe = 4, 10, 10, 10, 10, 10 |
---|
| 120 | flag_dqwak = 4, 5, 10, 10, 10, 10 |
---|
| 121 | flag_ftime_th = 4, 10, 10, 5, 10, 10 |
---|
[650] | 122 | # |
---|
[988] | 123 | # Parametres niveau de sorties differents fichiers |
---|
[650] | 124 | # |
---|
| 125 | ### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" |
---|
| 126 | # - lev_histhf=0 => pas de sorties histhf.nc |
---|
| 127 | # - lev_histhf=2 => defaut |
---|
| 128 | # - lev_histhf=3 => variables sur niveaux standards |
---|
[988] | 129 | # - lev_histhf=4 => champs 3d niveaux modele dans le fichier histhf.nc |
---|
| 130 | lev_histhf=5 |
---|
[650] | 131 | ### lev_histday=0-5, niveau de sortie fichier "histday.nc" |
---|
| 132 | # - lev_histday=0 => pas de sorties lev_histday.nc |
---|
| 133 | # - lev_histday=2 => defaut |
---|
| 134 | # - lev_histday=3 => + champs 3D => F. Lott |
---|
| 135 | # - lev_histday=4 => + champs sous-surfaces |
---|
| 136 | # - lev_histday=5 => + champs F. Aires |
---|
[988] | 137 | lev_histday=5 |
---|
[650] | 138 | ### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" |
---|
| 139 | # - lev_histmth=0 => pas de sorties lev_histmth.nc |
---|
| 140 | # - lev_histmth=2 => defaut |
---|
| 141 | # - lev_histmth=3 => albedo, rugosite sous-surfaces |
---|
| 142 | # - lev_histmth=4 => champs tendances 3d |
---|
[988] | 143 | lev_histmth=5 |
---|
| 144 | ### parametres pour hist*NMC.nc |
---|
| 145 | #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier |
---|
| 146 | # histdayNMC: 8 ou 17 |
---|
| 147 | lev_histdayNMC=8 |
---|
[1003] | 148 | ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence |
---|
| 149 | ok_histNMC = y, n, n |
---|
[650] | 150 | ### ecrit_day = frequence ecriture fichier histday.nc en jours |
---|
| 151 | ecrit_day=1. |
---|
| 152 | ### ecrit_mth = frequence ecriture fichier histmth.nc en jours |
---|
[697] | 153 | ecrit_mth=_ecrit-mth_ |
---|
[650] | 154 | ### freqin_isccp = frequence input en secondes du simulateur ISCCP |
---|
| 155 | freq_ISCCP=10800. |
---|
| 156 | ### freqout_isccp = frequence output en jours du simulateur ISCCP |
---|
[697] | 157 | ecrit_ISCCP=_ecrit-ISCCP_ |
---|
[650] | 158 | ### niveau du diagnostique de conservation d energie |
---|
| 159 | if_ebil=0 |
---|
| 160 | # |
---|
| 161 | # parametres KE |
---|
| 162 | # |
---|
| 163 | ### epmax = Efficacite precipitation maximale |
---|
[866] | 164 | epmax = .999 |
---|
[650] | 165 | ### ok_adj_ema = ?? pas utilise |
---|
| 166 | ok_adj_ema = n |
---|
| 167 | ### iflag_clw Flag calcul eau liquide |
---|
| 168 | # - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l) |
---|
| 169 | # - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l) |
---|
| 170 | # - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip |
---|
| 171 | iflag_clw = 1 |
---|
| 172 | # |
---|
| 173 | # parametres nuages |
---|
| 174 | # |
---|
| 175 | ### cld_lc_lsc contenu en eau liquide des nuages large-scale (fisrtilp) |
---|
[866] | 176 | cld_lc_lsc = 4.16e-4 |
---|
[650] | 177 | ### cld_lc_con contenu en eau liquide des nuages convectifs (fisrtilp) |
---|
[866] | 178 | cld_lc_con = 4.16e-4 |
---|
[686] | 179 | ### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) |
---|
[650] | 180 | cld_tau_lsc = 1800. |
---|
[686] | 181 | ### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) |
---|
[650] | 182 | cld_tau_con = 1800. |
---|
| 183 | ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) |
---|
[866] | 184 | ffallv_lsc = 0.5 |
---|
[650] | 185 | ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) |
---|
[866] | 186 | ffallv_con = 0.5 |
---|
[686] | 187 | ### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?) |
---|
[650] | 188 | coef_eva = 2.e-5 |
---|
| 189 | ### reevap_ice reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) |
---|
| 190 | reevap_ice = y |
---|
| 191 | ### iflag_cldcon flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) |
---|
| 192 | # - iflag_cldcon<=-1 diag. rain_Tiedtke |
---|
| 193 | # - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k) |
---|
| 194 | # - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc) |
---|
| 195 | # - iflag_cldcon=3, ratqs=ratqss |
---|
| 196 | iflag_cldcon = 3 |
---|
| 197 | ### iflag_pdf : flag calcul distribution sous-maille de l eau et des nuages |
---|
| 198 | # - iflag_pdf=0, version ratqs, |
---|
| 199 | # - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir |
---|
| 200 | # - des PDFS de Sandrine Bony |
---|
| 201 | iflag_pdf = 1 |
---|
| 202 | ### fact_cldcon constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) |
---|
| 203 | fact_cldcon = 1. |
---|
| 204 | ### facttemps= facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3) |
---|
[814] | 205 | facttemps = 0. |
---|
[650] | 206 | ## ok_newmicro =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) |
---|
| 207 | ok_newmicro = y |
---|
| 208 | ### iflag_ratqs=0 correspond a la version IPCC AR4 |
---|
| 209 | iflag_ratqs=0 |
---|
| 210 | ### ratqsbas ratqs en bas si iflag_cldcon=1 |
---|
| 211 | ratqsbas = 0.005 |
---|
| 212 | ### ratqshaut ratqs en haut pour ratqss "stables" |
---|
| 213 | ratqshaut = 0.33 |
---|
| 214 | ### rad_froid rayon cristaux des nuages de glace "froids" |
---|
| 215 | rad_froid = 35 |
---|
| 216 | ### rad_chau1 rayon goutelettes d eau chauds", en haut: k=4-klev |
---|
| 217 | rad_chau1=12 |
---|
| 218 | ### rad_chau2 rayon goutelettes d eau chauds", en bas: k=1-3 |
---|
| 219 | rad_chau2=11 |
---|
[686] | 220 | # |
---|
| 221 | # Coefficient et parametres sur les drags |
---|
| 222 | # |
---|
| 223 | f_cdrag_ter=1. |
---|
| 224 | f_cdrag_stable=1. |
---|
| 225 | f_cdrag_oce=0.8 |
---|
| 226 | f_rugoro=0. |
---|
| 227 | ### ksta_ter coef.diffusion minimale sur terre/sic/lic |
---|
[650] | 228 | ksta_ter=1.e-7 |
---|
| 229 | ### cdmmax = cdrag maximum pour le moment |
---|
| 230 | cdmmax = 2.5E-3 |
---|
| 231 | ### cdhmax = cdrag maximum pour l energie |
---|
| 232 | cdhmax = 2.0E-3 |
---|
| 233 | # |
---|
| 234 | # Parametres "orbitaux/ ere geologique" |
---|
| 235 | # |
---|
| 236 | ### R_ecc = Excentricite |
---|
| 237 | R_ecc = 0.016715 |
---|
| 238 | ### R_peri = Equinoxe |
---|
| 239 | R_peri = 102.7 |
---|
| 240 | ### R_incl = Inclinaison |
---|
| 241 | R_incl = 23.441 |
---|
| 242 | ### solaire = Constante solaire |
---|
[1026] | 243 | solaire = 1365.6537 |
---|
[650] | 244 | # |
---|
| 245 | # Taux gaz a effet de serre |
---|
| 246 | # |
---|
| 247 | ### co2_ppm = taux CO2 en ppm |
---|
[1040] | 248 | co2_ppm = 0.36886500E+03 |
---|
[650] | 249 | ### CH4_ppb = taux CH4 en ppb |
---|
[1040] | 250 | CH4_ppb = 0.17510225E+04 |
---|
[650] | 251 | ### N2O_ppb = taux N2O en ppb |
---|
[1040] | 252 | N2O_ppb = 0.31585000E+03 |
---|
[650] | 253 | ### CFC11_ppt = taux CFC11 en ppt |
---|
[1040] | 254 | CFC11_ppt = 5.18015181E+01 |
---|
[650] | 255 | ### CFC12_ppt = taux CFC12 en ppt |
---|
[1040] | 256 | CFC12_ppt = 0.99862742E+03 |
---|
[650] | 257 | # |
---|
| 258 | # Parametres effets directs/indirects des "aerosols" |
---|
| 259 | # |
---|
| 260 | ### ok_ade=y/n flag Aerosol direct effect |
---|
[762] | 261 | ok_ade=y |
---|
[650] | 262 | ### ok_aie=y/n flag Aerosol indirect effect |
---|
[762] | 263 | ok_aie=y |
---|
[814] | 264 | ### aer_type = Aerosol variation type : actuel / preind / scenario / annuel |
---|
[988] | 265 | aer_type=annuel |
---|
[762] | 266 | ### type of coupled aerosol =1 (default) =2 => bc only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol |
---|
| 267 | flag_aerosol=6 |
---|
[650] | 268 | ### bl95_b0 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
---|
| 269 | bl95_b0=1.7 |
---|
| 270 | ### bl95_b1 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
---|
| 271 | bl95_b1=0.2 |
---|
| 272 | # |
---|
[814] | 273 | # Parametre de lecture de l'ozone |
---|
| 274 | # |
---|
| 275 | # Allowed values are 0, 1 and 2 |
---|
| 276 | # 0: do not read an ozone climatology |
---|
| 277 | # 1: read a single ozone climatology that will be used day and night |
---|
| 278 | # 2: read two ozone climatologies, the average day and night climatology and the daylight climatology |
---|
[825] | 279 | read_climoz=2 |
---|
[814] | 280 | # |
---|
| 281 | # Parametres simulateur COSP (CFMIP Observational Simulator Package) |
---|
| 282 | # |
---|
| 283 | ### ok_cosp=y/n flag simulateur COSP |
---|
| 284 | ok_cosp=n |
---|
[1000] | 285 | ### freq_COSP = frequence d'appel de COSP en secondes |
---|
[814] | 286 | freq_COSP=10800. |
---|
[1000] | 287 | ### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc |
---|
| 288 | ok_mensuelCOSP=y |
---|
| 289 | ### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc |
---|
| 290 | ok_journeCOSP=y |
---|
| 291 | ### ok_hfCOSP=y/n, ecrire sorties haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc |
---|
| 292 | ok_hfCOSP=n |
---|
[814] | 293 | # |
---|
[650] | 294 | # Parametres simulateur ISCCP |
---|
| 295 | # |
---|
| 296 | ### ok_isccp=y/n flag simulateur ISCCP |
---|
| 297 | ok_isccp=n |
---|
| 298 | ### top_height = flag choix calcul nuages par le simulateur en utilisant |
---|
| 299 | # - les donnees IR et/ou VIS et l algorithme ISCCP-D1 |
---|
| 300 | # - top_height = 1 -> algo IR-VIS |
---|
| 301 | # - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" |
---|
| 302 | # - top_height = 3 -> algo IR |
---|
| 303 | top_height = 1 |
---|
| 304 | ### overlap = Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP |
---|
| 305 | # - overlap=1 Max overlap |
---|
| 306 | # - overlap=2 Random overlap |
---|
| 307 | # - overlap=3 Max/Random overlap |
---|
| 308 | overlap = 3 |
---|