source: CONFIG/LMDZOR/branches/LMDZOR_v4/LMDZOR/PARAM/physiq.def_L39 @ 1203

Last change on this file since 1203 was 1203, checked in by jgipsl, 14 years ago

Dev config LMDZOR_v4 :

Towards possibility to run LMDZOR experience as climate or AMIP
interannuel experience. All PARAM/*def and COMP/*driver are the same. To
switch change lmdz.card and orchidee.card.

  • Modified following parameters from lmdz.card and orchidee.card : VEGET_UPDATE, LAND_COVER_CHANGE, ok_guide, pmagic, aer_type
  • Added orchidee.card_AMIP, lmdz.card_AMIP : copied from AMIP/COMP
  • Uniformized lmdz.driver and orchidee.driver with functionality from AMIP experience
  • orchidee.driver : added function SRF_PeriodStart, desactaivate river routing for high resolution (280x280x39)
  • lmdz.driver : added lecture of GHG forcing files if they exist (if added in lmdz.card), desactivated fft filter in case of resolution

280x280x39 which contains zoom, added test if Bands file exist
(LMDZ_Bands_file_name)

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