source: CONFIG/LMDZREPR/branches/LMDZREPR_v1/EXP/PARAM/physiq.def_L39_AP @ 1599

Last change on this file since 1599 was 1599, checked in by jgipsl, 13 years ago

Added experience directory EXP for running LMDZ-REPROBUS. EXP is under work...
Input files for REPROBUS are not yet available at CCRT. To run this experience you need libIGCM/trunk with a modification in libIGCM_config.ksh not yet available at svn.

File size: 7.4 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
25# variable modife avec lmdz.card
26pmagic=999.9
27#
28#
29# Parametres fichiers de sortie
30#
31# Noms des fichiers
32phys_out_filenames=      histmth histday histhf  histhf3h histhf3hm histstn
33# Sortir ou non les fichiers
34phys_out_filekeys=       y       y       y       y       y          n
35# Niveaux de sorties
36phys_out_filelevels=     5       5       5       5       5          5
37### Type de fichier : global (n) ou stations (y)
38phys_out_filestations =  n       n       n       n       n          y
39### Frequences des sorties
40phys_out_filetimesteps = 1.mth,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
41
42
43# more variables for level 5 : vitw and pres
44# Vertical wind
45flag_vitw         =  2, 3, 5, 5, 1
46name_vitw         =  vitw
47# Air pressure
48flag_pres         =  2, 3, 5, 5, 1
49name_pres         =  pres
50#
51### parametres pour hist*NMC.nc
52#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
53#                histdayNMC: 8 ou 17
54lev_histdayNMC=8
55### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
56ok_histNMC = y, n, n
57
58### freqin_isccp = frequence input en secondes du simulateur ISCCP
59freq_ISCCP=10800.
60### freqout_isccp = frequence output en jours du simulateur ISCCP
61ecrit_ISCCP=_ecrit-ISCCP_
62### niveau du diagnostique de conservation d energie
63if_ebil=0
64#
65# parametres KE
66#
67### epmax = Efficacite precipitation maximale
68epmax = .999
69### ok_adj_ema = ?? pas utilise
70ok_adj_ema = n
71### iflag_clw Flag calcul eau liquide
72# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
73# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
74# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
75iflag_clw = 1
76#
77# parametres nuages
78#
79### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp)
80cld_lc_lsc = 4.16e-4
81### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp)
82cld_lc_con = 4.16e-4
83### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp)
84cld_tau_lsc = 1800.
85### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp)
86cld_tau_con = 1800.
87### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp)
88ffallv_lsc = 0.5
89### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp)
90ffallv_con = 0.5
91### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
92coef_eva = 2.e-5
93### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp)
94reevap_ice = y
95### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq)
96# - iflag_cldcon<=-1 diag. rain_Tiedtke
97# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
98# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
99# - iflag_cldcon=3,   ratqs=ratqss
100iflag_cldcon = 3
101### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages
102# - iflag_pdf=0, version ratqs,
103# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
104# -              des PDFS de Sandrine Bony
105iflag_pdf = 1
106### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F)
107fact_cldcon = 1.
108### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
109facttemps = 0.
110## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages)
111ok_newmicro = y
112### iflag_ratqs=0 correspond a la version IPCC AR4
113iflag_ratqs=0
114### ratqsbas     ratqs en bas si iflag_cldcon=1
115ratqsbas = 0.005
116### ratqshaut    ratqs en haut pour ratqss "stables" 
117ratqshaut = 0.33
118### rad_froid    rayon cristaux des nuages de glace "froids"
119rad_froid = 35
120### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev
121rad_chau1=12
122### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3 
123rad_chau2=11
124#
125# Coefficient et parametres sur les drags
126#
127f_cdrag_ter=1.
128f_cdrag_stable=1.
129f_cdrag_oce=0.8
130f_rugoro=0.
131### ksta_ter      coef.diffusion minimale sur terre/sic/lic
132ksta_ter=1.e-7
133### cdmmax =     cdrag maximum pour le moment
134cdmmax = 2.5E-3
135### cdhmax =     cdrag maximum pour l energie
136cdhmax = 2.0E-3
137#
138# Parametres "orbitaux/ ere geologique"
139#
140### R_ecc =      Excentricite
141R_ecc = 0.016715
142### R_peri =     Equinoxe
143R_peri = 102.7
144### R_incl =     Inclinaison
145R_incl = 23.441
146### solaire =    Constante solaire
147solaire = 1366.0896
148#
149# Taux gaz a effet de serre
150#
151### co2_ppm =    taux CO2 en ppm
152co2_ppm = 0.36886500E+03
153### CH4_ppb =    taux CH4 en ppb
154CH4_ppb = 0.17510225E+04
155### N2O_ppb =    taux N2O en ppb
156N2O_ppb = 0.31585000E+03
157### CFC11_ppt =  taux CFC11 en ppt
158CFC11_ppt = 5.18015181E+01
159### CFC12_ppt =  taux CFC12 en ppt
160CFC12_ppt = 0.99862742E+03
161#
162# Parametres effets directs/indirects des "aerosols"
163#
164### ok_ade=y/n   flag Aerosol direct effect
165ok_ade=y
166### ok_aie=y/n   flag Aerosol indirect effect
167ok_aie=y
168### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
169# variable modife avec lmdz.card
170aer_type=xxx
171###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
172flag_aerosol=6
173### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
174bl95_b0=1.7
175### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
176bl95_b1=0.2
177#
178# Parametre de lecture de l'ozone
179#
180# Allowed values are 0, 1 and 2
181# 0: do not read an ozone climatology
182# 1: read a single ozone climatology that will be used day and night
183# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
184read_climoz=2
185#
186# Parametres simulateur COSP (CFMIP Observational Simulator Package)
187#
188### ok_cosp=y/n flag simulateur COSP
189ok_cosp=n
190### freq_COSP = frequence d'appel de COSP en secondes
191freq_COSP=10800.
192### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
193ok_mensuelCOSP=y
194### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
195ok_journeCOSP=y
196### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
197ok_hfCOSP=n
198#
199# Parametres simulateur ISCCP
200#
201### ok_isccp=y/n flag simulateur ISCCP
202ok_isccp=n
203### top_height = flag choix calcul nuages par le simulateur en utilisant
204# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
205# - top_height = 1 -> algo IR-VIS
206# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
207# - top_height = 3 -> algo IR
208top_height = 1
209### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
210# - overlap=1    Max overlap
211# - overlap=2    Random overlap
212# - overlap=3    Max/Random overlap
213overlap = 3
Note: See TracBrowser for help on using the repository browser.