##############################################################################
# PARAMETRES ET CLES DE CONFIGURATION LMDZ
##############################################################################
# D: Valeur par default
#
# Controle des sorties
#*********************
OK_journe=y
OK_mensuel=y
ok_hf=n
OK_instan=n
#ok_regdyn : y/n calcul/non des regymes dynamiques sur regions pre-definies
ok_regdyn=n
# - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
ok_histNMC = n, n, n
#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
#                histdayNMC: 8 ou 17
lev_histdayNMC=8
# A transformer en la ligne ci-dessous quand ce sera compatible avec libIGCM
phys_out_regfkey=n       n       n       n       n      n
phys_out_filekeys=n       n       n       n       n      n    
phys_out_filenames=      histmth histday histhf  histhf2h histLES histstn
phys_out_filelevels=     10      4       4       0       1
phys_out_filetypes=      ave(X)  ave(X)  ave(X) ave(X) ave(X) inst(X)
phys_out_filetimesteps=  1mth   1day    6hr     2hr     2hr   1TS
flag_t2m = 1, 1, 1, 0, 1
flag_u10m = 1, 1, 1, 0, 1  
flag_v10m = 1, 1, 1, 0, 1
flag_flat = 5, 1, 10, 0, 1 
flag_sens = 1, 1, 10, 0, 1
flag_topl = 1, 1, 10, 0, 1
flag_pluc = 1, 1, 1, 0, 1
flag_plul = 1, 1, 1, 0, 1
flag_precip = 1, 1, 1, 0, 1
flag_cape = 1, 10, 10, 0, 1
flag_cin = 1, 1, 1, 0, 1
flag_fbase = 1, 10, 10, 0, 1
flag_ale_wk = 1, 1, 1, 0, 1
flag_ale_bl = 1, 1, 1, 0, 1
flag_alp_wk = 1, 1, 1, 0, 1
flag_ale = 1, 1, 1, 0, 1
flag_alp = 1, 1, 1, 0, 1
flag_alp_bl = 1, 1, 1, 0, 1
flag_dtphy = 2, 10, 10, 0, 1
flag_dqphy = 2, 10, 10, 0, 1

flag_temp = 2, 1, 10, 0, 1
flag_ovap = 2, 10, 10, 0, 1
flag_rhum = 2, 1, 10, 0, 1
flag_dtthe = 2, 10, 10, 0, 1
flag_dtcon = 2, 10, 10, 0, 1
flag_dtwak = 2, 10, 10, 0, 1
flag_dqthe = 2, 10, 10, 0, 1
flag_dqcon = 2, 10, 10, 0, 1
flag_dqwak = 2, 10, 10, 0, 1

flag_cldm = 1, 1, 10, 0, 10
flag_cldl = 1, 1, 10, 0, 10
flag_cldh = 1, 1, 10, 0, 10
flag_tops = 1, 1, 10, 0, 10
flag_tops0 = 1, 1, 10, 0, 10
flag_topl0 = 1, 1, 10, 0, 10
flag_sols = 1, 1, 10, 0, 10
flag_soll = 1, 1, 10, 0, 10
flag_tsol = 1, 1, 10, 0, 10

flag_s2 = 1, 1, 10, 0, 10

flag_n2 = 1, 1, 10, 0, 10
flag_zmax_th = 1, 1, 10, 0, 10
flag_zlcl_th = 1, 1, 10, 0, 10
flag_ale_bl_stat = 1, 1, 10, 0, 10
falg_e_bl_trig = 1, 1, 10, 0, 10
flag_pg_roba_notrig = 1, 1, 10, 0, 10

# Sorties pour Marine
flag_cldemi = 1, 1, 10, 10, 10
flag_iwcon = 1, 1, 10, 10, 10
flag_rneb  = 1, 1, 10, 10, 10
flag_ref_ice = 1, 1, 10, 10, 10

# Sandrine
flag_u850 = 1, 1, 10, 10, 10
flag_v850 = 1, 1, 10, 10, 10
flag_t850 = 1, 1, 10, 10, 10
flag_u300 = 1, 1, 10, 10, 10
flag_v300 = 1, 1, 10, 10, 10
flag_t300 = 1, 1, 10, 10, 10
flag_u500 = 1, 1, 10, 10, 10
flag_v500 = 1, 1, 10, 10, 10
flag_t500 = 1, 1, 10, 10, 10
flag_u200 = 1, 1, 10, 10, 10
flag_v200 = 1, 1, 10, 10, 10
flag_t200 = 1, 1, 10, 10, 10

#
#Couplage avec autres modules
#******************************
#option de couplage pour ocean (D:force)
type_ocean=force
version_ocean=nemo
#avec ou sans orchidee (D:n)
VEGET=KEY_VEGET
#type_run=AMIP, ENSP, clim (D:AMIP)
type_run=CLIM
##  Soil Model  ou non               
soil_model=y

#Code de transfert radiatif
#**************************
# ajout de ce flag pour activer/desactiver le rayonnement (MPL)
# 0 : pas de rayonnement. 1: on active le rayonnement (D=1)
iflag_radia=1
## Nombre appels des routines de rayonnements ( par jour)                 
nbapp_rad=24

# Parametres "orbitaux/ ere geologique"
#######################################
#Parametres orbitaux
#excentricite
R_ecc=0
#equinoxe
R_peri=102.7
#inclinaison
R_incl=23.441
#constante solaire
solaire=1366.0896
##  Facteur additif pour l albedo
pmagic=0.
#
# Taux gaz a effet de serre
#######################################
#concentration des gaz (co2 ch4   n2o  cfc11 cfc12)
#                  (D: 348. 1650. 306.  280.  484.)
co2_ppm=0.36886500E+03
#RCO2=co2_ppm * 1.0e-06  * 44.011/28.97= 5.286789092164308E-04
CH4_ppb=0.17510225E+04
#RCH4=1.65E-06* 16.043/28.97= 9.137366240938903E-07
N2O_ppb=0.31585000E+03
#RN2O=306.E-09* 44.013/28.97= 4.648939592682085E-07
CFC11_ppt=5.18015181E+01
#RCFC11=280.E-12* 137.3686/28.97= 1.327690990680013E-09
CFC12_ppt=0.99862742E+03
#RCFC12=484.E-12* 120.9140/28.97= 2.020102726958923E-09
#
# Parametres effets directs/indirects des "aerosols"
#######################################
ok_cdnc=n
### 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=0
### 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=READ_CLIMOZ
#
# 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
#

# Specifique aquaplanete
########################

ok_all_xml=y
#qsol0=100.
flag_aerosol=0
ok_ade=.false.
ok_aie=.false.
iflag_ice_thermo=0