## $Id$ # # PARAMETRES ANCIENNEMENT DANS gcm.def ## Cycle diurne ou non cycle_diurne=y ## Soil Model ou non soil_model=y ## Choix ou non de New oliq new_oliq=y ## Activation ou non de la parametrisation de Hines pour la strato ok_hines=y ## Orodr ou non pour l orographie ok_orodr=y ## Orolf ou non pour l orographie ok_orolf=y ## Si = .T. , lecture du fichier limit avec la bonne annee ok_limitvrai=n ## Nombre d'appels des routines de rayonnements ( par jour) nbapp_rad=24 ## Flag pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect iflag_con=30 ## Facteur multiplication des precip convectives dans KE cvl_corr=1.0 ## Facteur additif pour l'albedo pmagic=0.0 # # # # Parametres fichiers de sortie # ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc OK_journe=n ### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc OK_mensuel=y ### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la physique) OK_instan=n ### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc ok_hf=n # Noms des fichiers phys_out_filenames= histmth histday histhf histhf3h histhf3hm histstn # Sortir ou non les fichiers phys_out_filekeys= y y y y y n # Niveaux de sorties phys_out_filelevels= 5 5 5 5 5 5 ### Type de fichier : global (n) ou stations (y) phys_out_filestations = n n n n n y ### Frequences des sorties phys_out_filetimesteps = 1.mth, 1.day, 0.25day, 0.125day, 0.125day, 1800.s ### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott flag_t100 = 1, 7, 7, 5, 10, 10 flag_u100 = 1, 7, 7, 5, 10, 10 flag_v100 = 1, 7, 7, 5, 10, 10 flag_w100 = 1, 7, 7, 5, 10, 10 flag_t50 = 1, 7, 7, 5, 10, 10 flag_u50 = 1, 7, 7, 5, 10, 10 flag_v50 = 1, 7, 7, 5, 10, 10 flag_w50 = 1, 7, 7, 5, 10, 10 flag_t10 = 1, 7, 7, 5, 10, 10 flag_u10 = 1, 7, 7, 5, 10, 10 flag_v10 = 1, 7, 7, 5, 10, 10 flag_w10 = 1, 7, 7, 5, 10, 10 # more variables for level 5 histday/month: vitw and pres # Vertical wind flag_vitw = 2, 3, 7, 6, 6, 6 # Air pressure flag_pres = 2, 3, 7, 6, 6, 6 # w500 for histday flag_w500 = 1, 5, 3, 10, 10, 10 # t700 for histday flag_t700 = 1, 5, 3, 10, 10, 10 # Cloud fraction flag_rneb = 2, 5, 10, 5, 10, 10 # Aga sortie avec niveau 2 flag_Aga = 2, 10, 10, 10, 10, 10 # niveaux maximum 5 mensuel pour variables CFMIP flag_rlu = 5, 10, 10, 10, 10, 10 flag_rsu = 5, 10, 10, 10, 10, 10 flag_rld = 5, 10, 10, 10, 10, 10 flag_rsd = 5, 10, 10, 10, 10, 10 flag_rlucs = 5, 10, 10, 10, 10, 10 flag_rsucs = 5, 10, 10, 10, 10, 10 flag_rldcs = 5, 10, 10, 10, 10, 10 flag_rsdcs = 5, 10, 10, 10, 10, 10 flag_tnt = 5, 10, 10, 10, 10, 10 flag_dtdyn = 4, 10, 10, 10, 10, 10 flag_dtphy = 4, 10, 10, 10, 10, 10 flag_dtlschr = 5, 10, 10, 10, 10, 10 flag_tntr = 5, 10, 10, 10, 10, 10 flag_tntc = 5, 10, 10, 10, 10, 10 flag_tntscpbl = 5, 10, 10, 10, 10, 10 flag_tnhus = 5, 10, 10, 10, 10, 10 flag_dqdyn = 4, 10, 10, 10, 10, 10 flag_tnhusc = 5, 10, 10, 10, 10, 10 flag_dqvdf = 5, 10, 10, 10, 10, 10 flag_tnhusscpbl = 5, 10, 10, 10, 10, 10 flag_dqphy = 2, 10, 10, 10, 10, 10 flag_evu = 5, 10, 10, 10, 10, 10 flag_kz = 4, 10, 10, 10, 10, 10 ### # more variables for t tendency monthly ### flag_dtvdf = 5, 10, 10, 10, 10, 10 flag_dteva = 5, 10, 10, 10, 10, 10 flag_dtcon = 4, 10, 10, 10, 10, 10 flag_dtajs = 4, 10, 10, 10, 10, 10 flag_dtlsc = 4, 10, 10, 10, 10, 10 flag_dtoro = 4, 10, 10, 10, 10, 10 flag_dtlif = 4, 10, 10, 10, 10, 10 flag_dthin = 4, 10, 10, 10, 10, 10 flag_dtec = 4, 10, 10, 10, 10, 10 ### more variables for q tendency monthly flag_dqeva = 5, 10, 10, 10, 10, 10 flag_dqcon = 4, 10, 10, 10, 10, 10 flag_dqajs = 4, 10, 10, 10, 10, 10 flag_dqlsc = 5, 10, 10, 10, 10, 10 ### nouvelle physique monthly flag_dtthe = 4, 10, 10, 10, 10, 10 flag_dtwak = 4, 5, 10, 10, 10, 10 flag_dqthe = 4, 10, 10, 10, 10, 10 flag_dqwak = 4, 5, 10, 10, 10, 10 flag_ftime_th = 4, 10, 10, 5, 10, 10 # # Parametres niveau de sorties differents fichiers # ### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" # - lev_histhf=0 => pas de sorties histhf.nc # - lev_histhf=2 => defaut # - lev_histhf=3 => variables sur niveaux standards # - lev_histhf=4 => champs 3d niveaux modele dans le fichier histhf.nc lev_histhf=5 ### lev_histday=0-5, niveau de sortie fichier "histday.nc" # - lev_histday=0 => pas de sorties lev_histday.nc # - lev_histday=2 => defaut # - lev_histday=3 => + champs 3D => F. Lott # - lev_histday=4 => + champs sous-surfaces # - lev_histday=5 => + champs F. Aires lev_histday=5 ### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" # - lev_histmth=0 => pas de sorties lev_histmth.nc # - lev_histmth=2 => defaut # - lev_histmth=3 => albedo, rugosite sous-surfaces # - lev_histmth=4 => champs tendances 3d lev_histmth=5 ### 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 ### ecrit_day = frequence ecriture fichier histday.nc en jours ecrit_day=1. ### ecrit_mth = frequence ecriture fichier histmth.nc en jours ecrit_mth=_ecrit-mth_ ### freqin_isccp = frequence input en secondes du simulateur ISCCP freq_ISCCP=10800. ### freqout_isccp = frequence output en jours du simulateur ISCCP ecrit_ISCCP=_ecrit-ISCCP_ ### niveau du diagnostique de conservation d energie if_ebil=0 # # parametres KE # ### epmax = Efficacite precipitation maximale epmax = .999 ### ok_adj_ema = ?? pas utilise ok_adj_ema = n ### iflag_clw Flag calcul eau liquide # - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l) # - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l) # - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip iflag_clw = 1 # # parametres nuages # ### cld_lc_lsc contenu en eau liquide des nuages large-scale (fisrtilp) cld_lc_lsc = 4.16e-4 ### cld_lc_con contenu en eau liquide des nuages convectifs (fisrtilp) cld_lc_con = 4.16e-4 ### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) cld_tau_lsc = 1800. ### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) cld_tau_con = 1800. ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) ffallv_lsc = 0.5 ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) ffallv_con = 0.5 ### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?) coef_eva = 2.e-5 ### reevap_ice reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) reevap_ice = y ### iflag_cldcon flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) # - iflag_cldcon<=-1 diag. rain_Tiedtke # - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k) # - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc) # - iflag_cldcon=3, ratqs=ratqss iflag_cldcon = 3 ### iflag_pdf : flag calcul distribution sous-maille de l eau et des nuages # - iflag_pdf=0, version ratqs, # - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir # - des PDFS de Sandrine Bony iflag_pdf = 1 ### fact_cldcon constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) fact_cldcon = 1. ### facttemps= facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3) facttemps = 0. ## ok_newmicro =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) ok_newmicro = y ### iflag_ratqs=0 correspond a la version IPCC AR4 iflag_ratqs=0 ### ratqsbas ratqs en bas si iflag_cldcon=1 ratqsbas = 0.005 ### ratqshaut ratqs en haut pour ratqss "stables" ratqshaut = 0.33 ### rad_froid rayon cristaux des nuages de glace "froids" rad_froid = 35 ### rad_chau1 rayon goutelettes d eau chauds", en haut: k=4-klev rad_chau1=12 ### rad_chau2 rayon goutelettes d eau chauds", en bas: k=1-3 rad_chau2=11 # # Coefficient et parametres sur les drags # f_cdrag_ter=1. f_cdrag_stable=1. f_cdrag_oce=0.8 f_rugoro=0. ### ksta_ter coef.diffusion minimale sur terre/sic/lic ksta_ter=1.e-7 ### cdmmax = cdrag maximum pour le moment cdmmax = 2.5E-3 ### cdhmax = cdrag maximum pour l energie cdhmax = 2.0E-3 # # Parametres "orbitaux/ ere geologique" # ### R_ecc = Excentricite R_ecc = 0.016715 ### R_peri = Equinoxe R_peri = 102.7 ### R_incl = Inclinaison R_incl = 23.441 ### solaire = Constante solaire solaire = 1365.6537 # # Taux gaz a effet de serre # ### co2_ppm = taux CO2 en ppm co2_ppm = 0.36886500E+03 ### CH4_ppb = taux CH4 en ppb CH4_ppb = 0.17510225E+04 ### N2O_ppb = taux N2O en ppb N2O_ppb = 0.31585000E+03 ### CFC11_ppt = taux CFC11 en ppt CFC11_ppt = 5.18015181E+01 ### CFC12_ppt = taux CFC12 en ppt CFC12_ppt = 0.99862742E+03 # # Parametres effets directs/indirects des "aerosols" # ### ok_ade=y/n flag Aerosol direct effect ok_ade=y ### ok_aie=y/n flag Aerosol indirect effect ok_aie=y ### aer_type = Aerosol variation type : actuel / preind / scenario / annuel aer_type=annuel ### 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=2 # # 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 ### 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