Ignore:
Timestamp:
Jan 20, 2010, 4:31:14 PM (14 years ago)
Author:
Laurent Fairhead
Message:

Work on the config files (*.def):

  • the parameters used by default are those of a reference simulation adjusted for radiative equilibrium for the 96x94x39 resolution (CLOUD09)
  • information on the coupling used (ocean and land) is transferred from the physiq.def file to the run.def file to ease the transfer of parameters files between forced and coupled simulations

Travail sur les fichiers de configuration (les*.def):

  • on prend les paramètres d'une simulation de référence ajusté pour un équilibre radiatif (CLOUD09) par défaut
  • les informations sur la configuration du couplage (océan et végétation) sont déplacées de physiq.def à run.def (pour un échange plus aisé des fichiers de paramétrisations entre les configurations forcées et couplées)
File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4V5.0-dev/physiq.def

    r1279 r1300  
    1 #
    21## $Id$
    32#
    4 #
    5 # Automatically generated make config: don t edit
    6 #
    7 type_ocean=force
    8 # avec ou sans orchidee
    9 VEGET=n
    10 #type_run = AMIP, ENSP, clim
    11 type_run=AMIP
    12 #
    13 # Controle des sorties
    14 # sorties moyennees tous les jours  dans histday.nc
    15 OK_journe=y
    16 # sorties moyennees tous les mois  dans histmth.nc
    17 OK_mensuel=y
    18 # sorties moyennees toutes les 6 ou bien 3h dans histhf.nc
    19 ok_hf=n
    20 # sorties moyennees tous les pas de temps de la physique dans histins.nc
    21 OK_instan=n
    22 #
    23 ecrit_mth=30.
    24 ecrit_day=1.
    25 ecrit_hf=0.25
    26 #
    27 #niveau de sortie "hf" lev_histhf
    28 lev_histhf=4
    29 #niveau de sortie "day" lev_histday
    30 lev_histday=5
    31 #niveau de sortie "mth" lev_histmth
    32 lev_histmth=4
    33 
    34 # parametres KE
    35 if_ebil=0
    36 epmax = .99
    37 ok_adj_ema = n
    38 iflag_clw = 1
    39 #
    40 # parametres nuages
    41 cld_lc_lsc = 2.6e-4
    42 cld_lc_con = 2.6e-4
    43 cld_tau_lsc = 3600.
    44 cld_tau_con = 3600.
    45 ffallv_lsc = 1.
    46 ffallv_con = 1.
    47 coef_eva = 2.e-5
    48 reevap_ice = y
    49 iflag_cldcon = 3
    50 iflag_pdf = 1
    51 fact_cldcon = 1.
    52 facttemps = 1.e-4
    53 ok_newmicro = y
    54 iflag_ratqs=0
    55 ratqsbas = 0.005
    56 ratqshaut = 0.33
    57 rad_froid = 35
    58 rad_chau1=12
    59 rad_chau2=11
    60 ksta_ter=1.e-7
    61 ksta=1.e-10
    62 #ok_kzmin : calcul Kzmin dans la CL de surface
    63 ok_kzmin=y
    64 #
    65 # parametres climatique
    66 R_ecc = 0.016715
    67 R_peri = 102.7
    68 R_incl = 23.441
    69 solaire = 1365.
    70 co2_ppm = 348.
    71 #RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97
    72 #RCO2 = 348. * 1.0e-06  * 44.011/28.97
    73 #RCO2 =   5.286789092164308E-04
    74 #RCO2 = 425.43e-06
    75 CH4_ppb = 1650.
    76 #RCH4 = 1.65E-06* 16.043/28.97
    77 #RCH4 =   9.137366240938903E-07
    78 N2O_ppb = 306.
    79 #RN2O = 306.E-09* 44.013/28.97
    80 #RN2O =    4.648939592682085E-07
    81 CFC11_ppt = 280.
    82 #RCFC11 = 280.E-12* 137.3686/28.97
    83 #RCFC11 =    1.327690990680013E-09
    84 CFC12_ppt = 484.
    85 #RCFC12 = 484.E-12* 120.9140/28.97
    86 #RCFC12 =    2.020102726958923E-09
    87 #
    88 # effets direct et indirect des aerosols
    89 ok_ade=n
    90 ok_aie=n
    91 bl95_b0=1.7
    92 bl95_b1=0.2
    93 #
    94 # parametres simulateur ISCCP
    95 #ok_isccp : y/n avec/sans simulateur ISCCP
    96 ok_isccp=n
    97 #top_height = 1 ou 3
    98 top_height = 1
    99 #overlap = 1, 2 ou 3
    100 overlap = 3
    101 #cdmmax
    102 cdmmax = 2.5E-3
    103 #cdhmax
    104 cdhmax = 2.0E-3
    105 #
    106 #ok_regdyn : y/n calcul/non des regymes dynamiques sur regions pre-definies
    107 ok_regdyn=y
    108 #
    109 # Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE nouvelle physique, 30 KE IPCC)
    110 iflag_con=30
    111 #
    112 # activation thermiques wake, ...
    113 iflag_thermals = 0
    114 nsplit_thermals =1
    115 tau_thermals=1800.
    116 iflag_pbl = 1
    117 iflag_coupl=0
    118 iflag_wake=0
    119 iflag_clos=0
    120 iflag_mix=1
    121 qqa1=0.
    122 qqa2=1.
    123 ## frequence (en  jours ) de l'ecriture du fichier histphy               
    124 ecritphy=30
     3# PARAMETRES ANCIENNEMENT DANS gcm.def
    1254##  Cycle diurne  ou non                 
    1265cycle_diurne=y
    … …  
    1298##  Choix ou non  de  New oliq               
    1309new_oliq=y
     10##  Activation ou non de la parametrisation de Hines pour la strato
     11ok_hines=n
    13112##  Orodr  ou  non   pour l orographie             
    13213ok_orodr=y
    … …  
    13617ok_limitvrai=n
    13718## Nombre  d'appels des routines de rayonnements ( par jour)                 
    138 nbapp_rad=12
     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
    13922## Facteur multiplication des precip convectives dans KE
    14023cvl_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 TracChangeset for help on using the changeset viewer.