## $Id: physiq.def 1300 2010-01-20 15:31:14Z fairhead $
#
# 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=n
##  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.008
#
#
#
# Parametres fichiers de sortie
#
### type_run = type run par rapport aux fichiers et variables de sortie
# - type_run = CLIM/ENSP (=1)
# - type_run = AMIP/CFMI (=2)
type_run=AMIP
### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
OK_journe=n
### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc
OK_mensuel=y
### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la  physique)
OK_instan=n
### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc
ok_hf=n
#
# Parametres niveau de sorties differents fichiers 
#
### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" 
# - lev_histhf=0 => pas de sorties histhf.nc
# - lev_histhf=2 => defaut
# - lev_histhf=3 => variables sur niveaux standards
# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc
lev_histhf=2
### lev_histday=0-5, niveau de sortie fichier "histday.nc"
# - lev_histday=0 => pas de sorties lev_histday.nc
# - lev_histday=2 => defaut
# - lev_histday=3 => + champs 3D => F. Lott
# - lev_histday=4 => + champs sous-surfaces
# - lev_histday=5 => + champs F. Aires
lev_histday=2
### lev_histmth=0-4, niveau de sortie fichier "histmth.nc"
# - lev_histmth=0 => pas de sorties lev_histmth.nc
# - lev_histmth=2 => defaut
# - lev_histmth=3 => albedo, rugosite sous-surfaces
# - lev_histmth=4 => champs tendances 3d
lev_histmth=2
### ecrit_hf = frequence ecriture fichier histhf.nc en jours
ecrit_hf=0.250
### ecrit_day = frequence ecriture fichier histday.nc en jours
ecrit_day=1.
### ecrit_mth = frequence ecriture fichier histmth.nc en jours
ecrit_mth=30
### freqin_isccp = frequence input en secondes du simulateur ISCCP
freq_ISCCP=10800.
### freqout_isccp = frequence output en jours du simulateur ISCCP
ecrit_ISCCP=30
### 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 = 1361.
#
# Taux gaz a effet de serre
#
### co2_ppm =    taux CO2 en ppm
co2_ppm = 367.
### CH4_ppb =    taux CH4 en ppb
CH4_ppb = 1760.
### N2O_ppb =    taux N2O en ppb
N2O_ppb = 316.
### CFC11_ppt =  taux CFC11 en ppt
CFC11_ppt = 741.2
### CFC12_ppt =  taux CFC12 en ppt
CFC12_ppt = 535.
#
# Parametres effets directs/indirects des "aerosols"
#
### ok_ade=y/n   flag Aerosol direct effect
ok_ade=n
### ok_aie=y/n   flag Aerosol indirect effect
ok_aie=n
### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
aer_type=actuel
###  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=0
#
# 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.
#
# 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
