1 | ## $Id$ |
---|
2 | # |
---|
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 |
---|
22 | ## Facteur multiplication des precip convectives dans KE |
---|
23 | cvl_corr=1.0 |
---|
24 | ## Facteur additif pour l'albedo |
---|
25 | pmagic=0.0 |
---|
26 | # |
---|
27 | # |
---|
28 | # |
---|
29 | # Parametres fichiers de sortie |
---|
30 | # |
---|
31 | ### type_run = type run par rapport aux fichiers et variables de sortie |
---|
32 | # - type_run = CLIM/ENSP (=1) |
---|
33 | # - type_run = AMIP/CFMI (=2) |
---|
34 | type_run=CLIM |
---|
35 | ### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc |
---|
36 | OK_journe=n |
---|
37 | ### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc |
---|
38 | OK_mensuel=y |
---|
39 | ### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la physique) |
---|
40 | OK_instan=n |
---|
41 | ### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc |
---|
42 | ok_hf=n |
---|
43 | # |
---|
44 | # more variables for level 5 : vitw and pres |
---|
45 | # Vertical wind |
---|
46 | flag_vitw = 2, 3, 7, 6, 6 |
---|
47 | name_vitw = vitw |
---|
48 | # Air pressure |
---|
49 | flag_pres = 2, 3, 7, 6, 6 |
---|
50 | name_pres = pres |
---|
51 | # w500 for histday |
---|
52 | flag_w500 = 1, 5, 3, 10, 10 |
---|
53 | name_w500 = w500 |
---|
54 | # t700 for histday |
---|
55 | flag_t700 = 1, 5, 3, 10, 10 |
---|
56 | name_t700 = t700 |
---|
57 | # Cloud fraction |
---|
58 | flag_rneb = 2, 5, 10, 5, 10 |
---|
59 | name_rneb = rneb |
---|
60 | # Parametres niveau de sorties differents fichiers |
---|
61 | # |
---|
62 | ### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" |
---|
63 | # - lev_histhf=0 => pas de sorties histhf.nc |
---|
64 | # - lev_histhf=2 => defaut |
---|
65 | # - lev_histhf=3 => variables sur niveaux standards |
---|
66 | # - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc |
---|
67 | lev_histhf=2 |
---|
68 | ### lev_histday=0-5, niveau de sortie fichier "histday.nc" |
---|
69 | # - lev_histday=0 => pas de sorties lev_histday.nc |
---|
70 | # - lev_histday=2 => defaut |
---|
71 | # - lev_histday=3 => + champs 3D => F. Lott |
---|
72 | # - lev_histday=4 => + champs sous-surfaces |
---|
73 | # - lev_histday=5 => + champs F. Aires |
---|
74 | lev_histday=5 |
---|
75 | ### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" |
---|
76 | # - lev_histmth=0 => pas de sorties lev_histmth.nc |
---|
77 | # - lev_histmth=2 => defaut |
---|
78 | # - lev_histmth=3 => albedo, rugosite sous-surfaces |
---|
79 | # - lev_histmth=4 => champs tendances 3d |
---|
80 | lev_histmth=5 |
---|
81 | ### Aga sortie avec niveau 2 |
---|
82 | flag_Aga = 2, 10, 10, 10, 10 |
---|
83 | ### parametres pour hist*NMC.nc |
---|
84 | #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier |
---|
85 | # histdayNMC: 8 ou 17 |
---|
86 | lev_histdayNMC=8 |
---|
87 | ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence |
---|
88 | ok_histNMC = y, n, n |
---|
89 | ### ecrit_ins = frequence ecriture fichier histhf3h.nc en jours |
---|
90 | ecrit_ins=0.125 |
---|
91 | ### ecrit_hf = frequence ecriture fichier histhf.nc en jours |
---|
92 | ecrit_hf=0.250 |
---|
93 | ### ecrit_day = frequence ecriture fichier histday.nc en jours |
---|
94 | ecrit_day=1. |
---|
95 | ### ecrit_mth = frequence ecriture fichier histmth.nc en jours |
---|
96 | ecrit_mth=_ecrit-mth_ |
---|
97 | ### freqin_isccp = frequence input en secondes du simulateur ISCCP |
---|
98 | freq_ISCCP=10800. |
---|
99 | ### freqout_isccp = frequence output en jours du simulateur ISCCP |
---|
100 | ecrit_ISCCP=_ecrit-ISCCP_ |
---|
101 | ### niveau du diagnostique de conservation d energie |
---|
102 | if_ebil=0 |
---|
103 | # |
---|
104 | # parametres KE |
---|
105 | # |
---|
106 | ### epmax = Efficacite precipitation maximale |
---|
107 | epmax = .999 |
---|
108 | ### ok_adj_ema = ?? pas utilise |
---|
109 | ok_adj_ema = n |
---|
110 | ### iflag_clw Flag calcul eau liquide |
---|
111 | # - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l) |
---|
112 | # - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l) |
---|
113 | # - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip |
---|
114 | iflag_clw = 1 |
---|
115 | # |
---|
116 | # parametres nuages |
---|
117 | # |
---|
118 | ### cld_lc_lsc contenu en eau liquide des nuages large-scale (fisrtilp) |
---|
119 | cld_lc_lsc = 4.16e-4 |
---|
120 | ### cld_lc_con contenu en eau liquide des nuages convectifs (fisrtilp) |
---|
121 | cld_lc_con = 4.16e-4 |
---|
122 | ### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) |
---|
123 | cld_tau_lsc = 1800. |
---|
124 | ### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) |
---|
125 | cld_tau_con = 1800. |
---|
126 | ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) |
---|
127 | ffallv_lsc = 0.5 |
---|
128 | ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) |
---|
129 | ffallv_con = 0.5 |
---|
130 | ### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?) |
---|
131 | coef_eva = 2.e-5 |
---|
132 | ### reevap_ice reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) |
---|
133 | reevap_ice = y |
---|
134 | ### iflag_cldcon flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) |
---|
135 | # - iflag_cldcon<=-1 diag. rain_Tiedtke |
---|
136 | # - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k) |
---|
137 | # - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc) |
---|
138 | # - iflag_cldcon=3, ratqs=ratqss |
---|
139 | iflag_cldcon = 3 |
---|
140 | ### iflag_pdf : flag calcul distribution sous-maille de l eau et des nuages |
---|
141 | # - iflag_pdf=0, version ratqs, |
---|
142 | # - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir |
---|
143 | # - des PDFS de Sandrine Bony |
---|
144 | iflag_pdf = 1 |
---|
145 | ### fact_cldcon constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) |
---|
146 | fact_cldcon = 1. |
---|
147 | ### facttemps= facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3) |
---|
148 | facttemps = 0. |
---|
149 | ## ok_newmicro =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) |
---|
150 | ok_newmicro = y |
---|
151 | ### iflag_ratqs=0 correspond a la version IPCC AR4 |
---|
152 | iflag_ratqs=0 |
---|
153 | ### ratqsbas ratqs en bas si iflag_cldcon=1 |
---|
154 | ratqsbas = 0.005 |
---|
155 | ### ratqshaut ratqs en haut pour ratqss "stables" |
---|
156 | ratqshaut = 0.33 |
---|
157 | ### rad_froid rayon cristaux des nuages de glace "froids" |
---|
158 | rad_froid = 35 |
---|
159 | ### rad_chau1 rayon goutelettes d eau chauds", en haut: k=4-klev |
---|
160 | rad_chau1=12 |
---|
161 | ### rad_chau2 rayon goutelettes d eau chauds", en bas: k=1-3 |
---|
162 | rad_chau2=11 |
---|
163 | # |
---|
164 | # Coefficient et parametres sur les drags |
---|
165 | # |
---|
166 | f_cdrag_ter=1. |
---|
167 | f_cdrag_stable=1. |
---|
168 | f_cdrag_oce=0.8 |
---|
169 | f_rugoro=0. |
---|
170 | ### ksta_ter coef.diffusion minimale sur terre/sic/lic |
---|
171 | ksta_ter=1.e-7 |
---|
172 | ### cdmmax = cdrag maximum pour le moment |
---|
173 | cdmmax = 2.5E-3 |
---|
174 | ### cdhmax = cdrag maximum pour l energie |
---|
175 | cdhmax = 2.0E-3 |
---|
176 | # |
---|
177 | # Parametres "orbitaux/ ere geologique" |
---|
178 | # |
---|
179 | ### R_ecc = Excentricite |
---|
180 | R_ecc = 0.016715 |
---|
181 | ### R_peri = Equinoxe |
---|
182 | R_peri = 102.7 |
---|
183 | ### R_incl = Inclinaison |
---|
184 | R_incl = 23.441 |
---|
185 | ### solaire = Constante solaire |
---|
186 | solaire = 1365.6537 |
---|
187 | # |
---|
188 | # Taux gaz a effet de serre |
---|
189 | # |
---|
190 | ### co2_ppm = taux CO2 en ppm |
---|
191 | co2_ppm = 0.28472500E+03 |
---|
192 | ### CH4_ppb = taux CH4 en ppb |
---|
193 | CH4_ppb = 0.79097924E+03 |
---|
194 | ### N2O_ppb = taux N2O en ppb |
---|
195 | N2O_ppb = 0.27542506E+03 |
---|
196 | ### CFC11_ppt = taux CFC11 en ppt |
---|
197 | CFC11_ppt = 0. |
---|
198 | ### CFC12_ppt = taux CFC12 en ppt |
---|
199 | CFC12_ppt = 0. |
---|
200 | # |
---|
201 | # Parametres effets directs/indirects des "aerosols" |
---|
202 | # |
---|
203 | ### ok_ade=y/n flag Aerosol direct effect |
---|
204 | ok_ade=y |
---|
205 | ### ok_aie=y/n flag Aerosol indirect effect |
---|
206 | ok_aie=y |
---|
207 | ### aer_type = Aerosol variation type : actuel / preind / scenario / annuel |
---|
208 | aer_type=preind |
---|
209 | ### type of coupled aerosol =1 (default) =2 => bc only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol |
---|
210 | flag_aerosol=6 |
---|
211 | ### bl95_b0 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
---|
212 | bl95_b0=1.7 |
---|
213 | ### bl95_b1 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
---|
214 | bl95_b1=0.2 |
---|
215 | # |
---|
216 | # Parametre de lecture de l'ozone |
---|
217 | # |
---|
218 | # Allowed values are 0, 1 and 2 |
---|
219 | # 0: do not read an ozone climatology |
---|
220 | # 1: read a single ozone climatology that will be used day and night |
---|
221 | # 2: read two ozone climatologies, the average day and night climatology and the daylight climatology |
---|
222 | read_climoz=2 |
---|
223 | # |
---|
224 | # Parametres simulateur COSP (CFMIP Observational Simulator Package) |
---|
225 | # |
---|
226 | ### ok_cosp=y/n flag simulateur COSP |
---|
227 | ok_cosp=n |
---|
228 | ### freq_COSP = frequence d'appel de COSP en secondes |
---|
229 | freq_COSP=10800. |
---|
230 | ### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc |
---|
231 | ok_mensuelCOSP=y |
---|
232 | ### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc |
---|
233 | ok_journeCOSP=y |
---|
234 | ### ok_hfCOSP=y/n, ecrire sorties haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc |
---|
235 | ok_hfCOSP=n |
---|
236 | # |
---|
237 | # Parametres simulateur ISCCP |
---|
238 | # |
---|
239 | ### ok_isccp=y/n flag simulateur ISCCP |
---|
240 | ok_isccp=n |
---|
241 | ### top_height = flag choix calcul nuages par le simulateur en utilisant |
---|
242 | # - les donnees IR et/ou VIS et l algorithme ISCCP-D1 |
---|
243 | # - top_height = 1 -> algo IR-VIS |
---|
244 | # - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" |
---|
245 | # - top_height = 3 -> algo IR |
---|
246 | top_height = 1 |
---|
247 | ### overlap = Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP |
---|
248 | # - overlap=1 Max overlap |
---|
249 | # - overlap=2 Random overlap |
---|
250 | # - overlap=3 Max/Random overlap |
---|
251 | overlap = 3 |
---|