source: CONFIG/IPSLCM/IPSLCM5-2/CpiControl/PARAM/physiq.def_L39 @ 1539

## $Id$
#
# PARAMETRES ANCIENNEMENT DANS gcm.def
##  Cycle diurne  ou non                 
cycle_diurne=y
##  Soil Model  ou non               
soil_model=y
##  Choix ou non  de  New oliq               
new_oliq=y
##  Activation ou non de la parametrisation de Hines pour la strato
ok_hines=y
##  Orodr  ou  non   pour l orographie              
ok_orodr=y
##  Orolf  ou  non   pour l orographie              
ok_orolf=y
##   Si = .T. ,  lecture du fichier limit avec la bonne annee             
ok_limitvrai=n
## Nombre  d'appels des routines de rayonnements ( par jour)                 
nbapp_rad=24
##  Flag  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
iflag_con=30
## Facteur multiplication des precip convectives dans KE
cvl_corr=1.0
##  Facteur additif pour l'albedo
pmagic=0.0
#
#
#
# Parametres fichiers de sortie
#
# Noms des fichiers
phys_out_filenames=      histmth histday histhf   histhf3h  histhf3hm histstn
# Sortir ou non les fichiers
phys_out_filekeys=       n       n       n        n         n         n
# Niveaux de sorties
phys_out_filelevels=     5       5       2        5         5         5
### Type de fichier : global (n) ou stations (y)
phys_out_filestations =  n       n       n        n         n         y
### Frequences des sorties
phys_out_filetimesteps = 1.mth,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
#
#
### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott
flag_t100       = 1, 7, 7, 5, 10, 10
flag_u100       = 1, 7, 7, 5, 10, 10
flag_v100       = 1, 7, 7, 5, 10, 10
flag_w100       = 1, 7, 7, 5, 10, 10
flag_t50        = 1, 7, 7, 5, 10, 10
flag_u50        = 1, 7, 7, 5, 10, 10
flag_v50        = 1, 7, 7, 5, 10, 10
flag_w50        = 1, 7, 7, 5, 10, 10
flag_t10        = 1, 7, 7, 5, 10, 10
flag_u10        = 1, 7, 7, 5, 10, 10
flag_v10        = 1, 7, 7, 5, 10, 10
flag_w10        = 1, 7, 7, 5, 10, 10
# more variables for level 5 histday/month: vitw and pres
# Vertical wind
flag_vitw         =  2, 3, 7, 6, 6, 6
# Air pressure
flag_pres         =  2, 3, 7, 6, 6, 6
# w500 for histday
flag_w500         =  1, 5, 3, 10, 10, 10
# t700 for histday
flag_t700         =  1, 5, 3, 10, 10, 10
# Cloud fraction
flag_rneb         =  2, 5, 10, 5, 10, 10
# Aga sortie avec niveau 2
flag_Aga          =  2, 10, 10, 10, 10, 10
# niveaux maximum 5 mensuel pour variables CFMIP
flag_rlu          =  5, 10, 10, 10, 10, 10
flag_rsu          =  5, 10, 10, 10, 10, 10
flag_rld          =  5, 10, 10, 10, 10, 10
flag_rsd          =  5, 10, 10, 10, 10, 10
flag_rlucs        =  5, 10, 10, 10, 10, 10
flag_rsucs        =  5, 10, 10, 10, 10, 10
flag_rldcs        =  5, 10, 10, 10, 10, 10
flag_rsdcs        =  5, 10, 10, 10, 10, 10
flag_tnt          =  5, 10, 10, 10, 10, 10
flag_dtdyn        =  4, 10, 10, 10, 10, 10
flag_dtphy        =  4, 10, 10, 10, 10, 10
flag_dtlschr      =  5, 10, 10, 10, 10, 10
flag_tntr         =  5, 10, 10, 10, 10, 10
flag_tntc         =  5, 10, 10, 10, 10, 10
flag_tntscpbl     =  5, 10, 10, 10, 10, 10
flag_tnhus        =  5, 10, 10, 10, 10, 10
flag_dqdyn        =  4, 10, 10, 10, 10, 10
flag_tnhusc       =  5, 10, 10, 10, 10, 10
flag_dqvdf        =  5, 10, 10, 10, 10, 10
flag_tnhusscpbl   =  5, 10, 10, 10, 10, 10
flag_dqphy        =  2, 10, 10, 10, 10, 10
flag_evu          =  5, 10, 10, 10, 10, 10
flag_kz           =  4, 10, 10, 10, 10, 10
###
# more variables for t tendency monthly
###
flag_dtvdf        =  5, 10, 10, 10, 10, 10
flag_dteva        =  5, 10, 10, 10, 10, 10
flag_dtcon        =  4, 10, 10, 10, 10, 10
flag_dtajs        =  4, 10, 10, 10, 10, 10
flag_dtlsc        =  4, 10, 10, 10, 10, 10
flag_dtoro        =  4, 10, 10, 10, 10, 10
flag_dtlif        =  4, 10, 10, 10, 10, 10
flag_dthin        =  4, 10, 10, 10, 10, 10
flag_dtec         =  4, 10, 10, 10, 10, 10
### more variables for q tendency monthly
flag_dqeva        =  5, 10, 10, 10, 10, 10
flag_dqcon        =  4, 10, 10, 10, 10, 10
flag_dqajs        =  4, 10, 10, 10, 10, 10
flag_dqlsc        =  5, 10, 10, 10, 10, 10
### nouvelle physique monthly
flag_dtthe        =  4, 10, 10, 10, 10, 10
flag_dtwak        =  4, 5, 10, 10, 10, 10
flag_dqthe        =  4, 10, 10, 10, 10, 10
flag_dqwak        =  4, 5, 10, 10, 10, 10
flag_ftime_th     =  4, 10, 10, 5, 10, 10
#
# Parametres niveau de sorties differents fichiers
#
### parametres pour hist*NMC.nc
#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
#                histdayNMC: 8 ou 17
lev_histdayNMC=8
### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
ok_histNMC = y, n, n
### freqin_isccp = frequence input en secondes du simulateur ISCCP
freq_ISCCP=10800.
### freqout_isccp = frequence output en jours du simulateur ISCCP
ecrit_ISCCP=_ecrit-ISCCP_
### niveau du diagnostique de conservation d energie
if_ebil=0
#
# parametres KE
#
### epmax = Efficacite precipitation maximale
epmax = .999
### ok_adj_ema = ?? pas utilise
ok_adj_ema = n
### iflag_clw Flag calcul eau liquide
# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
iflag_clw = 1
# 
# parametres nuages
#
### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp) 
cld_lc_lsc = 4.16e-4
### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp) 
cld_lc_con = 4.16e-4
### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) 
cld_tau_lsc = 1800.
### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) 
cld_tau_con = 1800.
### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) 
ffallv_lsc = 0.5
### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) 
ffallv_con = 0.5
### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
coef_eva = 2.e-5
### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) 
reevap_ice = y
### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) 
# - iflag_cldcon<=-1 diag. rain_Tiedtke 
# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
# - iflag_cldcon=3,   ratqs=ratqss
iflag_cldcon = 3
### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages 
# - iflag_pdf=0, version ratqs, 
# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
# -              des PDFS de Sandrine Bony 
iflag_pdf = 1
### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) 
fact_cldcon = 1.
### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
facttemps = 0.
## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) 
ok_newmicro = y
### iflag_ratqs=0 correspond a la version IPCC AR4 
iflag_ratqs=0
### ratqsbas     ratqs en bas si iflag_cldcon=1 
ratqsbas = 0.005
### ratqshaut    ratqs en haut pour ratqss "stables"  
ratqshaut = 0.33
### rad_froid    rayon cristaux des nuages de glace "froids" 
rad_froid = 35
### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev 
rad_chau1=12
### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3  
rad_chau2=11
# 
# Coefficient et parametres sur les drags
#
f_cdrag_ter=1.
f_cdrag_stable=1.
f_cdrag_oce=0.8
f_rugoro=0.
### ksta_ter      coef.diffusion minimale sur terre/sic/lic 
ksta_ter=1.e-7
### cdmmax =     cdrag maximum pour le moment
cdmmax = 2.5E-3
### cdhmax =     cdrag maximum pour l energie
cdhmax = 2.0E-3
#
# Parametres "orbitaux/ ere geologique"
#
### R_ecc =      Excentricite
R_ecc = 0.016715
### R_peri =     Equinoxe
R_peri = 102.7
### R_incl =     Inclinaison
R_incl = 23.441
### solaire =    Constante solaire
solaire = 1365.6537
#
# Taux gaz a effet de serre
#
### co2_ppm =    taux CO2 en ppm
co2_ppm = 0.28472500E+03
### CH4_ppb =    taux CH4 en ppb
CH4_ppb = 0.79097924E+03
### N2O_ppb =    taux N2O en ppb
N2O_ppb = 0.27542506E+03
### CFC11_ppt =  taux CFC11 en ppt
CFC11_ppt = 0.
### CFC12_ppt =  taux CFC12 en ppt
CFC12_ppt = 0.
#
# Parametres effets directs/indirects des "aerosols"
#
### ok_ade=y/n   flag Aerosol direct effect
ok_ade=y
### ok_aie=y/n   flag Aerosol indirect effect
ok_aie=y
### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
aer_type=preind
###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
flag_aerosol=6 
### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
bl95_b0=1.7
### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
bl95_b1=0.2
#
# Parametre de lecture de l'ozone
#
# Allowed values are 0, 1 and 2
# 0: do not read an ozone climatology
# 1: read a single ozone climatology that will be used day and night
# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
read_climoz=2
#
# Parametres simulateur COSP (CFMIP Observational Simulator Package)
#
### ok_cosp=y/n flag simulateur COSP
ok_cosp=n
### freq_COSP = frequence d'appel de COSP en secondes
freq_COSP=10800.
### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
ok_mensuelCOSP=y
### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
ok_journeCOSP=y
### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
ok_hfCOSP=n
#
# Parametres simulateur ISCCP
#
### ok_isccp=y/n flag simulateur ISCCP
ok_isccp=n
### top_height = flag choix calcul nuages par le simulateur en utilisant 
# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
# - top_height = 1 -> algo IR-VIS
# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
# - top_height = 3 -> algo IR
top_height = 1
### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
# - overlap=1    Max overlap
# - overlap=2    Random overlap
# - overlap=3    Max/Random overlap
overlap = 3