Changeset 411 in lmdz_wrf for branches/LMDZ_WRFmeas


Ignore:
Timestamp:
May 13, 2015, 4:59:30 PM (10 years ago)
Author:
lfita
Message:

Getting a more 'state-of-the-art' LMDZ configuration for high resolution

Location:
branches/LMDZ_WRFmeas/WRFV3/run
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • branches/LMDZ_WRFmeas/WRFV3/run/config.def

    r1 r411  
    1414ok_regdyn=y
    1515# A transformer en la ligne ci-dessous quand ce sera compatible avec libIGCM
    16 phys_out_filekeys=       y       y       n       y       n
    17 phys_out_filenames=      histmth histday histhf  histins histLES
    18 phys_out_filetimesteps=  5day   1day    1hr     3hr     6hr
    19 phys_out_filelevels=     10      10      0       10       4
    20 phys_out_filetypes=      ave(X)  ave(X)  ave(X)  inst(X) inst(X)
     16phys_out_filekeys=       y       y       n       y       n        n
     17phys_out_filenames=      histmth histday histhf  histins histLES  histstations
     18phys_out_filetimesteps=  5day    1day    1hr     1hr     6hr      15mn
     19phys_out_filelevels=     10      10      0       10      4        10
     20phys_out_filetypes=      ave(X)  ave(X)  ave(X)  inst(X) inst(X)  inst(x)
    2121
    2222#

  • branches/LMDZ_WRFmeas/WRFV3/run/gcm.def

    r1 r411  
    11## $Id: gcm.def 1403 2010-07-01 09:02:53Z fairhead $
    22## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )     
    3 day_step=144
     3day_step=8640
    44## periode pour le pas Matsuno (en pas)
    5 iperiod=6
     5iperiod=5
    66## periode de la dissipation (en pas)
    7 idissip=6
     7idissip=5
    88## choix de l'operateur de dissipation (star ou  non star )
    99lstardis=y
    … …  
    1515niterh=2
    1616## temps de dissipation des plus petites long.d ondes pour u,v (gradiv) 
    17 tetagdiv=18000.
     17tetagdiv=300.
    1818## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
    19 tetagrot=18000.
     19tetagrot=300.
    2020## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
    21 tetatemp=18000.
     21tetatemp=300.
    2222## period appel a la dissipation (multiple iperiod, 0 si calcul automatique)
    2323dissip_period=0
    … …  
    3030## 1: avec physique (e.g. physique phylmd)
    3131## 2: avec rappel newtonien dans la dynamique                                         
    32 iflag_phys=76
     32iflag_phys=1
    3333## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ?
    3434## (sans fichiers de demarrage, initialisation des champs par iniacademic
    3535##  dans la dynamique)
    36 read_start=n
     36read_start=y
    3737## periode de la physique (en pas dynamiques, n'a de sens que si iflag_phys=1)                                       
    38 iphysiq=1
     38iphysiq=20
    3939# Nombre de sous pas de temps pour la physique
    4040nsplit_phys=1
    4141##  Avec ou sans strato
    42 ok_strato=n
     42ok_strato=y
    4343#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
    44 iflag_top_bound=0
     44iflag_top_bound=2
    4545#  Coefficient pour la couche eponge (valeur derniere couche)
    4646tau_top_bound=5.e-5

  • branches/LMDZ_WRFmeas/WRFV3/run/orchidee.def

    r1 r411  
    11#
    2 #
    3 # Parameter file for LMDZ4OR_v2 configuration
    4 # See comments : http://forge.ipsl.jussieu.fr/orchidee/
    5 #
    6 STOMATE_OK_CO2=TRUE
    7 # STOMATE_OK_STOMATE is not set
    8 # STOMATE_OK_DGVM is not set
    9 # STOMATE_WATCHOUT is not set
    10 SECHIBA_restart_in=start_sech.nc
    11 SECHIBA_rest_out=restart_sech.nc
    12 SECHIBA_reset_time=y
    13 #
    14 OUTPUT_FILE=sechiba_out.nc
    15 WRITE_STEP=2592000
    16 SECHIBA_HISTLEVEL=5
    17 #
    18 SECHIBA_HISTFILE2 = FALSE
    19 SECHIBA_OUTPUT_FILE2 = sechiba_out_2.nc
    20 WRITE_STEP2 = 86400.0
    21 SECHIBA_HISTLEVEL2 = 1
    22 #
    23 STOMATE_OUTPUT_FILE=stomate_history.nc
    24 STOMATE_HIST_DT=10.
    25 STOMATE_HISTLEVEL=0
    26 SECHIBA_DAY=0.0
    27 SECHIBA_ZCANOP=0.5
    28 DT_SLOW=86400.
    29 # IMPOSE_VEG is not set
    30 VEGETATION_FILE=carteveg5km.nc
    31 # VEGETATION_FILE=pft_new.nc
    32 DIFFUCO_LEAFCI=233.
    33 CONDVEG_SNOWA=default
    34 # IMPOSE_AZE is not set
    35 SOILALB_FILE=soils_param.nc
    36 SOILTYPE_FILE=soils_param.nc
    37 ENERBIL_TSURF=280.
    38 HYDROL_SNOW=0.0
    39 HYDROL_SNOWAGE=0.0
    40 HYDROL_SNOWICE=0.0
    41 HYDROL_SNOWICEAGE=0.0
    42 HYDROL_HDRY=1.0
    43 HYDROL_HUMR=1.0
    44 HYDROL_BQSB=default
    45 HYDROL_GQSB=0.0
    46 HYDROL_DSG=0.0
    47 HYDROL_DSP=default
    48 HYDROL_QSV=0.0
    49 HYDROL_OK_HDIFF=n
    50 HYDROL_TAU_HDIFF=1800.
    51 THERMOSOIL_TPRO=280.
    52 RIVER_ROUTING=y
    53 ROUTING_FILE=routing.nc
    54 LAI_MAP=y
    55 LAI_FILE=lai2D.nc
    56 SECHIBA_QSINT=0.02
    57 ALB_BARE_MODEL = FALSE
    58 PERCENT_THROUGHFALL_PFT = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30.
    59 RVEG_PFT = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5
    60 CDRAG_FROM_GCM = .TRUE.
    61 #LAND_USE=y
    62 #VEGET_YEAR=0
    63 #VEGET_UPDATE=1Y
     2#**************************************************************************
     3#                    Parameter file list for ORCHIDEE
     4#**************************************************************************
     5
     6# Input and output
     7#**************************************************************************
     8# Restart the time from the GCM.
     9# default = n
     10SECHIBA_reset_time = y
     11
     12# Name of restart to read for initial conditions
     13# default = NONE
     14SECHIBA_restart_in = NONE
     15
     16# Name of restart to read for initial conditions of STOMATE
     17# default = NONE
     18STOMATE_RESTART_FILEIN = NONE
     19
     20# Use XIOS for writing diagnostics file
     21# defulat = n
     22XIOS_ORCHIDEE_OK = n
     23
     24# Flag to activate sechiba_out_2.nc history file for SECHIBA
     25# default  = FALSE
     26SECHIBA_HISTFILE2 = y
     27
     28# SECHIBA history output level (0..10)
     29# default = 5
     30SECHIBA_HISTLEVEL = 5
     31
     32# SECHIBA history 2 output level (0..10)
     33# default = 1
     34SECHIBA_HISTLEVEL2 = 5
     35
     36# STOMATE history output level (0..10)
     37# default = 10
     38STOMATE_HISTLEVEL = 5
     39
     40# Writefrequency in seconds in sechiba_history.nc
     41# default = 86400.0
     42WRITE_STEP = 86400.0
     43
     44# Writefrequency in seconds sechiba_out_2.nc
     45# default = 1800.0
     46WRITE_STEP2 = 3600.0
     47
     48# Writefrequency in days in stomate_history.nc
     49# default = 10.
     50STOMATE_HIST_DT = 1
     51
     52# Writefrequency in days or -1 for monthly output in stomate_ipcc_history.nc
     53# default = 0.
     54STOMATE_IPCC_HIST_DT = -1
     55
     56
     57# Hydrology parameters
     58#**************************************************************************
     59# Activate the multi-layer diffusion scheme adapted from CWRR
     60# by default the Choisnel hydrology is used.
     61# default = n
     62HYDROL_CWRR = n
     63
     64# Total depth of soil reservoir
     65# default = 2./4. depending on HYDROL_CWRR
     66HYDROL_SOIL_DEPTH = 4.
     67
     68# Root profile
     69HYDROL_HUMCSTE = 5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1.
     70
     71# Activate river routing
     72# default = n
     73RIVER_ROUTING = n
     74
     75
     76#**************************************************************************
     77# Activate Stomate component
     78# default = n
     79STOMATE_OK_STOMATE = y
     80
     81# Activate calculations of CO2 according to Farqhuar and Ball
     82# default = n
     83STOMATE_OK_CO2 = y
     84
     85# Use cdrag coefficient from gcm
     86CDRAG_FROM_GCM = y
     87
     88# Interception reservoir coefficient.
     89# default = 0.1
     90SECHIBA_QSINT = 0.02
     91
     92# Parmeters related to vegetation map
     93#**************************************************************************
     94# Read a land_use vegetation map
     95# default = n
     96LAND_USE = y
     97
     98# Year of the land_use vegetation map readed
     99# year off the pft map
     100# default = 282
     101VEGET_YEAR = 1
     102
     103# The parameter is used to bypass veget_year count
     104# and reinitialize it with VEGET_YEAR parameter.
     105# Then it is possible to change LAND USE file.
     106# If LAND_USE
     107# default = y
     108VEGET_REINIT = n
     109
     110# Update vegetation frequency (since 2.0 version)
     111# The veget datas will be update at this period if LAND_USE
     112# default = 1Y
     113VEGET_UPDATE = 0Y
     114
     115# treat land use modifications
     116# With this variable, you can use a Land Use map
     117# to simulate anthropic modifications such as   
     118# deforestation.                               
     119# If LAND_USE
     120# default = y
     121LAND_COVER_CHANGE = n
     122

  • branches/LMDZ_WRFmeas/WRFV3/run/physiq.def

    r1 r411  
    1111f_cdrag_stable=1.
    1212f_cdrag_ter=1.
    13 f_cdrag_oce=0.8
     13f_cdrag_oce=0.7
    1414#cdmmax
    1515cdmmax=2.5E-3
    … …  
    3232#******************
    3333#seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
    34 cld_lc_lsc=4.16e-4
    35 cld_lc_con=4.16e-4
     34cld_lc_lsc=6e-4
     35cld_lc_con=6e-4
    3636#constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.)
    3737cld_tau_lsc=1800.
    3838cld_tau_con=1800.
    3939#facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
    40 ffallv_lsc=0.5
    41 ffallv_con=0.5
     40ffallv_lsc=1.35
     41ffallv_con=1.35
    4242#coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
    43 coef_eva=2e-5
     43coef_eva=1e-4
    4444reevap_ice=y
    4545#calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)
    46 iflag_cldcon=3
     46iflag_cldcon=6
    4747fact_cldcon=1.
    4848facttemps=0.
    … …  
    5151#calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
    5252ok_newmicro=y
    53 iflag_ratqs=0
    54 ratqsbas=0.005
    55 ratqshaut=0.33
     53iflag_ratqs=2
     54ratqsbas=0.002
     55ratqshaut=0.25
    5656#rayon effectif pour les nuages de glace (D:35)
    5757rad_froid=35
    … …  
    6262new_oliq=y
    6363#
     64rei_min=20.
     65rei_max=61.29
     66#
     67t_glace_min=258.
     68t_glace_max=273.13
     69#
    6470# Flags Convection
    6571#*****************
    … …  
    6975if_ebil=0
    7076#efficacite de precepitation maxim (D:.993)
    71 epmax=0.999
     77epmax=0.997
    7278#ajustement convectif sec au debut de la convection (D: n 0)
    7379ok_adj_ema=n
    74 iflag_clw=1
    75 #fermeture convective (D:1, 1:AR4, 2:ALE and ALP)
    76 iflag_clos=1
    77 #loi de melange a l entrainement  (D:1, 0: plate, 1=AR4:PDF)
     80iflag_clw=0
     81#fermeture convective (D:1, 0:AR4, 2:ALE and ALP)
     82iflag_clos=2
     83#loi de melange a l entrainement  (D:1, 0:AR4=PDF plate, 1=PDF)
    7884iflag_mix=1
    7985#poids des PDFs plate et en cloche (D: 1 0)
    80 qqa1=0.
    81 qqa2=1.
     86qqa1=1.
     87qqa2=0.
    8288# Facteur multiplication des precip convectives dans KE
    8389cvl_corr=1.0
     90#
     91Fmax=0.65
    8492#
    8593# Flags Couche limite
    … …  
    9098nsplit_thermals=1
    9199tau_thermals=0.
    92 iflag_thermals_ed=0
     100iflag_thermals_ed=10
    93101iflag_thermals_optflux=0
    94102#shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada)
    … …  
    99107#ok_kzmin : calcul Kzmin dans la CL de surface
    100108ok_kzmin=y
    101 #couplage ou non avec la convection (0:AR4 et 1:nouvelle physique)
    102 iflag_coupl=0
     109#couplage ou non avec la convection
     110iflag_coupl=5 (0:AR4 et 1:nouvelle physique)
    103111seuil_inversion=-0.08
    104 
     112#rendement sur la puissance fournie par les thermiques a la convection
     113alp_bl_k=0.5
     114#
    105115# Flags wakes
    106116#************
    107117#activation ou non des wakes (D:0, 0:AR4, 1:nouvelle physique)
    108118iflag_wake=0
    109 alp_offset=-0.2
     119alp_offset=0.
    110120
     121### nrlmd le 10/04/2012: Declenchement ###
     122#Déclelenchement par couche-limite (0: deterministe, 1: stochastique)
     123iflag_trig_bl=0
     124#Surface critique au LCL
     125s_trig=1.2e7
     126#Temps de decorrelation entre scenes de cumulus
     127tau_trig_shallow=600
     128#Temps de decorrelation entre scenes de cumulonimbus
     129tau_trig_deep=1800
     130### fin nrlmd le 10/04/2012
     131
     132### nrlmd le 10/04/2012: Fermeture ###
     133#Fermeture par couche-limite (0: deterministe, 1: statistique, 2: statistique+convergence grande échelle)
     134iflag_clos_bl=0
     135### fin nrlmd le 10/04/2012
     136

  • branches/LMDZ_WRFmeas/WRFV3/run/run.def

    r1 r411  
    1212ip_ebil_phy=1
    1313ip_ebil_dyn=1
    14 calend=earth_366d
     14calend=earth_365d
    1515## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
    1616dayref=13
    … …  
    3030periodav= 5.
    3131# Values of solarlong0 are all defaults
    32 #solarlong0 = 0
     32solarlong0 = 0
    3333## activation du calcul d equilibrage de charge
    3434adjust=n
Note: See TracChangeset for help on using the changeset viewer.