############################################################################## # PARAMETRES ET CLES DE CONFIGURATION LMDZ # Cette parie sera bientot extrait dans un fichier config.def ############################################################################## # D: Valeur par default # # Controle des sorties #********************* #ok_regdyn : y/n calcul/non des regymes dynamiques sur regions pre-definies ok_regdyn=y # A transformer en la ligne ci-dessous quand ce sera compatible avec libIGCM phys_out_filekeys= y y y y y n phys_out_filenames= histmth histday histhf histhf3h histhf3hm histstn phys_out_filelevels= 5 5 5 5 5 5 phys_out_filetypes= ave(X) ave(X) ave(X) inst(X) inst(X) phys_out_filetimesteps= 1.mth, 1.day, 0.25day, 0.125day, 0.125day, 1800.s ### parametres pour hist*NMC.nc #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier # histdayNMC: 8 ou 17 lev_histdayNMC=8 ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence ok_histNMC = y, n, n # #Couplage avec autres modules #****************************** ## Soil Model ou non soil_model=y ## Cycle diurne ou non cycle_diurne=y #Code de transfert radiatif #************************** # ajout de ce flag pour activer/desactiver le rayonnement (MPL) # 0 : pas de rayonnement. 1: on active le rayonnement (D=1) iflag_radia=1 ## Nombre appels des routines de rayonnements ( par jour) nbapp_rad=24 # Parametres "orbitaux/ ere geologique" ####################################### #Parametres orbitaux #excentricite R_ecc=0.016715 #equinoxe R_peri=102.7 #inclinaison R_incl=23.441 #constante solaire solaire = 1366.0896 ## Facteur additif pour l albedo pmagic=999.9 # # Taux gaz a effet de serre ####################################### #concentration des gaz (co2 ch4 n2o cfc11 cfc12) # (D: 348. 1650. 306. 280. 484.) co2_ppm = 0.36886500E+03 #RCO2=co2_ppm * 1.0e-06 * 44.011/28.97= 5.286789092164308E-04 CH4_ppb = 0.17510225E+04 #RCH4=1.65E-06* 16.043/28.97= 9.137366240938903E-07 N2O_ppb = 0.31585000E+03 #RN2O=306.E-09* 44.013/28.97= 4.648939592682085E-07 CFC11_ppt = 5.18015181E+01 #RCFC11=280.E-12* 137.3686/28.97= 1.327690990680013E-09 CFC12_ppt = 0.99862742E+03 #RCFC12=484.E-12* 120.9140/28.97= 2.020102726958923E-09 # # Parametres effets directs/indirects des "aerosols" ####################################### ### ok_ade=y/n flag Aerosol direct effect ok_ade=n ### ok_aie=y/n flag Aerosol indirect effect ok_aie=n ### aer_type = Aerosol variation type : actuel / preind / scenario / annuel aer_type=xxx ### type of coupled aerosol =1 (default) =2 => bc only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol flag_aerosol=6 ### bl95_b0 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) bl95_b0=1.7 ### bl95_b1 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) bl95_b1=0.2 # # Parametre de lecture de l'ozone ####################################### # Allowed values are 0, 1 and 2 # 0: do not read an ozone climatology # 1: read a single ozone climatology that will be used day and night # 2: read two ozone climatologies, the average day and night climatology and the daylight climatology read_climoz=0 # # Parametres simulateur COSP (CFMIP Observational Simulator Package) ####################################### ### ok_cosp=y/n flag simulateur COSP ok_cosp=n ### freq_COSP = frequence d'appel de COSP en secondes freq_COSP=10800. ### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc ok_mensuelCOSP=y ### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc ok_journeCOSP=y ### ok_hfCOSP=y/n, ecrire sorties haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc ok_hfCOSP=n # # Parametres simulateur ISCCP ####################################### ### ok_isccp=y/n flag simulateur ISCCP ok_isccp=n ### freqin_isccp = frequence input en secondes du simulateur ISCCP freq_ISCCP=10800. ### top_height = flag choix calcul nuages par le simulateur en utilisant # - les donnees IR et/ou VIS et l algorithme ISCCP-D1 # - top_height = 1 -> algo IR-VIS # - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" # - top_height = 3 -> algo IR top_height = 1 ### overlap = Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP # - overlap=1 Max overlap # - overlap=2 Random overlap # - overlap=3 Max/Random overlap overlap = 3 # # Specifique aquaplanete ######################## #qsol0=100. ############################################################################## # PARAMETRES ET CLES DE LA PHYSIQUE LMDZ # Cette partie constitue la partie : physiq.def_L39_NPv3.0 ############################################################################## # D: Valeur par default # #Parametres Orographiques et cdrags #********************************** ## Si=.T. , lecture du fichier limit avec la bonne annee ok_limitvrai=n #Cdrags f_cdrag_stable=1. f_cdrag_ter=1. f_cdrag_oce=0.7 #cdmmax cdmmax=2.5E-3 #cdhmax cdhmax=2.0E-3 ## Orodr ou non pour l orographie ok_orodr=y ## Orolf ou non pour l orographie ok_orolf=y #Rugoro f_rugoro=0. # # Rayonnement #********************* # activation du nouveau code de rayonnement RRTM # 0 : Ancien code et 1 : RRTM (D=0) iflag_rrtm=0 # ## Activation ou non de la parametrisation de Hines pour la strato #****************************************************************** ok_hines=y # # Parametres nuages #****************** #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4) cld_lc_lsc=6.e-4 cld_lc_con=6.e-4 #constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.) cld_tau_lsc=1800. cld_tau_con=1800. #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1) ffallv_lsc=1.35 ffallv_con=1.35 #coefficient sur la reevaporation de la pluie (D: 2.e-5 n) coef_eva=1e-4 reevap_ice=y #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) iflag_cldcon=6 fact_cldcon=1. facttemps=0. #calcul eau condensee et fraction nuageuse a partir des PDF (D:0, 0:version avec ratqs sinon nouvelles PDFs) iflag_pdf=1 #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3) ok_newmicro=y iflag_ratqs=2 ratqsbas=0.002 ratqshaut=0.25 #rayon effectif pour les nuages de glace (D:35) rad_froid=35 #rayons effectifs pour les nuages eau liq (tailles des gouttes d eau nuageuse, D: 13 9) rad_chau1=12 rad_chau2=11 ## Choix ou non de New oliq new_oliq=y # rei_min=20. rei_max=61.29 # t_glace_min=258. t_glace_max=273.13 # # Flags Convection #***************** #flag pour la convection (D:2, 1:LMD, 2:Tiedtke, 3:KE nouvelle physique, 30:KE IPCC) iflag_con=3 #niveau de sorties des diagnos sur la conservation energie if_ebil=0 #efficacite de precepitation maxim (D:.993) epmax=0.997 #ajustement convectif sec au debut de la convection (D: n 0) ok_adj_ema=n iflag_clw=0 #fermeture convective (D:1, 0:AR4, 2:ALE and ALP) iflag_clos=2 #loi de melange a l entrainement (D:1, 0:AR4=PDF plate, 1=PDF) iflag_mix=1 #poids des PDFs plate et en cloche (D: 1 0) qqa1=1. qqa2=0. # Facteur multiplication des precip convectives dans KE cvl_corr=1.0 # Fmax=0.65 # # Flags Couche limite #******************** #shema de CL des thermiques (D:0, 0:ajustement sec,=>1:versions thermiques) iflag_thermals=15 #pas de temps du spliting pour les thermiques nsplit_thermals=1 tau_thermals=0. iflag_thermals_ed=10 iflag_thermals_optflux=0 #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) iflag_pbl=8 #diffusion turbulente ksta_ter=1.e-7 ksta=1.e-10 #ok_kzmin : calcul Kzmin dans la CL de surface ok_kzmin=y #couplage ou non avec la convection iflag_coupl=5 (0:AR4 et 1:nouvelle physique) seuil_inversion=-0.08 #rendement sur la puissance fournie par les thermiques a la convection alp_bl_k=0.5 # # Flags wakes #************ #activation ou non des wakes (D:0, 0:AR4, 1:nouvelle physique) iflag_wake=1 alp_offset=0. nsplit_phys=2