| 1 | ############################################################################## |
|---|
| 2 | # PARAMETRES ET CLES DE CONFIGURATION LMDZ |
|---|
| 3 | # Cette parie sera bientot extrait dans un fichier config.def |
|---|
| 4 | ############################################################################## |
|---|
| 5 | # D: Valeur par default |
|---|
| 6 | # |
|---|
| 7 | # Controle des sorties |
|---|
| 8 | #********************* |
|---|
| 9 | #ok_regdyn : y/n calcul/non des regymes dynamiques sur regions pre-definies |
|---|
| 10 | ok_regdyn=y |
|---|
| 11 | # A transformer en la ligne ci-dessous quand ce sera compatible avec libIGCM |
|---|
| 12 | phys_out_filekeys= y y y y y n |
|---|
| 13 | phys_out_filenames= histmth histday histhf histhf3h histhf3hm histstn |
|---|
| 14 | phys_out_filelevels= 5 5 5 5 5 5 |
|---|
| 15 | phys_out_filetypes= ave(X) ave(X) ave(X) inst(X) inst(X) |
|---|
| 16 | phys_out_filetimesteps= 1.mth, 1.day, 0.25day, 0.125day, 0.125day, 1800.s |
|---|
| 17 | ### parametres pour hist*NMC.nc |
|---|
| 18 | #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier |
|---|
| 19 | # histdayNMC: 8 ou 17 |
|---|
| 20 | lev_histdayNMC=8 |
|---|
| 21 | ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence |
|---|
| 22 | ok_histNMC = y, n, n |
|---|
| 23 | # |
|---|
| 24 | #Couplage avec autres modules |
|---|
| 25 | #****************************** |
|---|
| 26 | ## Soil Model ou non |
|---|
| 27 | soil_model=y |
|---|
| 28 | ## Cycle diurne ou non |
|---|
| 29 | cycle_diurne=y |
|---|
| 30 | |
|---|
| 31 | #Code de transfert radiatif |
|---|
| 32 | #************************** |
|---|
| 33 | # ajout de ce flag pour activer/desactiver le rayonnement (MPL) |
|---|
| 34 | # 0 : pas de rayonnement. 1: on active le rayonnement (D=1) |
|---|
| 35 | iflag_radia=1 |
|---|
| 36 | ## Nombre appels des routines de rayonnements ( par jour) |
|---|
| 37 | nbapp_rad=24 |
|---|
| 38 | |
|---|
| 39 | # Parametres "orbitaux/ ere geologique" |
|---|
| 40 | ####################################### |
|---|
| 41 | #Parametres orbitaux |
|---|
| 42 | #excentricite |
|---|
| 43 | R_ecc=0.016715 |
|---|
| 44 | #equinoxe |
|---|
| 45 | R_peri=102.7 |
|---|
| 46 | #inclinaison |
|---|
| 47 | R_incl=23.441 |
|---|
| 48 | #constante solaire |
|---|
| 49 | solaire = 1366.0896 |
|---|
| 50 | ## Facteur additif pour l albedo |
|---|
| 51 | pmagic=999.9 |
|---|
| 52 | # |
|---|
| 53 | # Taux gaz a effet de serre |
|---|
| 54 | ####################################### |
|---|
| 55 | #concentration des gaz (co2 ch4 n2o cfc11 cfc12) |
|---|
| 56 | # (D: 348. 1650. 306. 280. 484.) |
|---|
| 57 | co2_ppm = 0.36886500E+03 |
|---|
| 58 | #RCO2=co2_ppm * 1.0e-06 * 44.011/28.97= 5.286789092164308E-04 |
|---|
| 59 | CH4_ppb = 0.17510225E+04 |
|---|
| 60 | #RCH4=1.65E-06* 16.043/28.97= 9.137366240938903E-07 |
|---|
| 61 | N2O_ppb = 0.31585000E+03 |
|---|
| 62 | #RN2O=306.E-09* 44.013/28.97= 4.648939592682085E-07 |
|---|
| 63 | CFC11_ppt = 5.18015181E+01 |
|---|
| 64 | #RCFC11=280.E-12* 137.3686/28.97= 1.327690990680013E-09 |
|---|
| 65 | CFC12_ppt = 0.99862742E+03 |
|---|
| 66 | #RCFC12=484.E-12* 120.9140/28.97= 2.020102726958923E-09 |
|---|
| 67 | # |
|---|
| 68 | # Parametres effets directs/indirects des "aerosols" |
|---|
| 69 | ####################################### |
|---|
| 70 | ### ok_ade=y/n flag Aerosol direct effect |
|---|
| 71 | ok_ade=y |
|---|
| 72 | ### ok_aie=y/n flag Aerosol indirect effect |
|---|
| 73 | ok_aie=y |
|---|
| 74 | ### aer_type = Aerosol variation type : actuel / preind / scenario / annuel |
|---|
| 75 | aer_type=xxx |
|---|
| 76 | ### type of coupled aerosol =1 (default) =2 => bc only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol |
|---|
| 77 | flag_aerosol=6 |
|---|
| 78 | ### bl95_b0 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
|---|
| 79 | bl95_b0=1.7 |
|---|
| 80 | ### bl95_b1 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
|---|
| 81 | bl95_b1=0.2 |
|---|
| 82 | # |
|---|
| 83 | # Parametre de lecture de l'ozone |
|---|
| 84 | ####################################### |
|---|
| 85 | # Allowed values are 0, 1 and 2 |
|---|
| 86 | # 0: do not read an ozone climatology |
|---|
| 87 | # 1: read a single ozone climatology that will be used day and night |
|---|
| 88 | # 2: read two ozone climatologies, the average day and night climatology and the daylight climatology |
|---|
| 89 | read_climoz=2 |
|---|
| 90 | # |
|---|
| 91 | # Parametres simulateur COSP (CFMIP Observational Simulator Package) |
|---|
| 92 | ####################################### |
|---|
| 93 | ### ok_cosp=y/n flag simulateur COSP |
|---|
| 94 | ok_cosp=n |
|---|
| 95 | ### freq_COSP = frequence d'appel de COSP en secondes |
|---|
| 96 | freq_COSP=10800. |
|---|
| 97 | ### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc |
|---|
| 98 | ok_mensuelCOSP=y |
|---|
| 99 | ### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc |
|---|
| 100 | ok_journeCOSP=y |
|---|
| 101 | ### ok_hfCOSP=y/n, ecrire sorties haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc |
|---|
| 102 | ok_hfCOSP=n |
|---|
| 103 | # |
|---|
| 104 | # Parametres simulateur ISCCP |
|---|
| 105 | ####################################### |
|---|
| 106 | ### ok_isccp=y/n flag simulateur ISCCP |
|---|
| 107 | ok_isccp=n |
|---|
| 108 | ### freqin_isccp = frequence input en secondes du simulateur ISCCP |
|---|
| 109 | freq_ISCCP=10800. |
|---|
| 110 | ### top_height = flag choix calcul nuages par le simulateur en utilisant |
|---|
| 111 | # - les donnees IR et/ou VIS et l algorithme ISCCP-D1 |
|---|
| 112 | # - top_height = 1 -> algo IR-VIS |
|---|
| 113 | # - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" |
|---|
| 114 | # - top_height = 3 -> algo IR |
|---|
| 115 | top_height = 1 |
|---|
| 116 | ### overlap = Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP |
|---|
| 117 | # - overlap=1 Max overlap |
|---|
| 118 | # - overlap=2 Random overlap |
|---|
| 119 | # - overlap=3 Max/Random overlap |
|---|
| 120 | overlap = 3 |
|---|
| 121 | # |
|---|
| 122 | |
|---|
| 123 | # Specifique aquaplanete |
|---|
| 124 | ######################## |
|---|
| 125 | #qsol0=100. |
|---|
| 126 | ############################################################################## |
|---|
| 127 | # PARAMETRES ET CLES DE LA PHYSIQUE LMDZ |
|---|
| 128 | # Cette partie constitue la partie : physiq.def_L39_NPv3.0 |
|---|
| 129 | ############################################################################## |
|---|
| 130 | # D: Valeur par default |
|---|
| 131 | # |
|---|
| 132 | #Parametres Orographiques et cdrags |
|---|
| 133 | #********************************** |
|---|
| 134 | ## Si=.T. , lecture du fichier limit avec la bonne annee |
|---|
| 135 | ok_limitvrai=n |
|---|
| 136 | #Cdrags |
|---|
| 137 | f_cdrag_stable=1. |
|---|
| 138 | f_cdrag_ter=1. |
|---|
| 139 | f_cdrag_oce=0.7 |
|---|
| 140 | #cdmmax |
|---|
| 141 | cdmmax=2.5E-3 |
|---|
| 142 | #cdhmax |
|---|
| 143 | cdhmax=2.0E-3 |
|---|
| 144 | ## Orodr ou non pour l orographie |
|---|
| 145 | ok_orodr=y |
|---|
| 146 | ## Orolf ou non pour l orographie |
|---|
| 147 | ok_orolf=y |
|---|
| 148 | #Rugoro |
|---|
| 149 | f_rugoro=0. |
|---|
| 150 | # |
|---|
| 151 | # Rayonnement |
|---|
| 152 | #********************* |
|---|
| 153 | # activation du nouveau code de rayonnement RRTM |
|---|
| 154 | # 0 : Ancien code et 1 : RRTM (D=0) |
|---|
| 155 | iflag_rrtm=0 |
|---|
| 156 | # |
|---|
| 157 | ## Activation ou non de la parametrisation de Hines pour la strato |
|---|
| 158 | #****************************************************************** |
|---|
| 159 | ok_hines=y |
|---|
| 160 | # |
|---|
| 161 | # Parametres nuages |
|---|
| 162 | #****************** |
|---|
| 163 | #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4) |
|---|
| 164 | cld_lc_lsc=6.e-4 |
|---|
| 165 | cld_lc_con=6.e-4 |
|---|
| 166 | #constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.) |
|---|
| 167 | cld_tau_lsc=1800. |
|---|
| 168 | cld_tau_con=1800. |
|---|
| 169 | #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1) |
|---|
| 170 | ffallv_lsc=1.35 |
|---|
| 171 | ffallv_con=1.35 |
|---|
| 172 | #coefficient sur la reevaporation de la pluie (D: 2.e-5 n) |
|---|
| 173 | coef_eva=1e-4 |
|---|
| 174 | reevap_ice=y |
|---|
| 175 | #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) |
|---|
| 176 | iflag_cldcon=6 |
|---|
| 177 | fact_cldcon=1. |
|---|
| 178 | facttemps=0. |
|---|
| 179 | #calcul eau condensee et fraction nuageuse a partir des PDF (D:0, 0:version avec ratqs sinon nouvelles PDFs) |
|---|
| 180 | iflag_pdf=1 |
|---|
| 181 | #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3) |
|---|
| 182 | ok_newmicro=y |
|---|
| 183 | iflag_ratqs=2 |
|---|
| 184 | ratqsbas=0.002 |
|---|
| 185 | ratqshaut=0.25 |
|---|
| 186 | #rayon effectif pour les nuages de glace (D:35) |
|---|
| 187 | rad_froid=35 |
|---|
| 188 | #rayons effectifs pour les nuages eau liq (tailles des gouttes d eau nuageuse, D: 13 9) |
|---|
| 189 | rad_chau1=12 |
|---|
| 190 | rad_chau2=11 |
|---|
| 191 | ## Choix ou non de New oliq |
|---|
| 192 | new_oliq=y |
|---|
| 193 | # |
|---|
| 194 | rei_min=20. |
|---|
| 195 | rei_max=61.29 |
|---|
| 196 | # |
|---|
| 197 | t_glace_min=258. |
|---|
| 198 | t_glace_max=273.13 |
|---|
| 199 | # |
|---|
| 200 | # Flags Convection |
|---|
| 201 | #***************** |
|---|
| 202 | #flag pour la convection (D:2, 1:LMD, 2:Tiedtke, 3:KE nouvelle physique, 30:KE IPCC) |
|---|
| 203 | iflag_con=3 |
|---|
| 204 | #niveau de sorties des diagnos sur la conservation energie |
|---|
| 205 | if_ebil=0 |
|---|
| 206 | #efficacite de precepitation maxim (D:.993) |
|---|
| 207 | epmax=0.997 |
|---|
| 208 | #ajustement convectif sec au debut de la convection (D: n 0) |
|---|
| 209 | ok_adj_ema=n |
|---|
| 210 | iflag_clw=0 |
|---|
| 211 | #fermeture convective (D:1, 0:AR4, 2:ALE and ALP) |
|---|
| 212 | iflag_clos=2 |
|---|
| 213 | #loi de melange a l entrainement (D:1, 0:AR4=PDF plate, 1=PDF) |
|---|
| 214 | iflag_mix=1 |
|---|
| 215 | #poids des PDFs plate et en cloche (D: 1 0) |
|---|
| 216 | qqa1=1. |
|---|
| 217 | qqa2=0. |
|---|
| 218 | # Facteur multiplication des precip convectives dans KE |
|---|
| 219 | cvl_corr=1.0 |
|---|
| 220 | # |
|---|
| 221 | Fmax=0.65 |
|---|
| 222 | # |
|---|
| 223 | # Flags Couche limite |
|---|
| 224 | #******************** |
|---|
| 225 | #shema de CL des thermiques (D:0, 0:ajustement sec,=>1:versions thermiques) |
|---|
| 226 | iflag_thermals=15 |
|---|
| 227 | #pas de temps du spliting pour les thermiques |
|---|
| 228 | nsplit_thermals=1 |
|---|
| 229 | tau_thermals=0. |
|---|
| 230 | iflag_thermals_ed=10 |
|---|
| 231 | iflag_thermals_optflux=0 |
|---|
| 232 | #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) |
|---|
| 233 | iflag_pbl=8 |
|---|
| 234 | #diffusion turbulente |
|---|
| 235 | ksta_ter=1.e-7 |
|---|
| 236 | ksta=1.e-10 |
|---|
| 237 | #ok_kzmin : calcul Kzmin dans la CL de surface |
|---|
| 238 | ok_kzmin=y |
|---|
| 239 | #couplage ou non avec la convection |
|---|
| 240 | iflag_coupl=5 (0:AR4 et 1:nouvelle physique) |
|---|
| 241 | seuil_inversion=-0.08 |
|---|
| 242 | #rendement sur la puissance fournie par les thermiques a la convection |
|---|
| 243 | alp_bl_k=0.5 |
|---|
| 244 | # |
|---|
| 245 | # Flags wakes |
|---|
| 246 | #************ |
|---|
| 247 | #activation ou non des wakes (D:0, 0:AR4, 1:nouvelle physique) |
|---|
| 248 | iflag_wake=1 |
|---|
| 249 | alp_offset=0. |
|---|
| 250 | |
|---|
| 251 | nsplit_phys=2 |
|---|
| 252 | |
|---|
