source: CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXPBIOS/PARAM/orchidee.def @ 738

Last change on this file since 738 was 738, checked in by mafoipsl, 15 years ago

Commit again to keep history : Adapt EXPBIOS.
EXPBIOS is a new configuration to activate BIO-CHEMICAL processes (STOMATE) with Land Use 1860 map.

See test in http://dods.extra.cea.fr/data/p86manci/IPSLCM5/CM5STO/

File size: 39.5 KB
Line 
1#
2#**************************************************************************
3#                    Namelist for ORCHIDEE
4#**************************************************************************
5#
6#
7#**************************************************************************
8#          OPTIONS NOT SET
9#**************************************************************************
10#
11#
12#**************************************************************************
13#          Management of display in the run of ORCHIDEE
14#**************************************************************************
15
16# Model chatting level
17# level of online diagnostics in STOMATE (0-4)
18# With this variable, you can determine how much online information STOMATE
19#  gives during the run. 0 means virtually no info.
20BAVARD = 1
21# default = 1
22
23# Flag for debug information
24# This option allows to switch on the output of debug
25#         information without recompiling the code.
26DEBUG_INFO = n
27#default = n
28
29# ORCHIDEE will print more messages
30# This flag permits to print more debug messages in the run.
31LONGPRINT = n
32#default = n
33
34#---------------------------------------------------------------------
35
36# Should the output follow the ALMA convention
37# If this logical flag is set to true the model
38#  will output all its data according to the ALMA
39#  convention. It is the recommended way to write
40#  data out of ORCHIDEE.
41ALMA_OUTPUT = n
42# default = n
43
44# To reset the time coming from SECHIBA restart file
45# This option allows the model to override the time
46#  found in the restart file of SECHIBA with the time
47#  of the first call. That is the restart time of the GCM.
48SECHIBA_reset_time = n
49# default = n
50
51#**************************************************************************
52#          Files : incoming / forcing / restart /output
53#**************************************************************************
54# Ancillary files :
55#---------------------------------------------------------------------
56
57# Name of file from which the vegetation map is to be read
58# If !IMPOSE_VEG
59# If LAND_USE
60#   default = pft_new.nc
61#   The name of the file to be opened to read a vegetation
62#   map (in pft) is to be given here.
63# If !LAND_USE
64#   default = ../surfmap/carteveg5km.nc
65#   The name of the file to be opened to read the vegetation
66#   map is to be given here. Usualy SECHIBA runs with a 5kmx5km
67#   map which is derived from the IGBP one. We assume that we have
68#   a classification in 87 types. This is Olson modified by Viovy.
69VEGETATION_FILE = PFTmap.20C3M.nc
70
71
72# Name of file from which the bare soil albedo
73# If !IMPOSE_AZE
74# The name of the file to be opened to read the soil types from
75#  which we derive then the bare soil albedos. This file is 1x1
76#  deg and based on the soil colors defined by Wilson and Henderson-Seller.
77SOILALB_FILE = soils_param.nc
78# default = ../surfmap/soils_param.nc
79
80# Name of file from which soil types are read
81# If !IMPOSE_VEG
82# The name of the file to be opened to read the soil types.
83#  The data from this file is then interpolated to the grid of
84#  of the model. The aim is to get fractions for sand loam and
85#  clay in each grid box. This information is used for soil hydrology
86#  and respiration.
87SOILTYPE_FILE = soils_param.nc
88# default = ../surfmap/soils_param.nc
89
90# Name of file from which the reference
91# The name of the file to be opened to read
92#  temperature is read
93#  the reference surface temperature.
94#  The data from this file is then interpolated
95#  to the grid of the model.
96#  The aim is to get a reference temperature either
97#  to initialize the corresponding prognostic model
98#  variable correctly (ok_dgvm = TRUE) or to impose it
99#  as boundary condition (ok_dgvm = FALSE)
100REFTEMP_FILE = reftemp.nc
101# default = reftemp.nc
102
103# Forcing file name
104# Name of file containing the forcing data
105# This is the name of the file which should be opened
106# for reading the forcing data of the dim0 model.
107# The format of the file has to be netCDF and COADS
108# compliant. Cabauw.nc, islscp_for.nc, WG_cru.nc
109FORCING_FILE = islscp_for.nc
110# default = islscp_for.nc
111
112# Input and output restart file for the driver
113#---------------------------------------------------------------------
114
115# Name of restart to READ for initial conditions
116# This is the name of the file which will be opened
117#  to extract the initial values of all prognostic
118#  values of the model. This has to be a netCDF file.
119#  Not truly COADS compliant. NONE will mean that
120#  no restart file is to be expected.
121RESTART_FILEIN = NONE
122# default = NONE
123
124# Name of restart files to be created by the driver
125# This variable give the  name for
126#  the restart file. The restart software within
127#  IOIPSL will add .nc if needed
128RESTART_FILEOUT = driver_rest_out.nc
129# default = driver_rest_out.nc
130
131
132# Input and output restart file for SECHIBA :
133#---------------------------------------------------------------------
134
135# Name of restart to READ for initial conditions
136# This is the name of the file which will be opened
137#  to extract the initial values of all prognostic
138#  values of the model. This has to be a netCDF file.
139#  Not truly COADS compliant. NONE will mean that
140#  no restart file is to be expected.
141SECHIBA_restart_in = NONE
142# default = NONE
143
144# Name of restart files to be created by SECHIBA
145# This variable give the name for the restart files.
146#  The restart software within IOIPSL will add .nc if needed.
147SECHIBA_rest_out = sechiba_rest_out.nc
148# default = sechiba_rest_out.nc
149
150# Input and output restart file for STOMATE :
151#---------------------------------------------------------------------
152
153# Name of restart to READ for initial conditions of STOMATE
154# If STOMATE_OK_STOMATE || STOMATE_WATCHOUT
155# This is the name of the file which will be opened of STOMATE
156#   to extract the initial values of all prognostic values of STOMATE.
157STOMATE_RESTART_FILEIN = NONE
158# default = NONE
159
160# Name of restart files to be created by STOMATE
161# If STOMATE_OK_STOMATE || STOMATE_WATCHOUT
162# This is the name of the file which will be opened
163#        to write the final values of all prognostic values
164#        of STOMATE.
165STOMATE_RESTART_FILEOUT = stomate_rest_out.nc
166# default = stomate_restart.nc
167
168# Forcing files for TESTSTOMATE and FORCESOIL
169#---------------------------------------------------------------------
170
171# Name of STOMATE's forcing file
172# Name that will be given to STOMATE's offline forcing file
173STOMATE_FORCING_NAME = stomate_forcing.nc
174#default = NONE
175
176# Size of STOMATE forcing data in memory (MB)
177# This variable determines how many
178#  forcing states will be kept in memory.
179#  Must be a compromise between memory
180#  use and frequeny of disk access.
181STOMATE_FORCING_MEMSIZE = 50
182# default = 50
183
184# Name of STOMATE's carbon forcing file
185# Name that will be given to STOMATE's carbon offline forcing file
186STOMATE_CFORCING_NAME = stomate_Cforcing.nc
187# default = NONE
188
189
190# Produced forcing file name (SECHIBA puis STOMATE) :
191#---------------------------------------------------------------------
192
193# ORCHIDEE will write out its forcing to a file
194# This flag allows to write to a file all the variables
195#  which are used to force the land-surface. The file
196#  has exactly the same format than a normal off-line forcing
197#  and thus this forcing can be used for forcing ORCHIDEE.
198ORCHIDEE_WATCHOUT = n
199# default = n
200
201# Filenane for the ORCHIDEE forcing file
202# If ORCHIDEE_WATCHOUT
203# This is the name of the file in which the
204#  forcing used here will be written for later use.
205WATCHOUT_FILE = orchidee_watchout.nc
206# default = orchidee_watchout.nc
207
208# ORCHIDEE will write out with this frequency
209# If ORCHIDEE_WATCHOUT
210# This flag indicates the frequency of the write of the variables.
211DT_WATCHOUT = 1800
212# default = dt
213
214# STOMATE does minimum service
215# set to TRUE if you want STOMATE to read
216#  and write its start files and keep track
217#  of longer-term biometeorological variables.
218#  This is useful if OK_STOMATE is not set,
219#  but if you intend to activate STOMATE later.
220#  In that case, this run can serve as a
221#  spinup for longer-term biometeorological
222#  variables.
223STOMATE_WATCHOUT = n
224# default = n
225
226# Output file name (SECHIBA and STOMATE) :
227#---------------------------------------------------------------------
228# Name of file in which the output is going
229# This file is going to be created by the model
230#  to be written
231#  and will contain the output from the model.
232#  This file is a truly COADS compliant netCDF file.
233#  It will be generated by the hist software from
234#  the IOIPSL package.
235OUTPUT_FILE = sechiba_history.nc
236# default = cabauw_out.nc
237
238# Flag to switch on histfile 2 for SECHIBA (hi-frequency ?)
239# This Flag switch on the second SECHIBA writing for hi (or low)
240#  frequency writing. This second output is optional and not written
241#  by default.
242SECHIBA_HISTFILE2 = FALSE
243# default  = FALSE
244
245# Name of file in which the output number 2 is going
246#   to be written
247# If SECHIBA_HISTFILE2
248# This file is going to be created by the model
249#   and will contain the output 2 from the model.
250SECHIBA_OUTPUT_FILE2 = sechiba_out_2.nc
251# default  = sechiba_out_2.nc
252
253# Name of file in which STOMATE's output is going to be written
254# This file is going to be created by the model
255#  and will contain the output from the model.
256#  This file is a truly COADS compliant netCDF file.
257#  It will be generated by the hist software from
258#  the IOIPSL package.
259STOMATE_OUTPUT_FILE = stomate_history.nc
260# default = stomate_history.nc
261
262# Write levels for outputs files (number of variables) :
263#---------------------------------------------------------------------
264
265# SECHIBA history output level (0..10)
266# Chooses the list of variables in the history file.
267#  Values between 0: nothing is written; 10: everything is
268#  written are available More details can be found on the web under documentation.
269#  web under documentation.
270SECHIBA_HISTLEVEL = 5
271# default = 5
272
273# SECHIBA history 2 output level (0..10)
274# If SECHIBA_HISTFILE2
275# Chooses the list of variables in the history file.
276#   Values between 0: nothing is written; 10: everything is
277#   written are available More details can be found on the web under documentation.
278#   web under documentation.
279# First level contains all ORCHIDEE outputs.
280SECHIBA_HISTLEVEL2 = 1
281# default = 1
282
283# STOMATE history output level (0..10)
284#  0: nothing is written; 10: everything is written
285STOMATE_HISTLEVEL = 10
286# default = 10
287
288# Write frequency for output files (SECHIBA in seconds et
289# STOMATE in days) :
290#---------------------------------------------------------------------
291# Frequency in seconds at which to WRITE output
292# This variables gives the frequency the output of
293#  the model should be written into the netCDF file.
294#  It does not affect the frequency at which the
295#  operations such as averaging are done.
296WRITE_STEP = 86400.0
297# default = 86400.0
298
299# Frequency in seconds at which to WRITE output
300# If SECHIBA_HISTFILE2
301# This variables gives the frequency the output 2 of
302#   the model should be written into the netCDF file.
303#   It does not affect the frequency at which the
304#   operations such as averaging are done.
305#   That is IF the coding of the calls to histdef
306#   are correct !
307WRITE_STEP2 = 1800.0
308# default = 1800.0
309
310# STOMATE history time step (d)
311# Time step of the STOMATE history file
312# Care : this variable must be higher than DT_SLOW
313STOMATE_HIST_DT = 10.
314# default = 10.
315
316#**************************************************************************
317#                             Area location
318#**************************************************************************
319#  The model will use the smalest regions from
320#  region specified here and the one of the forcing file.
321
322# Western limit of region
323# Western limit of the region we are
324#  interested in. Between -180 and +180 degrees
325LIMIT_WEST = -180.
326# default = -180.
327
328# Eastern limit of region
329# Eastern limit of the region we are
330#  interested in. Between -180 and +180 degrees
331LIMIT_EAST = 180.
332# default = 180.
333
334# Northern limit of region
335# Northern limit of the region we are
336#  interested in. Between +90 and -90 degrees
337LIMIT_NORTH = 90.
338# default = 90.
339
340# Southern limit of region
341# Southern limit of the region we are
342#  interested in. Between 90 and -90 degrees
343LIMIT_SOUTH = -90.
344# default = -90.
345
346##**************************************************************************
347#                       Simulation parameters
348#**************************************************************************
349
350# method of forcing
351# A method is proposed by which the first atmospheric
352#   level is not directly forced by observations but
353#   relaxed with a time constant towards observations.
354#   For the moment the methods tends to smooth too much
355#   the diurnal cycle and introduces a time shift.
356#   A more sophisticated method is needed.
357RELAXATION = n
358# default = n
359
360# Time constant of the relaxation layer RELAXATION
361# The time constant associated to the atmospheric
362#  conditions which are going to be computed
363#  in the relaxed layer. To avoid too much
364#  damping the value should be larger than 1000.
365RELAX_A = 1000.
366# default = 1000.0
367
368# Height at which T and Q are given
369# The atmospheric variables (temperature and specific
370#  humidity) are measured at a specific level.
371#  The height of this level is needed to compute
372#  correctly the turbulent transfer coefficients.
373#  Look at the description of the forcing
374#  DATA for the correct value.
375HEIGHT_LEV1 = 2.0
376# default  = 2.0
377
378# Height at which the wind is given
379# The height at which wind is needed to compute
380#  correctly the turbulent transfer coefficients.
381HEIGHT_LEVW = 10.0
382# default  = 10.0
383
384#---------------------------------------------------------------------
385# Weather generator or not :
386#---------------------------------------------------------------------
387
388# Allow weather generator to create data.
389# This flag allows the forcing-reader to generate
390#  synthetic data if the data in the file is too sparse
391#  and the temporal resolution would not be enough to
392#  run the model.
393ALLOW_WEATHERGEN = n
394# default = n
395
396# North-South Resolution
397# If ALLOW_WEATHERGEN
398# North-South Resolution of the region we are
399#  interested in. In degrees
400MERID_RES = 2.
401# default = 2.
402
403# East-West Resolution
404# If ALLOW_WEATHERGEN
405# East-West Resolution of the region we are
406#  interested in. In degrees
407ZONAL_RES = 2.
408# default = 2.
409
410# Use prescribed values
411# If ALLOW_WEATHERGEN
412# If this is set to 1, the weather generator
413#   uses the monthly mean values for daily means.
414#   If it is set to 0, the weather generator
415#   uses statistical relationships to derive daily
416#   values from monthly means.
417IPPREC = 0
418# default = 0
419
420# Interpolation  or not IF split is larger than 1
421# Choose IF you wish to interpolate linearly or not.
422NO_INTER = y
423INTER_LIN = n
424# default :
425#  NO_INTER = y
426#  INTER_LIN = n
427
428# Exact monthly precipitation
429# If ALLOW_WEATHERGEN
430# If this is set to y, the weather generator
431#   will generate pseudo-random precipitations
432#   whose monthly mean is exactly the prescribed one.
433#   In this case, the daily precipitation (for rainy
434#   days) is constant (that is, some days have 0 precip,
435#   the other days have precip = Precip_month/n_precip,
436#   where n_precip is the prescribed number of rainy days
437#   per month).
438WEATHGEN_PRECIP_EXACT = n
439# default = n
440
441# Calling frequency of weather generator (s)
442# Determines how often the weather generator
443#  is called (time step in s). Should be equal
444#  to or larger than Sechiba's time step (say,
445#  up to 6 times Sechiba's time step or so).
446DT_WEATHGEN = 1800.
447# default = 1800.
448
449# Conserve net radiation in the forcing
450# When the interpolation is used the net radiation
451#  provided by the forcing is not conserved anymore.
452#  This should be avoided and thus this option should
453#  be TRUE (y).
454#  This option is not used for short-wave if the
455#  time-step of the forcing is longer than an hour.
456#  It does not make sense to try and reconstruct
457#  a diurnal cycle and at the same time conserve the
458#  incoming solar radiation.
459NETRAD_CONS = y
460# default = y
461
462# Write weather from generator into a forcing file
463# This flag makes the weather generator dump its
464#  generated weather into a forcing file which can
465#  then be used to get the same forcing on different
466#  machines. This only works correctly if there is
467#  a restart file (otherwise the forcing at the first
468#  time step is slightly wrong).
469DUMP_WEATHER = n
470# default = n
471
472# Name of the file that contains
473#  the weather from generator
474# If DUMP_WEATHER
475DUMP_WEATHER_FILE = weather_dump.nc
476# default = 'weather_dump.nc'
477
478# Dump weather data on gathered grid
479# If 'y', the weather data are gathered
480#  for all land points.
481# If DUMP_WEATHER
482DUMP_WEATHER_GATHERED = y
483# default = y
484
485
486# Read Orbital Parameters
487
488# Eccentricity Effect
489# Use prescribed values
490# IF ALLOW_WEATHERGEN
491ECCENTRICITY = 0.016724
492# default = 0.016724
493
494# Longitude of perihelie
495# Use prescribed values
496# If ALLOW_WEATHERGEN
497PERIHELIE = 102.04
498# default = 102.04
499
500# Use prescribed values
501# If ALLOW_WEATHERGEN
502OBLIQUITY = 23.446
503# default = 23.446
504
505#**************************************************************************
506# length of simulation :
507#---------------------------------------------------------------------
508# Length of the integration in time.
509# Length of integration. By default the entire length
510#        of the forcing is used. The FORMAT of this date can
511#        be either of the following :
512# n   : time step n within the forcing file
513# nS  : n seconds after the first time-step in the file
514# nD  : n days after the first time-step
515# nM  : n month after the first time-step (year of 365 days)
516# nY  : n years after the first time-step (year of 365 days)
517#        Or combinations :
518# nYmM: n years and m month
519TIME_LENGTH = default
520# default = depend on the time length and the number of time step in forcing file
521#         = itau_len = itau_fin-itau_dep
522
523
524# split time step :
525#---------------------------------------------------------------------
526
527# Splits the timestep imposed by the forcing
528# With this value the time step of the forcing
529#  will be devided. In principle this can be run
530#  in explicit mode but it is strongly suggested
531#  to use the implicit method so that the
532#  atmospheric forcing has a smooth evolution.
533SPLIT_DT = 12
534# default = 12
535
536#  Time in the forcing file at which the model is started.
537#  This time give the point in time at which the model
538#  should be started.
539#  If exists, the date of the restart file is use.
540#  The FORMAT of this date can be either of the following :
541#  n   : time step n within the forcing file
542#  nS  : n seconds after the first time-step in the file
543#  nD  : n days after the first time-step
544#  nM  : n month after the first time-step (year of 365 days)
545#  nY  : n years after the first time-step (year of 365 days)
546#      Or combinations :
547#  nYmM: n years and m month
548TIME_SKIP = 0
549# default = 0
550
551# Number of time steps per year for carbon spinup.
552FORCESOIL_STEP_PER_YEAR = 12
553# default = 12
554
555# Number of years saved for carbon spinup.
556FORCESOIL_NB_YEAR = 1
557# default = 1
558
559# Spread the precipitation.
560# Spread the precipitaiton over n steps of the splited forcing time step.
561#  This is ONLY applied if the forcing time step has been splited (SPLIT_DT).
562#  If the value indicated is greater than SPLIT_DT, SPLIT_DT is used for it.
563SPRED_PREC = 1
564# default = 1
565
566
567
568#---------------------------------------------------------------------
569# Parametrization :
570#---------------------------------------------------------------------
571
572# Activate STOMATE?
573# set to TRUE if STOMATE is to be activated
574STOMATE_OK_STOMATE = n
575# default = n
576
577# Activate DGVM?
578# set to TRUE if Dynamic Vegetation DGVM is to be activated
579STOMATE_OK_DGVM = n
580# default = n
581
582# Activate CO2?
583# set to TRUE if photosynthesis is to be activated
584STOMATE_OK_CO2 = n
585# default = n
586
587# Flag to force the value of atmospheric CO2 for vegetation.
588# If this flag is set to true, the ATM_CO2 parameter is used
589#  to prescribe the atmospheric CO2.
590# This Flag is only use in couple mode.
591FORCE_CO2_VEG = FALSE
592# default = FALSE
593
594# Value for atm CO2.
595# If FORCE_CO2_VEG (in not forced mode)
596# Value to prescribe the atm CO2.
597#  For pre-industrial simulations, the value is 286.2 .
598#  348. for 1990 year.
599ATM_CO2 = 350.
600# default = 350.
601
602
603# Index of grid point for online diagnostics
604# This is the index of the grid point which will be used for online diagnostics.
605STOMATE_DIAGPT = 1
606# default = 1
607
608# constant tree mortality
609# If yes, then a constant mortality is applied to trees.
610#  Otherwise, mortality is a function of the trees'
611#  vigour (as in LPJ).
612LPJ_GAP_CONST_MORT = y
613# default = y
614
615# no fire allowed
616# With this variable, you can allow or not
617#  the estimation of CO2 lost by fire
618FIRE_DISABLE = n
619# default = n
620
621
622# parameters describing the surface (vegetation + soil) :
623#---------------------------------------------------------------------
624#
625# Should the vegetation be prescribed
626# This flag allows the user to impose a vegetation distribution
627#  and its characterisitcs. It is espacially interesting for 0D
628#  simulations. On the globe it does not make too much sense as
629#  it imposes the same vegetation everywhere
630IMPOSE_VEG = n
631# default = n
632
633# Flag to use old "interpolation" of vegetation map.
634# IF NOT IMPOSE_VEG and NOT LAND_USE
635#  If you want to recover the old (ie orchidee_1_2 branch)
636#   "interpolation" of vegetation map.
637SLOWPROC_VEGET_OLD_INTERPOL = n
638# default = n
639
640# Vegetation distribution within the mesh (0-dim mode)
641# If IMPOSE_VEG
642# The fraction of vegetation is read from the restart file. If
643#  it is not found there we will use the values provided here.
644SECHIBA_VEG__01 = 0.2
645SECHIBA_VEG__02 = 0.0
646SECHIBA_VEG__03 = 0.0
647SECHIBA_VEG__04 = 0.0
648SECHIBA_VEG__05 = 0.0
649SECHIBA_VEG__06 = 0.0
650SECHIBA_VEG__07 = 0.0
651SECHIBA_VEG__08 = 0.0
652SECHIBA_VEG__09 = 0.0
653SECHIBA_VEG__10 = 0.8
654SECHIBA_VEG__11 = 0.0
655SECHIBA_VEG__12 = 0.0
656SECHIBA_VEG__13 = 0.0
657# default = 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8, 0.0, 0.0, 0.0
658
659# Maximum vegetation distribution within the mesh (0-dim mode)
660# If IMPOSE_VEG
661# The fraction of vegetation is read from the restart file. If
662#  it is not found there we will use the values provided here.
663SECHIBA_VEGMAX__01 = 0.2
664SECHIBA_VEGMAX__02 = 0.0
665SECHIBA_VEGMAX__03 = 0.0
666SECHIBA_VEGMAX__04 = 0.0
667SECHIBA_VEGMAX__05 = 0.0
668SECHIBA_VEGMAX__06 = 0.0
669SECHIBA_VEGMAX__07 = 0.0
670SECHIBA_VEGMAX__08 = 0.0
671SECHIBA_VEGMAX__09 = 0.0
672SECHIBA_VEGMAX__10 = 0.8
673SECHIBA_VEGMAX__11 = 0.0
674SECHIBA_VEGMAX__12 = 0.0
675SECHIBA_VEGMAX__13 = 0.0
676# default = 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8, 0.0, 0.0, 0.0
677
678# LAI for all vegetation types (0-dim mode)
679# If IMPOSE_VEG
680# The maximum LAI used in the 0dim mode. The values should be found
681#  in the restart file. The new values of LAI will be computed anyway
682#  at the end of the current day. The need for this variable is caused
683#  by the fact that the model may stop during a day and thus we have not
684#  yet been through the routines which compute the new surface conditions.
685SECHIBA_LAI__01 = 0.
686SECHIBA_LAI__02 = 8.
687SECHIBA_LAI__03 = 8.
688SECHIBA_LAI__04 = 4.
689SECHIBA_LAI__05 = 4.5
690SECHIBA_LAI__06 = 4.5
691SECHIBA_LAI__07 = 4.
692SECHIBA_LAI__08 = 4.5
693SECHIBA_LAI__09 = 4.
694SECHIBA_LAI__10 = 2.
695SECHIBA_LAI__11 = 2.
696SECHIBA_LAI__12 = 2.
697SECHIBA_LAI__13 = 2.
698# default = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2.
699
700# Height for all vegetation types (m)
701# If IMPOSE_VEG
702# The height used in the 0dim mode. The values should be found
703#  in the restart file. The new values of height will be computed anyway
704#  at the end of the current day. The need for this variable is caused
705#  by the fact that the model may stop during a day and thus we have not
706#  yet been through the routines which compute the new surface conditions.
707SLOWPROC_HEIGHT__01 = 0.
708SLOWPROC_HEIGHT__02 = 50.
709SLOWPROC_HEIGHT__03 = 50.
710SLOWPROC_HEIGHT__04 = 30.
711SLOWPROC_HEIGHT__05 = 30.
712SLOWPROC_HEIGHT__06 = 30.
713SLOWPROC_HEIGHT__07 = 20.
714SLOWPROC_HEIGHT__08 = 20.
715SLOWPROC_HEIGHT__09 = 20.
716SLOWPROC_HEIGHT__10 = .2
717SLOWPROC_HEIGHT__11 = .2
718SLOWPROC_HEIGHT__12 = .4
719SLOWPROC_HEIGHT__13 = .4
720# default = 0., 30., 30., 20., 20., 20., 15., 15., 15., .5, .6, 1.0, 1.0
721
722
723# Fraction of the 3 soil types (0-dim mode)
724# If IMPOSE_VEG
725# Determines the fraction for the 3 soil types
726#  in the mesh in the following order : sand loam and clay.
727SOIL_FRACTIONS__01 = 0.28
728SOIL_FRACTIONS__02 = 0.52
729SOIL_FRACTIONS__03 = 0.20
730# default = 0.28, 0.52, 0.20
731
732# Temperature used for the initial guess of LAI
733# If there is no LAI in the restart file, we may need
734#  a temperature that is used to guess the initial LAI.
735SLOWPROC_LAI_TEMPDIAG = 280.
736# default = 280.
737
738# Soil level (m) used for canopy development
739# If STOMATE is not activated.
740# The temperature at this soil depth is used to determine the LAI when
741#   STOMATE is not activated.
742SECHIBA_ZCANOP = 0.5
743# default = 0.5
744
745# Fraction of other surface types within the mesh (0-dim mode)
746# If IMPOSE_VEG
747# The fraction of ice, lakes, etc. is read from the restart file. If
748#  it is not found there we will use the values provided here.
749#  For the moment, there is only ice.
750SECHIBA_FRAC_NOBIO = 0.0
751# default = 0.0
752
753# Fraction of the clay fraction (0-dim mode)
754# If IMPOSE_VEG
755# Determines the fraction of clay in the grid box.
756CLAY_FRACTION = 0.2
757# default = 0.2
758
759# Should the surface parameters be prescribed
760# This flag allows the user to impose the surface parameters
761#  (Albedo Roughness and Emissivity). It is espacially interesting for 0D
762#  simulations. On the globe it does not make too much sense as
763#  it imposes the same vegetation everywhere
764IMPOSE_AZE = n
765# default = n
766
767# Emissivity of the surface for LW radiation
768# If IMPOSE_AZE
769# The surface emissivity used for compution the LE emission
770#  of the surface in a 0-dim version. Values range between
771#  0.97 and 1.. The GCM uses 0.98.
772CONDVEG_EMIS = 1.0
773# default = 1.0
774
775# SW visible albedo for the surface
776# If IMPOSE_AZE
777# Surface albedo in visible wavelengths to be used
778#  on the point if a 0-dim version of SECHIBA is used.
779#  Look at the description of the forcing data for
780#  the correct value.
781CONDVEG_ALBVIS = 0.25
782# default = 0.25
783
784# SW near infrared albedo for the surface
785# If IMPOSE_AZE
786# Surface albedo in near infrared wavelengths to be used
787#  on the point if a 0-dim version of SECHIBA is used.
788#  Look at the description of the forcing data for
789#  the correct value.
790CONDVEG_ALBNIR = 0.25
791# default = 0.25
792
793# Average method for z0
794# If this flag is set to true (y) then the neutral Cdrag
795#  is averaged instead of the log(z0). This should be
796#  the prefered option. We still wish to keep the other
797#  option so we can come back if needed. If this is
798#  desired then one should set Z0CDRAG_AVE = n
799Z0CDRAG_AVE = y
800# default = y
801
802# Surface roughness (m)
803# If IMPOSE_AZE
804# Surface rougness to be used on the point if a 0-dim version
805#  of SECHIBA is used. Look at the description of the forcing 
806#  data for the correct value.
807CONDVEG_Z0 = 0.15
808# default = 0.15_stnd
809
810# Height to be added to the height of the first level (m)
811# If IMPOSE_AZE
812# ORCHIDEE assumes that the atmospheric level height is counted
813#  from the zero wind level. Thus to take into account the roughness
814#  of tall vegetation we need to correct this by a certain fraction
815#  of the vegetation height. This is called the roughness height in
816#  ORCHIDEE talk.
817ROUGHHEIGHT = 0.0
818# default = 0.0
819
820# The snow albedo used by SECHIBA
821# This option allows the user to impose a snow albedo.
822#  Default behaviour is to use the model of snow albedo
823#  developed by Chalita (1993).
824CONDVEG_SNOWA = default
825# default = use the model of snow albedo developed by Chalita
826
827# Switch bare soil albedo dependent (if TRUE) on soil wetness
828# If TRUE, the model for bare soil albedo is the old formulation.
829#  Then it depend on the soil dry or wetness. If FALSE, it is the
830#  new computation that is taken, it is only function of soil color.
831ALB_BARE_MODEL = FALSE
832# default = FALSE
833
834# Initial snow mass if not found in restart
835# The initial value of snow mass if its value is not found
836#   in the restart file. This should only be used if the model is
837#   started without a restart file.
838HYDROL_SNOW = 0.0
839# default = 0.0
840
841
842# Initial snow age if not found in restart
843# The initial value of snow age if its value is not found
844#  in the restart file. This should only be used if the model is
845#  started without a restart file.
846HYDROL_SNOWAGE = 0.0
847# default = 0.0
848
849# Initial snow amount on ice, lakes, etc. if not found in restart
850# The initial value of snow if its value is not found
851#  in the restart file. This should only be used if the model is
852#  started without a restart file.
853HYDROL_SNOW_NOBIO = 0.0
854# default = 0.0
855
856# Initial snow age on ice, lakes, etc. if not found in restart
857# The initial value of snow age if its value is not found
858#  in the restart file. This should only be used if the model is
859#  started without a restart file.
860HYDROL_SNOW_NOBIO_AGE = 0.0
861# default = 0.0
862
863# Initial dry soil height if not found in restart for ORCHIDEE_1.3 to 1.5 Tags only.
864# The initial value of dry soil height if its value is not found
865# in the restart file. This should only be used if the model is
866# started without a restart file.
867HYDROL_HDRY = 0.0
868# default = 0.0
869
870# Initial soil moisture stress if not found in restart
871# The initial value of soil moisture stress if its value is not found
872#  in the restart file. This should only be used if the model is
873#  started without a restart file.
874HYDROL_HUMR = 1.0
875# default = 1.0
876
877# Total depth of soil reservoir
878HYDROL_SOIL_DEPTH = 2.
879# default = 2.
880
881# Initial restart deep soil moisture if not found in restart
882# The initial value of deep soil moisture if its value is not found
883#  in the restart file. This should only be used if the model is
884#  started without a restart file. Default behaviour is a saturated soil.
885HYDROL_BQSB = default
886# default = Maximum quantity of water (Kg/M3) * Total depth of soil reservoir = 150. * 2
887
888# Initial upper soil moisture if not found in restart
889# The initial value of upper soil moisture if its value is not found
890#  in the restart file. This should only be used if the model is
891#  started without a restart file.
892HYDROL_GQSB = 0.0
893# default = 0.0
894
895# Initial upper reservoir depth if not found in restart
896# The initial value of upper reservoir depth if its value is not found
897#  in the restart file. This should only be used if the model is
898#  started without a restart file.
899HYDROL_DSG = 0.0
900# default = 0.0
901
902# Initial dry soil above upper reservoir if not found in restart
903# The initial value of dry soil above upper reservoir if its value
904#  in the restart file. This should only be used if the model is
905#  started without a restart file. The default behaviour
906#  is to compute it from the variables above. Should be OK most of
907#  the time.
908HYDROL_DSP = default
909# default = Total depth of soil reservoir - HYDROL_BQSB / Maximum quantity of water (Kg/M3) = 0.0
910
911# Initial water on canopy if not found in restart
912# The initial value of moisture on canopy if its value
913#  in the restart file. This should only be used if the model is
914#  started without a restart file.
915HYDROL_QSV = 0.0
916# default = 0.0
917
918# Soil moisture on each soil tile and levels
919# The initial value of mc if its value is not found
920#  in the restart file. This should only be used if the model is
921#  started without a restart file.
922HYDROL_MOISTURE_CONTENT = 0.3
923# default = 0.3
924
925# US_NVM_NSTM_NSLM
926# The initial value of us (relative moisture) if its value is not found
927#  in the restart file. This should only be used if the model is
928#  started without a restart file.
929US_INIT = 0.0
930# default = 0.0
931
932# Coefficient for free drainage at bottom
933# The initial value of free drainage if its value is not found
934#  in the restart file. This should only be used if the model is
935#  started without a restart file.
936FREE_DRAIN_COEF = 1.0, 1.0, 1.0
937# default = 1.0, 1.0, 1.0
938
939# Bare soil evap on each soil if not found in restart
940# The initial value of bare soils evap if its value is not found
941#  in the restart file. This should only be used if the model is
942#  started without a restart file.
943EVAPNU_SOIL = 0.0
944# default = 0.0
945
946
947# Initial temperature if not found in restart
948# The initial value of surface temperature if its value is not found
949#  in the restart file. This should only be used if the model is
950#  started without a restart file.
951ENERBIL_TSURF = 280.
952# default = 280.
953
954# Initial Soil Potential Evaporation
955# The initial value of soil potential evaporation if its value
956#  is not found in the restart file. This should only be used if
957#  the model is started without a restart file.
958ENERBIL_EVAPOT = 0.0
959# default = 0.0
960
961# Initial soil temperature profile if not found in restart
962# The initial value of the temperature profile in the soil if
963#   its value is not found in the restart file. This should only
964#   be used if the model is started without a restart file. Here
965#   we only require one value as we will assume a constant
966#   throughout the column.
967THERMOSOIL_TPRO = 280.
968# default = 280.
969
970# Initial leaf CO2 level if not found in restart
971# The initial value of leaf_ci if its value is not found
972#  in the restart file. This should only be used if the model is
973#  started without a restart file.
974DIFFUCO_LEAFCI = 233.
975# default = 233.
976
977
978# Keep cdrag coefficient from gcm.
979# Set to .TRUE. if you want q_cdrag coming from GCM.
980#  Keep cdrag coefficient from gcm for latent and sensible heat fluxes.
981#  TRUE if q_cdrag on initialization is non zero (FALSE for off-line runs).
982CDRAG_FROM_GCM = y
983# default =  IF q_cdrag == 0 ldq_cdrag_from_gcm = .FALSE. ELSE .TRUE.
984
985
986# Artificial parameter to increase or decrease canopy resistance
987# Add from Nathalie - the 28 of March 2006 - advice from Fred Hourdin
988# By PFT.
989RVEG_PFT = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5
990# default = 1.
991
992
993# Interception reservoir coefficient.
994# Transforms leaf area index into size of interception reservoir
995#  for slowproc_derivvar or stomate.
996SECHIBA_QSINT = 0.02
997# default = 0.1
998
999#**************************************************************************
1000# LAI
1001#**************************************************************************
1002
1003# Read the LAI map
1004# It is possible to read a 12 month LAI map which will
1005#  then be interpolated to daily values as needed.
1006#  If n => type_of_lai (constant_veg.f90)
1007#     - mean    : lai(ji,jv) = undemi * (llaimax(jv) + llaimin(jv))
1008#     - inter   : llaimin(jv) + tempfunc(stempdiag(ji,lcanop)) * (llaimax(jv) - llaimin(jv))
1009LAI_MAP = n
1010# default = n
1011
1012# Name of file from which the vegetation map is to be read
1013# If LAI_MAP
1014# The name of the file to be opened to read the LAI
1015#  map is to be given here. Usualy SECHIBA runs with a 5kmx5km
1016#  map which is derived from a Nicolas VIOVY one.
1017LAI_FILE = lai2D.nc
1018# default = ../surfmap/lai2D.nc
1019
1020# Flag to use old "interpolation" of LAI
1021# If LAI_MAP
1022# If you want to recover the old (ie orchidee_1_2 branch)
1023# "interpolation" of LAI map.
1024SLOWPROC_LAI_OLD_INTERPOL = n
1025# default = n
1026
1027#**************************************************************************
1028# LAND_USE
1029#**************************************************************************
1030
1031# Read a land_use vegetation map
1032# pft values are needed, max time axis is 293
1033LAND_USE = y
1034# default = n
1035
1036# Year of the land_use vegetation map readed
1037# year off the pft map
1038# If LAND_USE (11 = 1860 - 1850 +1 for PFTmap.20C3M.nc)
1039VEGET_YEAR = 11
1040# default = 282
1041
1042# Update vegetation frequency (since 2.0 version)
1043# The veget datas will be update each this time step.
1044# If LAND_USE
1045VEGET_UPDATE = 0Y
1046# default = 1Y
1047
1048# treat land use modifications
1049# With this variable, you can use a Land Use map
1050# to simulate anthropic modifications such as   
1051# deforestation.                               
1052# If LAND_USE
1053LAND_COVER_CHANGE = n
1054# default = y
1055
1056#**************************************************************************
1057
1058# agriculture allowed?
1059# With this variable, you can determine
1060#  whether agriculture is allowed
1061AGRICULTURE = y
1062# default = y
1063
1064# Harvert model for agricol PFTs.
1065# Compute harvest above ground biomass for agriculture.
1066# Change daily turnover.
1067HARVEST_AGRI = y
1068# default = y
1069
1070# herbivores allowed?
1071# With this variable, you can activate herbivores
1072HERBIVORES = n
1073# default = n
1074
1075# treat expansion of PFTs across a grid cell?
1076# With this variable, you can determine
1077#  whether we treat expansion of PFTs across a
1078#  grid cell.
1079TREAT_EXPANSION = n
1080# default = n
1081
1082#**************************************************************************
1083
1084# Time within the day simulated
1085# This is the time spent simulating the current day. This variable is
1086#  prognostic as it will trigger all the computations which are
1087#  only done once a day.
1088SECHIBA_DAY = 0.0
1089# default = 0.0
1090
1091# Time step of STOMATE and other slow processes
1092# Time step (s) of regular update of vegetation
1093#  cover, LAI etc. This is also the time step
1094#  of STOMATE.
1095DT_SLOW = 86400.
1096# default = un_jour = 86400.
1097
1098#**************************************************************************
1099
1100# Allows to switch on the multilayer hydrology of CWRR
1101# This flag allows the user to decide if the vertical
1102#  hydrology should be treated using the multi-layer
1103#  diffusion scheme adapted from CWRR by Patricia de Rosnay.
1104#  by default the Choisnel hydrology is used.
1105HYDROL_CWRR = n
1106# default = n
1107
1108# do horizontal diffusion?
1109# If TRUE, then water can diffuse horizontally between
1110#  the PFTs' water reservoirs.
1111HYDROL_OK_HDIFF = n
1112# default = n
1113 
1114
1115# time scale (s) for horizontal diffusion of water
1116# If HYDROL_OK_HDIFF
1117# Defines how fast diffusion occurs horizontally between
1118#  the individual PFTs' water reservoirs. If infinite, no
1119#  diffusion.
1120HYDROL_TAU_HDIFF = 1800.
1121# default = 86400.
1122
1123# Percent of precip that is not intercepted by the canopy (only for TAG 1.6).
1124# During one rainfall event, PERCENT_THROUGHFALL% of the incident rainfall
1125#  will get directly to the ground without being intercepted.
1126PERCENT_THROUGHFALL = 30.
1127# default = 30.
1128
1129# Percent by PFT of precip that is not intercepted by the canopy (since TAG 1.8).
1130# During one rainfall event, PERCENT_THROUGHFALL_PFT% of the incident rainfall
1131#  will get directly to the ground without being intercepted, for each PFT..
1132PERCENT_THROUGHFALL_PFT = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30.
1133# default = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30.
1134
1135
1136# Decides if we route the water or not
1137# This flag allows the user to decide if the runoff
1138#  and drainage should be routed to the ocean
1139#  and to downstream grid boxes.
1140RIVER_ROUTING = n
1141# default = n
1142
1143# Name of file which contains the routing information
1144# The file provided here should allow the routing module to
1145#  read the high resolution grid of basins and the flow direction
1146#  from one mesh to the other.
1147ROUTING_FILE = routing.nc
1148# default = routing.nc
1149
1150# Time step of the routing scheme
1151# If RIVER_ROUTING
1152# This values gives the time step in seconds of the routing scheme.
1153#   It should be multiple of the main time step of ORCHIDEE. One day
1154#   is a good value.
1155ROUTING_TIMESTEP = 86400
1156# default = 86400
1157
1158# Number of rivers
1159# If RIVER_ROUTING
1160# This parameter chooses the number of largest river basins
1161#  which should be treated as independently as rivers and not
1162#  flow into the oceans as diffusion coastal flow.
1163ROUTING_RIVERS = 50
1164# default = 50
1165
1166# Should we compute an irrigation flux
1167# This parameters allows the user to ask the model
1168#  to compute an irigation flux. This performed for the
1169#  on very simple hypothesis. The idea is to have a good
1170#  map of irrigated areas and a simple function which estimates
1171#  the need to irrigate.
1172DO_IRRIGATION = n
1173# default = n
1174
1175# Name of file which contains the map of irrigated areas
1176# If IRRIGATE
1177# The name of the file to be opened to read the field
1178#  with the area in m^2 of the area irrigated within each
1179#  0.5 0.5 deg grid box. The map currently used is the one
1180#  developed by the Center for Environmental Systems Research
1181#  in Kassel (1995).
1182IRRIGATION_FILE = irrigated.nc
1183# default = irrigated.nc
1184
1185# Should we include floodplains
1186# This parameters allows the user to ask the model
1187#  to take into account the flood plains and return
1188#  the water into the soil moisture. It then can go
1189#  back to the atmopshere. This tried to simulate
1190#  internal deltas of rivers.
1191DO_FLOODPLAINS = n
1192# default = n
1193
1194#**************************************************************************
Note: See TracBrowser for help on using the repository browser.