Ignore:
Timestamp:
Jun 17, 2022, 4:24:49 PM (2 years ago)
Author:
lguez
Message:

Sync latest trunk changes to branch LMDZ-ECRAD.

Location:
LMDZ6/branches/LMDZ-ECRAD
Files:
11 edited
19 copied

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/LMDZ-ECRAD

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/context_lmdz.xml

    r3721 r4171  
    118118        <axis axis_ref="plev" />
    119119     </grid>
     120
     121     <grid id="grid_out_plev17">
     122        <domain id="dom_out" />
     123          <axis long_name="pressure" n_glo="17" name="plev_std" positive="down" standard_name="air_pressure" unit="Pa" value="(0,16)[ 100000. 92500. 85000. 70000. 60000. 50000. 40000. 30000. 25000. 20000. 15000. 10000. 7000. 5000. 3000. 2000. 1000. ]"><interpolate_axis coordinate="pres" order="1" type="polynomial" />
     124        </axis>
     125     </grid>
     126
    120127     <grid id="grid_out_plev">
    121128        <domain domain_ref="dom_out" />

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/context_lmdz_sans_cosp.xml

    r3721 r4171  
    107107     </grid>
    108108
     109     <grid id="grid_out_plev17">
     110        <domain id="dom_out" />
     111          <axis axis_ref="pres" long_name="pressure" n_glo="17" name="plev_std" positive="down" standard_name="air_pressure" unit="Pa" value="(0,16)[ 100000. 92500. 85000. 70000. 60000. 50000. 40000. 30000. 25000. 20000. 15000. 10000. 7000. 5000. 3000. 2000. 1000. ]"><interpolate_axis coordinate="pres" order="1" type="polynomial" />
     112        </axis>
     113     </grid>
     114
    109115     <grid id="grid_out_plev">
    110116        <domain domain_ref="dom_out" />

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/field_def_lmdz.xml

    r3791 r4171  
    5656        <field id="t2m_oce"    long_name="Temp 2m oce"    unit="K" />
    5757        <field id="t2m_sic"    long_name="Temp 2m sic"    unit="K" />
     58        <field id="nt2mout"    long_name="nt2m hors intervalle, calc.complet" unit="-" />
     59        <field id="nq2mout"    long_name="nq2m hors intervalle, calc. complet" unit="-" />
     60        <field id="nu2mout"    long_name="nu2m hors intervalle, calc. complet" unit="-" />
     61            <field id="nt2moutfg"  long_name="nt2m hors intervalle, calc. complet/fg" unit="-" />
     62            <field id="nq2moutfg"  long_name="nq2m hors intervalle, calc. complet/fg" unit="-" />
     63            <field id="nu2moutfg"  long_name="nu2m hors intervalle, calc. complet/fg" unit="-" />
    5864        <field id="t2m_probsup25"   field_ref="t2m" long_name="Prob. t2m exceeds 25 degC" unit="-"> t2m &gt; 298.15 </field>
    5965        <field id="t2m_probsup28"   field_ref="t2m" long_name="Prob. t2m exceeds 28 degC" unit="-"> t2m &gt; 301.15 </field>
    … …  
    204210        <field id="tsol_oce"    long_name="Temperature oce"    unit="K" />
    205211        <field id="tsol_sic"    long_name="Temperature sic"    unit="K" />
    206         <field id="evappot_ter"    long_name="Temperature ter"    unit="K" />
    207         <field id="evappot_lic"    long_name="Temperature lic"    unit="K" />
    208         <field id="evappot_oce"    long_name="Temperature oce"    unit="K" />
    209         <field id="evappot_sic"    long_name="Temperature sic"    unit="K" />
     212        <field id="evappot_ter"    long_name="Potential evaporation ter"    unit="kg/(m2*s)" />
     213        <field id="evappot_lic"    long_name="Potential evaporation lic"    unit="kg/(m2*s)" />
     214        <field id="evappot_oce"    long_name="Potential evaporation oce"    unit="kg/(m2*s)" />
     215        <field id="evappot_sic"    long_name="Potential evaporation sic"    unit="kg/(m2*s)" />
    210216        <field id="sens_ter"    long_name="Sensible heat flux ter"    unit="W/m2" />
    211217        <field id="sens_lic"    long_name="Sensible heat flux lic"    unit="W/m2" />
    … …  
    520526        <field id="flx_co2_land"   long_name="CO2 flux from land"             unit="kg CO2/m2/s" />    <!-- Added OB -->
    521527        <field id="flx_co2_ocean"  long_name="CO2 flux from ocean"            unit="kg CO2/m2/s" />    <!-- Added OB -->
     528        <field id="flx_co2_ocean_cor"  long_name="correction of the CO2 flux from ocean"          unit="kg CO2/m2/s" />    <!-- Added PC -->
     529        <field id="flx_co2_land_cor"   long_name="correction of the CO2 flux from land"           unit="kg CO2/m2/s" />    <!-- Added PC -->
    522530        <field id="flx_co2_ff"     long_name="CO2 flux from ff"               unit="kg CO2/m2/s" />    <!-- Added OB -->
    523531        <field id="flx_co2_bb"     long_name="CO2 flux from bb"               unit="kg CO2/m2/s" />    <!-- Added OB -->
     532        <field id="delta_SST" long_name="ocean-air interface temperature minus bulk SST" unit="K" detect_missing_value=".true." />
     533        <field id="delta_sal"
     534               long_name="ocean-air interface salinity minus bulk salinity"
     535               unit="ppt" detect_missing_value=".true." />
     536        <field id="dS_ns" long_name="delta salinity near surface" unit="ppt" detect_missing_value=".true." />
     537        <field id="dT_ns" long_name="sub-skin temperature minus foundation temperature" unit="K" detect_missing_value=".true." />
     538        <field id="dTer"               long_name="interface temperature minus sub-skin temperature"            unit="K" detect_missing_value=".true." />
     539        <field id="dSer"
     540               long_name="salinity variation in the diffusive microlayer"
     541               unit="ppt" detect_missing_value=".true." />
     542        <field id="tkt" long_name="thickness of thermal microlayer"            unit="m"                detect_missing_value=".true." />
     543        <field id="tks" long_name="thickness of salinity microlayer"           unit="m"                detect_missing_value=".true." />
     544        <field id="taur" long_name="momentum flux due to rain"         unit="Pa"               detect_missing_value=".true." />
     545        <field id="SSS" long_name="bulk sea-surface salinity" unit="ppt" detect_missing_value=".true." />
    524546    </field_group>
    525547
    … …  
    736758        <field id="cool_volc"    long_name="LW cooling rate from volcano"    unit="K/s" />
    737759        <field id="pvap"    long_name="pvap intermediary variable" unit="-">pres*ovap*461.5 / (287.04*(1.+ (10.9491/18.0153)*ovap)) </field>
     760        <field id="oclr"    long_name="Clear sky total water"    unit="kg/kg" />
     761        <field id="ocld"    long_name="Cloudy sky total water"   unit="kg/kg" />
     762        <field id="oss"     long_name="ISSR total water"    unit="kg/kg" />
     763        <field id="ovc"     long_name="In-cloud vapor"    unit="kg/kg" />
     764        <field id="rnebclr"    long_name="Clear sky fraction"    unit="-" />
     765        <field id="rnebss"     long_name="ISSR fraction"     unit="-" />
     766        <field id="rnebseri"   long_name="Cloud fraction"    unit="-" />
     767        <field id="gammass"    long_name="Supersaturation ratio"    unit="-" />
     768        <field id="N1ss"       long_name="N1ss"    unit="-" />
     769        <field id="N2ss"       long_name="N2ss"    unit="-" />
     770        <field id="drnebsub"   long_name="Cloud fraction change because of sublimation"    unit="s-1" />
     771        <field id="drnebcon"   long_name="Cloud fraction change because of condensation"   unit="s-1" />
     772        <field id="drnebtur"   long_name="Cloud fraction change because of turbulence"     unit="s-1" />
     773        <field id="drnebavi"   long_name="Cloud fraction change because of aviation"       unit="s-1" />
     774        <field id="qsatl"      long_name="Saturation with respect to liquid water"    unit="-" />
     775        <field id="qsats"      long_name="Saturation with respect to liquid ice"      unit="-" />
     776        <field id="flightm"    long_name="Aviation meters flown"      unit="m/s" />
     777        <field id="flighth2o"  long_name="Aviation H2O emitted"       unit="kg H2O/s" />
     778        <field id="fcontrP"    long_name="Gridbox fraction with potential persistent contrails"     unit="-" />
     779        <field id="fcontrN"    long_name="Gridbox fraction with potential non-persistent contrails"     unit="-" />
     780        <field id="qcontr"     long_name="Contrail qcontr"     unit="Pa" />
     781        <field id="qcontr2"     long_name="Contrail qcontr2"     unit="kg/kg" />
     782        <field id="Tcontr"     long_name="Contrail Tcontr"     unit="K" />
    738783    </field_group>
    739784

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/file_def_histLES_lmdz.xml

    r3507 r4171  
    22<file_definition>
    33    <file_group id="defile">
    4         <file id="histLES" name="XhistLES" output_freq="1h" output_level="11" enabled="false" compression_level="0">
     4        <file id="histLES" name="XhistLES" output_freq="1h" output_level="11" enabled=".FALSE." compression_level="0">
    55           
    66            <!-- VARS 1D -->

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/file_def_histday_lmdz.xml

    r3516 r4171  
    638638                <field field_ref="l_mixmin_sic" level="10" />
    639639                <field field_ref="ozone_daylight" level="10" />
     640                <field field_ref="oclr" level="10" />
     641                <field field_ref="ocld" level="10" />
     642                <field field_ref="oss" level="10" />
     643                <field field_ref="ovc" level="10" />
     644                <field field_ref="rnebclr" level="10" />
     645                <field field_ref="rnebss" level="10" />
     646                <field field_ref="rnebseri" level="10" />
     647                <field field_ref="gammass" level="10" />
     648                <field field_ref="N1ss" level="10" />
     649                <field field_ref="N2ss" level="10" />
     650                <field field_ref="drnebsub" level="10" />
     651                <field field_ref="drnebcon" level="10" />
     652                <field field_ref="drnebtur" level="10" />
     653                <field field_ref="drnebavi" level="10" />
     654                <field field_ref="qsatl" level="10" />
     655                <field field_ref="qsats" level="10" />
     656                <field field_ref="flightm" level="10" />
     657                <field field_ref="flighth2o" level="10" />
     658                <field field_ref="fcontrP" level="10" />
     659                <field field_ref="fcontrN" level="10" />
     660                <field field_ref="qcontr" level="10" />
     661                <field field_ref="qcontr2" level="10" />
     662                <field field_ref="Tcontr" level="10" />
    640663            </field_group>
    641664           

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/file_def_histhf_lmdz.xml

    r3622 r4171  
    485485                <field field_ref="treedrg_oce" level="10" />
    486486                <field field_ref="treedrg_sic" level="10" />
    487 
    488487                <field field_ref="cldtau" level="10" />
    489488                <field field_ref="cldemi" level="10" />
    … …  
    500499                <field field_ref="lwcon" level="10" />
    501500                <field field_ref="iwcon" level="10" />
    502                 <field field_ref="temp" level="10" />
     501                <field field_ref="temp" level="4" />
    503502                <field field_ref="theta" level="10" />
    504503                <field field_ref="ovap" level="10" />
    … …  
    654653                <field field_ref="rsdcs4co2" level="10" />
    655654                <field field_ref="rldcs4co2" level="10" />
     655                <field field_ref="oclr" level="10" />
     656                <field field_ref="ocld" level="10" />
     657                <field field_ref="oss" level="10" />
     658                <field field_ref="ovc" level="10" />
     659                <field field_ref="rnebclr" level="10" />
     660                <field field_ref="rnebss" level="10" operation="instant" />
     661                <field field_ref="rnebseri" level="10" />
     662                <field field_ref="gammass" level="10" />
     663                <field field_ref="N1ss" level="10" />
     664                <field field_ref="N2ss" level="10" />
     665                <field field_ref="drnebsub" level="10" />
     666                <field field_ref="drnebcon" level="10" />
     667                <field field_ref="drnebtur" level="10" />
     668                <field field_ref="drnebavi" level="10" />
     669                <field field_ref="qsatl" level="10" />
     670                <field field_ref="qsats" level="10" />
     671                <field field_ref="fcontrP" level="10" />
     672                <field field_ref="fcontrN" level="10" />
     673                <field field_ref="qcontr" level="10" operation="instant" />
     674                <field field_ref="qcontr2" level="10" operation="instant" />
     675                <field field_ref="Tcontr" level="10" operation="instant" />
    656676              </field_group>
    657677

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/file_def_histins_lmdz.xml

    r3507 r4171  
    22<file_definition>
    33    <file_group id="defile">
    4         <file id="histins" name="Xhistins" output_freq="1ts" output_level="11" enabled="false" compression_level="0">
     4        <file id="histins" name="Xhistins" output_freq="1ts" output_level="11" enabled=".FALSE." compression_level="0">
    55           
    66            <!-- VARS 1D -->
    … …  
    374374                <field field_ref="lcc" level="10" />
    375375                <field field_ref="wvapp" level="10" />
    376                 <field field_ref="ozone_daylight" level="10" />
    377376                <field field_ref="albe_ter" level="10" />
    378377                <field field_ref="albe_lic" level="10" />
    … …  
    626625                <field field_ref="rsdcs4co2" level="10" />
    627626                <field field_ref="rldcs4co2" level="10" />
     627                <field field_ref="ozone_daylight" level="10" />
    628628            </field_group>
    629629          </field_group>

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/file_def_histmth_lmdz.xml

    r3507 r4171  
    280280                <field field_ref="colO3_strat" level="2" />     <!-- Added ThL -->
    281281                <field field_ref="colO3_trop"  level="2" />     <!-- Added ThL -->
     282                <field field_ref="flx_co2_ocean"  level="1" />
     283                <field field_ref="flx_co2_land"   level="1" />
     284                <field field_ref="flx_co2_ocean_cor" level="1" /> <!-- Added PC -->
     285                <field field_ref="flx_co2_land_cor"  level="1" /> <!-- Added PC -->
     286                <field field_ref="flx_co2_ff"  level="1" />
     287                <field field_ref="flx_co2_bb"  level="1" />
    282288
    283289                <field_group operation="average" detect_missing_value=".true.">
    … …  
    689695            </field_group>
    690696
     697<!--
     698                <field field_ref="psbg" level="5" />
     699                <field field_ref="tro3" level="5" />
     700                <field field_ref="tro3_daylight" level="5" />
     701-->
     702
     703            <!-- VARS 3D -->
     704            <field_group operation="average" grid_ref="grid_out_plev17">
     705                <field field_ref="temp" name="ta" grid_ref="grid_out_plev17" level="4" />
     706                <field field_ref="geop" name="zg" grid_ref="grid_out_plev17" level="4" />
     707                <field field_ref="ovap" name="hus" grid_ref="grid_out_plev17" level="4" />
     708                <field field_ref="rhum" name="hur" grid_ref="grid_out_plev17" level="4" />
     709                <field field_ref="vitu" name="ua" grid_ref="grid_out_plev17" level="4" />
     710                <field field_ref="vitv" name="va" grid_ref="grid_out_plev17" level="4" />
     711                <field field_ref="vitw" name="wap" grid_ref="grid_out_plev17" level="4" />
     712<!--
     713                <field field_ref="tro3" name="tro3" grid_ref="grid_out_plev17" level="4" />
     714                <field field_ref="tro3_daylight" name="tro3_daylight" grid_ref="grid_out_plev17" level="4" />
     715                <field field_ref="tnondef" name="psbg" grid_ref="grid_out_plev17" level="4" />
     716                <field field_ref="vitu_vitv" name="uv" grid_ref="grid_out_plev17" level="4" />
     717                <field field_ref="vitv_ovap" name="vq" grid_ref="grid_out_plev17" level="4" />
     718                <field field_ref="vitv_temp" name="vT" grid_ref="grid_out_plev17" level="4" />
     719                <field field_ref="vitw_ovap" name="wq" grid_ref="grid_out_plev17" level="4" />
     720                <field field_ref="vitv_geop" name="vphi" grid_ref="grid_out_plev17" level="4" />
     721                <field field_ref="vitw_temp" name="wT" grid_ref="grid_out_plev17" level="4" />
     722                <field field_ref="vtiu_vitu" name="u2" grid_ref="grid_out_plev17" level="4" />
     723                <field field_ref="vitv_vitv" name="v2" grid_ref="grid_out_plev17" level="4" />
     724                <field field_ref="temp_temp" name="T2" grid_ref="grid_out_plev17" level="4" />
     725-->
     726            </field_group>
     727
    691728            <!-- VARS 3D -->
    692729            <field_group operation="average" grid_ref="grid_out_spectband">

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/namelist_ecrad

    r3904 r4171  
    4646do_lw_side_emissivity   = true,
    4747do_3d_lw_multilayer_effects = true,
    48 directory_name          = "/d3/idelkadi/Ecrad/ecrad/data",         ! Location of configuration files
     48directory_name          = "data",         ! Location of configuration files
    4949! directory_name                = "../../data",         ! Location of configuration files
    5050

  • LMDZ6/branches/LMDZ-ECRAD/DefLists/physiq.def_NPv6.1

    r3621 r4171  
    11#
    2 # $Id: physiq.def_NPv6.0.15_S 3593 2018-02-20 10:07:02Z jgipsl $
     2# $Id: physiq.def_NPv6.1 4020 2018-09-05 17:28:38Z jgipsl $
    33#
    44
    … …  
    163163inertie_lic=2000
    164164inertie_sic=2000
    165     # tau_gl = constante de rappel de la temperature a la surface de la glace (en nombre de jours)
    166     # Utilise celon le forcage : tau_gl=2 pour Amip et tau_gl=1 pour Hadisst2
    167 tau_gl= 2.
     165
     166
    168167          # ???
    169168iflag_albedo=1
    170169
    171           # Frequence appel convection. Nombre appels par jour
     170          # Frequence appel rayonnement. Nombre appels par jour
    172171nbapp_cv=48
    173172
    … …  
    369368ptcrit=500
    370369
    371           # Precipitation ejection from adiabatic ascents
    372           # cvflag_prec_eject: when TRUE ejection is activated. Default: False
    373           # ejectliq : 0. =>liquid precipitation is carried up in the adiab updraughts
    374           #            1. =>liquid precipitation is ejected from adiab updraughts
    375           # ejectice   0. =>solid precipitation is carried up in the adiab updraughts
    376           #            1. =>solid precipitation is ejected from adiab updraughts
    377           # ejectliq may take only 0. and 1. values
    378           # ejectice may take any value between 0. and 1.
    379           #   Default : cvflag_prec_eject=n, ejectliq=0., ejectice=0.
    380 cvflag_prec_eject=n
    381 ejectliq=0.
    382 ejectice=0.
    383 
    384 
    385           # Check whether parcel level temperature and specific humidity
    386           # are reasonable and demand that LCL be high enough
    387 cv_flag_feed=2
    388 
    389           # save 20% CPU
    390 ok_optim_yield=y
    391 
    392           # Convective clouds
    393 iflag_cld_cv=0
    394 
    395           # To get back to bug where convective output variables
    396           # were set to 0 at the start of the physics
    397 ok_bug_cv_trac=n
    398 
    399            # solving of the temperature equation for the adiabatic ascent
    400            # 0 -> AJ scheme in two loops ; ice fraction given by envt. temp.
    401            # 1 -> JYG scheme in one loop ; ice fraction given by envt. temp.
    402            # 2 -> JYG scheme in one loop ; ice fraction given by adiab. asc. temp.
    403            # Default: 0
    404            # 2 is costly and ill tested; experts only.
    405 icvflag_Tpa=0
    406 
    407            # taking into account the dependence of qv on qt at saturation
    408            # Default: n
    409 qsat_depends_on_qt=n
    410 
    411 ##########################################################
    412 # Fermeture convective et couplages
    413 ##########################################################
    414 
    415           # Appel separe des thermiques
    416 iflag_pbl_split=10
    417 
    418           #rendement sur la puissance fournie par les thermiques a la convection
    419 alp_bl_k=0.5
    420 
    421           #Constant value to be added to Alp (D: 0.)
    422 alp_offset=0
    423 
    424           #Convection must reach a level at which T<t_top_max to be active
    425 t_top_max=1000
    426 
    427           ### nrlmd le 10/04/2012: Declenchement ###
    428           #Declenchement par couche-limite (0: deterministe, 1: stochastique)
    429 iflag_trig_bl=1
    430 
    431           #Surface critique au LCL
    432 s_trig=12000000
    433 
    434           #Temps de decorrelation entre scenes de cumulus
    435 tau_trig_shallow=1200
    436 
    437           #Temps de decorrelation entre scenes de cumulonimbus
    438 tau_trig_deep=1200
    439 
    440           # Pour la stabilisation
    441 random_notrig_max=0.99
    442 
    443           # Temperature a la quelle on fond les descentes glacees
    444 tmax_fonte_cv=278.15
    445 
    446           #Fermeture par couche-limite (Nicolas Rochetin)
    447           #(0: deterministe, 1: statistique, 2: statistique+convergence
    448           # grande echelle)
    449 iflag_clos_bl=1
    450 
    451           #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
    452           # def: 0.33
    453 stark=0.33
    454 
    455           #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided
    456           # by wakes, def: 0.25
    457 alpk=0.25
    458 
    459 #
    460 ##########################################################
    461 # Couche limite convective / modele du thermique
    462 ##########################################################
    463 
    464           # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique)
    465 iflag_thermals=18
    466 
    467           # pas de temps du spliting pour les thermiques
    468           # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
    469 nsplit_thermals=1
    470 
    471           # tau_thermals pour avoir une constante de temps sur les thermiques.
    472           # Non valide
    473 tau_thermals=0
    474 
    475           # Flag controlant le detrainement et l entrainement
    476 iflag_thermals_ed=8
    477 
    478           # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer
    479           # le detrainement (A. Jam)
    480 fact_thermals_ed_dz=0.07
    481 
    482           # Une option de plus sur les thermiques.
    483 iflag_thermals_optflux=0
    484 
    485           #couplage ou non avec la convection
    486 iflag_coupl=5
    487 
    488           # Seuil pour desactiver les thermiques en cas de forte inversion
    489           # en sommet de couche limite. Actif si iflag_thermals impair.
    490 seuil_inversion=-0.08
    491 
    492           # Reactivation du calcul d un zmax continu pour les thermiques
    493           # reactivation d une fermeture humide des thermiques
    494           # iflag_thermals_closure=2 (D=1)
    495 iflag_thermals_closure=2
    496 
    497 
    498 ##########################################################
    499 # Flags wakes
    500 ##########################################################
    501 
    502           #Wake scheme switch (D:0, 0:AR4, 1:New Physics)
    503 iflag_wake=1
    504 
    505           #multiplicative factor of the damping by gravity waves, def: 4.
    506 coefgw=4
    507 
    508           #wake density = number of wake centers per m2, def: 8.E-12
    509           #wdens_ref=8.E-12
    510 wdens_ref_o=1e-09
    511 wdens_ref_l=8e-12
    512 
    513           # Ajustement convectif prealable au calcul des poches
    514 ok_adjwk=y
    515 
    516370          # Prevent some crashes
    517371          # Filter out bad wakes
    … …  
    520374iflag_alp_wk_cond=1
    521375
     376
     377          # Check whether parcel level temperature and specific humidity
     378          # are reasonable and demand that LCL be high enough
     379cv_flag_feed=2
     380
     381          # save 20% CPU
     382ok_optim_yield=y
     383
     384          # Convective clouds
     385iflag_cld_cv=0
     386
     387          # To get back to bug where convective output variables
     388          # were set to 0 at the start of the physics
     389#libIGCM version : ok_bug_cv_trac=_AUTO_: DEFAULT = n
     390ok_bug_cv_trac=n
     391
     392##########################################################
     393# Fermeture convective et couplages
     394##########################################################
     395
     396          # Appel separe des thermiques
     397iflag_pbl_split=10
     398
     399          #rendement sur la puissance fournie par les thermiques a la convection
     400alp_bl_k=0.5
     401
     402          #Constant value to be added to Alp (D: 0.)
     403alp_offset=0
     404
     405          #Convection must reach a level at which T<t_top_max to be active
     406t_top_max=1000
     407
     408          ### nrlmd le 10/04/2012: Declenchement ###
     409          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
     410iflag_trig_bl=1
     411
     412          #Surface critique au LCL
     413s_trig=12000000
     414
     415          #Temps de decorrelation entre scenes de cumulus
     416tau_trig_shallow=1200
     417
     418          #Temps de decorrelation entre scenes de cumulonimbus
     419tau_trig_deep=1200
     420
     421          # Pour la stabilisation
     422random_notrig_max=0.99
     423
     424          # Temperature a la quelle on fond les descentes glacees
     425tmax_fonte_cv=278.15
     426
     427          #Fermeture par couche-limite (Nicolas Rochetin)
     428          #(0: deterministe, 1: statistique, 2: statistique+convergence
     429          # grande echelle)
     430iflag_clos_bl=1
     431
     432          #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
     433          # def: 0.33
     434stark=0.33
     435
     436          #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided
     437          # by wakes, def: 0.25
     438alpk=0.25
     439
     440#
     441##########################################################
     442# Couche limite convective / modele du thermique
     443##########################################################
     444
     445          # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique)
     446iflag_thermals=18
     447
     448          # pas de temps du spliting pour les thermiques
     449          # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
     450nsplit_thermals=1
     451
     452          # tau_thermals pour avoir une constante de temps sur les thermiques.
     453          # Non valide
     454tau_thermals=0
     455
     456          # Flag controlant le detrainement et l entrainement
     457iflag_thermals_ed=8
     458
     459          # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer
     460          # le detrainement (A. Jam)
     461fact_thermals_ed_dz=0.07
     462
     463          # Une option de plus sur les thermiques.
     464iflag_thermals_optflux=0
     465
     466          #couplage ou non avec la convection
     467iflag_coupl=5
     468
     469          # Seuil pour desactiver les thermiques en cas de forte inversion
     470          # en sommet de couche limite. Actif si iflag_thermals impair.
     471seuil_inversion=-0.08
     472
     473          # Reactivation du calcul d un zmax continu pour les thermiques
     474          # reactivation d une fermeture humide des thermiques
     475          # iflag_thermals_closure=2 (D=1)
     476iflag_thermals_closure=2
     477
     478
     479##########################################################
     480# Flags wakes
     481##########################################################
     482
     483          #Wake scheme switch (D:0, 0:AR4, 1:New Physics)
     484iflag_wake=1
     485
     486          #multiplicative factor of the damping by gravity waves, def: 4.
     487coefgw=4
     488
     489          #wake density = number of wake centers per m2, def: 8.E-12
     490          #wdens_ref=8.E-12
     491wdens_ref_o=1e-09
     492wdens_ref_l=8e-12
     493
     494          # Ajustement convectif prealable au calcul des poches
     495ok_adjwk=y
    522496
    523497
    … …  
    557531
    558532cdnc_min=10.
    559 
    560 
Note: See TracChangeset for help on using the changeset viewer.