source: CMIP5/aqua4xCO2/v3cf.aqua4xCO2/PARAM/physiq.def_L39 @ 1481

Last change on this file since 1481 was 1481, checked in by musat, 13 years ago

v3cf.aqua4xCO2 IPSL CMIP5/CMIP2 simulation registration.

File size: 13.9 KB
Line 
1## $Id$
2#
3# PARAMETRES ANCIENNEMENT DANS gcm.def
4##  Cycle diurne  ou non                 
5cycle_diurne=y
6##  Soil Model  ou non               
7soil_model=y
8##  Choix ou non  de  New oliq               
9new_oliq=y
10##  Activation ou non de la parametrisation de Hines pour la strato
11ok_hines=y
12##  Orodr  ou  non   pour l orographie             
13ok_orodr=y
14##  Orolf  ou  non   pour l orographie             
15ok_orolf=y
16##   Si = .T. ,  lecture du fichier limit avec la bonne annee             
17ok_limitvrai=n
18## Nombre  d'appels des routines de rayonnements ( par jour)                 
19nbapp_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
21iflag_con=30
22## Facteur multiplication des precip convectives dans KE
23cvl_corr=1.0
24##  Facteur additif pour l'albedo
25pmagic=0.0
26#
27#
28#
29# Parametres fichiers de sortie
30#
31### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
32OK_journe=n
33### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc
34OK_mensuel=y
35### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la  physique)
36OK_instan=n
37### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc
38ok_hf=n
39# Noms des fichiers
40phys_out_filenames=      histmth histday histhf histhf3h histhf3hm histstn
41# Sortir ou non les fichiers
42phys_out_filekeys=       y       y       y       y       y        y
43# Niveaux de sorties
44phys_out_filelevels=     5       5       5       5       5        5
45### Type de fichier : global (n) ou stations (y)
46phys_out_filestations =  n       n       n       n       n        y
47### Frequences des sorties
48phys_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
50flag_t100       = 1, 7, 7, 5, 10, 10
51flag_u100       = 1, 7, 7, 5, 10, 10
52flag_v100       = 1, 7, 7, 5, 10, 10
53flag_w100       = 1, 7, 7, 5, 10, 10
54flag_t50        = 1, 7, 7, 5, 10, 10
55flag_u50        = 1, 7, 7, 5, 10, 10
56flag_v50        = 1, 7, 7, 5, 10, 10
57flag_w50        = 1, 7, 7, 5, 10, 10
58flag_t10        = 1, 7, 7, 5, 10, 10
59flag_u10        = 1, 7, 7, 5, 10, 10
60flag_v10        = 1, 7, 7, 5, 10, 10
61flag_w10        = 1, 7, 7, 5, 10, 10
62### niveaux maximum 5 pour variables CFMIP
63flag_t2m        =  1, 1, 1, 5, 10, 5
64flag_t2m_min    =  1, 1, 10, 5, 10, 10
65flag_t2m_max    =  1, 1, 10, 5, 10, 10
66flag_slp        =  1, 1, 1, 5, 10, 5
67flag_psol       =  1, 1, 1, 5, 10, 5
68flag_u10m       =  1, 1, 1, 5, 10, 5
69flag_v10m       =  1, 1, 1, 5, 10, 5
70flag_wind10m    =  1, 1, 1, 5, 10, 5
71flag_rh2m       =  5, 5, 10, 5, 10, 5
72flag_q2m        =  1, 1, 1, 5, 10, 5
73flag_precip     =  1, 1, 1, 5, 5, 5
74flag_snow       =  1, 1, 10, 5, 5, 5
75flag_pluc       =  1, 1, 1, 5, 5, 5
76flag_evap       =  1, 1, 10, 5, 10, 5
77flag_taux       =   1, 10, 10, 5, 10, 5
78flag_tauy       =  1, 10, 10, 5, 10, 5
79flag_flat       =  5, 1, 10, 5, 5, 5
80flag_sens       =  1, 1, 10, 5, 5, 5
81flag_LWdnSFC    =  1, 4, 10, 5, 5, 5
82flag_LWupSFCclr =  1, 4, 10, 5, 5, 5
83flag_SWdnSFC    =  1, 1, 10, 5, 5, 5
84flag_SWupSFC    =  1, 4, 10, 5, 5, 5
85flag_SWdnSFCclr =  1, 4, 10, 5, 5, 5
86flag_SWupSFCclr =  1, 4, 10, 5, 5, 5
87flag_LWdnSFCclr =  1, 4, 10, 5, 5, 5
88flag_SWdnTOA    =  1, 4, 10, 5, 5, 5
89flag_SWupTOA    =  1, 4, 10, 5, 5, 5
90flag_topl       =  1, 1, 10, 5, 5, 5
91flag_topl0      =  1, 5, 10, 5, 5, 5
92flag_SWupTOAclr =  1, 4, 10, 5, 5, 5
93flag_prw        =  1, 1, 10, 5, 5, 5
94flag_cldt       =  1, 1, 2, 5, 5, 5
95flag_lwp        =  1, 5, 10, 5, 5, 5
96flag_iwp        =  1, 5, 10, 5, 5, 5
97flag_nettop     =  1, 4, 10, 5, 5, 5
98flag_pbase      =  1, 5, 10, 5, 5, 5
99flag_ptop       = 1, 5, 10, 5, 5, 5
100flag_ftime_con  = 4, 10, 10, 5, 5, 5
101flag_zfull      =  2, 3, 10, 5, 5, 5
102flag_zhalf      = 2, 3, 10, 5, 5, 5
103flag_pres       =  2, 3, 10, 5, 5, 5
104flag_paprs      = 2, 3, 10, 5, 5, 5
105flag_temp       = 2, 3, 4, 5, 10, 5
106flag_h2o        =  4, 10, 10, 5, 10, 10
107flag_lwcon      =  2, 5, 10, 5, 10, 5
108flag_iwcon      =  2, 5, 10, 5, 10, 5
109flag_pr_con_l   =  2, 10, 10, 5, 10, 10
110flag_pr_lsc_l   =  2, 10, 10, 5, 10, 10
111flag_pr_con_i   =  2, 10, 10, 5, 10, 10
112flag_pr_lsc_i   =  2, 10, 10, 5, 10, 10
113flag_cldtau     =  10, 5, 10, 5, 10, 10
114flag_cldemi     =  10, 5, 10, 5, 10, 10
115flag_rlu        =  5, 10, 10, 10, 10, 5
116flag_rsu        =  5, 10, 10, 10, 10, 5
117flag_rld        =  5, 10, 10, 10, 10, 5
118flag_rsd        =  5, 10, 10, 10, 10, 5
119flag_rlucs      =  5, 10, 10, 10, 10, 5
120flag_rsucs      =  5, 10, 10, 10, 10, 5
121flag_rldcs      =  5, 10, 10, 10, 10, 5
122flag_rsdcs      =  5, 10, 10, 10, 10, 5
123flag_tnt        =  5, 10, 10, 10, 10, 5
124flag_dtdyn      = 4, 10, 10, 10, 10, 5
125flag_dtphy      = 4, 10, 10, 10, 10, 5
126flag_dtlschr    = 5, 10, 10, 10, 10, 5
127flag_tntr       = 5, 10, 10, 10, 10, 5
128flag_tntc       = 5, 10, 10, 10, 10, 5
129flag_tntscpbl   = 5, 10, 10, 10, 10, 5
130flag_rhum       = 2, 5, 10, 10, 10, 5
131flag_ovap       = 2, 3, 4, 10, 10, 5
132flag_tnhus      = 5, 10, 10, 10, 10, 5
133flag_dqdyn      = 4, 10, 10, 10, 10, 5
134flag_tnhusc     =  5, 10, 10, 10, 10, 5
135flag_dqvdf      = 5, 10, 10, 10, 10, 5
136flag_tnhusscpbl =  5, 10, 10, 10, 10, 5
137flag_dqphy      = 2, 10, 10, 10, 10, 5
138flag_evu        =  5, 10, 10, 10, 10, 5
139flag_kz         = 4, 10, 10, 10, 10, 5
140flag_rneb       = 2, 5, 10, 5, 10, 5
141flag_mc         = 4, 5, 10, 10, 10, 5
142flag_ref_liq    = 4, 10, 10, 10, 10, 5
143flag_ref_ice    = 4, 10, 10, 10, 10, 5
144flag_rsut4co2   = 10, 10, 10, 10, 10, 10
145flag_rlut4co2   = 10, 10, 10, 10, 10, 10
146flag_rsutcs4co2 = 10, 10, 10, 10, 10, 10
147flag_rlutcs4co2 = 10, 10, 10, 10, 10, 10
148flag_rlu4co2    = 10, 10, 10, 10, 10, 10
149flag_rsu4co2    = 10, 10, 10, 10, 10, 10
150flag_rld4co2    = 10, 10, 10, 10, 10, 10
151flag_rsd4co2    = 10, 10, 10, 10, 10, 10
152flag_rlucs4co2  = 10, 10, 10, 10, 10, 10
153flag_rsucs4co2  = 10, 10, 10, 10, 10, 10
154flag_rldcs4co2  = 10, 10, 10, 10, 10, 10
155flag_rsdcs4co2  = 10, 10, 10, 10, 10, 10
156flag_w500       = 1, 5, 7, 10, 10, 10
157flag_t700       = 1, 5, 7, 10, 10, 10
158flag_LWupSFC    = 1, 4, 10, 10, 5, 5
159flag_vitu       = 2, 3, 4, 6, 10, 5
160flag_vitv       = 2, 3, 4, 6, 10, 5
161flag_vitw       = 2, 3, 5, 5, 5, 5
162flag_geop       = 2, 3, 10, 10, 10, 5
163flag_tsol       =  1, 1, 1, 5, 5, 5
164###
165# more variables for t tendency
166###
167flag_dtvdf      = 5, 10, 10, 10, 10, 5
168flag_dteva      = 5, 10, 10, 10, 10, 5
169flag_dtcon      = 4, 10, 10, 10, 10, 5
170flag_dtajs      = 4, 10, 10, 10, 10, 5
171flag_dtlsc      = 4, 10, 10, 10, 10, 5
172flag_dtoro      =  4, 10, 10, 10, 10, 5
173flag_dtlif      =  4, 10, 10, 10, 10, 5
174flag_dthin      =  4, 10, 10, 10, 10, 5
175flag_dtec       =  4, 10, 10, 10, 10, 5
176### more variables for q tendency
177flag_dqeva      = 5, 10, 10, 10, 10, 5
178flag_dqcon      = 4, 10, 10, 10, 10, 5
179flag_dqajs      = 4, 10, 10, 10, 10, 5
180flag_dqlsc      = 5, 10, 10, 10, 10, 5
181### nouvelle physique
182flag_dtthe      = 4, 10, 10, 10, 10, 5
183flag_dtwak      = 4, 5, 10, 10, 10, 5
184flag_dqthe      = 4, 10, 10, 10, 10, 5
185flag_dqwak      = 4, 5, 10, 10, 10, 5
186flag_ftime_th   = 4, 10, 10, 5, 10, 5
187#
188# Parametres niveau de sorties differents fichiers
189#
190### lev_histhf=0-4, niveau de sortie fichier "histhf.nc"
191# - lev_histhf=0 => pas de sorties histhf.nc
192# - lev_histhf=2 => defaut
193# - lev_histhf=3 => variables sur niveaux standards
194# - lev_histhf=4 => champs 3d niveaux modele dans le fichier histhf.nc
195lev_histhf=5
196### lev_histday=0-5, niveau de sortie fichier "histday.nc"
197# - lev_histday=0 => pas de sorties lev_histday.nc
198# - lev_histday=2 => defaut
199# - lev_histday=3 => + champs 3D => F. Lott
200# - lev_histday=4 => + champs sous-surfaces
201# - lev_histday=5 => + champs F. Aires
202lev_histday=5
203### lev_histmth=0-4, niveau de sortie fichier "histmth.nc"
204# - lev_histmth=0 => pas de sorties lev_histmth.nc
205# - lev_histmth=2 => defaut
206# - lev_histmth=3 => albedo, rugosite sous-surfaces
207# - lev_histmth=4 => champs tendances 3d
208lev_histmth=5
209### parametres pour hist*NMC.nc
210#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
211#                histdayNMC: 8 ou 17
212lev_histdayNMC=8
213### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
214ok_histNMC = y, n, n
215### ecrit_LES = frequence ecriture fichier histhf3hm.nc en jours
216### ecrit_ins = frequence ecriture fichier histhf5h.nc en jours
217ecrit_ins=0.125
218### ecrit_hf = frequence ecriture fichier histhf.nc en jours
219ecrit_hf=0.250
220### ecrit_day = frequence ecriture fichier histday.nc en jours
221ecrit_day=1.
222### ecrit_mth = frequence ecriture fichier histmth.nc en jours
223ecrit_mth=_ecrit-mth_
224### freqin_isccp = frequence input en secondes du simulateur ISCCP
225freq_ISCCP=10800.
226### freqout_isccp = frequence output en jours du simulateur ISCCP
227ecrit_ISCCP=_ecrit-ISCCP_
228### niveau du diagnostique de conservation d energie
229if_ebil=0
230#
231# parametres KE
232#
233### epmax = Efficacite precipitation maximale
234epmax = .999
235### ok_adj_ema = ?? pas utilise
236ok_adj_ema = n
237### iflag_clw Flag calcul eau liquide
238# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
239# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
240# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
241iflag_clw = 1
242#
243# parametres nuages
244#
245### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp)
246cld_lc_lsc = 4.16e-4
247### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp)
248cld_lc_con = 4.16e-4
249### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp)
250cld_tau_lsc = 1800.
251### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp)
252cld_tau_con = 1800.
253### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp)
254ffallv_lsc = 0.5
255### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp)
256ffallv_con = 0.5
257### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
258coef_eva = 2.e-5
259### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp)
260reevap_ice = y
261### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq)
262# - iflag_cldcon<=-1 diag. rain_Tiedtke
263# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
264# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
265# - iflag_cldcon=3,   ratqs=ratqss
266iflag_cldcon = 3
267### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages
268# - iflag_pdf=0, version ratqs,
269# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
270# -              des PDFS de Sandrine Bony
271iflag_pdf = 1
272### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F)
273fact_cldcon = 1.
274### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
275facttemps = 0.
276## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages)
277ok_newmicro = y
278### iflag_ratqs=0 correspond a la version IPCC AR4
279iflag_ratqs=0
280### ratqsbas     ratqs en bas si iflag_cldcon=1
281ratqsbas = 0.005
282### ratqshaut    ratqs en haut pour ratqss "stables" 
283ratqshaut = 0.33
284### rad_froid    rayon cristaux des nuages de glace "froids"
285rad_froid = 35
286### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev
287rad_chau1=12
288### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3 
289rad_chau2=11
290#
291# Coefficient et parametres sur les drags
292#
293f_cdrag_ter=1.
294f_cdrag_stable=1.
295f_cdrag_oce=0.8
296f_rugoro=0.
297### ksta_ter      coef.diffusion minimale sur terre/sic/lic
298ksta_ter=1.e-7
299### cdmmax =     cdrag maximum pour le moment
300cdmmax = 2.5E-3
301### cdhmax =     cdrag maximum pour l energie
302cdhmax = 2.0E-3
303#
304# Parametres "orbitaux/ ere geologique"
305#
306### R_ecc =      Excentricite
307R_ecc = 0.
308### R_peri =     Equinoxe
309R_peri = 102.7
310### R_incl =     Inclinaison
311R_incl = 0.
312### solarlong0 = longitude solaire
313solarlong0 = 0.
314### solaire =    Constante solaire
315solaire = 1365.
316#
317# Taux gaz a effet de serre
318#
319### co2_ppm =    taux CO2 en ppm
320co2_ppm = 1.392E+03
321### CH4_ppb =    taux CH4 en ppb
322CH4_ppb = 0.1650E+04
323### N2O_ppb =    taux N2O en ppb
324N2O_ppb = 0.306E+03
325### CFC11_ppt =  taux CFC11 en ppt
326CFC11_ppt = 0.1E+01
327### CFC12_ppt =  taux CFC12 en ppt
328CFC12_ppt = 0.
329#
330# Taux gaz a effet de serre - 4xco2
331#
332### co2_ppm_per =    taux CO2 en ppm
333co2_ppm_per = 1.392E+03
334#
335# Parametres effets directs/indirects des "aerosols"
336#
337### ok_ade=y/n   flag Aerosol direct effect
338ok_ade=n
339### ok_aie=y/n   flag Aerosol indirect effect
340ok_aie=n
341### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
342aer_type=annuel
343###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
344flag_aerosol=6
345### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
346bl95_b0=1.7
347### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
348bl95_b1=0.2
349#
350# Parametre de lecture de l'ozone
351#
352# Allowed values are 0, 1 and 2
353# 0: do not read an ozone climatology
354# 1: read a single ozone climatology that will be used day and night
355# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
356read_climoz=1
357#
358# Parametres simulateur COSP (CFMIP Observational Simulator Package)
359#
360### ok_cosp=y/n flag simulateur COSP
361ok_cosp=n
362### freq_COSP = frequence d'appel de COSP en secondes
363freq_COSP=10800.
364### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
365ok_mensuelCOSP=y
366### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
367ok_journeCOSP=y
368### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
369ok_hfCOSP=n
370#
371# Parametres simulateur ISCCP
372#
373### ok_isccp=y/n flag simulateur ISCCP
374ok_isccp=n
375### top_height = flag choix calcul nuages par le simulateur en utilisant
376# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
377# - top_height = 1 -> algo IR-VIS
378# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
379# - top_height = 3 -> algo IR
380top_height = 1
381### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
382# - overlap=1    Max overlap
383# - overlap=2    Random overlap
384# - overlap=3    Max/Random overlap
385overlap = 3
386#
387qsol0=100.
Note: See TracBrowser for help on using the repository browser.