source: LMDZ4/trunk/physiq.def @ 2061

Last change on this file since 2061 was 1403, checked in by Laurent Fairhead, 14 years ago

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

  • IPSLCM5A configuration (comparison between trunk and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent
  • MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

  • IPSLCM5A configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent
  • MH07 configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files


Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

  • pour une config IPSLCM5A (comparaison tronc/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux (à part sur RN et Pb)
    • fichiers histoire égaux
  • pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux
    • fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

  • pour une config IPSLCM5A:
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc
  • pour une config MH07
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 7.5 KB
Line 
1## $Id: physiq.def 1403 2010-07-01 09:02:53Z fairhead $
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=n
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.008
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)
34type_run=AMIP
35### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
36OK_journe=n
37### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc
38OK_mensuel=y
39### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la  physique)
40OK_instan=n
41### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc
42ok_hf=n
43#
44# Parametres niveau de sorties differents fichiers
45#
46### lev_histhf=0-4, niveau de sortie fichier "histhf.nc"
47# - lev_histhf=0 => pas de sorties histhf.nc
48# - lev_histhf=2 => defaut
49# - lev_histhf=3 => variables sur niveaux standards
50# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc
51lev_histhf=2
52### lev_histday=0-5, niveau de sortie fichier "histday.nc"
53# - lev_histday=0 => pas de sorties lev_histday.nc
54# - lev_histday=2 => defaut
55# - lev_histday=3 => + champs 3D => F. Lott
56# - lev_histday=4 => + champs sous-surfaces
57# - lev_histday=5 => + champs F. Aires
58lev_histday=2
59### lev_histmth=0-4, niveau de sortie fichier "histmth.nc"
60# - lev_histmth=0 => pas de sorties lev_histmth.nc
61# - lev_histmth=2 => defaut
62# - lev_histmth=3 => albedo, rugosite sous-surfaces
63# - lev_histmth=4 => champs tendances 3d
64lev_histmth=2
65### ecrit_hf = frequence ecriture fichier histhf.nc en jours
66ecrit_hf=0.250
67### ecrit_day = frequence ecriture fichier histday.nc en jours
68ecrit_day=1.
69### ecrit_mth = frequence ecriture fichier histmth.nc en jours
70ecrit_mth=30
71### freqin_isccp = frequence input en secondes du simulateur ISCCP
72freq_ISCCP=10800.
73### freqout_isccp = frequence output en jours du simulateur ISCCP
74ecrit_ISCCP=30
75### niveau du diagnostique de conservation d energie
76if_ebil=0
77#
78# parametres KE
79#
80### epmax = Efficacite precipitation maximale
81epmax = .999
82### ok_adj_ema = ?? pas utilise
83ok_adj_ema = n
84### iflag_clw Flag calcul eau liquide
85# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
86# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
87# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
88iflag_clw = 1
89#
90# parametres nuages
91#
92### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp)
93cld_lc_lsc = 4.16e-4
94### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp)
95cld_lc_con = 4.16e-4
96### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp)
97cld_tau_lsc = 1800.
98### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp)
99cld_tau_con = 1800.
100### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp)
101ffallv_lsc = 0.5
102### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp)
103ffallv_con = 0.5
104### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
105coef_eva = 2.e-5
106### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp)
107reevap_ice = y
108### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq)
109# - iflag_cldcon<=-1 diag. rain_Tiedtke
110# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
111# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
112# - iflag_cldcon=3,   ratqs=ratqss
113iflag_cldcon = 3
114### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages
115# - iflag_pdf=0, version ratqs,
116# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
117# -              des PDFS de Sandrine Bony
118iflag_pdf = 1
119### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F)
120fact_cldcon = 1.
121### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
122facttemps = 0.
123## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages)
124ok_newmicro = y
125### iflag_ratqs=0 correspond a la version IPCC AR4
126iflag_ratqs=0
127### ratqsbas     ratqs en bas si iflag_cldcon=1
128ratqsbas = 0.005
129### ratqshaut    ratqs en haut pour ratqss "stables" 
130ratqshaut = 0.33
131### rad_froid    rayon cristaux des nuages de glace "froids"
132rad_froid = 35
133### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev
134rad_chau1=12
135### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3 
136rad_chau2=11
137#
138# Coefficient et parametres sur les drags
139#
140f_cdrag_ter=1.
141f_cdrag_stable=1.
142f_cdrag_oce=0.8
143f_rugoro=0.
144### ksta_ter      coef.diffusion minimale sur terre/sic/lic
145ksta_ter=1.e-7
146### cdmmax =     cdrag maximum pour le moment
147cdmmax = 2.5E-3
148### cdhmax =     cdrag maximum pour l energie
149cdhmax = 2.0E-3
150#
151# Parametres "orbitaux/ ere geologique"
152#
153### R_ecc =      Excentricite
154R_ecc = 0.016715
155### R_peri =     Equinoxe
156R_peri = 102.7
157### R_incl =     Inclinaison
158R_incl = 23.441
159### solaire =    Constante solaire
160solaire = 1361.
161#
162# Taux gaz a effet de serre
163#
164### co2_ppm =    taux CO2 en ppm
165co2_ppm = 367.
166### CH4_ppb =    taux CH4 en ppb
167CH4_ppb = 1760.
168### N2O_ppb =    taux N2O en ppb
169N2O_ppb = 316.
170### CFC11_ppt =  taux CFC11 en ppt
171CFC11_ppt = 741.2
172### CFC12_ppt =  taux CFC12 en ppt
173CFC12_ppt = 535.
174#
175# Parametres effets directs/indirects des "aerosols"
176#
177### ok_ade=y/n   flag Aerosol direct effect
178ok_ade=n
179### ok_aie=y/n   flag Aerosol indirect effect
180ok_aie=n
181### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
182aer_type=actuel
183###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
184flag_aerosol=6
185### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
186bl95_b0=1.7
187### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
188bl95_b1=0.2
189#
190# Parametre de lecture de l'ozone
191#
192# Allowed values are 0, 1 and 2
193# 0: do not read an ozone climatology
194# 1: read a single ozone climatology that will be used day and night
195# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
196read_climoz=0
197#
198# Parametres simulateur COSP (CFMIP Observational Simulator Package)
199#
200### ok_cosp=y/n flag simulateur COSP
201ok_cosp=n
202## freq_COSP = frequence d'appel de COSP en secondes
203freq_COSP=10800.
204#
205# Parametres simulateur ISCCP
206#
207### ok_isccp=y/n flag simulateur ISCCP
208ok_isccp=n
209### top_height = flag choix calcul nuages par le simulateur en utilisant
210# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
211# - top_height = 1 -> algo IR-VIS
212# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
213# - top_height = 3 -> algo IR
214top_height = 1
215### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
216# - overlap=1    Max overlap
217# - overlap=2    Random overlap
218# - overlap=3    Max/Random overlap
219overlap = 3
Note: See TracBrowser for help on using the repository browser.